We develop and implement a Stochastic Discrete Event Simulation(SDES)algorithm to model the housing re-covery trajectory after an extreme event.The algorithm models discrete events and their underlying uncertainties i...We develop and implement a Stochastic Discrete Event Simulation(SDES)algorithm to model the housing re-covery trajectory after an extreme event.The algorithm models discrete events and their underlying uncertainties in each construction phase.Specifically,the algorithm is developed for the Government Assisted Owner Driven(GAOD)reconstruction system to simulate long-term recovery trajectory.SDES,as a flexible modeling approach,can simulate any housing recovery scenario that follows phased reconstruction.The 2015 M 7.8 Gorkha earthquake sequence in Nepal is considered the extreme event,with 796,245 buildings requiring reconstruction.We present some recovery trajectories from severely hit,crisis hit,and earthquake hit parishes,comparing them with the actual reconstruction progress.We also assess quality and improvement of reconstructed buildings using seismic fragility functions,compared to pre-earthquake constructions.Housing recovery uncertainties are dissected in relation to reconstruction pace.We conclude that the vast majority of the reconstructed buildings followed the Build Back Better(BBB)approach and missed the opportunity to pursue the Build Back Resilient(BBR)approach due to multifaceted challenges ranging from unclear policies to economic constraints.We critically assess the GAOD vs Owner Driven(OD)recovery framework and conclude that insurance-supported and technically assisted OD approach could be the most suitable model for post extreme event housing recovery.展开更多
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.展开更多
Storm-enhanced density(SED)and the tongue of ionization(TOI)are key ionospheric storm-time structures whose rapid evolution and fine-scale variability remain challenging to capture with conventional empirical high-lat...Storm-enhanced density(SED)and the tongue of ionization(TOI)are key ionospheric storm-time structures whose rapid evolution and fine-scale variability remain challenging to capture with conventional empirical high-latitude drivers.In this study,we examine the May 10–11,2024,superstorm using the Thermosphere–Ionosphere–Electrodynamics General Circulation Model(TIEGCM)with observation-constrained high-latitude forcing.Auroral precipitation parameters(energy flux and mean energy)are assimilated from a Defense Meteorological Satellite Program(DMSP)Special Sensor Ultraviolet Spectrographic Imager(SSUSI)using a multi-resolution Gaussian process(Lattice Kriging)approach,whereas high-latitude convection potentials are derived by assimilating Super Dual Auroral Radar Network(SuperDARN)observations with the Thomas and Shepherd(2018)model(TS18).For comparison,an additional simulation is performed using empirical models for both convection and auroral forcing.The results show that during the main phase of the May 10 storm,the data-driven simulation provides a more realistic depiction of the SED source region than does the empirical model run by capturing its rapid intensification more clearly and reproducing its spatial location and structural features with higher fidelity.These improvements lead to a more accurate representation of its poleward extension into the polar cap that develops into the TOI.Above the ionospheric F2 peak over the SED source region,SuperDARN-constrained potentials generate stronger and more localized E×B drifts that dominate plasma uplift and drive its transport into the polar cap,although neutral winds and downward ambipolar diffusion partially offset these effects.Below the F2 peak,neutral winds and photochemical processes play a major role in shaping the spatial extent and intensity of the SED and TOI.These results highlight the role of observation-constrained high-latitude drivers in representing ionosphere–thermosphere responses during extreme storms and suggest their relevance for improving physical interpretation and model performance.展开更多
Under the background of achieving energy efficiency and carbon reduction,the use of Continuous Casting-Direct Hot Charging Rolling(CC-DHCR)technology was explored at the Hot Rolling Plant of Baosteel,China.As temperat...Under the background of achieving energy efficiency and carbon reduction,the use of Continuous Casting-Direct Hot Charging Rolling(CC-DHCR)technology was explored at the Hot Rolling Plant of Baosteel,China.As temperature variation and temperature uniformity directly affect the rolling quality of the billet,numerical simulations of the billet temperature profile changes in the CC-DHCR were conducted at the start of the industrial trial,and the billet temperature distribution and cross-section temperature difference during its transportation and throughout the heating process were analyzed.According to numerical simulation calculations,the average temperature of billet discharged from the heating furnace was 1150℃,which after subsequent controlled cooling met the final rolling temperature requirements of(880±30)℃for this kind of steel.The maximum temperature difference of the billet discharged from the furnace was within 35 K,which meets the billet heating uniformity requirements of the rolling process.The simulation results were compared with experimental results,and the rationality of the simulation was validated.In addition,the industrial trial billet was rolled,and the rolling quality was good during tracking.In this CC-DHCR industrial trial,the fuel consumption was 28.6 kgce/t(kilogram coal equivalent per ton),a reduction of 33.7%compared with the current traditional process,and CO 2 emissions were reduced by 38.09 kg/t.展开更多
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.展开更多
文摘We develop and implement a Stochastic Discrete Event Simulation(SDES)algorithm to model the housing re-covery trajectory after an extreme event.The algorithm models discrete events and their underlying uncertainties in each construction phase.Specifically,the algorithm is developed for the Government Assisted Owner Driven(GAOD)reconstruction system to simulate long-term recovery trajectory.SDES,as a flexible modeling approach,can simulate any housing recovery scenario that follows phased reconstruction.The 2015 M 7.8 Gorkha earthquake sequence in Nepal is considered the extreme event,with 796,245 buildings requiring reconstruction.We present some recovery trajectories from severely hit,crisis hit,and earthquake hit parishes,comparing them with the actual reconstruction progress.We also assess quality and improvement of reconstructed buildings using seismic fragility functions,compared to pre-earthquake constructions.Housing recovery uncertainties are dissected in relation to reconstruction pace.We conclude that the vast majority of the reconstructed buildings followed the Build Back Better(BBB)approach and missed the opportunity to pursue the Build Back Resilient(BBR)approach due to multifaceted challenges ranging from unclear policies to economic constraints.We critically assess the GAOD vs Owner Driven(OD)recovery framework and conclude that insurance-supported and technically assisted OD approach could be the most suitable model for post extreme event housing recovery.
基金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.
基金The Shandong Provincial Natural Science Foundation(Grant No.ZR2022JQ18)supported this worksupported by the National Natural Science Foundation of China(NNFSC)Youth Program(Grant No.42304168)+1 种基金supported by the National Key R&D Program of China(Grant No.2022YFF0504400)the NNSFC(Grant Nos.42188101 and 42174210)。
文摘Storm-enhanced density(SED)and the tongue of ionization(TOI)are key ionospheric storm-time structures whose rapid evolution and fine-scale variability remain challenging to capture with conventional empirical high-latitude drivers.In this study,we examine the May 10–11,2024,superstorm using the Thermosphere–Ionosphere–Electrodynamics General Circulation Model(TIEGCM)with observation-constrained high-latitude forcing.Auroral precipitation parameters(energy flux and mean energy)are assimilated from a Defense Meteorological Satellite Program(DMSP)Special Sensor Ultraviolet Spectrographic Imager(SSUSI)using a multi-resolution Gaussian process(Lattice Kriging)approach,whereas high-latitude convection potentials are derived by assimilating Super Dual Auroral Radar Network(SuperDARN)observations with the Thomas and Shepherd(2018)model(TS18).For comparison,an additional simulation is performed using empirical models for both convection and auroral forcing.The results show that during the main phase of the May 10 storm,the data-driven simulation provides a more realistic depiction of the SED source region than does the empirical model run by capturing its rapid intensification more clearly and reproducing its spatial location and structural features with higher fidelity.These improvements lead to a more accurate representation of its poleward extension into the polar cap that develops into the TOI.Above the ionospheric F2 peak over the SED source region,SuperDARN-constrained potentials generate stronger and more localized E×B drifts that dominate plasma uplift and drive its transport into the polar cap,although neutral winds and downward ambipolar diffusion partially offset these effects.Below the F2 peak,neutral winds and photochemical processes play a major role in shaping the spatial extent and intensity of the SED and TOI.These results highlight the role of observation-constrained high-latitude drivers in representing ionosphere–thermosphere responses during extreme storms and suggest their relevance for improving physical interpretation and model performance.
文摘Under the background of achieving energy efficiency and carbon reduction,the use of Continuous Casting-Direct Hot Charging Rolling(CC-DHCR)technology was explored at the Hot Rolling Plant of Baosteel,China.As temperature variation and temperature uniformity directly affect the rolling quality of the billet,numerical simulations of the billet temperature profile changes in the CC-DHCR were conducted at the start of the industrial trial,and the billet temperature distribution and cross-section temperature difference during its transportation and throughout the heating process were analyzed.According to numerical simulation calculations,the average temperature of billet discharged from the heating furnace was 1150℃,which after subsequent controlled cooling met the final rolling temperature requirements of(880±30)℃for this kind of steel.The maximum temperature difference of the billet discharged from the furnace was within 35 K,which meets the billet heating uniformity requirements of the rolling process.The simulation results were compared with experimental results,and the rationality of the simulation was validated.In addition,the industrial trial billet was rolled,and the rolling quality was good during tracking.In this CC-DHCR industrial trial,the fuel consumption was 28.6 kgce/t(kilogram coal equivalent per ton),a reduction of 33.7%compared with the current traditional process,and CO 2 emissions were reduced by 38.09 kg/t.
文摘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.
