The stability of slopes and tunnels is controlled by rock discontinuities,and the rock discontinuities roughness and the sliding direction play a signifcant role in shear failure.However,three-dimensional roughness ev...The stability of slopes and tunnels is controlled by rock discontinuities,and the rock discontinuities roughness and the sliding direction play a signifcant role in shear failure.However,three-dimensional roughness evaluation considering shear directions is scare,and the internal shear fracturing processes,micromechanical mechanisms and failure precursor of rock discontinuities are not well understood.Therefore,this study proposes a novel roughness evaluation index to quantitatively analyze the anisotropic characteristics of rock discontinuities.In conjunction with shear tests,a novel 3D-GBM modelling method considering the micromineral constituent and particle size distribution characteristics of granite as well as the geometric shape of discontinuities was realized.The strength,macro and micro-fracture characteristics,visual anisotropic shear evolution process and microfailure mechanism of granite discontinuities at diferent roughness and shear direction were investigated.Finally,the spatial and temporal evolutions of AE parameter b-value and magnitude M were further analyzed to reveal the shear fracture precursor of granite discontinuities.展开更多
Deepwater drilling riser is the key equipment connecting the subsea wellhead and floating drilling platform.Due to complex marine environment,vortex-induced vibration(ViV)will be generated on riser,which will induce f...Deepwater drilling riser is the key equipment connecting the subsea wellhead and floating drilling platform.Due to complex marine environment,vortex-induced vibration(ViV)will be generated on riser,which will induce fatigue failure and even cause unpredictable drilling accidents.Therefore,it is important to study the ViV characteristics of deepwater drilling riser and reveal the main controlling factors for ensuring the safe and efficient operation of deepwater drilling engineering.In this paper,the ViV of deepwater drilling riser is numerically simulated in time domain based on the discrete vortex method(DvM).A hydrodynamic analysis model and governing equation of VIV is proposed with considering the effect of riser motion using DVM and slice method,where the governing equation is solved by Runge-Kutta method.Model validation is performed,which verified the correctness and accuracy of the mechanical model and the solution method.On this basis,the influence of the number of control points,current velocity,riser outer diameter,shear flow and top tension on the ViV characteristics of deepwater drilling risers are discussed in detail.The results show that with the increase of current velocity,the vibration amplitude of deepwater drilling riser decreases obviously,while the vibration frequency increases gradually.However,if the outer diameter of riser increases,the vibration amplitude increases,while the vibration frequency decreases gradually.The top tension also has great influence on the VIV of riser.When the top tension is 1.25 G,the VIV is suppressed to a certain extent.This study has guiding significance for optimal design and engineering control of deepwater drilling riser.展开更多
Cylindrical cross sections are critical components in offshore structures, including jacket platform legs, pipelines, mooring lines, and risers. Thesecylindrical structures are subjected to vortex-induced vibrations (...Cylindrical cross sections are critical components in offshore structures, including jacket platform legs, pipelines, mooring lines, and risers. Thesecylindrical structures are subjected to vortex-induced vibrations (VIV) due to strong ocean currents, where vortices generated during fluid flowresult in significant vibrations in crossflow and in-flow directions. Such vibrations can lead to severe damage to platforms, cables, and risersystems. Consequently, mitigating VIV caused by vortex-induced forces is important. This study investigates the hydrodynamic performance offive strake models relative to a bare cylinder at moderate Reynolds numbers. The models encompass one conventional continuous helical strake(HS) and four helical discrete strake (HDS) with varying segment spacing between the fins. The hydrodynamic performance, specifically liftand drag force coefficients, was computed using a Reynolds averaged Navier –Stokes-based CFD solver and validated with experimentalmeasurements. The conventional HS suppresses 95% of the lift force but increases the drag force by up to a maximum of 48% in measurements.The HDS suppress the lift force by 70%–88% and increase the drag force by 15%–30%, which is less than the increase observed with the HS.Flow visualization showed that HS and HDS cylinders mitigate vortex-induced forces by altering the vortex-shedding pattern along the length ofthe cylinder. The HDS achieves a reduction in drag compared with the conventional continuous HS. The segment spacing is found to significantlyimpact the reduction in vortex-induced forces.展开更多
We propose a symplectic partitioned Runge-Kutta （SPRK） method with eighth-order spatial accuracy based on the extended Hamiltonian system of the acoustic waveequation. Known as the eighth-order NSPRK method, this te...We propose a symplectic partitioned Runge-Kutta （SPRK） method with eighth-order spatial accuracy based on the extended Hamiltonian system of the acoustic waveequation. Known as the eighth-order NSPRK method, this technique uses an eighth-orderaccurate nearly analytic discrete （NAD） operator to discretize high-order spatial differentialoperators and employs a second-order SPRK method to discretize temporal derivatives.The stability criteria and numerical dispersion relations of the eighth-order NSPRK methodare given by a semi-analytical method and are tested by numerical experiments. We alsoshow the differences of the numerical dispersions between the eighth-order NSPRK methodand conventional numerical methods such as the fourth-order NSPRK method, the eighth-order Lax-Wendroff correction （LWC） method and the eighth-order staggered-grid （SG）method. The result shows that the ability of the eighth-order NSPRK method to suppress thenumerical dispersion is obviously superior to that of the conventional numerical methods. Inthe same computational environment, to eliminate visible numerical dispersions, the eighth-order NSPRK is approximately 2.5 times faster than the fourth-order NSPRK and 3.4 timesfaster than the fourth-order SPRK, and the memory requirement is only approximately47.17% of the fourth-order NSPRK method and 49.41% of the fourth-order SPRK method,which indicates the highest computational efficiency. Modeling examples for the two-layermodels such as the heterogeneous and Marmousi models show that the wavefields generatedby the eighth-order NSPRK method are very clear with no visible numerical dispersion.These numerical experiments illustrate that the eighth-order NSPRK method can effectivelysuppress numerical dispersion when coarse grids are adopted. Therefore, this methodcan greatly decrease computer memory requirement and accelerate the forward modelingproductivity. In general, the eighth-order NSPRK method has tremendous potential value forseismic exploration and seismology research.展开更多
According to the innate characteristic of four types of furnace, the copper flash continuous smelting （CFCS） furnace can be considered a synthetic reactor of two relatively independent processes： flash matte smelti...According to the innate characteristic of four types of furnace, the copper flash continuous smelting （CFCS） furnace can be considered a synthetic reactor of two relatively independent processes： flash matte smelting process （FMSP） and copper continuous converting process （CCCP）. Then, the CFCS thermodynamic model was proposed by establishing the multi-phase equilibrium model of FMSP and the local-equilibrium model of CCCP, respectively, and by combining them through the smelting intermediates. Subsequently, the influences of the furnace structures were investigated using the model on the formation of blister copper, the Fe3O4 behavior, the copper loss in slag and the copper recovery rate. The results show that the type D furnace, with double flues and a slag partition wall, is an ideal CFCS reactor compared with the other three types furnaces. For CFCS, it is effective to design a partition wall in the furnace to make FMSP and CCCP perform in two relatively independent zones, respectively, and to make smelting gas and converting gas discharge from respective flues.展开更多
This work is focused on the development of computational algorithms to create a simulator for solving the heat transfer during the continuous casting process of steel. The temperatures and the solid shell thickness pr...This work is focused on the development of computational algorithms to create a simulator for solving the heat transfer during the continuous casting process of steel. The temperatures and the solid shell thickness profiles were calculated and displayed on the screen for a billet through a defined continuous casting plant （CCP）. The algorithms developed to calculate billet temperatures, involve the solutions of the corresponding equations for the heat removal conditions such as radiation, forced convection, and conduction according to the billet posi- tion through the CCP. This is done by a simultaneous comparison with the kinematics model previously developed. A finite difference method known as Cmnk-Nicholson is applied to solve the two-dimensional computational array （2D model）. Enthalpy （HIJ） and temperature （TIL） in every node are updated at each step time. The routines to display the results have been developed using a graphical user interface （GUI） in the programming language C＋＋. Finally, the results obtained are compared with those of industrial trials for the surface temperature of three steel casters with different plant configurations in different casting conditions.展开更多
The aim of this study was to simulate the solidification process of beam blank continuous casting, and then find the reasons for the typical defects of the beam blank. A two-dimensional transient coupled finite elemen...The aim of this study was to simulate the solidification process of beam blank continuous casting, and then find the reasons for the typical defects of the beam blank. A two-dimensional transient coupled finite element model has been developed to compute the temperature and stress profile in beam blank continuous casting. The enthalpy method was used in the heat conduction equation. The thermo-mechanical property in the mushy zone was taken into consideration in this calculation. It is shown that at the mold exit the thickness of the shell had its maximum value at the flange tip and its minimum value at the fillet. The temperature had a great fluctuation on the surface of the beam blank in the secondary cooling zone. At the unbending point, the surface temperature of the web was in the brittleness temperature range under the present condition. To ensure the quality, it is necessary to weaken the intensity of secondary cooling. At the mold exit the equivalent stress and strain have higher values at the flange tip and at the web. From the spray 1 to the unbending point, the maximum values of stress and strain gradually moved to the internal section of the flange tip and the web. However, whenever, there were bigger stress and strain values near the flange tip and the web than in the other parts, it must be very easy to generate cracks at those positions. Now, online verification of this simulation has been developed, which has proved to be very useful and efficient to instruct the practical production of beam blank continuous casting.展开更多
The steady-state temperature field of horizontal core-filling continuous casting （HCFC） for producing copper cladding aluminum rods was simulated by finite element method to investigate the effects of key processing...The steady-state temperature field of horizontal core-filling continuous casting （HCFC） for producing copper cladding aluminum rods was simulated by finite element method to investigate the effects of key processing parameters on the positions of solid-liquid interfaces （SLIs） of copper and aluminum. It is found that mandrel tube length and mean withdrawing speed have significant effects on the SLI positions of both copper and aluminum. Aluminum casting temperature （TAI） （1003-1123 K） and secondary cooling water flux （600-900 L.h-1） have little effect on the SLI of copper but cause the SLI of aluminum to move 2-4 mm. When TA1 is in a range of 1043-1123 K, the liquid aluminum can fill continuously into the pre-solidified copper tube. Based on the numerical simulation, reasonable processing parameters were determined.展开更多
A three-dimensional heat transfer model for continuous steel slab casting has been developed with realistic spray cooling patterns and a coupled microsegregation solidification model that calculates the solidification...A three-dimensional heat transfer model for continuous steel slab casting has been developed with realistic spray cooling patterns and a coupled microsegregation solidification model that calculates the solidification path for multi-component steels. Temperature and composition dependent properties are implemented in a database for 15 chemical species. Considerable effort is made to accurately model the spray cooling heat transfer. Each spray nozzle position and distribution is considered, including variations of the spray patterns with flow rate, and spray overlap. Nozzle type, layout, nozzle-to-slab distance, and spray span and flux are variable. Natural convection, thermal radiation and contact cooling of individual rolls are computed. The present model provides more comprehensive information and realistic slab surface temperatures than results from a model using the 'averaged' treatment of boundary conditions. Cooling operating conditions and parameters of individual spray nozzles can be analyzed to optimize nozzle spray distribution, improve product quality, and troubleshoot issues such as nozzle clogging that may arise during production. One spray cooling correlation is used for the entire machine, achieving as good or better agreement with surface temperature measurements than was found previously for the model using an 'averaged' treatment of boundary conditions and using three machine-segment-dependent correlations.展开更多
Piled embankments are widely used in highway and railway engineering due to their economy and efficiency inovercoming several issues encountered in constructing embankments over weak soils. Soil arching, caused by the...Piled embankments are widely used in highway and railway engineering due to their economy and efficiency inovercoming several issues encountered in constructing embankments over weak soils. Soil arching, caused by the pile-subsoilrelative displacement （△s）, plays an important role in reducing the embankment load falling on weak soil, however, the funda-mental characteristics （e.g., formation and features） of soil arching remain poorly understood. In this study, a series of discreteelement method （DEM） modellings are performed to study the formation and features of soil arching with the variation of As inpiled embankments with or without geosynthetic reinforcement. Firstly, calibration for the modelling parameters is carried out bycomparing the DEM results with the experimental data obtained from the existing literature. Secondly, the analysis of the macro-and micro-behaviours is performed in detail. Finally, a parametric study is conducted in an effort to identify the influences of threekey factors on soil arching： the friction coefficient of the embankment fill （f）, the embankment height （h）, and the pile clear spacing（s-a）. Numerical results indicate that △s is a key factor governing the formation and features of soil arching in embankments. Tobe specific, soil arching gradually evolves from two inclined shear planes at a small △s to a hemispherical arch at a relatively largeAs. Then, with a continuous increase in △s, the soil arching height gradually increases and finally approaches a constant value of0.8（s-a） （i.e., the maximum soil arching height）. For a given case, the higher the soil arching height, the greater the degree of soilarching effect. The parametric study shows that the friction coefficient of the embankment fill has a negligible influence on theformation and features of soil arching. However, embankment height is a key factor governing the formation and features of soilarching. In addition, pile clear spacing has a significant effect on the formation of soil arching, but not on its features.展开更多
A new method called mixed Lagrangian and Eulerian （MiLE） method was used to simulate the continuous casting process in a mold of free-cutting steel 38MnVS.The simulation results are basically in agreement with exper...A new method called mixed Lagrangian and Eulerian （MiLE） method was used to simulate the continuous casting process in a mold of free-cutting steel 38MnVS.The simulation results are basically in agreement with experimental data in the literature,achieving the three-dimensional visualization of temperature distribution,melt flow,shell thickness,and stress distribution of blooms in a mold.It is shown that the flow velocity of steel melt becomes smaller gradually as the casting proceeds.When the flow reaches a certain depth,two types of flow patterns can be observed in the upper zone of the mold.The first flow pattern is to flow downwards,and the second one is to flow upwards to the meniscus.The corner temperature is higher,and the thickness is thinner than those in the mid-face.The effective stress in the corner area is much bigger than that in the mid-face,indicating that the corner area is the dangerous zone of cracking.展开更多
Mold is the heart of the continuous casting machine. Heat transfer and solidification in a water- cooled mold are the most important factors during the continuous casting of steel. For studying the temperature distrib...Mold is the heart of the continuous casting machine. Heat transfer and solidification in a water- cooled mold are the most important factors during the continuous casting of steel. For studying the temperature distribution of a mold wall, a simulated apparatus of mold was designed and experiments were performed by it. The measured results indicated that the mold wall temperature approaches the temperature of cooling-water. An equivalent thermal-conductivity coefficient was proposed and deduced on the basis of the conclusion of the experiments. This coefficient was applied to solve the heat transfer between the melt and cooling water, and to characterize the heat transfer capacity of the mold. By this equivalent thermal-conductivity coefficient, it is very easy and convenient to numerically simulate the solidification process of continuous casting. And the calculation results are in agreement with the experiments. The effects of custing speed and water flow rate on the mold temperature field were also discussed.展开更多
The magnetic,heat transfer and flow phenomenon occurring in the continuous casting process under the mold electromagnetic stirring was further analyzed by solving the 3-D electromagnetic field mathematical model and f...The magnetic,heat transfer and flow phenomenon occurring in the continuous casting process under the mold electromagnetic stirring was further analyzed by solving the 3-D electromagnetic field mathematical model and flow solidification model with finite element method and finite volume method,respectively.The results indicate that the solidified shell thickness located in the effective stirring region fluctuates because of the unsteady scouring under the mold electromagnetic stirring.The maximum rotational velocity is a key parameter to the solidification of the billet when controlling the stirring intensity.When the rotational velocity reaches 0.32 m/s,the mush zone enlarges significantly and the solidification rate is further accelerated.The number of vortexes in the lower recirculation zone is not only two and depends on the stirring parameters.Besides,the secondary flow is closely associated with the solidification.Compared with the results of the model ignoring the influence of solidification on the flow of molten steel,the flow pattern within the lower recirculation region changes dramatically,and thus a coupling analysis of the flow,heat transfer,and solidification is essential when simulating the electromagnetic continuous casting process.展开更多
Heat transfer and deformation of initial solidification shell in soft contact continuous casting mold under high frequency electromagnetic field were analyzed using numerical simulation method; the relative electromag...Heat transfer and deformation of initial solidification shell in soft contact continuous casting mold under high frequency electromagnetic field were analyzed using numerical simulation method; the relative electromagnetic parameters were obtained from the previous studies. Owing to the induction heating of a high frequency electromagnetic field （20 kHz）, the thickness of initial solidification shell decreases, and the temperature of strand surface and slit copper mold increases when compared with the case without the electromagnetic filed. The viscosity of flux de- creases because of the induction heating of the high frequency electromagnetic field, and the dimension of the flux channel increases with electromagnetic pressure; thus, the deformation behavior of initial solidification shell was different before and after the action of high frequency electromagnetic field. Furthermore, the abatement mechanism of oscillation marks under high frequency electromagnetic field was explained.展开更多
A 3D/2D hybrid multi-physical-field mathematical model,which takes into consideration the thermosolutal buoyance,was developed to predict the macrosegregation of gear steel 20CrMnTi continuously cast by a curved bille...A 3D/2D hybrid multi-physical-field mathematical model,which takes into consideration the thermosolutal buoyance,was developed to predict the macrosegregation of gear steel 20CrMnTi continuously cast by a curved billet caster with size of 160 mm×160 mm,and investigated the effect of final electromagnetic stirring(F-EMS)on the fluid flow,heat transfer and solute distribution in the liquid core of continuously cast steel.The results show that the application of F-EMS eliminates the effect of thermosolutal buoyancy on the asymmetric distribution of carbon concentration in the cross section of billet and accelerates the final solidification of resident molten steel in the liquid core of strand,but promotes the negative carbon segregation near the billet center.When the gear steel 20CrMnTi is cast at the temperature of 1803 K and speed of 1.7 m/min,the solidification end advances forward from 9.84 to 9.72 m,and center carbon segregation ratio of billet decreases from 1.24 to 1.17 with the increase in current density of F-EMS from 0 to 350 A.展开更多
The discrete particle method was used to simulate the distribution of gas holdup in a gas-liquid standard Rushton stirred tank. The gas phase was treated as a large number of bubbles and their trajectories were tracke...The discrete particle method was used to simulate the distribution of gas holdup in a gas-liquid standard Rushton stirred tank. The gas phase was treated as a large number of bubbles and their trajectories were tracked with the results of motion equations. The two-way approach was performed to couple the interphase momentum exchange. The turbulent dispersion of bubbles with a size distribution was modeled using a stochastic tracking model, and the added mass force was involved to account for the effect of bubble acceleration on the surrounding fluid. The predicted gas holdup distribution showed that this method could give reasonable prediction comparable to the reported experimental data when the effect of turbulence was took into account in modification for drag coefficient.展开更多
The unsteady turbulent flow during the continuous casting of steel is important, because it influences critical phenomena that affect steel quality. Unsteady three-dimensional flow in the mold region of the liquid poo...The unsteady turbulent flow during the continuous casting of steel is important, because it influences critical phenomena that affect steel quality. Unsteady three-dimensional flow in the mold region of the liquid pool during continuous casting of steel slabs has been computed using realistic geometries starting from the submerged inlet nozzle to the mold. The cassette filter function was used to deal with unsteady Navier-Stokes equation, and then the turbulent flow in the thin slab CCM was simulated with the large eddy simulation method combined with the Smagorinsky sub-grid scale model in this paper. And the model was verified by the Particle Image Velocimetry (PIV) experimental results which was got from a relate scientific literature. In this thesis, by means of LES, the flow characteristics in the thin slab CCM were acquired, such as the vortex distribution, the formation of the large eddy coherent structures, development, shedding and fracture process. In the same time, the turbulent asymmetric distribution was revealed even the nozzle in the centre position. Interactions between the two halves cause large velocity fluctuations near the meniscus. And the vortex is located at the low velocity side adjacent to the SEN. Along with the unsteady time development, the unsteady turbulent large vortex structures of the liquid steel in the CCM presented periodic bias flow distribution, and the period is about 20 seconds.展开更多
In order to continuously simulate multi-pass plate rolling process,a 3-D elastic hollow-roll model was proposed and an auto mesh-refining module with data passing was developed and integrated with FE software,Marc.The...In order to continuously simulate multi-pass plate rolling process,a 3-D elastic hollow-roll model was proposed and an auto mesh-refining module with data passing was developed and integrated with FE software,Marc.The hollow-roll model has equivalent stiffness of bending resistance and deformation to the real solid and much less meshes,so the computational time is greatly reduced.Based on these,the factors influencing plate profile,such as the roll-bending force,initial crown,thermal crown and heat transfer during rolling and inter-pass cooling can be taken into account in the simulation.The auto mesh-refining module with data passing can automatically refine and re-number elements and transfer the nodal and elemental results to the new meshes.Furthermore,the 3-D modeling routine is parametrically developed and can be run independently of Marc pre-processing program.A seven-pass industrial hot rolling process was continuously simulated to validate the accuracy of model.By comparison of the calculated results with the industrial measured data,the rolling force,temperature and plate profile are in good accordance with the measured ones.展开更多
The discrete element method (DEM) was used to simulate the flow characteristic and strength characteristic of the conditioned sands in the earth pressure balance (EPB) tunneling. In the laboratory the conditioned sand...The discrete element method (DEM) was used to simulate the flow characteristic and strength characteristic of the conditioned sands in the earth pressure balance (EPB) tunneling. In the laboratory the conditioned sands were reproduced and the slump test and the direct shear test of the conditioned sands were implemented. A DEM equivalent model that can simulate the macro mechanical characteristic of the conditioned sands was proposed,and the corresponding numerical models of the slump test and the shear test were established. By selecting proper DEM model parameters,the errors of the slump values between the simulation results and the test results are in the range of 10.3%-14.3%,and the error of the curves between the shear displacement and the shear stress calculated with the DEM simulation is 4.68%-16.5% compared with that of the laboratory direct shear test. This illustrates that the proposed DEM equivalent model can approximately simulate the mechanical characteristics of the conditioned sands,which provides the basis for further simulation of the interaction between the conditioned soil and the chamber pressure system of the EPB machine.展开更多
基金Financial support to complete this study was provided in part by National Natural Science Foundation of China(52109119)the Open Research Fund of State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin(China Institute of Water Resources and Hydropower Research)(IWHR-SKL-202202)+1 种基金the China Postdoctoral Science Foundation Project(2022M723408)the Systematic Project of Guangxi Key Laboratory of Disaster Prevention and Engineering Safety(2020ZDK007)。
文摘The stability of slopes and tunnels is controlled by rock discontinuities,and the rock discontinuities roughness and the sliding direction play a signifcant role in shear failure.However,three-dimensional roughness evaluation considering shear directions is scare,and the internal shear fracturing processes,micromechanical mechanisms and failure precursor of rock discontinuities are not well understood.Therefore,this study proposes a novel roughness evaluation index to quantitatively analyze the anisotropic characteristics of rock discontinuities.In conjunction with shear tests,a novel 3D-GBM modelling method considering the micromineral constituent and particle size distribution characteristics of granite as well as the geometric shape of discontinuities was realized.The strength,macro and micro-fracture characteristics,visual anisotropic shear evolution process and microfailure mechanism of granite discontinuities at diferent roughness and shear direction were investigated.Finally,the spatial and temporal evolutions of AE parameter b-value and magnitude M were further analyzed to reveal the shear fracture precursor of granite discontinuities.
基金the financial support from National Key R&D Program of China(Grant number:2024YFC2815100)Natural Science Foundation of China(Grant number:52322110)Beijing Nova Program(Grant number:20230484341).
文摘Deepwater drilling riser is the key equipment connecting the subsea wellhead and floating drilling platform.Due to complex marine environment,vortex-induced vibration(ViV)will be generated on riser,which will induce fatigue failure and even cause unpredictable drilling accidents.Therefore,it is important to study the ViV characteristics of deepwater drilling riser and reveal the main controlling factors for ensuring the safe and efficient operation of deepwater drilling engineering.In this paper,the ViV of deepwater drilling riser is numerically simulated in time domain based on the discrete vortex method(DvM).A hydrodynamic analysis model and governing equation of VIV is proposed with considering the effect of riser motion using DVM and slice method,where the governing equation is solved by Runge-Kutta method.Model validation is performed,which verified the correctness and accuracy of the mechanical model and the solution method.On this basis,the influence of the number of control points,current velocity,riser outer diameter,shear flow and top tension on the ViV characteristics of deepwater drilling risers are discussed in detail.The results show that with the increase of current velocity,the vibration amplitude of deepwater drilling riser decreases obviously,while the vibration frequency increases gradually.However,if the outer diameter of riser increases,the vibration amplitude increases,while the vibration frequency decreases gradually.The top tension also has great influence on the VIV of riser.When the top tension is 1.25 G,the VIV is suppressed to a certain extent.This study has guiding significance for optimal design and engineering control of deepwater drilling riser.
文摘Cylindrical cross sections are critical components in offshore structures, including jacket platform legs, pipelines, mooring lines, and risers. Thesecylindrical structures are subjected to vortex-induced vibrations (VIV) due to strong ocean currents, where vortices generated during fluid flowresult in significant vibrations in crossflow and in-flow directions. Such vibrations can lead to severe damage to platforms, cables, and risersystems. Consequently, mitigating VIV caused by vortex-induced forces is important. This study investigates the hydrodynamic performance offive strake models relative to a bare cylinder at moderate Reynolds numbers. The models encompass one conventional continuous helical strake(HS) and four helical discrete strake (HDS) with varying segment spacing between the fins. The hydrodynamic performance, specifically liftand drag force coefficients, was computed using a Reynolds averaged Navier –Stokes-based CFD solver and validated with experimentalmeasurements. The conventional HS suppresses 95% of the lift force but increases the drag force by up to a maximum of 48% in measurements.The HDS suppress the lift force by 70%–88% and increase the drag force by 15%–30%, which is less than the increase observed with the HS.Flow visualization showed that HS and HDS cylinders mitigate vortex-induced forces by altering the vortex-shedding pattern along the length ofthe cylinder. The HDS achieves a reduction in drag compared with the conventional continuous HS. The segment spacing is found to significantlyimpact the reduction in vortex-induced forces.
基金This research was supported by the National Natural Science Foundation of China （Nos. 41230210 and 41204074）, the Science Foundation of the Education Department of Yunnan Province （No. 2013Z152）, and Statoil Company （Contract No. 4502502663）.
文摘We propose a symplectic partitioned Runge-Kutta （SPRK） method with eighth-order spatial accuracy based on the extended Hamiltonian system of the acoustic waveequation. Known as the eighth-order NSPRK method, this technique uses an eighth-orderaccurate nearly analytic discrete （NAD） operator to discretize high-order spatial differentialoperators and employs a second-order SPRK method to discretize temporal derivatives.The stability criteria and numerical dispersion relations of the eighth-order NSPRK methodare given by a semi-analytical method and are tested by numerical experiments. We alsoshow the differences of the numerical dispersions between the eighth-order NSPRK methodand conventional numerical methods such as the fourth-order NSPRK method, the eighth-order Lax-Wendroff correction （LWC） method and the eighth-order staggered-grid （SG）method. The result shows that the ability of the eighth-order NSPRK method to suppress thenumerical dispersion is obviously superior to that of the conventional numerical methods. Inthe same computational environment, to eliminate visible numerical dispersions, the eighth-order NSPRK is approximately 2.5 times faster than the fourth-order NSPRK and 3.4 timesfaster than the fourth-order SPRK, and the memory requirement is only approximately47.17% of the fourth-order NSPRK method and 49.41% of the fourth-order SPRK method,which indicates the highest computational efficiency. Modeling examples for the two-layermodels such as the heterogeneous and Marmousi models show that the wavefields generatedby the eighth-order NSPRK method are very clear with no visible numerical dispersion.These numerical experiments illustrate that the eighth-order NSPRK method can effectivelysuppress numerical dispersion when coarse grids are adopted. Therefore, this methodcan greatly decrease computer memory requirement and accelerate the forward modelingproductivity. In general, the eighth-order NSPRK method has tremendous potential value forseismic exploration and seismology research.
基金Project (50904027) supported by the National Natural Science Foundation of ChinaProject (2013BAB03B05) supported by the National Key Technology R&D Program of China+1 种基金Project (20133BCB23018) supported by the Foundation for Young Scientist(Jinggang Star)of Jiangxi Province,ChinaProject (2012ZBAB206002) supported by the Natural Science Foundation of Jiangxi Province,China
文摘According to the innate characteristic of four types of furnace, the copper flash continuous smelting （CFCS） furnace can be considered a synthetic reactor of two relatively independent processes： flash matte smelting process （FMSP） and copper continuous converting process （CCCP）. Then, the CFCS thermodynamic model was proposed by establishing the multi-phase equilibrium model of FMSP and the local-equilibrium model of CCCP, respectively, and by combining them through the smelting intermediates. Subsequently, the influences of the furnace structures were investigated using the model on the formation of blister copper, the Fe3O4 behavior, the copper loss in slag and the copper recovery rate. The results show that the type D furnace, with double flues and a slag partition wall, is an ideal CFCS reactor compared with the other three types furnaces. For CFCS, it is effective to design a partition wall in the furnace to make FMSP and CCCP perform in two relatively independent zones, respectively, and to make smelting gas and converting gas discharge from respective flues.
文摘This work is focused on the development of computational algorithms to create a simulator for solving the heat transfer during the continuous casting process of steel. The temperatures and the solid shell thickness profiles were calculated and displayed on the screen for a billet through a defined continuous casting plant （CCP）. The algorithms developed to calculate billet temperatures, involve the solutions of the corresponding equations for the heat removal conditions such as radiation, forced convection, and conduction according to the billet posi- tion through the CCP. This is done by a simultaneous comparison with the kinematics model previously developed. A finite difference method known as Cmnk-Nicholson is applied to solve the two-dimensional computational array （2D model）. Enthalpy （HIJ） and temperature （TIL） in every node are updated at each step time. The routines to display the results have been developed using a graphical user interface （GUI） in the programming language C＋＋. Finally, the results obtained are compared with those of industrial trials for the surface temperature of three steel casters with different plant configurations in different casting conditions.
基金supported by the Hebei Provincial Natural Science Foundation of China(No.E2007000591).
文摘The aim of this study was to simulate the solidification process of beam blank continuous casting, and then find the reasons for the typical defects of the beam blank. A two-dimensional transient coupled finite element model has been developed to compute the temperature and stress profile in beam blank continuous casting. The enthalpy method was used in the heat conduction equation. The thermo-mechanical property in the mushy zone was taken into consideration in this calculation. It is shown that at the mold exit the thickness of the shell had its maximum value at the flange tip and its minimum value at the fillet. The temperature had a great fluctuation on the surface of the beam blank in the secondary cooling zone. At the unbending point, the surface temperature of the web was in the brittleness temperature range under the present condition. To ensure the quality, it is necessary to weaken the intensity of secondary cooling. At the mold exit the equivalent stress and strain have higher values at the flange tip and at the web. From the spray 1 to the unbending point, the maximum values of stress and strain gradually moved to the internal section of the flange tip and the web. However, whenever, there were bigger stress and strain values near the flange tip and the web than in the other parts, it must be very easy to generate cracks at those positions. Now, online verification of this simulation has been developed, which has proved to be very useful and efficient to instruct the practical production of beam blank continuous casting.
基金financially supported by the National High Technology Research and Development Program of China (No. 2013AA030706 and No. 2009AA03Z532)the Fundamental Research Funds for the Central Universities of China (No. FRF-TP-12-146A)
文摘The steady-state temperature field of horizontal core-filling continuous casting （HCFC） for producing copper cladding aluminum rods was simulated by finite element method to investigate the effects of key processing parameters on the positions of solid-liquid interfaces （SLIs） of copper and aluminum. It is found that mandrel tube length and mean withdrawing speed have significant effects on the SLI positions of both copper and aluminum. Aluminum casting temperature （TAI） （1003-1123 K） and secondary cooling water flux （600-900 L.h-1） have little effect on the SLI of copper but cause the SLI of aluminum to move 2-4 mm. When TA1 is in a range of 1043-1123 K, the liquid aluminum can fill continuously into the pre-solidified copper tube. Based on the numerical simulation, reasonable processing parameters were determined.
基金The authors are grateful to IPSCO Inc. for supporting this research program and for providing surface temperature measure-ments, and measurements of the spray nozzle flow characteris-tics. Without the nozzle flow pattern measurements made by Dr. L.K.Chia
文摘A three-dimensional heat transfer model for continuous steel slab casting has been developed with realistic spray cooling patterns and a coupled microsegregation solidification model that calculates the solidification path for multi-component steels. Temperature and composition dependent properties are implemented in a database for 15 chemical species. Considerable effort is made to accurately model the spray cooling heat transfer. Each spray nozzle position and distribution is considered, including variations of the spray patterns with flow rate, and spray overlap. Nozzle type, layout, nozzle-to-slab distance, and spray span and flux are variable. Natural convection, thermal radiation and contact cooling of individual rolls are computed. The present model provides more comprehensive information and realistic slab surface temperatures than results from a model using the 'averaged' treatment of boundary conditions. Cooling operating conditions and parameters of individual spray nozzles can be analyzed to optimize nozzle spray distribution, improve product quality, and troubleshoot issues such as nozzle clogging that may arise during production. One spray cooling correlation is used for the entire machine, achieving as good or better agreement with surface temperature measurements than was found previously for the model using an 'averaged' treatment of boundary conditions and using three machine-segment-dependent correlations.
基金supported by the National Key Research and Development Program of China(2016YFC0800208)the National Natural Science Foundation of China(Nos.51278216,51478201,51308241,and 51608316)
文摘Piled embankments are widely used in highway and railway engineering due to their economy and efficiency inovercoming several issues encountered in constructing embankments over weak soils. Soil arching, caused by the pile-subsoilrelative displacement （△s）, plays an important role in reducing the embankment load falling on weak soil, however, the funda-mental characteristics （e.g., formation and features） of soil arching remain poorly understood. In this study, a series of discreteelement method （DEM） modellings are performed to study the formation and features of soil arching with the variation of As inpiled embankments with or without geosynthetic reinforcement. Firstly, calibration for the modelling parameters is carried out bycomparing the DEM results with the experimental data obtained from the existing literature. Secondly, the analysis of the macro-and micro-behaviours is performed in detail. Finally, a parametric study is conducted in an effort to identify the influences of threekey factors on soil arching： the friction coefficient of the embankment fill （f）, the embankment height （h）, and the pile clear spacing（s-a）. Numerical results indicate that △s is a key factor governing the formation and features of soil arching in embankments. Tobe specific, soil arching gradually evolves from two inclined shear planes at a small △s to a hemispherical arch at a relatively largeAs. Then, with a continuous increase in △s, the soil arching height gradually increases and finally approaches a constant value of0.8（s-a） （i.e., the maximum soil arching height）. For a given case, the higher the soil arching height, the greater the degree of soilarching effect. The parametric study shows that the friction coefficient of the embankment fill has a negligible influence on theformation and features of soil arching. However, embankment height is a key factor governing the formation and features of soilarching. In addition, pile clear spacing has a significant effect on the formation of soil arching, but not on its features.
基金supported by the National Natural Science Foundation of China (No.50874007)
文摘A new method called mixed Lagrangian and Eulerian （MiLE） method was used to simulate the continuous casting process in a mold of free-cutting steel 38MnVS.The simulation results are basically in agreement with experimental data in the literature,achieving the three-dimensional visualization of temperature distribution,melt flow,shell thickness,and stress distribution of blooms in a mold.It is shown that the flow velocity of steel melt becomes smaller gradually as the casting proceeds.When the flow reaches a certain depth,two types of flow patterns can be observed in the upper zone of the mold.The first flow pattern is to flow downwards,and the second one is to flow upwards to the meniscus.The corner temperature is higher,and the thickness is thinner than those in the mid-face.The effective stress in the corner area is much bigger than that in the mid-face,indicating that the corner area is the dangerous zone of cracking.
基金the National Natural Science Foundation of China （No. 599995442）.
文摘Mold is the heart of the continuous casting machine. Heat transfer and solidification in a water- cooled mold are the most important factors during the continuous casting of steel. For studying the temperature distribution of a mold wall, a simulated apparatus of mold was designed and experiments were performed by it. The measured results indicated that the mold wall temperature approaches the temperature of cooling-water. An equivalent thermal-conductivity coefficient was proposed and deduced on the basis of the conclusion of the experiments. This coefficient was applied to solve the heat transfer between the melt and cooling water, and to characterize the heat transfer capacity of the mold. By this equivalent thermal-conductivity coefficient, it is very easy and convenient to numerically simulate the solidification process of continuous casting. And the calculation results are in agreement with the experiments. The effects of custing speed and water flow rate on the mold temperature field were also discussed.
文摘The magnetic,heat transfer and flow phenomenon occurring in the continuous casting process under the mold electromagnetic stirring was further analyzed by solving the 3-D electromagnetic field mathematical model and flow solidification model with finite element method and finite volume method,respectively.The results indicate that the solidified shell thickness located in the effective stirring region fluctuates because of the unsteady scouring under the mold electromagnetic stirring.The maximum rotational velocity is a key parameter to the solidification of the billet when controlling the stirring intensity.When the rotational velocity reaches 0.32 m/s,the mush zone enlarges significantly and the solidification rate is further accelerated.The number of vortexes in the lower recirculation zone is not only two and depends on the stirring parameters.Besides,the secondary flow is closely associated with the solidification.Compared with the results of the model ignoring the influence of solidification on the flow of molten steel,the flow pattern within the lower recirculation region changes dramatically,and thus a coupling analysis of the flow,heat transfer,and solidification is essential when simulating the electromagnetic continuous casting process.
基金Item Sponsored by National Natural Science Foundation of China (59734080)
文摘Heat transfer and deformation of initial solidification shell in soft contact continuous casting mold under high frequency electromagnetic field were analyzed using numerical simulation method; the relative electromagnetic parameters were obtained from the previous studies. Owing to the induction heating of a high frequency electromagnetic field （20 kHz）, the thickness of initial solidification shell decreases, and the temperature of strand surface and slit copper mold increases when compared with the case without the electromagnetic filed. The viscosity of flux de- creases because of the induction heating of the high frequency electromagnetic field, and the dimension of the flux channel increases with electromagnetic pressure; thus, the deformation behavior of initial solidification shell was different before and after the action of high frequency electromagnetic field. Furthermore, the abatement mechanism of oscillation marks under high frequency electromagnetic field was explained.
基金The authors gratefully acknowledge the financial support of the National Key Research and Development Plan(Nos.2017YFB0304100 and 2016YFB0300105)National Natural Science Foundation of China(Nos.51674072,51704151,and 51804067)Fundamental Research Funds for the Central Universities(Nos.N182504014,N170708020,and N172503013).
文摘A 3D/2D hybrid multi-physical-field mathematical model,which takes into consideration the thermosolutal buoyance,was developed to predict the macrosegregation of gear steel 20CrMnTi continuously cast by a curved billet caster with size of 160 mm×160 mm,and investigated the effect of final electromagnetic stirring(F-EMS)on the fluid flow,heat transfer and solute distribution in the liquid core of continuously cast steel.The results show that the application of F-EMS eliminates the effect of thermosolutal buoyancy on the asymmetric distribution of carbon concentration in the cross section of billet and accelerates the final solidification of resident molten steel in the liquid core of strand,but promotes the negative carbon segregation near the billet center.When the gear steel 20CrMnTi is cast at the temperature of 1803 K and speed of 1.7 m/min,the solidification end advances forward from 9.84 to 9.72 m,and center carbon segregation ratio of billet decreases from 1.24 to 1.17 with the increase in current density of F-EMS from 0 to 350 A.
基金Supported by the National Natural Science Foundation of China （No.20776121） and the Specialized Research Fund for the Doctoral Program of Higher Education （No.20050530001）, and the Scientific Research Fund of Hunan Provincial Education Department （No.07C765）.
文摘The discrete particle method was used to simulate the distribution of gas holdup in a gas-liquid standard Rushton stirred tank. The gas phase was treated as a large number of bubbles and their trajectories were tracked with the results of motion equations. The two-way approach was performed to couple the interphase momentum exchange. The turbulent dispersion of bubbles with a size distribution was modeled using a stochastic tracking model, and the added mass force was involved to account for the effect of bubble acceleration on the surrounding fluid. The predicted gas holdup distribution showed that this method could give reasonable prediction comparable to the reported experimental data when the effect of turbulence was took into account in modification for drag coefficient.
文摘The unsteady turbulent flow during the continuous casting of steel is important, because it influences critical phenomena that affect steel quality. Unsteady three-dimensional flow in the mold region of the liquid pool during continuous casting of steel slabs has been computed using realistic geometries starting from the submerged inlet nozzle to the mold. The cassette filter function was used to deal with unsteady Navier-Stokes equation, and then the turbulent flow in the thin slab CCM was simulated with the large eddy simulation method combined with the Smagorinsky sub-grid scale model in this paper. And the model was verified by the Particle Image Velocimetry (PIV) experimental results which was got from a relate scientific literature. In this thesis, by means of LES, the flow characteristics in the thin slab CCM were acquired, such as the vortex distribution, the formation of the large eddy coherent structures, development, shedding and fracture process. In the same time, the turbulent asymmetric distribution was revealed even the nozzle in the centre position. Interactions between the two halves cause large velocity fluctuations near the meniscus. And the vortex is located at the low velocity side adjacent to the SEN. Along with the unsteady time development, the unsteady turbulent large vortex structures of the liquid steel in the CCM presented periodic bias flow distribution, and the period is about 20 seconds.
基金Project(20050248007) supported by the Specialized Research Fund for the Doctoral Program of Higher Education of China
文摘In order to continuously simulate multi-pass plate rolling process,a 3-D elastic hollow-roll model was proposed and an auto mesh-refining module with data passing was developed and integrated with FE software,Marc.The hollow-roll model has equivalent stiffness of bending resistance and deformation to the real solid and much less meshes,so the computational time is greatly reduced.Based on these,the factors influencing plate profile,such as the roll-bending force,initial crown,thermal crown and heat transfer during rolling and inter-pass cooling can be taken into account in the simulation.The auto mesh-refining module with data passing can automatically refine and re-number elements and transfer the nodal and elemental results to the new meshes.Furthermore,the 3-D modeling routine is parametrically developed and can be run independently of Marc pre-processing program.A seven-pass industrial hot rolling process was continuously simulated to validate the accuracy of model.By comparison of the calculated results with the industrial measured data,the rolling force,temperature and plate profile are in good accordance with the measured ones.
基金Project (2007CB714006) supported by the National Basic Research Program of China
文摘The discrete element method (DEM) was used to simulate the flow characteristic and strength characteristic of the conditioned sands in the earth pressure balance (EPB) tunneling. In the laboratory the conditioned sands were reproduced and the slump test and the direct shear test of the conditioned sands were implemented. A DEM equivalent model that can simulate the macro mechanical characteristic of the conditioned sands was proposed,and the corresponding numerical models of the slump test and the shear test were established. By selecting proper DEM model parameters,the errors of the slump values between the simulation results and the test results are in the range of 10.3%-14.3%,and the error of the curves between the shear displacement and the shear stress calculated with the DEM simulation is 4.68%-16.5% compared with that of the laboratory direct shear test. This illustrates that the proposed DEM equivalent model can approximately simulate the mechanical characteristics of the conditioned sands,which provides the basis for further simulation of the interaction between the conditioned soil and the chamber pressure system of the EPB machine.
