Numerical Simulation of Fire-Smoke Diffusion Caused by Vehicles in a Tunnel

Urban tunnels are generally narrow and fire smoke can hardly diffuse.In the present study,numerical simulation is used to analyze the diffusion of high temperature smoke produced by fire inside a specific tunnel(the Kaiyuan tunnel).The results are compared with similar data relating to other tests to determine the validity of the numerical method.Moreover,the critical velocity obtained by numerical simulation of 5 MW,20 MW,and 50 MW fires in curved and linear sections of the considered tunnel is compared with the values obtained using empirical formulas.The results show that,for the tunnel ventilation design,it is necessary to consider the fan pressurization at different sections and the fan pressurization should be higher at curved sections than that at linear sections.The safety of personnel escaping under different critical velocity values in the linear section has also been considered.On the basis of our findings,if only relying on natural ventilation,people can escape safely for the case of small fires,whereas for medium and large fires,it is necessary to turn on mechanical ventilation in time(and in order to avoid the danger caused by rapid diffusion of smoke,the timing of mechanical ventilation should be carefully tuned).

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.

Numerical simulation on the multiphase flow and reoxidation of the molten steel in a two-strand tundish during ladle change

A 3D mathematical model was proposed to investigate the molten steel–slag–air multiphase flow in a two-strand slab continuous casting(CC)tundish during ladle change.The study focused on the exposure of the molten steel and the subsequent reoxidation occurrence.The exposure of the molten steel was calculated using the coupled realizable k–εmodel and volume of fluid(VOF)model.The diffusion of dissolved oxygen was determined by solving the user-defined scalar(UDS)equation.Moreover,the user-defined function(UDF)was used to describe the source term in the UDS equation and determine the oxidation rate and oxidation position.The effect of the refilling speed on the molten steel exposure and dissolved oxygen content was also discussed.Increasing the refilling speed during ladle change reduced the refilling time and the exposure duration of the molten steel.However,the elevated refilling speed enlarged the slag eyes and increased the average dissolved oxygen content within the tundish,thereby exacerbating the reoxidation phenomenon.In addition,the time required for the molten steel with a high dissolved oxygen content to exit the tundish varied with the refilling speed.When the inlet speed was 3.0 m·s^(-1)during ladle change,the molten steel with a high dissolved oxygen content exited the outlet in a short period,reaching a maximum dissolved oxygen content of 0.000525wt%.Conversely,when the inlet speed was 1.8 m·s^(-1),the maximum dissolved oxygen content was 0.000382wt%.The refilling speed during the ladle change process must be appropriately decreased to minimize reoxidation effects and enhance the steel product quality.

In-situ measurement via the flow-through method and numerical simulations for radon exhalation during measurements of the radon exhalation rate

Small-scale measurements of the radon exhalation rate using the flow-through and closed-loop methods were conducted on the surface of a uranium tailing pond to better understand the differences between the two methods.An abnormal radon exhalation behavior was observed,leading to computational fluid dynamics(CFD)-based simulations in which dynamic radon migration in a porous medium and accumulation chamber was considered.Based on the in-situ experimental and numerical simulation results,variations in the radon exhalation rate subject to permeability,flow rate,and insertion depth were quantified and analyzed.The in-situ radon exhalation rates measured using the flow-through method were higher than those measured using the closed-loop method,which could be explained by the negative pressure difference between the inside and outside of the chamber during the measurements.The consistency of the variations in the radon exhalation rate between the experiments and simulations suggests the reliability of CFD-based techniques in obtaining the dynamic evolution of transient radon exhalation rates for diffusion and convection at the porous medium-air interface.The synergistic effects of the three factors(insertion depth,flow rate,and permeability)on the negative pressure difference and measured exhalation rate were quantified,and multivariate regression models were established,with positive correlations in most cases;the exhalation rate decreased with increasing insertion depth at a permeability of 1×10^(−11) m^(2).CFD-based simulations can provide theoretical guidance for improving the flow-through method and thus achieve accurate measurements.

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 of flow field deposition and erosion characteristics around bridge-road transition section

