Multi-physical fields distribution in billet during helical electromagnetic stirring:A numerical simulation research

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.

Explosion resistance performance of reinforced concrete box girder coated with polyurea:Model test and numerical simulation

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.

Analysis of debris flow control effect and hazard assessment in Xinqiao Gully,Wenchuan M_(s)8.0 earthquake area based on numerical simulation

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.

Numerical simulation on sand sedimentation and erosion characteristics around HDPE sheet sand barrier under different wind angles

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.

Numerical study on the cavity dynamics for vertical water entries of twin spheres

In this study, a three dimensional(3D) numerical model of six-degrees-of-freedom(6DOF) is applied to simulate the water entries of twin spheres side-by-side at different lateral distances and time intervals.The turbulence structure is described using the shear-stress transport k-ω(SST k-ω) model, and the volume of fluid(VOF) method is used to track the complex air-liquid interface. The motion of spheres during water entry is simulated using an independent overset grid. The numerical model is verified by comparing the cavity evolution results from simulations and experiments. Numerical results reveal that the time interval between the twin water entries evidently affects cavity expansion and contraction behaviors in the radial direction. However, this influence is significantly weakened by increasing the lateral distance between the two spheres. In synchronous water entries, pressure is reduced on the midline of two cavities during surface closure, which is directly related to the cavity volume. The evolution of vortexes inside the two cavities is analyzed using a velocity vector field, which is affected by the lateral distance and time interval of water entries.

Investigations on High-Speed Flash Boiling Atomization of Fuel Based on Numerical Simulations

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.

Flow field, sedimentation, and erosion characteristics around folded linear HDPE sheet sand fence: Numerical simulation study

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.

Numerical Analysis of the Influence of the Impinging Distance on the Scouring Efficiency of Submerged Jets

Water jet technology is widely used in submerged buried pipes as a non-traditional trenching process,often invol-ving a complex sediment response.An important adjustable and influential engineering variable in this technol-ogy is represented by the impinging distance.In this study,the FLOW-3D software was used to simulate the jet scouring of sand beds in a submerged environment.In particular,four sets of experimental conditions were con-sidered to discern the relationship between the maximum scour depth and mass and the impinging distance.As shown by the results,a critical impinging distance h0 exists by which the static scour depth can be maximized;the scour mass ratio between dynamic and static conditions decreases as the impinging distance increases.Moreover,the profile contours are similar when the erosion parameter Ec is in the range 0.35<Ec<2.Empirical equations applicable for predicting the jet trenching contour under both dynamic and static scour modes are also provided in this study.

Numerical simulations for radon migration and exhalation behavior during measuring radon exhalation rate with closed-loop method

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.

Numerical Study on the Effect of Vortex Generators on the Aerodynamic Drag of a High-Speed Train

A relatively high aerodynamic drag is an important factor that hinders the further acceleration of high-speed trains.Using the shear stress transport(SST)k-ωturbulence model,the effect of various vortex generator types on the aerodynamic characteristics of an ICE2(Inter-city Electricity)train has been investigated.The results indi-cate that the vortex generators with wider triangle,trapezoid,and micro-ramp arranged on the surface of the tail car can significantly change the distribution of surface pressure and affect the vorticity intensity in the wake.This alteration effectively reduces the resistance of the tail car.Meanwhile,the micro-ramp vortex generator with its convergent structure at the rear exhibits enhancedflow-guiding capabilities,resulting in a 15.4%reduction in the drag of the tail car.

Numerical Simulation of Slurry Diffusion in Fractured Rocks Considering a Time-Varying Viscosity

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.

Numerical Simulations of the Flow Field around a Cylindrical Lightning Rod

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.

Numerical simulation of the influence factors of NO_(x) emission on a W-type radiant tube

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%.

The Effects of Thickness and Location of PCMon the Building’s Passive Temperature-Control–A Numerical Study

