A 3D mathematical model was proposed to investigate the molten steel–slag–air multiphase flow in a two-strand slab continuous casting(CC)tundish during ladle change.The study focused on the exposure of the molten st...A 3D mathematical model was proposed to investigate the molten steel–slag–air multiphase flow in a two-strand slab continuous casting(CC)tundish during ladle change.The study focused on the exposure of the molten steel and the subsequent reoxidation occurrence.The exposure of the molten steel was calculated using the coupled realizable k–εmodel and volume of fluid(VOF)model.The diffusion of dissolved oxygen was determined by solving the user-defined scalar(UDS)equation.Moreover,the user-defined function(UDF)was used to describe the source term in the UDS equation and determine the oxidation rate and oxidation position.The effect of the refilling speed on the molten steel exposure and dissolved oxygen content was also discussed.Increasing the refilling speed during ladle change reduced the refilling time and the exposure duration of the molten steel.However,the elevated refilling speed enlarged the slag eyes and increased the average dissolved oxygen content within the tundish,thereby exacerbating the reoxidation phenomenon.In addition,the time required for the molten steel with a high dissolved oxygen content to exit the tundish varied with the refilling speed.When the inlet speed was 3.0 m·s^(-1)during ladle change,the molten steel with a high dissolved oxygen content exited the outlet in a short period,reaching a maximum dissolved oxygen content of 0.000525wt%.Conversely,when the inlet speed was 1.8 m·s^(-1),the maximum dissolved oxygen content was 0.000382wt%.The refilling speed during the ladle change process must be appropriately decreased to minimize reoxidation effects and enhance the steel product quality.展开更多
Electromagnetic stirring is one of the widely applied techniques to modify the quality of casting billets.Different from conventional rotate stirring,the helical stirring is more professional in assisting multi-dimens...Electromagnetic stirring is one of the widely applied techniques to modify the quality of casting billets.Different from conventional rotate stirring,the helical stirring is more professional in assisting multi-dimensional flow of molten metal and eliminating solidification defects.In this study,the single-winding helical stirring(SWHS)was introduced,offering advantages such as smaller volume and lower electromagnetic shielding compared to traditional helical stirring methods.Following a comprehensive numerical simulation,the stirring parameters of SWHS were adjusted to yoke inclination angle of 43°and frequency of 12 Hz.The higher electromagnetic force and flow velocity in drawing direction,as well as the lower temperature gradient induced by the SWHS,are positive factors for homogeneous solidification of billet.The experimental results on Al-8%Si alloy and 0.4%C-1.1%Mn steel demonstrate that compared to rotate stirring,the SWHS process can induce better billet quality and is more effective in accelerating the equiaxed expansion and reducing element segregation.The SWHS process can enhance the equiaxed ratio of the billet by 58.3%and reduce segregation degree of carbon element by 10.97%.Consequently,SWHS holds great promise as a potential approach for improving the quality of continuous casting billets.展开更多
Accurate measurements of the radon exhalation rate help identify and evaluate radon risk regions in the environment.Among these measurement methods,the closed-loop method is frequently used.However,traditional experim...Accurate measurements of the radon exhalation rate help identify and evaluate radon risk regions in the environment.Among these measurement methods,the closed-loop method is frequently used.However,traditional experiments are insufficient or cannot analyze the radon migration and exhalation patterns at the gas–solid interface in the accumulation chamber.The CFD-based technique was applied to predict the radon concentration distribution in a limited space,allowing radon accumulation and exhalation inside the chamber intuitively and visually.In this study,three radon exhalation rates were defined,and two structural ventilation tubes were designed for the chamber.The consistency of the simulated results with the variation in the radon exhalation rate in a previous experiment or analytical solution was verified.The effects of the vent tube structure and flow rate on the radon uniformity in the chamber;permeability,insertion depth,and flow rate on the radon exhalation rate and the effective diffusion coefficient on back-diffusion were investigated.Based on the results,increasing the inser-tion depth from 1 to 5 cm decreased the effective decay constant by 19.55%,whereas the curve-fitted radon exhalation rate decreased(lower than the initial value)as the deviation from the initial value increased by approximately 7%.Increasing the effective diffusion coefficient from 2.77×10^(-7) to 7.77×10^(-6) m^(2) s^(-1) made the deviation expand from 2.14 to 15.96%.The conclusion is that an increased insertion depth helps reduce leakage in the chamber,subject to notable back-diffusion,and that the closed-loop method is reasonably used for porous media with a low effective diffusion coefficient in view of the back-diffusion effect.The CFD-based simulation is expected to provide guidance for the optimization of the radon exhalation rate measurement method and,thus,the accurate measurement of the radon exhalation rate.展开更多
Wind-sand flow generates erosion and deposition around obstacles such as bridges and roadbeds, resulting in sand damage and endangering railway systems in sandy regions. Previous studies have mainly focused on the flo...Wind-sand flow generates erosion and deposition around obstacles such as bridges and roadbeds, resulting in sand damage and endangering railway systems in sandy regions. Previous studies have mainly focused on the flow field around roadbeds, overlooking detailed examinations of sand particle erosion and deposition patterns near bridges and roadbeds. This study employs numerical simulations to analyze the influence of varying heights and wind speeds on sand deposition and erosion characteristics at different locations: the bridge-road transition section(side piers), middle piers, and roadbeds. The results show that the side piers, experience greater accumulation than the middle piers. Similarly, the leeward side of the roadbed witnesses more deposition compared to the windward side. Another finding reveals a reduced sand deposition length as the vertical profile, in alignment with the wind direction, moves further from the bridge abutments at the same clearance height. As wind speeds rise, there’s a decline in sand deposition and a marked increase in erosion around the side piers, middle piers and roadbeds. In conclusion, a bridge clearance that’s too low can cause intense sand damage near the side piers, while an extremely high roadbed may lead to extensive surface sand deposition. Hence, railway bridges in areas prone to sandy winds should strike a balance in clearance height. This research provides valuable guidelines for determining the most suitable bridge and roadbed heights in regions affected by wind and sand.展开更多
For the safety of railroad operations,sand barriers are utilized to mitigate wind-sand disaster effects.These disasters,characterized by multi-directional wind patterns,result in diverse angles among the barriers.In t...For the safety of railroad operations,sand barriers are utilized to mitigate wind-sand disaster effects.These disasters,characterized by multi-directional wind patterns,result in diverse angles among the barriers.In this study,using numerical simulations,we examined the behavior of High Density Polyethylene(HDPE)sheet sand barriers under different wind angles,focusing on flow field distribution,windproof efficiency,and sedimentation erosion dynamics.This study discovered that at a steady wind speed,airflow velocity varies as the angle between the airflow and the HDPE barrier changes.Specifically,a 90°angle results in the widest low-speed airflow area on the barrier’s downwind side.If the airflow is not perpendicular to the barrier,it prompts a lateral airflow movement which decreases as the angle expands.The windproof efficiency correlates directly with this angle but inversely with the wind’s speed.Notably,with a wind angle of 90°,wind speed drops by 81%.The minimum wind speed is found at 5.1H(the sand barrier height)on the barrier’s downwind side.As the angle grows,the barrier’s windproof efficiency improves,extending its protective reach.Sedimentation is most prominent on the barrier’s downwind side,as the wind angle shifts from 30°to 90°,the sand sedimentation area on the barrier’s downwind side enlarges by 14.8H.As the angle grows,sedimentation intensifies,eventually overtakes the forward erosion and enlarges the sedimentation area.展开更多
To study the anti-explosion protection effect of polyurea coating on reinforced concrete box girder,two segmental girder specimens were made at a scale of 1:3,numbered as G(without polyurea coating)and PCG(with polyur...To study the anti-explosion protection effect of polyurea coating on reinforced concrete box girder,two segmental girder specimens were made at a scale of 1:3,numbered as G(without polyurea coating)and PCG(with polyurea coating).The failure characteristics and dynamic responses of the specimens were compared through conducting explosion tests.The reliability of the numerical simulation using LS-DYNA software was verified by the test results.The effects of different scaled distances,reinforcement ratios,concrete strengths,coating thicknesses and ranges of polyurea were studied.The results show that the polyurea coating can effectively enhance the anti-explosion performance of the girder.The top plate of middle chamber in specimen G forms an elliptical penetrating hole,while that in specimen PCG only shows a very slight local dent.The peak vertical displacement and residual displacement of PCG decrease by 74.8% and 73.7%,respectively,compared with those of specimen G.For the TNT explosion with small equivalent,the polyurea coating has a more significant protective effect on reducing the size of fracture.With the increase of TNT equivalent,the protective effect of polyurea on reducing girder displacement becomes more significant.The optimal reinforcement ratio,concrete strength,thickness and range of polyurea coating were also drawn.展开更多
Flash boiling atomization(FBA)is a promising approach for enhancing spray atomization,which can generate a fine and more evenly distributed spray by increasing the fuel injection temperature or reducing the ambient pr...Flash boiling atomization(FBA)is a promising approach for enhancing spray atomization,which can generate a fine and more evenly distributed spray by increasing the fuel injection temperature or reducing the ambient pressure.However,when the outlet speed of the nozzle exceeds 400 m/s,investigating high-speed flash boiling atomization(HFBA)becomes quite challenging.This difficulty arises fromthe involvement ofmany complex physical processes and the requirement for a very fine mesh in numerical simulations.In this study,an HFBA model for gasoline direct injection(GDI)is established.This model incorporates primary and secondary atomization,as well as vaporization and boilingmodels,to describe the development process of the flash boiling spray.Compared to lowspeed FBA,these physical processes significantly impact HFBA.In this model,the Eulerian description is utilized for modeling the gas,and the Lagrangian description is applied to model the droplets,which effectively captures the movement of the droplets and avoids excessive mesh in the Eulerian coordinates.Under various conditions,numerical solutions of the Sauter mean diameter(SMD)for GDI show good agreement with experimental data,validating the proposed model’s performance.Simulations based on this HFBA model investigate the influences of fuel injection temperature and ambient pressure on the atomization process.Numerical analyses of the velocity field,temperature field,vapor mass fraction distribution,particle size distribution,and spray penetration length under different superheat degrees reveal that high injection temperature or low ambient pressure significantly affects the formation of small and dispersed droplet distribution.This effect is conducive to the refinement of spray particles and enhances atomization.展开更多
Wind and sand hazards are serious in the Milan Gobi area of the Xinjiang section of the Korla Railway. In order to ensure the safe operation of railroads, there is a need for wind and sand protection in heavily sandy ...Wind and sand hazards are serious in the Milan Gobi area of the Xinjiang section of the Korla Railway. In order to ensure the safe operation of railroads, there is a need for wind and sand protection in heavily sandy areas. The wind and sand flow in the region is notably bi-directional. To shield railroads from sand, a unique sand fence made of folded linear high-density polyethylene(HDPE) is used, aligning with the principle that the dominant wind direction is perpendicular to the fence. This study employed field observations and numerical simulations to investigate the effectiveness of these HDPE sand fences in altering flow field distribution and offering protection. It also explored how these fences affect the deposition and erosion of sand particles. Findings revealed a significant reduction in wind speed near the fence corner;the minimum horizontal wind speed on the leeward side of the first sand fence(LSF) decreased dramatically from 3 m/s to 0.64 m/s. The vortex area on the LSF markedly impacted horizontal wind speeds. Within the LSF, sand deposition was a primary occurrence. As wind speeds increased, the deposition zone shrank, whereas the positive erosion zone expanded. Close to the folded corners of the HDPE sand fence, there was a notable shift from the positive erosion zone to a deposition zone. Field tests and numerical simulations confirmed the high windproof efficiency(WE) and sand resistance efficiency(SE) in the HDPE sand fence. Folded linear HDPE sheet sand fence can effectively slow down the incoming flow and reduce the sand content, thus achieving good wind and sand protection. This study provides essential theoretical guidance for the design and improvement of wind and sand protection systems in railroad engineering.展开更多
Magnetohydrodynamic(MHD)induction pumps are contactless pumps able to withstand harsh environments.The rate of fluid flow through the pump directly affects the efficiency and stability of the device.To explore the inf...Magnetohydrodynamic(MHD)induction pumps are contactless pumps able to withstand harsh environments.The rate of fluid flow through the pump directly affects the efficiency and stability of the device.To explore the influence of induction pump settings on the related delivery speed,in this study,a numerical model for coupled electromagnetic and flow field effects is introduced and used to simulate liquid metal lithium flow in the induction pump.The effects of current intensity,frequency,coil turns and coil winding size on the velocity of the working fluid are analyzed.It is shown that the first three parameters have a significant impact,while changes in the coil turns have a negligible influence.The maximum increase in working fluid velocity within the pump for the parameter combination investigated in this paper is approximately 618%.As the frequency is increased from 20 to 60 Hz,the maximum increase in the mean flow rate of the working fluid is approximately 241%.These research findings are intended to support the design and optimization of these devices.展开更多
To analyze the effects of a time-varying viscosity on the penetration length of grouting,in this study cement slur-ries with varying water-cement ratios have been investigated using the Bingham’sfluidflow equation and ...To analyze the effects of a time-varying viscosity on the penetration length of grouting,in this study cement slur-ries with varying water-cement ratios have been investigated using the Bingham’sfluidflow equation and a dis-crete element method.Afluid-solid coupling numerical model has been introduced accordingly,and its accuracy has been validated through comparison of theoretical and numerical solutions.For different fracture forms(a single fracture,a branch fracture,and a fracture network),the influence of the time-varying viscosity on the slurry length range has been investigated,considering the change in the fracture aperture.The results show that under different fracture forms and the same grouting process conditions,the influence of the time-varying viscosity on the seepage length is 0.350 m.展开更多
As an important lightning protection device in substations,lightning rods are susceptible to vibration and potential structural damage under wind loads.In order to understand their vibration mechanism,it is necessary ...As an important lightning protection device in substations,lightning rods are susceptible to vibration and potential structural damage under wind loads.In order to understand their vibration mechanism,it is necessary to conduct flow analysis.In this study,numerical simulations of the flow field around a 330 kV cylindrical lightning rod with different diameters were performed using the SST k-ωmodel.The flow patterns in different segments of the lightning rod at the same reference wind speed(wind speed at a height of 10 m)and the flow patterns in the same segment at different reference wind speeds were investigated.The variations of lift coefficient,drag coefficient,and vorticity distribution were obtained.The results showed that vortex shedding phenomena occurred in all segments of the lightning rod,and the strength of vortex shedding increased with decreasing diameter.The vorticity magnitude and the root mean square magnitudes of the lift coefficient and drag coefficient also increased accordingly.The time history curves of the lift coefficient and drag coefficient on the surface of the lightning rod exhibited sinusoidal patterns with a single dominant frequency.For the same segment,as the wind speed increased in a certain range,the root mean square values of the lift coefficient and drag coefficient decreased,while their dominant frequencies increased.Moreover,there was a proportional relationship between the dominant frequencies of the lift coefficient and drag coefficient.The findings of this study can provide valuable insights for the refined design of lightning rods with similar structures.展开更多
The radiant tube burner was modeled and analyzed by the numerical simulation method to investigate the influence factors and rules of NO_(x) emissions in a W-type radiant tube.These factors,which include air preheatin...The radiant tube burner was modeled and analyzed by the numerical simulation method to investigate the influence factors and rules of NO_(x) emissions in a W-type radiant tube.These factors,which include air preheating temperature,excess air coefficient,and fuel gas composition,were modified to study their effects on NO_(x) emissions under varying working conditions.Simulation results were compared with the theoretical calculation value based on chemical reaction equilibrium theory and the onsite experimental value to verify the simulation accuracy.The results show that NO_(x) emissions rise with increasing air preheating temperatures.NO_(x) production increases to an extreme value and then decreases during the oxygen-poor to oxygen-enriched process with the rise of the excess air coefficient.Enhancing the proportion of coke oven gas in the fuel gas raises the combustion temperature as well as the NO_(x) discharge.Both the thermal efficiency and NO_(x) emissions should be balanced.Therefore,the recommended values based on the simulation results are as follows:the air preheating temperature should not exceed 400℃,the excess air coefficient should be between 1.1 and 1.2,and the volume fraction of the coke oven gas should not exceed 30%.展开更多
In order to study the residual stress distribution law of welded joints of arch ribs of large-span steel pipe concrete arch bridges,numerical simulation of temperature,stress and strain fields based on ABAQUS for weld...In order to study the residual stress distribution law of welded joints of arch ribs of large-span steel pipe concrete arch bridges,numerical simulation of temperature,stress and strain fields based on ABAQUS for welded joints of arch-ribbed steel tubes using 7-,8-and 9-layer welds is carried out and its accuracy is demonstrated.The steel pipe welding temperature changes,residual stress distribution,different processes residual stress changes in the law,the prediction of post-weld residual stress distribution and deformation are studied in this paper.The results show that the temperature field values and test results are more consistent with the accuracy of numerical simulation of welding,the welding process is mainly in the form of heat transfer;Residual high stresses are predominantly distributed in the Fusion zone(FZ)and Heat-affected zone(HAZ),with residual stress levels tending to decrease from the center of the weld along the axial path,the maximum stress appears in the FZ and HAZ junction;The number of welding layers has an effect on the residual stress distribution,the number of welding layers increases,the residual stress tends to decrease,while the FZ and HAZ high stress area range shrinks;Increasing the number of plies will increase the amount of residual distortion.展开更多
Xinqiao Gully is located in the area of the 2008 Wenchuan M_(s)8.0 earthquake in Sichuan province,China.Based on the investigation of the 2023"6-26"Xinqiao Gully debris flow event,this study assessed the eff...Xinqiao Gully is located in the area of the 2008 Wenchuan M_(s)8.0 earthquake in Sichuan province,China.Based on the investigation of the 2023"6-26"Xinqiao Gully debris flow event,this study assessed the effectiveness of the debris flow control project and evaluated the debris flow hazards.Through field investigation and numerical simulation methods,the indicators of flow intensity reduction rate and storage capacity fullness were proposed to quantify the effectiveness of the engineering measures in the debris flow event.The simulation results show that the debris flow control project reduced the flow intensity by41.05%to 64.61%.The storage capacity of the dam decreases gradually from upstream to the mouth of the gully,thus effectively intercepting and controlling the debris flow.By evaluating the debris flow of different recurrence intervals,further measures are recommended for managing debris flow events.展开更多
Cemented paste backfill(CPB)technology is a green mining method used to control underground goaves and tailings ponds.The curing process of CPB in the stope is the product of a thermo-hydro-mechanical-chemical multi-f...Cemented paste backfill(CPB)technology is a green mining method used to control underground goaves and tailings ponds.The curing process of CPB in the stope is the product of a thermo-hydro-mechanical-chemical multi-field performance interaction.At present,research on the multi-field performance of CPB mainly includes indoor similar simulation experiments,in-situ multi-field performance monitoring experiments,multi-field performance coupling model construction of CPB,and numerical simulation of the multi-field performance of CPB.Because it is hard to study the in-situ multi-field performance of CPB in the real stope,most current research on in-situ multi-field performance adopts the numerical simulation method.By simulating the conditions of CPB in the real stope(e.g.,maintenance environment,stope geometry,drainage conditions,and barricade and backfilling rates),the multi-field performance of CPB is further studied.This paper summarizes the mathematical models employed in the numerical simulation and lists the engineering application cases of numerical simulation in the in-situ multi-field performance of CPB.Finally,it proposes that the multi-field performance of CPB needs to strengthen the theoretical study of multi-field performance,form the strength design criterion based on the multi-field performance of CPB,perform a full-range numerical simulation of the multi-field performance of CPB,develop a pre-warning technology for the CPB safety of CPB,develop automatic and wireless sensors for the multi-field performance monitoring of CPB,and realize the application and popularization of CPB monitoring technology.展开更多
Based on numerical simulations,this study highlights the sedimentation and erosion problems around a sand barrier through the relationship between the shear stress of the surface around the sand barrier and the critic...Based on numerical simulations,this study highlights the sedimentation and erosion problems around a sand barrier through the relationship between the shear stress of the surface around the sand barrier and the critical shear stress of sand grains.The numerical simulation results were verified using data measured by the wind tunnel test.The results showed that when the porosity was the same,the size and position of the vortex on the leeward side of the sand barrier were related to the inlet wind speed.As the wind speed increased,the vortex volume increased and the positions of the separation and reattachment points moved toward the leeward side.When the porosity of the sand barrier was 30%,the strength of the acceleration zone above the sand barrier was the highest,and the strength of the acceleration zone was negatively correlated with the porosity.Sand erosion and sedimentation distance were related to wind speed.With an increase in wind speed,the sand grain forward erosion or reverse erosion areas on the leeward side of the sand barrier gradually replaced the sedimentation area.With an increase in porosity,the sand sedimentation distance on the leeward side of the sand barrier gradually shortened,and the sand erosion area gradually disappeared.The sand sedimentation distance on the leeward side of the sand barrier with 30%porosity was the longest.The numerical simulation results were in good agreement with the wind tunnel test results.Based on the sand erosion and sedimentation results of the numerical simulation and wind tunnel test,when the porosity was 30%,the protection effect of the High Density Polyethylene(HDPE)board sand barrier was best.展开更多
The dendrite growth behavior of high-strength steel during slab continuous casting with a traveling-wave magnetic field was studied in this paper. The morphology of the solidification structure and composition distrib...The dendrite growth behavior of high-strength steel during slab continuous casting with a traveling-wave magnetic field was studied in this paper. The morphology of the solidification structure and composition distribution were analyzed. Results showed that the columnar crystals could deflect and break when the traveling-wave magnetic field had low current intensity. With the increase in current intensity, the secondary dendrite arm spacing and solute permeability decreased, and the columnar crystal transformed into an equiaxed crystal. The electromagnetic force caused by the traveling-wave magnetic field changed the temperature gradient and velocity magnitude and promoted the breaking and fusing of dendrites. Dendrite compactness and composition uniformity were arranged in descending order as follows:columnar-toequiaxed transition (high current intensity), columnar crystal zone (low current intensity), columnar-to-equiaxed transition (low current intensity), and equiaxed crystal zone (high current intensity). Verified numerical simulation results combined with the boundary layer theory of solidification front and dendrite breaking–fusing model revealed the dendrite deflection mechanism and growth process. When thermal stress is not considered, and no narrow segment can be found in the dendrite, the velocity magnitude on the solidification front of liquid steel can reach up to 0.041 m/s before the dendrites break.展开更多
Piles in a group experience additional displacements in soil due to pile-to-pile interactions apart from those resulting from the external loading.The effect of these interactions determined assuming soil as an elasti...Piles in a group experience additional displacements in soil due to pile-to-pile interactions apart from those resulting from the external loading.The effect of these interactions determined assuming soil as an elastic and/or viscoelastic material on pile head impedance functions of the pile group is studied by relating the group stiffness to the static stiffness of a single pile.However,the prevailing elastic solutions may misestimate the resulting pile group response due to the lack of consideration for either soil(material)and/or soil-pile interface nonlinearities.It is well established that soil behaves nonlinearly under moderate-to-high loading amplitudes,and besides,the soil-pile interface nonlinearity can exist even at small loading amplitudes.This study addresses the effects of these nonlinearities on the vertical impedance functions of a 3×3-pile group using numerical methods by direct analyses and superposition using pile-to-pile interaction factors.The numerical results are validated using scaled model tests under 1 g conditions.The results highlight the overestimation of pile-to-pile interactions in the pile group when assuming elastic soil conditions.The cases either by direct analyses or superposition approach involving soil and soil-pile interface nonlinearities agree well with the experimental pile group responses under close-to-elastic and nonlinear conditions.展开更多
A large amount of data can partly assure good fitting quality for the trained neural networks.When the quantity of experimental or on-site monitoring data is commonly insufficient and the quality is difficult to contr...A large amount of data can partly assure good fitting quality for the trained neural networks.When the quantity of experimental or on-site monitoring data is commonly insufficient and the quality is difficult to control in engineering practice,numerical simulations can provide a large amount of controlled high quality data.Once the neural networks are trained by such data,they can be used for predicting the properties/responses of the engineering objects instantly,saving the further computing efforts of simulation tools.Correspondingly,a strategy for efficiently transferring the input and output data used and obtained in numerical simulations to neural networks is desirable for engineers and programmers.In this work,we proposed a simple image representation strategy of numerical simulations,where the input and output data are all represented by images.The temporal and spatial information is kept and the data are greatly compressed.In addition,the results are readable for not only computers but also human resources.Some examples are given,indicating the effectiveness of the proposed strategy.展开更多
Thermo-hydro-mechanical-chemical(THMC)interactions are prevalent during CO_(2)geological sequestration(CGS).In this study,a sequential coupling THMC numerical simulation program was constructed,which can be used to ex...Thermo-hydro-mechanical-chemical(THMC)interactions are prevalent during CO_(2)geological sequestration(CGS).In this study,a sequential coupling THMC numerical simulation program was constructed,which can be used to explore the following issues of CGS:fluid and heat flow,solute transport;stresses,displacements and rock failures related to geo-mechanical effects;equilibrium and kinetic chemical reactions;chemical damage to mechanical properties of the rock.Then,the coupling program was applied to the Ordos CGS Project to study the formation response under the multifield interaction caused by CO_(2)injection.The simulation results show that the mechanical process dominates the short CO_(2)injection period.Specifically,the formation’s permeability near the injection well increases by 43%,due to the reduction of effective stress,which significantly promotes the lateral migration of CO_(2).When the injection rate exceeds 0.15 million tons per year,the cohesion of the reservoir rock is not enough to resist the shear force inside the rock and rock failure may occur.During the subsequent long-term sequestration period(200 years),the influence of mineral reactions gradually increases.Due to calcite dissolution,the shear modulus of caprock is predicted to decrease by 7.6%,which will to some extent increase the risk of rock failure.展开更多
基金the National Key R&D Program(No.2023YFB3709900)the National Natural Science Foundation of China(Nos.U22A20171 and 52104343)the High Steel Central(HSC)at North China University of Science and Technology and Yanshan Univ ersity,China。
文摘A 3D mathematical model was proposed to investigate the molten steel–slag–air multiphase flow in a two-strand slab continuous casting(CC)tundish during ladle change.The study focused on the exposure of the molten steel and the subsequent reoxidation occurrence.The exposure of the molten steel was calculated using the coupled realizable k–εmodel and volume of fluid(VOF)model.The diffusion of dissolved oxygen was determined by solving the user-defined scalar(UDS)equation.Moreover,the user-defined function(UDF)was used to describe the source term in the UDS equation and determine the oxidation rate and oxidation position.The effect of the refilling speed on the molten steel exposure and dissolved oxygen content was also discussed.Increasing the refilling speed during ladle change reduced the refilling time and the exposure duration of the molten steel.However,the elevated refilling speed enlarged the slag eyes and increased the average dissolved oxygen content within the tundish,thereby exacerbating the reoxidation phenomenon.In addition,the time required for the molten steel with a high dissolved oxygen content to exit the tundish varied with the refilling speed.When the inlet speed was 3.0 m·s^(-1)during ladle change,the molten steel with a high dissolved oxygen content exited the outlet in a short period,reaching a maximum dissolved oxygen content of 0.000525wt%.Conversely,when the inlet speed was 1.8 m·s^(-1),the maximum dissolved oxygen content was 0.000382wt%.The refilling speed during the ladle change process must be appropriately decreased to minimize reoxidation effects and enhance the steel product quality.
基金financially supported by the National Key R&D Projects(No.2021YFB3702000)the Regional Company Projects in Ansteel Beijing Research Institute(No.2022BJB07GF&No.2022BJB-13GF)。
文摘Electromagnetic stirring is one of the widely applied techniques to modify the quality of casting billets.Different from conventional rotate stirring,the helical stirring is more professional in assisting multi-dimensional flow of molten metal and eliminating solidification defects.In this study,the single-winding helical stirring(SWHS)was introduced,offering advantages such as smaller volume and lower electromagnetic shielding compared to traditional helical stirring methods.Following a comprehensive numerical simulation,the stirring parameters of SWHS were adjusted to yoke inclination angle of 43°and frequency of 12 Hz.The higher electromagnetic force and flow velocity in drawing direction,as well as the lower temperature gradient induced by the SWHS,are positive factors for homogeneous solidification of billet.The experimental results on Al-8%Si alloy and 0.4%C-1.1%Mn steel demonstrate that compared to rotate stirring,the SWHS process can induce better billet quality and is more effective in accelerating the equiaxed expansion and reducing element segregation.The SWHS process can enhance the equiaxed ratio of the billet by 58.3%and reduce segregation degree of carbon element by 10.97%.Consequently,SWHS holds great promise as a potential approach for improving the quality of continuous casting billets.
基金This work was supported by the National Natural Science Foundation of China(No.11575080)the National Natural Science Foundation of Hunan Province,China(No.2022JJ30482)the Hunan Provincial Innovation Foundation for Postgraduates(No.QL20220206).
文摘Accurate measurements of the radon exhalation rate help identify and evaluate radon risk regions in the environment.Among these measurement methods,the closed-loop method is frequently used.However,traditional experiments are insufficient or cannot analyze the radon migration and exhalation patterns at the gas–solid interface in the accumulation chamber.The CFD-based technique was applied to predict the radon concentration distribution in a limited space,allowing radon accumulation and exhalation inside the chamber intuitively and visually.In this study,three radon exhalation rates were defined,and two structural ventilation tubes were designed for the chamber.The consistency of the simulated results with the variation in the radon exhalation rate in a previous experiment or analytical solution was verified.The effects of the vent tube structure and flow rate on the radon uniformity in the chamber;permeability,insertion depth,and flow rate on the radon exhalation rate and the effective diffusion coefficient on back-diffusion were investigated.Based on the results,increasing the inser-tion depth from 1 to 5 cm decreased the effective decay constant by 19.55%,whereas the curve-fitted radon exhalation rate decreased(lower than the initial value)as the deviation from the initial value increased by approximately 7%.Increasing the effective diffusion coefficient from 2.77×10^(-7) to 7.77×10^(-6) m^(2) s^(-1) made the deviation expand from 2.14 to 15.96%.The conclusion is that an increased insertion depth helps reduce leakage in the chamber,subject to notable back-diffusion,and that the closed-loop method is reasonably used for porous media with a low effective diffusion coefficient in view of the back-diffusion effect.The CFD-based simulation is expected to provide guidance for the optimization of the radon exhalation rate measurement method and,thus,the accurate measurement of the radon exhalation rate.
基金financially supported by the fellowship of the China Postdoctoral Science Foundation (2021M703466)the Natural Science Foundation of Gansu Province, China (20JR10RA231)the Natural Science Foundation of Gansu Province, China (22JR5RA050)。
文摘Wind-sand flow generates erosion and deposition around obstacles such as bridges and roadbeds, resulting in sand damage and endangering railway systems in sandy regions. Previous studies have mainly focused on the flow field around roadbeds, overlooking detailed examinations of sand particle erosion and deposition patterns near bridges and roadbeds. This study employs numerical simulations to analyze the influence of varying heights and wind speeds on sand deposition and erosion characteristics at different locations: the bridge-road transition section(side piers), middle piers, and roadbeds. The results show that the side piers, experience greater accumulation than the middle piers. Similarly, the leeward side of the roadbed witnesses more deposition compared to the windward side. Another finding reveals a reduced sand deposition length as the vertical profile, in alignment with the wind direction, moves further from the bridge abutments at the same clearance height. As wind speeds rise, there’s a decline in sand deposition and a marked increase in erosion around the side piers, middle piers and roadbeds. In conclusion, a bridge clearance that’s too low can cause intense sand damage near the side piers, while an extremely high roadbed may lead to extensive surface sand deposition. Hence, railway bridges in areas prone to sandy winds should strike a balance in clearance height. This research provides valuable guidelines for determining the most suitable bridge and roadbed heights in regions affected by wind and sand.
基金financially supported by the Natural Science Foundation of Gansu Province,China(22JR5RA050,20JR10RA231)the fellowship of the China Postdoctoral Science Foundation(2021M703466)the Basic Research Innovation Group Project of Gansu Province,China(21JR7RA347).
文摘For the safety of railroad operations,sand barriers are utilized to mitigate wind-sand disaster effects.These disasters,characterized by multi-directional wind patterns,result in diverse angles among the barriers.In this study,using numerical simulations,we examined the behavior of High Density Polyethylene(HDPE)sheet sand barriers under different wind angles,focusing on flow field distribution,windproof efficiency,and sedimentation erosion dynamics.This study discovered that at a steady wind speed,airflow velocity varies as the angle between the airflow and the HDPE barrier changes.Specifically,a 90°angle results in the widest low-speed airflow area on the barrier’s downwind side.If the airflow is not perpendicular to the barrier,it prompts a lateral airflow movement which decreases as the angle expands.The windproof efficiency correlates directly with this angle but inversely with the wind’s speed.Notably,with a wind angle of 90°,wind speed drops by 81%.The minimum wind speed is found at 5.1H(the sand barrier height)on the barrier’s downwind side.As the angle grows,the barrier’s windproof efficiency improves,extending its protective reach.Sedimentation is most prominent on the barrier’s downwind side,as the wind angle shifts from 30°to 90°,the sand sedimentation area on the barrier’s downwind side enlarges by 14.8H.As the angle grows,sedimentation intensifies,eventually overtakes the forward erosion and enlarges the sedimentation area.
基金the Natural Science Foundation of Jiangsu Province(Grant No.BK20200494)China Postdoctoral Science Foundation(Grant No.2021M701725)+3 种基金Jiangsu Postdoctoral Research Funding Program(Grant No.2021K522C)Fundamental Research Funds for the Central Universities(Grant No.30919011246)National Natural Science Foundation of China(Grant No.52278188)Natural Science Foundation of Jiangsu Province(Grant No.BK20211196)。
文摘To study the anti-explosion protection effect of polyurea coating on reinforced concrete box girder,two segmental girder specimens were made at a scale of 1:3,numbered as G(without polyurea coating)and PCG(with polyurea coating).The failure characteristics and dynamic responses of the specimens were compared through conducting explosion tests.The reliability of the numerical simulation using LS-DYNA software was verified by the test results.The effects of different scaled distances,reinforcement ratios,concrete strengths,coating thicknesses and ranges of polyurea were studied.The results show that the polyurea coating can effectively enhance the anti-explosion performance of the girder.The top plate of middle chamber in specimen G forms an elliptical penetrating hole,while that in specimen PCG only shows a very slight local dent.The peak vertical displacement and residual displacement of PCG decrease by 74.8% and 73.7%,respectively,compared with those of specimen G.For the TNT explosion with small equivalent,the polyurea coating has a more significant protective effect on reducing the size of fracture.With the increase of TNT equivalent,the protective effect of polyurea on reducing girder displacement becomes more significant.The optimal reinforcement ratio,concrete strength,thickness and range of polyurea coating were also drawn.
基金supported by the National Natural Science Foundation of China(Project Nos.12272270,11972261).
文摘Flash boiling atomization(FBA)is a promising approach for enhancing spray atomization,which can generate a fine and more evenly distributed spray by increasing the fuel injection temperature or reducing the ambient pressure.However,when the outlet speed of the nozzle exceeds 400 m/s,investigating high-speed flash boiling atomization(HFBA)becomes quite challenging.This difficulty arises fromthe involvement ofmany complex physical processes and the requirement for a very fine mesh in numerical simulations.In this study,an HFBA model for gasoline direct injection(GDI)is established.This model incorporates primary and secondary atomization,as well as vaporization and boilingmodels,to describe the development process of the flash boiling spray.Compared to lowspeed FBA,these physical processes significantly impact HFBA.In this model,the Eulerian description is utilized for modeling the gas,and the Lagrangian description is applied to model the droplets,which effectively captures the movement of the droplets and avoids excessive mesh in the Eulerian coordinates.Under various conditions,numerical solutions of the Sauter mean diameter(SMD)for GDI show good agreement with experimental data,validating the proposed model’s performance.Simulations based on this HFBA model investigate the influences of fuel injection temperature and ambient pressure on the atomization process.Numerical analyses of the velocity field,temperature field,vapor mass fraction distribution,particle size distribution,and spray penetration length under different superheat degrees reveal that high injection temperature or low ambient pressure significantly affects the formation of small and dispersed droplet distribution.This effect is conducive to the refinement of spray particles and enhances atomization.
基金financially supported by the Chang Jiang Scholar and Innovation Team Development Plan of China (IRT_15R29)the Basic Research Innovation Group Project of Gansu Province, China (21JR7RA347)the Natural Science Foundation of Gansu Province, China (20JR10RA231)。
文摘Wind and sand hazards are serious in the Milan Gobi area of the Xinjiang section of the Korla Railway. In order to ensure the safe operation of railroads, there is a need for wind and sand protection in heavily sandy areas. The wind and sand flow in the region is notably bi-directional. To shield railroads from sand, a unique sand fence made of folded linear high-density polyethylene(HDPE) is used, aligning with the principle that the dominant wind direction is perpendicular to the fence. This study employed field observations and numerical simulations to investigate the effectiveness of these HDPE sand fences in altering flow field distribution and offering protection. It also explored how these fences affect the deposition and erosion of sand particles. Findings revealed a significant reduction in wind speed near the fence corner;the minimum horizontal wind speed on the leeward side of the first sand fence(LSF) decreased dramatically from 3 m/s to 0.64 m/s. The vortex area on the LSF markedly impacted horizontal wind speeds. Within the LSF, sand deposition was a primary occurrence. As wind speeds increased, the deposition zone shrank, whereas the positive erosion zone expanded. Close to the folded corners of the HDPE sand fence, there was a notable shift from the positive erosion zone to a deposition zone. Field tests and numerical simulations confirmed the high windproof efficiency(WE) and sand resistance efficiency(SE) in the HDPE sand fence. Folded linear HDPE sheet sand fence can effectively slow down the incoming flow and reduce the sand content, thus achieving good wind and sand protection. This study provides essential theoretical guidance for the design and improvement of wind and sand protection systems in railroad engineering.
文摘Magnetohydrodynamic(MHD)induction pumps are contactless pumps able to withstand harsh environments.The rate of fluid flow through the pump directly affects the efficiency and stability of the device.To explore the influence of induction pump settings on the related delivery speed,in this study,a numerical model for coupled electromagnetic and flow field effects is introduced and used to simulate liquid metal lithium flow in the induction pump.The effects of current intensity,frequency,coil turns and coil winding size on the velocity of the working fluid are analyzed.It is shown that the first three parameters have a significant impact,while changes in the coil turns have a negligible influence.The maximum increase in working fluid velocity within the pump for the parameter combination investigated in this paper is approximately 618%.As the frequency is increased from 20 to 60 Hz,the maximum increase in the mean flow rate of the working fluid is approximately 241%.These research findings are intended to support the design and optimization of these devices.
基金supported by the National Natural Science Foundation of China(Grant Numbers:U22A20234,42277170)the Key Research and Development Project of Hubei Province(Grant Number:2020BCB073).
文摘To analyze the effects of a time-varying viscosity on the penetration length of grouting,in this study cement slur-ries with varying water-cement ratios have been investigated using the Bingham’sfluidflow equation and a dis-crete element method.Afluid-solid coupling numerical model has been introduced accordingly,and its accuracy has been validated through comparison of theoretical and numerical solutions.For different fracture forms(a single fracture,a branch fracture,and a fracture network),the influence of the time-varying viscosity on the slurry length range has been investigated,considering the change in the fracture aperture.The results show that under different fracture forms and the same grouting process conditions,the influence of the time-varying viscosity on the seepage length is 0.350 m.
基金supported by State Grid Ningxia Electric Power Co.,Ltd.under Grant 5229CG220006Natural Science Foundation of Ningxia Province under Grant 2022AAC03629.
文摘As an important lightning protection device in substations,lightning rods are susceptible to vibration and potential structural damage under wind loads.In order to understand their vibration mechanism,it is necessary to conduct flow analysis.In this study,numerical simulations of the flow field around a 330 kV cylindrical lightning rod with different diameters were performed using the SST k-ωmodel.The flow patterns in different segments of the lightning rod at the same reference wind speed(wind speed at a height of 10 m)and the flow patterns in the same segment at different reference wind speeds were investigated.The variations of lift coefficient,drag coefficient,and vorticity distribution were obtained.The results showed that vortex shedding phenomena occurred in all segments of the lightning rod,and the strength of vortex shedding increased with decreasing diameter.The vorticity magnitude and the root mean square magnitudes of the lift coefficient and drag coefficient also increased accordingly.The time history curves of the lift coefficient and drag coefficient on the surface of the lightning rod exhibited sinusoidal patterns with a single dominant frequency.For the same segment,as the wind speed increased in a certain range,the root mean square values of the lift coefficient and drag coefficient decreased,while their dominant frequencies increased.Moreover,there was a proportional relationship between the dominant frequencies of the lift coefficient and drag coefficient.The findings of this study can provide valuable insights for the refined design of lightning rods with similar structures.
文摘The radiant tube burner was modeled and analyzed by the numerical simulation method to investigate the influence factors and rules of NO_(x) emissions in a W-type radiant tube.These factors,which include air preheating temperature,excess air coefficient,and fuel gas composition,were modified to study their effects on NO_(x) emissions under varying working conditions.Simulation results were compared with the theoretical calculation value based on chemical reaction equilibrium theory and the onsite experimental value to verify the simulation accuracy.The results show that NO_(x) emissions rise with increasing air preheating temperatures.NO_(x) production increases to an extreme value and then decreases during the oxygen-poor to oxygen-enriched process with the rise of the excess air coefficient.Enhancing the proportion of coke oven gas in the fuel gas raises the combustion temperature as well as the NO_(x) discharge.Both the thermal efficiency and NO_(x) emissions should be balanced.Therefore,the recommended values based on the simulation results are as follows:the air preheating temperature should not exceed 400℃,the excess air coefficient should be between 1.1 and 1.2,and the volume fraction of the coke oven gas should not exceed 30%.
基金Sponsored by the National Natural Science Foundation of China(Grant No.52268048)the Guangxi Key Technology Research and Development Program(Grant No.GUI-KEAB23026101)the Guangxi Science and Technology Major Special Project(Grant No.GUI-KEAA22068066).
文摘In order to study the residual stress distribution law of welded joints of arch ribs of large-span steel pipe concrete arch bridges,numerical simulation of temperature,stress and strain fields based on ABAQUS for welded joints of arch-ribbed steel tubes using 7-,8-and 9-layer welds is carried out and its accuracy is demonstrated.The steel pipe welding temperature changes,residual stress distribution,different processes residual stress changes in the law,the prediction of post-weld residual stress distribution and deformation are studied in this paper.The results show that the temperature field values and test results are more consistent with the accuracy of numerical simulation of welding,the welding process is mainly in the form of heat transfer;Residual high stresses are predominantly distributed in the Fusion zone(FZ)and Heat-affected zone(HAZ),with residual stress levels tending to decrease from the center of the weld along the axial path,the maximum stress appears in the FZ and HAZ junction;The number of welding layers has an effect on the residual stress distribution,the number of welding layers increases,the residual stress tends to decrease,while the FZ and HAZ high stress area range shrinks;Increasing the number of plies will increase the amount of residual distortion.
基金supported by the project of the China Geological Survey(No.DD20221746)the National Natural Science Foundation of China(Grant Nos.41101086)。
文摘Xinqiao Gully is located in the area of the 2008 Wenchuan M_(s)8.0 earthquake in Sichuan province,China.Based on the investigation of the 2023"6-26"Xinqiao Gully debris flow event,this study assessed the effectiveness of the debris flow control project and evaluated the debris flow hazards.Through field investigation and numerical simulation methods,the indicators of flow intensity reduction rate and storage capacity fullness were proposed to quantify the effectiveness of the engineering measures in the debris flow event.The simulation results show that the debris flow control project reduced the flow intensity by41.05%to 64.61%.The storage capacity of the dam decreases gradually from upstream to the mouth of the gully,thus effectively intercepting and controlling the debris flow.By evaluating the debris flow of different recurrence intervals,further measures are recommended for managing debris flow events.
基金financially supported by the National Natural Science Foundation of China(No.52130404)the State Key Laboratory for GeoMechanics and Deep Underground Engineering China University of Mining&Technology,Beijing(No.SKLGDUEK2127)+1 种基金the Fundamental Research Funds for the Central Universities(Nos.FRF-TP-19-002C2Z,FRF-IDRY-20-031)and the Fundamental Research Funds for the Central Universities and the Youth Teacher International Exchange&Growth Program(No.QNXM20220002)。
文摘Cemented paste backfill(CPB)technology is a green mining method used to control underground goaves and tailings ponds.The curing process of CPB in the stope is the product of a thermo-hydro-mechanical-chemical multi-field performance interaction.At present,research on the multi-field performance of CPB mainly includes indoor similar simulation experiments,in-situ multi-field performance monitoring experiments,multi-field performance coupling model construction of CPB,and numerical simulation of the multi-field performance of CPB.Because it is hard to study the in-situ multi-field performance of CPB in the real stope,most current research on in-situ multi-field performance adopts the numerical simulation method.By simulating the conditions of CPB in the real stope(e.g.,maintenance environment,stope geometry,drainage conditions,and barricade and backfilling rates),the multi-field performance of CPB is further studied.This paper summarizes the mathematical models employed in the numerical simulation and lists the engineering application cases of numerical simulation in the in-situ multi-field performance of CPB.Finally,it proposes that the multi-field performance of CPB needs to strengthen the theoretical study of multi-field performance,form the strength design criterion based on the multi-field performance of CPB,perform a full-range numerical simulation of the multi-field performance of CPB,develop a pre-warning technology for the CPB safety of CPB,develop automatic and wireless sensors for the multi-field performance monitoring of CPB,and realize the application and popularization of CPB monitoring technology.
基金financially supported by the fellowship of the China Postdoctoral Science Foundation(2021M703466)the Natural Science Foundation of Gansu Province,China(20JR10RA231)+1 种基金the Basic Research Innovation Group Project of Gansu Province,China(21JR7RA347)an Special Funds for Guiding Local Scientific and Technological Development by the Central Government(22ZY1QA005)。
文摘Based on numerical simulations,this study highlights the sedimentation and erosion problems around a sand barrier through the relationship between the shear stress of the surface around the sand barrier and the critical shear stress of sand grains.The numerical simulation results were verified using data measured by the wind tunnel test.The results showed that when the porosity was the same,the size and position of the vortex on the leeward side of the sand barrier were related to the inlet wind speed.As the wind speed increased,the vortex volume increased and the positions of the separation and reattachment points moved toward the leeward side.When the porosity of the sand barrier was 30%,the strength of the acceleration zone above the sand barrier was the highest,and the strength of the acceleration zone was negatively correlated with the porosity.Sand erosion and sedimentation distance were related to wind speed.With an increase in wind speed,the sand grain forward erosion or reverse erosion areas on the leeward side of the sand barrier gradually replaced the sedimentation area.With an increase in porosity,the sand sedimentation distance on the leeward side of the sand barrier gradually shortened,and the sand erosion area gradually disappeared.The sand sedimentation distance on the leeward side of the sand barrier with 30%porosity was the longest.The numerical simulation results were in good agreement with the wind tunnel test results.Based on the sand erosion and sedimentation results of the numerical simulation and wind tunnel test,when the porosity was 30%,the protection effect of the High Density Polyethylene(HDPE)board sand barrier was best.
基金financially supported by the National Natural Science Foundation of China (No.51774031)。
文摘The dendrite growth behavior of high-strength steel during slab continuous casting with a traveling-wave magnetic field was studied in this paper. The morphology of the solidification structure and composition distribution were analyzed. Results showed that the columnar crystals could deflect and break when the traveling-wave magnetic field had low current intensity. With the increase in current intensity, the secondary dendrite arm spacing and solute permeability decreased, and the columnar crystal transformed into an equiaxed crystal. The electromagnetic force caused by the traveling-wave magnetic field changed the temperature gradient and velocity magnitude and promoted the breaking and fusing of dendrites. Dendrite compactness and composition uniformity were arranged in descending order as follows:columnar-toequiaxed transition (high current intensity), columnar crystal zone (low current intensity), columnar-to-equiaxed transition (low current intensity), and equiaxed crystal zone (high current intensity). Verified numerical simulation results combined with the boundary layer theory of solidification front and dendrite breaking–fusing model revealed the dendrite deflection mechanism and growth process. When thermal stress is not considered, and no narrow segment can be found in the dendrite, the velocity magnitude on the solidification front of liquid steel can reach up to 0.041 m/s before the dendrites break.
文摘Piles in a group experience additional displacements in soil due to pile-to-pile interactions apart from those resulting from the external loading.The effect of these interactions determined assuming soil as an elastic and/or viscoelastic material on pile head impedance functions of the pile group is studied by relating the group stiffness to the static stiffness of a single pile.However,the prevailing elastic solutions may misestimate the resulting pile group response due to the lack of consideration for either soil(material)and/or soil-pile interface nonlinearities.It is well established that soil behaves nonlinearly under moderate-to-high loading amplitudes,and besides,the soil-pile interface nonlinearity can exist even at small loading amplitudes.This study addresses the effects of these nonlinearities on the vertical impedance functions of a 3×3-pile group using numerical methods by direct analyses and superposition using pile-to-pile interaction factors.The numerical results are validated using scaled model tests under 1 g conditions.The results highlight the overestimation of pile-to-pile interactions in the pile group when assuming elastic soil conditions.The cases either by direct analyses or superposition approach involving soil and soil-pile interface nonlinearities agree well with the experimental pile group responses under close-to-elastic and nonlinear conditions.
基金support from the National Natural Science Foundation of China(NSFC)(52178324).
文摘A large amount of data can partly assure good fitting quality for the trained neural networks.When the quantity of experimental or on-site monitoring data is commonly insufficient and the quality is difficult to control in engineering practice,numerical simulations can provide a large amount of controlled high quality data.Once the neural networks are trained by such data,they can be used for predicting the properties/responses of the engineering objects instantly,saving the further computing efforts of simulation tools.Correspondingly,a strategy for efficiently transferring the input and output data used and obtained in numerical simulations to neural networks is desirable for engineers and programmers.In this work,we proposed a simple image representation strategy of numerical simulations,where the input and output data are all represented by images.The temporal and spatial information is kept and the data are greatly compressed.In addition,the results are readable for not only computers but also human resources.Some examples are given,indicating the effectiveness of the proposed strategy.
基金jointly supported by the National Natural Science Foundation of China(Grant No.42141013)the CNPC Innovation Found(2021 D002-1102)the China Postdoctoral Science Foundation(No.2021M701379)。
文摘Thermo-hydro-mechanical-chemical(THMC)interactions are prevalent during CO_(2)geological sequestration(CGS).In this study,a sequential coupling THMC numerical simulation program was constructed,which can be used to explore the following issues of CGS:fluid and heat flow,solute transport;stresses,displacements and rock failures related to geo-mechanical effects;equilibrium and kinetic chemical reactions;chemical damage to mechanical properties of the rock.Then,the coupling program was applied to the Ordos CGS Project to study the formation response under the multifield interaction caused by CO_(2)injection.The simulation results show that the mechanical process dominates the short CO_(2)injection period.Specifically,the formation’s permeability near the injection well increases by 43%,due to the reduction of effective stress,which significantly promotes the lateral migration of CO_(2).When the injection rate exceeds 0.15 million tons per year,the cohesion of the reservoir rock is not enough to resist the shear force inside the rock and rock failure may occur.During the subsequent long-term sequestration period(200 years),the influence of mineral reactions gradually increases.Due to calcite dissolution,the shear modulus of caprock is predicted to decrease by 7.6%,which will to some extent increase the risk of rock failure.