Wind-sand flow generates erosion and deposition around obstacles such as bridges and roadbeds, resulting in sand damage and endangering railway systems in sandy regions. Previous studies have mainly focused on the flow field around roadbeds, overlooking detailed examinations of sand particle erosion and deposition patterns near bridges and roadbeds. This study employs numerical simulations to analyze the influence of varying heights and wind speeds on sand deposition and erosion characteristics at different locations: the bridge-road transition section(side piers), middle piers, and roadbeds. The results show that the side piers, experience greater accumulation than the middle piers. Similarly, the leeward side of the roadbed witnesses more deposition compared to the windward side. Another finding reveals a reduced sand deposition length as the vertical profile, in alignment with the wind direction, moves further from the bridge abutments at the same clearance height. As wind speeds rise, there’s a decline in sand deposition and a marked increase in erosion around the side piers, middle piers and roadbeds. In conclusion, a bridge clearance that’s too low can cause intense sand damage near the side piers, while an extremely high roadbed may lead to extensive surface sand deposition. Hence, railway bridges in areas prone to sandy winds should strike a balance in clearance height. This research provides valuable guidelines for determining the most suitable bridge and roadbed heights in regions affected by wind and sand.

Determination method of mesh size for numerical simulation of blast load in near-ground detonation

In order to improve the overall resilience of the urban infrastructures, it is required to conduct blast resistant design for important building structures in the city. For complex terrain in the city, it is recommended to determine the blast load on the structures via numerical simulation. Since the mesh size of the numerical model highly depends on the explosion scenario, there is no generally applicable approach for the mesh size selection. An efficient method to determine the mesh size of the numerical model of near-ground detonation based on explosion scenarios is proposed in this study. The effect of mesh size on the propagation of blast wave under different explosive weights was studied, and the correlations between the mesh size effect and the charge weight or the scaled distance was described. Based on the principle of the finite element method and Hopkinson-Cranz scaling law, a mesh size measurement unit related to the explosive weight was proposed as the criterion for determining the mesh size in the numerical simulation. Finally, the applicability of the method proposed in this paper was verified by comparing the results from numerical simulation and the explosion tests and was verified in AUTODYN.

Numerical simulation of melt flow and temperature field during DC casting 2024 aluminium alloy under different casting conditions

Casting speed,casting temperature and secondary cooling water flow rate are the main process parameters affecting the DC casting process.These parameters significantly influence the flow and temperature fields during casting,which are crucial for the quality of the ingot and can determine the success or failure of the casting operation.Numerical simulation,with the advantages of low cost,rapid execution,and visualized results,is an important method to study and optimize the DC casting process.In the present work,a simulation model of DC casting 2024 aluminum alloy was established,and the reliability of the model was verified.Then,the influence of casting parameters on flow field and temperature field was studied in detail by numerical simulation method.Results show that with the increase of casting speed,the melt flow becomes faster,the depths of slurry zone and mushy zone increase,and the variation of slurry zone depth is greater than that of mushy zone.With an increase in casting temperature,the melt flow rate increases,the depth of the slurry zone becomes shallower,and the depth of the mushy zone experiences only minor changes.The simulation results further indicate that the increase of the flow rate of the secondary cooling water slightly reduces the depths of both slurry and mushy zone.

Numerical simulation on directional solidification and heat treatment processes of turbine blades

Study on turbine blades is crucial due to their critical role in ensuring the efficient and reliable operation of aircraft engines.Nickel-based single crystal superalloys are extensively used in the hot manufacturing of turbine blades due to their exceptional high-temperature mechanical properties.The hot manufacturing of single crystal blades involves directional solidification and heat treatment.Experimental manufacturing of these blades is time-consuming,capital-intensive,and often insufficient to meet industrial demands.Numerical simulation techniques have gained widespread acceptance in blade manufacturing research due to their low energy consumption,high efficiency,and rapid turnaround time.This article introduces the modeling and simulation of hot manufacturing in single crystal blades.The discussion outlines the prevalent mathematical models employed in numerical simulations related to blade hot manufacturing.It encapsulates the advancements in research concerning macro to micro-level numerical simulation techniques for directional solidification and heat treatment processes.Furthermore,potential future trajectories for the numerical simulation of single crystal blade hot manufacturing are also discussed.