Building energy consumption and building carbon emissions both account for more than 20%of their total national values in China.Building employing phase change material(PCM)for passive temperature control shows a promising prospect in meeting the comfort demand and reducing energy consumption simultaneously.However,there is a lack of more detailed research on the interaction between the location and thickness of PCM and indoor natural convection,as well as indoor temperature distribution.In this study,the numerical model of a passive temperature-controlled building integrating the developed PCM module is established with the help of ANSYS.In which,the actual weather condition of Beijing city is set as the boundary conditions and the indoor natural convection is simulated with the consideration of radiation model.The effects of PCM’s thickness and location on the internal temperature field are analyzed and discussed.The results show that the room could maintain within the human comfort temperature range with the longest ratio of 94.10%and the shortest ratio of 51.04%as integrating PCM.In comparison,the value is only 26.70%without PCM.The room’s maximum temperature fluctuation can also be improved;it could be lowered by 64.4%compared to the normal condition.When the quantity of PCM is sufficient,further increasing the PCM amount results in a temperature fluctuation reduction of less than 0.1°C and does not increase the comfort time.Placing PCM on the wall induces an apparent variation in indoor temperature along the vertical direction.Conversely,placing PCM on the roof can lead to a heat transfer rate difference of up to seven times.The optimal placement of PCM depends on the difference between the environmental and phase change temperatures.If the difference is positive,placing PCM on the roof is more effective;conversely,the opposite holds.According to the results over the entire cycle,PCM application on vertical walls yields better performance.The significant difference in natural convection caused by the same thickness of PCM but different application positions,coupled with the influence of air movement on the melting and solidification of PCM,further impacts indoor temperature fluctuations and comfort.This study can provide guidance for the application location and thickness of PCM,especially for scenarios where temperature regulation is required at a specific time.

Numerical Modeling of Mass Transfer in the Interaction between River Biofilm and a Turbulent Boundary Layer

In this article dedicated to the modeling of vertical mass transfers between the biofilm and the bulk flow, we have, in the first instance, presented the methodology used, followed by the presentation of various results obtained through analyses conducted on velocity fields, different fluxes, and overall transfer coefficients. Due to numerical constraints (resolution of relevant spatial scales), we have restricted the analysis to low Schmidt numbers (Sc=0.1, Sc=1, and Sc=10) and a single roughness Reynolds number (Re*=150). The analysis of instantaneous concentration fields from various simulations revealed logarithmic concentration profiles above the canopy. In this zone, the concentration is relatively homogeneous for longer times. The analysis of results also showed that the contribution of molecular diffusion to the total flux depends on the Schmidt number. This contribution is negligible for Schmidt numbers Sc≥0.1, but nearly balances the turbulent flux for Sc=0.1. In the canopy, the local Sherwood number, given by the ratio of the total flux (within or above the canopy) to the molecular diffusion flux at the wall, also depends on the Schmidt number and varies significantly between the canopy and the region above. The exchange velocity, a purely hydrodynamic parameter, is independent of the Schmidt number and is on the order of 10% of in the present case. This study also reveals that nutrient absorption by organisms near the wall depends on the Schmidt number. Such absorption is facilitated by lower Schmidt numbers.

Experimental research and numerical simulation of the multi-field performance of cemented paste backfill:Review and future perspectives

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.

Numerical investigation of the mechanical behavior of the backfill–rock composite structure under triaxial compression

To ensure safe and economical backfill mining,the mechanical response of the backfill–rock interaction system needs to be understood.The numerical investigation of the mechanical behavior of backfill–rock composite structure(BRCS)under triaxial compression,which includes deformation,failure patterns,strength characteristics,and acoustic emission(AE)evolution,was proposed.The models used in the tests have one rough interface,two cement–iron tailings ratios(CTRs),four interface angles(IAs),and three confining pressures(CPs).Results showed that the deformation,strength characteristics,and failure patterns of BRCS under triaxial compression depend on IA,CP,and CTR.The stress–strain curves of BRCS under triaxial compression could be divided into five stages,namely,compaction,elasticity,yield,strain softening,and residual stress.The relevant AE counts have corresponding relationships with different stages.The triaxial compressive strengths of composites increase linearly with the increase of the CP.Furthermore,the CP stress strengthening effect occurs.When the IAs are45°and 60°,the failure areas of composites appear in the interface and backfill.When the IAs are 75°and 90°,the failure areas of composites appear in the backfill,interface,and rock.Moreover,the corresponding failure modes yield the combined shear failure.The research results provide the basis for further understanding of the stability of the BRCS.

Effects of sand sedimentation and wind erosion around sand barrier:Numerical simulation and wind tunnel test studies

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.

Vertical impedance functions of pile groups under low-to-high loading amplitudes:numerical simulations and experimental validation

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.

Numerical Simulation of Thermo-Hydro-Mechanical-Chemical Response Caused by CO_(2)Injection into Saline Geological Formations:A Case Study from the Ordos Project,China

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.