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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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%.展开更多
A transient 3D model was established to investigate the effect of spatial interaction of ultrasounds on the dual-frequency ultrasonic field in magnesium alloy melt.The effects of insertion depth and tip shape of the u...A transient 3D model was established to investigate the effect of spatial interaction of ultrasounds on the dual-frequency ultrasonic field in magnesium alloy melt.The effects of insertion depth and tip shape of the ultrasonic rods,input pressures and their ratio on the acoustic field distribution were discussed in detail.Additionally,the spacing,angle,and insertion depth of two ultrasonic rods significantly affect the interaction between distinct ultrasounds.As a result,various acoustic pressure distributions and cavitation regions are obtained.The spherical rods mitigate the longitudinal and transversal attenuation of acoustic pressure and expand the cavitation volume by 53.7%and 31.7%,respectively,compared to the plate and conical rods.Increasing the input pressure will enlarge the cavitation region but has no effect on the acoustic pressure distribution pattern.The acoustic pressure ratio significantly affects the pressure distribution and the cavitation region,and the best cavitation effect is obtained at the ratio of 2:1(P15:P20).展开更多
The Beishan pluton in Gansu of China was selected as the simulated model.The simulation results indicate that the formation of unloading joints in granite is mainly influenced by the unloading rate of confin-ing press...The Beishan pluton in Gansu of China was selected as the simulated model.The simulation results indicate that the formation of unloading joints in granite is mainly influenced by the unloading rate of confin-ing pressure.Among the rates tested,the slowest unloading rate 0.025 MPa/s is found to be most conducive to the development of unloading joints.Therefore,a slower unloading rate is favourable for the occurrence of unloading joints.A series of simulations with varying initial depths of uplift ranging from 900 m to 200 m were conducted.The results confirm that when the specimen rises to a depth of 550-500 m,the unloading joints begin to form.The uplift from a depth of 700-500 m,with variations in both vertical and lateral un-loading rates,was simulated.The generation of unloading joints exhibits a negative correlation with vertical unloading and no correlation with lateral unloading,indicating that the unloading joints are mainly controlled by the unloading of vertical pressure.Throughout the simulation process,the vertical joints exhibit irregular and unrealistic regularity,suggesting a more complex formation mechanism than that of the unloading joints.展开更多
The settling flux of biodeposition affects the environmental quality of cage culture areas and determines their environmental carrying capacity.Simple and effective simulation of the settling flux of biodeposition is ...The settling flux of biodeposition affects the environmental quality of cage culture areas and determines their environmental carrying capacity.Simple and effective simulation of the settling flux of biodeposition is extremely important for determining the spatial distribution of biodeposition.Theoretically,biodeposition in cage culture areas without specific emission rules can be simplified as point source pollution.Fluent is a fluid simulation software that can simulate the dispersion of particulate matter simply and efficiently.Based on the simplification of pollution sources and bays,the settling flux of biodeposition can be easily and effectively simulated by Fluent fluid software.In the present work,the feasibility of this method was evaluated by simulation of the settling flux of biodeposition in Maniao Bay,Hainan Province,China,and 20 sampling sites were selected for determining the settling fluxes.At sampling sites P1,P2,P3,P4,P5,Z1,Z2,Z3,Z4,A1,A2,A3,A4,B1,B2,C1,C2,C3 and C4,the measured settling fluxes of biodeposition were 26.02,15.78,10.77,58.16,6.57,72.17,12.37,12.11,106.64,150.96,22.59,11.41,18.03,7.90,19.23,7.06,11.84,5.19 and 2.57 g d^(−1)m^(−2),respectively.The simulated settling fluxes of biodeposition at the corresponding sites were 16.03,23.98,8.87,46.90,4.52,104.77,16.03,8.35,180.83,213.06,39.10,17.47,20.98,9.78,23.25,7.84,15.90,6.06 and 1.65 g d^(−1)m^(−2),respectively.There was a positive correlation between the simulated settling fluxes and measured ones(R=0.94,P=2.22×10^(−9)<0.05),which implies that the spatial differentiation of biodeposition flux was well simulated.Moreover,the posterior difference ratio of the simulation was 0.38,and the small error probability was 0.94,which means that the simulated results reached an acceptable level from the perspective of relative error.Thus,if nonpoint source pollution is simplified to point source pollution and open waters are simplified based on similarity theory,the setting flux of biodeposition in the open waters can be simply and effectively simulated by the fluid simulation software Fluent.展开更多
Asphaltene deposition is a significant problem during gas injection processes,as it can block the porous medium,the wellbore,and the involved facilities,significantly impacting reservoir productivity and ultimate oil re...Asphaltene deposition is a significant problem during gas injection processes,as it can block the porous medium,the wellbore,and the involved facilities,significantly impacting reservoir productivity and ultimate oil recovery.Only a few studies have investigated the numerical modeling of this potential effect in porous media.This study focuses on asphaltene deposition due to natural gas and CO_(2) injection.Predictions of the effect of gas injection on asphaltene deposition behavior have been made using a 3D numerical simulation model.The results indicate that the injection of natural gas exacerbates asphaltene deposition,leading to a significant reduction in permeability near the injection well and throughout the reservoir.This reduction in permeability strongly affects the ability of gas toflow through the reservoir,resulting in an improvement of the displacement front.The displacement effi-ciency of the injection gas process increases by up to 1.40%when gas is injected at 5500 psi,compared to the scenario where the asphaltene model is not considered.CO_(2) injection leads to a miscible process with crude oil,extracting light and intermediate components,which intensifies asphaltene precipitation and increases the viscosity of the remaining crude oil,ultimately reducing the recovery rate.展开更多
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.展开更多
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.展开更多
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 research on the multiphase flow characteristics of hydrate slurry is the key to implementing the risk prevention and control technology of hydrate slurry in deep-water oil and gas mixed transportation system.This ...The research on the multiphase flow characteristics of hydrate slurry is the key to implementing the risk prevention and control technology of hydrate slurry in deep-water oil and gas mixed transportation system.This paper established a geometric model based on the high-pressure hydrate slurry experimental loop.The model was used to carry out simulation research on the flow characteristics of gas-liquid-solid three-phase flow.The specific research is as follows:Firstly,the effects of factors such as slurry flow velocity,hydrate particle density,hydrate particle size,and hydrate volume fraction on the stratified smooth flow were specifically studied.Orthogonal test obtained particle size has the most influence on the particle concentration distribution.The slurry flow velocity is gradually increased based on stratified smooth flow.Various flow patterns were observed and their characteristics were analyzed.Secondly,increasing the slurry velocity to 2 m/s could achieve the slurry flow pattern of partial hydrate in the pipeline transition from stratified smooth flow to wavy flow.When the flow rate increases to 3 m/s,a violent wave forms throughout the entire loop.Based on wave flow,as the velocity increased to 4 m/s,and the flow pattern changed to slug flow.When the particle concentration was below 10%,the increase of the concentration would aggravate the slug flow trend;if the particle concentration was above 10%,the increase of the concentration would weaken the slug flow trend,the increase of particle density and liquid viscosity would weaken the tendency of slug flow.The relationship between the pressure drop gradients of several different flow patterns is:slug flow>wave flow>stratified smooth flow.展开更多
2G-NPR bolt (the 2nd generation Negative Poisson’s Ratio bolt) is a new type of bolt with high strength, high toughness and no yield platform. It has signifcant efects on improving the shear strength of jointed rock ...2G-NPR bolt (the 2nd generation Negative Poisson’s Ratio bolt) is a new type of bolt with high strength, high toughness and no yield platform. It has signifcant efects on improving the shear strength of jointed rock mass and controlling the stability of surrounding rock. To achieve an accurate simulation of bolted joint shear tests, we have studied a numerical simulation method that takes into account the 2G-NPR bolt's tensile–shear fracture criterion. Firstly, the indoor experimental study on the tensile–shear mechanical properties of 2G-NPR bolt is carried out to explore its mechanical properties under diferent tensile–shear angles, and the fracture criterion of 2G-NPR bolt considering the tensile–shear angle is established. Then, a three-dimensional numerical simulation method considering the tensile–shear mechanical constitutive and fracture criterion of 2G-NPR bolt, the elastoplastic mechanical behavior of surrounding rock and the damage and deterioration of grouting body is proposed. The feasibility and accuracy of the method are verifed by comparing with the indoor shear test results of 2G-NPR bolt anchorage joints. Finally, based on the numerical simulation results, the deformation and stress of the bolt, the distribution of the plastic zone of the rock mass, the stress distribution and the damage of the grouting body are analyzed in detail. The research results can provide a good reference value for the practical engineering application and shear mechanical performance analysis of 2G-NPR bolt.展开更多
基金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.
基金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.
基金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.
基金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.
基金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.
文摘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.
基金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.
基金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%.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.51974082 and 52274377)the Fundamental Research Funds for the Central Universities(Grant No.N2209001)the Programme of Introducing Talents of Discipline Innovation to Universities 2.0(the 111 Project 2.0 of China,Grant No.BP0719037)。
文摘A transient 3D model was established to investigate the effect of spatial interaction of ultrasounds on the dual-frequency ultrasonic field in magnesium alloy melt.The effects of insertion depth and tip shape of the ultrasonic rods,input pressures and their ratio on the acoustic field distribution were discussed in detail.Additionally,the spacing,angle,and insertion depth of two ultrasonic rods significantly affect the interaction between distinct ultrasounds.As a result,various acoustic pressure distributions and cavitation regions are obtained.The spherical rods mitigate the longitudinal and transversal attenuation of acoustic pressure and expand the cavitation volume by 53.7%and 31.7%,respectively,compared to the plate and conical rods.Increasing the input pressure will enlarge the cavitation region but has no effect on the acoustic pressure distribution pattern.The acoustic pressure ratio significantly affects the pressure distribution and the cavitation region,and the best cavitation effect is obtained at the ratio of 2:1(P15:P20).
文摘The Beishan pluton in Gansu of China was selected as the simulated model.The simulation results indicate that the formation of unloading joints in granite is mainly influenced by the unloading rate of confin-ing pressure.Among the rates tested,the slowest unloading rate 0.025 MPa/s is found to be most conducive to the development of unloading joints.Therefore,a slower unloading rate is favourable for the occurrence of unloading joints.A series of simulations with varying initial depths of uplift ranging from 900 m to 200 m were conducted.The results confirm that when the specimen rises to a depth of 550-500 m,the unloading joints begin to form.The uplift from a depth of 700-500 m,with variations in both vertical and lateral un-loading rates,was simulated.The generation of unloading joints exhibits a negative correlation with vertical unloading and no correlation with lateral unloading,indicating that the unloading joints are mainly controlled by the unloading of vertical pressure.Throughout the simulation process,the vertical joints exhibit irregular and unrealistic regularity,suggesting a more complex formation mechanism than that of the unloading joints.
基金support from the National Key Research and Development Program of China(No.2018YFD0900704)the National Natural Science Foundation of China(No.31972796).
文摘The settling flux of biodeposition affects the environmental quality of cage culture areas and determines their environmental carrying capacity.Simple and effective simulation of the settling flux of biodeposition is extremely important for determining the spatial distribution of biodeposition.Theoretically,biodeposition in cage culture areas without specific emission rules can be simplified as point source pollution.Fluent is a fluid simulation software that can simulate the dispersion of particulate matter simply and efficiently.Based on the simplification of pollution sources and bays,the settling flux of biodeposition can be easily and effectively simulated by Fluent fluid software.In the present work,the feasibility of this method was evaluated by simulation of the settling flux of biodeposition in Maniao Bay,Hainan Province,China,and 20 sampling sites were selected for determining the settling fluxes.At sampling sites P1,P2,P3,P4,P5,Z1,Z2,Z3,Z4,A1,A2,A3,A4,B1,B2,C1,C2,C3 and C4,the measured settling fluxes of biodeposition were 26.02,15.78,10.77,58.16,6.57,72.17,12.37,12.11,106.64,150.96,22.59,11.41,18.03,7.90,19.23,7.06,11.84,5.19 and 2.57 g d^(−1)m^(−2),respectively.The simulated settling fluxes of biodeposition at the corresponding sites were 16.03,23.98,8.87,46.90,4.52,104.77,16.03,8.35,180.83,213.06,39.10,17.47,20.98,9.78,23.25,7.84,15.90,6.06 and 1.65 g d^(−1)m^(−2),respectively.There was a positive correlation between the simulated settling fluxes and measured ones(R=0.94,P=2.22×10^(−9)<0.05),which implies that the spatial differentiation of biodeposition flux was well simulated.Moreover,the posterior difference ratio of the simulation was 0.38,and the small error probability was 0.94,which means that the simulated results reached an acceptable level from the perspective of relative error.Thus,if nonpoint source pollution is simplified to point source pollution and open waters are simplified based on similarity theory,the setting flux of biodeposition in the open waters can be simply and effectively simulated by the fluid simulation software Fluent.
基金funded by CNOOC Production Research Project(CCL2022SZPS0076).
文摘Asphaltene deposition is a significant problem during gas injection processes,as it can block the porous medium,the wellbore,and the involved facilities,significantly impacting reservoir productivity and ultimate oil recovery.Only a few studies have investigated the numerical modeling of this potential effect in porous media.This study focuses on asphaltene deposition due to natural gas and CO_(2) injection.Predictions of the effect of gas injection on asphaltene deposition behavior have been made using a 3D numerical simulation model.The results indicate that the injection of natural gas exacerbates asphaltene deposition,leading to a significant reduction in permeability near the injection well and throughout the reservoir.This reduction in permeability strongly affects the ability of gas toflow through the reservoir,resulting in an improvement of the displacement front.The displacement effi-ciency of the injection gas process increases by up to 1.40%when gas is injected at 5500 psi,compared to the scenario where the asphaltene model is not considered.CO_(2) injection leads to a miscible process with crude oil,extracting light and intermediate components,which intensifies asphaltene precipitation and increases the viscosity of the remaining crude oil,ultimately reducing the recovery rate.
基金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.
文摘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.
基金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.
基金supported by the National Natural Science Foundation of China(Grant No.52274061&52004039&51974037)China Postdoctoral Science Foundation(Grant No.2023T160717&2021M693908)+1 种基金The major project of universities affiliated to Jiangsu Province basic science(natural science)research(Grant No.21KJA440001)Jiangsu Qinglan Project,Changzhou Longcheng Talent Plan-Youth Science and Technology Talent Recruitment Project.
文摘The research on the multiphase flow characteristics of hydrate slurry is the key to implementing the risk prevention and control technology of hydrate slurry in deep-water oil and gas mixed transportation system.This paper established a geometric model based on the high-pressure hydrate slurry experimental loop.The model was used to carry out simulation research on the flow characteristics of gas-liquid-solid three-phase flow.The specific research is as follows:Firstly,the effects of factors such as slurry flow velocity,hydrate particle density,hydrate particle size,and hydrate volume fraction on the stratified smooth flow were specifically studied.Orthogonal test obtained particle size has the most influence on the particle concentration distribution.The slurry flow velocity is gradually increased based on stratified smooth flow.Various flow patterns were observed and their characteristics were analyzed.Secondly,increasing the slurry velocity to 2 m/s could achieve the slurry flow pattern of partial hydrate in the pipeline transition from stratified smooth flow to wavy flow.When the flow rate increases to 3 m/s,a violent wave forms throughout the entire loop.Based on wave flow,as the velocity increased to 4 m/s,and the flow pattern changed to slug flow.When the particle concentration was below 10%,the increase of the concentration would aggravate the slug flow trend;if the particle concentration was above 10%,the increase of the concentration would weaken the slug flow trend,the increase of particle density and liquid viscosity would weaken the tendency of slug flow.The relationship between the pressure drop gradients of several different flow patterns is:slug flow>wave flow>stratified smooth flow.
基金supported by the National Natural Science Foundation of China(NSFC)(41941018).
文摘2G-NPR bolt (the 2nd generation Negative Poisson’s Ratio bolt) is a new type of bolt with high strength, high toughness and no yield platform. It has signifcant efects on improving the shear strength of jointed rock mass and controlling the stability of surrounding rock. To achieve an accurate simulation of bolted joint shear tests, we have studied a numerical simulation method that takes into account the 2G-NPR bolt's tensile–shear fracture criterion. Firstly, the indoor experimental study on the tensile–shear mechanical properties of 2G-NPR bolt is carried out to explore its mechanical properties under diferent tensile–shear angles, and the fracture criterion of 2G-NPR bolt considering the tensile–shear angle is established. Then, a three-dimensional numerical simulation method considering the tensile–shear mechanical constitutive and fracture criterion of 2G-NPR bolt, the elastoplastic mechanical behavior of surrounding rock and the damage and deterioration of grouting body is proposed. The feasibility and accuracy of the method are verifed by comparing with the indoor shear test results of 2G-NPR bolt anchorage joints. Finally, based on the numerical simulation results, the deformation and stress of the bolt, the distribution of the plastic zone of the rock mass, the stress distribution and the damage of the grouting body are analyzed in detail. The research results can provide a good reference value for the practical engineering application and shear mechanical performance analysis of 2G-NPR bolt.