Numerical Simulation and Entropy Production Analysis of Centrifugal Pump with Various Viscosity

The fluid’s viscosity significantly affects the performance of a centrifugal pump.The entropy production method and leakage are employed to analyze the performance changes under various viscosities by numerical simulation and validated by experiments.The results showed that increasing viscosity reduces both the pump head and efficiency.In addition,the optimal operating point shifts to the left.Leakage is influenced by vortex distribution in the front chamber and boundary layer thickness in wear-ring clearance,leading to an initial increase and subsequent decrease in leakage with increasing viscosity.The total entropy production Spro,Total inside the pump rises with increasing viscosity.The different mechanisms dominate under varying conditions:Turbulent dissipation dominates at low viscosity.Under high-viscosity conditions,energy loss is primarily caused by direct dissipation Spro,D and wall entropy production Spro,W.This study provides a deeper and more objective understanding of the energy characteristics of centrifugal pumps handling fluids of various viscosity,potentially aiding in optimizing pump design and improving energy conversion efficiency.

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.

Numerical Simulation and Parallel Computing of Acoustic Wave Equation in Isotropic-Heterogeneous Media

In this paper,we consider the numerical implementation of the 2D wave equation in isotropic-heterogeneous media.The stability analysis of the scheme using the von Neumann stability method has been studied.We conducted a study on modeling the propagation of acoustic waves in a heterogeneous medium and performed numerical simulations in various heterogeneous media at different time steps.Developed parallel code using Compute Unified Device Architecture(CUDA)technology and tested on domains of various sizes.Performance analysis showed that our parallel approach showed significant speedup compared to sequential code on the Central Processing Unit(CPU).The proposed parallel visualization simulator can be an important tool for numerous wave control systems in engineering practice.

Integrated numerical simulation of hydraulic fracturing and production in shale gas well considering gas-water two-phase flow

Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.

Numerical simulation study on the mold strength of magnetic mold casting based on a coupled electromagnetic-structural method

The properties of the magnetic mold in magnetic mold casting directly determine the quality of the final cast parts.In this study,the magnetic mold properties in magnetic mold casting,were studied utilizing a coupled electromagnetic-structural method through numerical simulation.This study investigated key factors including equivalent stress,the distribution of tensile and compressive stresses,and the area ratio of tensile stress.It compared molds made entirely of magnetic materials with those made partially of magnetic materials.Simulation results indicate that as current increases from 4 A to 8 A,both the initial magnetic mold and the material-replaced magnetic mold initially show an increasing trend in equivalent stress,tensile-compressive stress,and the area ratio of tensile stress,peaking at 6 A before declining.After material replacement,the area ratio of tensile stress at 6 A decreases to 19.84%,representing a reduction of 29.72%.Magnetic molds comprising a combination of magnetic and non-magnetic materials exhibit sufficient strength and a reduced area ratio of tensile stress compared to those made entirely from magnetic materials.This study provides valuable insights for optimizing magnetic mold casting processes and offers practical guidance for advancing the application of magnetic molds.

High-Precision Flow Numerical Simulation and Productivity Evaluation of Shale Oil Considering Stress Sensitivity

Continental shale oil reservoirs,characterized by numerous bedding planes and micro-nano scale pores,feature significantly higher stress sensitivity compared to other types of reservoirs.However,research on suitable stress sensitivity characterization models is still limited.In this study,three commonly used stress sensitivity models for shale oil reservoirs were considered,and experiments on representative core samples were conducted.By fitting and comparing the data,the“exponential model”was identified as a characterization model that accurately represents stress sensitivity in continental shale oil reservoirs.To validate the accuracy of the model,a two-phase seepage mathematical model for shale oil reservoirs coupled with the exponential model was introduced.The model was discretely solved using the finite volume method,and its accuracy was verified through the commercial simulator CMG.The study evaluated the productivity of a typical horizontal well under different engineering,geological,and fracture conditions.The results indicate that considering stress sensitivity leads to a 13.57%reduction in production for the same matrix permeability.Additionally,as the fracture half-length and the number of fractures increase,and the bottomhole flowing pressure decreases,the reservoir stress sensitivity becomes higher.

Experimental and Three-Dimensional Numerical Simulation of Phenomena Induced by Submerged Oblique Jet Scouring

Scouring experiments were conducted using a three-dimensional laser scanning technology for angles of the jet spanning the interval from 0°to 30°,and the characteristics of the scour hole in equilibrium conditions were investigated accordingly.The results indicate that the optimal scouring effects occur when the jet angle is in the ranges between 15°and 20°.Moreover,the dimensionless profiles of the scour hole exhibit a high degree of similarity at different jet angles.Numerical simulations conducted using the Flow-3D software to investigate the bed shear stress along the jet impingement surface have shown that this stress is influenced by both the resultant force and the jet impingement surface area.It reaches its maximum value when the jet is vertical,decreases rapidly as the jet starts to tilt,then increases slightly,and decreases again significantly when the angle exceeds 20°.

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.

Coupled Numerical Simulation of Electromagnetic and Flow Fields in a Magnetohydrodynamic Induction Pump

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.

Numerical Simulation of Oil-Water Two-Phase Flow in Low Permeability Tight Reservoirs Based on Weighted Least Squares Meshless Method

In response to the complex characteristics of actual low-permeability tight reservoirs,this study develops a meshless-based numerical simulation method for oil-water two-phase flow in these reservoirs,considering complex boundary shapes.Utilizing radial basis function point interpolation,the method approximates shape functions for unknown functions within the nodal influence domain.The shape functions constructed by the aforementioned meshless interpolation method haveδ-function properties,which facilitate the handling of essential aspects like the controlled bottom-hole flow pressure in horizontal wells.Moreover,the meshless method offers greater flexibility and freedom compared to grid cell discretization,making it simpler to discretize complex geometries.A variational principle for the flow control equation group is introduced using a weighted least squares meshless method,and the pressure distribution is solved implicitly.Example results demonstrate that the computational outcomes of the meshless point cloud model,which has a relatively small degree of freedom,are in close agreement with those of the Discrete Fracture Model(DFM)employing refined grid partitioning,with pressure calculation accuracy exceeding 98.2%.Compared to high-resolution grid-based computational methods,the meshless method can achieve a better balance between computational efficiency and accuracy.Additionally,the impact of fracture half-length on the productivity of horizontal wells is discussed.The results indicate that increasing the fracture half-length is an effective strategy for enhancing production from the perspective of cumulative oil production.

Numerical Simulation of the Mechanical Stirring Process in a Tannin-Based Foaming Precursor Resin

Tannin foam is a new functional material.It can be widely applied to the automobile industry,construction industry,and packaging industry due to its wide range of raw materials,renewable,easily degraded,low cost and almost no pollution.Preparing tannin foam is a very complex process that includes high temperature,two phases,mechanical agitation,and phase change.To investigate the influence of the stirring velocity and paddle shape,simulation was calculated by making use of the volume of fluid(VOF)method and multiple reference frame(MRF)method in a three-dimensional flow field of tannin-based foaming precursor resin.The gas holdup and velocity magnitude were analysed with various conditions of mechanical velocities and paddle shape in the stirring flow field.The result shows the higher the velocity,the greater the disturbance and paddle shape between the eggbeater and the Rushton turbine,obviously the paddle shape of the eggbeater with a wider range of agitation,which can entrap more air into the tannin-based foaming precursor resin in a short time.Especially when the speed is 1500 rpm,the flow field of the Rushton turbine comes out of a ditch,which decreases the efficiency of mass transfer;there is less air to mix into the tannin-based foaming precursor resin,which causes unevenness.At the same time,the eggbeater shows the marvelous capability of hybrid as it has two vortexes and multiple cycles that make a difference from the Rushton turbine,which has only one vortex and two upper and lower loops;the structure makes the flow field more stable allowed evenness of flow field tannin-based foaming precursor resin.The results reveal that it is beneficial for tannin-based foaming precursor resin to use an eggbeater with a speed of 1500 rpm to reduce the consumption of resources while obtaining a uniform flow field.