The kinetic characteristics of plasma-assisted oxidative pyrolysis of ammonia are studied by using the global/fluid models hybrid solution method.Firstly,the stable products of plasma-assisted oxidative pyrolysis of a...The kinetic characteristics of plasma-assisted oxidative pyrolysis of ammonia are studied by using the global/fluid models hybrid solution method.Firstly,the stable products of plasma-assisted oxidative pyrolysis of ammonia are measured.The results show that the consumption of NH_(3)/O_(2)and the production of N_(2)/H_(2)change linearly with the increase of voltage,which indicates the decoupling of nonequilibrium molecular excitation and oxidative pyrolysis of ammonia at low temperatures.Secondly,the detailed reaction kinetics mechanism of ammonia oxidative pyrolysis stimulated by a nanosecond pulse voltage at low pressure and room temperature is established.Based on the reaction path analysis,the simplified mechanism is obtained.The detailed and simplified mechanism simulation results are compared with experimental data to verify the accuracy of the simplified mechanism.Finally,based on the simplified mechanism,the fluid model of ammonia oxidative pyrolysis stimulated by the nanosecond pulse plasma is established to study the pre-sheath/sheath behavior and the resultant consumption and formation of key species.The results show that the generation,development,and propagation of the pre-sheath have a great influence on the formation and consumption of species.The consumption of NH_(3)by the cathode pre-sheath is greater than that by the anode pre-sheath,but the opposite is true for OH and O(1S).However,within the sheath,almost all reactions do not occur.Further,by changing the parameters of nanosecond pulse power supply voltage,it is found that the electron number density,electron current density,and applied peak voltages are not the direct reasons for the structural changes of the sheath and pre-sheath.Furthermore,the discharge interval has little effect on the sheath structure and gas mixture breakdown.The research results of this paper not only help to understand the kinetic promotion of non-equilibrium excitation in the process of oxidative pyrolysis but also help to explore the influence of transport and chemical reaction kinetics on the oxidative pyrolysis of ammonia.展开更多
Configuring computational fluid dynamics(CFD)simulations typically demands extensive domain expertise,limiting broader access.Although large language models(LLMs)have advanced scientific computing,their use in automat...Configuring computational fluid dynamics(CFD)simulations typically demands extensive domain expertise,limiting broader access.Although large language models(LLMs)have advanced scientific computing,their use in automating CFD workflows is underdeveloped.We introduce a novel approach centered on domain-specific LLM adaptation.By fine-tuning Qwen2.5-7B-Instruct on NL2FOAM,our custom dataset of 28,716 natural language-to-OpenFOAM configuration pairs with chain-of-thought(CoT)annotations enables direct translation from natural language descriptions to executable CFD setups.A multi-agent system orchestrates the process,autonomously verifying inputs,generating configurations,running simulations,and correcting errors.Evaluation on a benchmark of 21 diverse flow cases demonstrates state-of-the-art performance,achieving 88.7%solution accuracy and 82.6%first-attempt success rate.This significantly outperforms larger general-purpose models such as Qwen2.5-72B-Instruct,DeepSeek-R1,and Llama3.3-70B-Instruct,while also requiring fewer correction iterations and maintaining high computational efficiency.The results highlight the critical role of domain-specific adaptation in deploying LLM assistants for complex engineering workflows.Our code and fine-tuned model have been deposited at https://github.com/YYgroup/AutoCFD.展开更多
This work investigates water-based micropolar hybrid nanofluid(MHNF) flow on an elongating variable porous sheet.Nanoparticles of diamond and copper have been used in the water to boost its thermal conductivity. The m...This work investigates water-based micropolar hybrid nanofluid(MHNF) flow on an elongating variable porous sheet.Nanoparticles of diamond and copper have been used in the water to boost its thermal conductivity. The motion of the fluid is taken as two-dimensional with the impact of a magnetic field in the normal direction. The variable, permeable, and stretching nature of sheet's surface sets the fluid into motion. Thermal and mass diffusions are controlled through the use of the Cattaneo–Christov flux model. A dataset is generated using MATLAB bvp4c package solver and employed to train an artificial neural network(ANN) based on the Levenberg–Marquardt back-propagation(LMBP) algorithm. It has been observed as an outcome of this study that the modeled problem achieves peak performance at epochs 637, 112, 4848, and 344 using ANN-LMBP. The linear velocity of the fluid weakens with progression in variable porous and magnetic factors.With an augmentation in magnetic factor, the micro-rotational velocity profiles are augmented on the domain 0 ≤ η < 1.5 due to the support of micro-rotations by Lorentz forces close to the sheet's surface, while they are suppressed on the domain 1.5 ≤ η < 6.0 due to opposing micro-rotations away from the sheet's surface. Thermal distributions are augmented with an upsurge in thermophoresis, Brownian motion, magnetic, and radiation factors, while they are suppressed with an upsurge in thermal relaxation parameter. Concentration profiles increase with an expansion in thermophoresis factor and are suppressed with an intensification of Brownian motion factor and solute relaxation factor. The absolute errors(AEs) are evaluated for all the four scenarios that fall within the range 10^(-3)–10^(-8) and are associated with the corresponding ANN configuration that demonstrates a fine degree of accuracy.展开更多
This paper investigates the capabilities of large language models(LLMs)to leverage,learn and create knowledge in solving computational fluid dynamics(CFD)problems through three categories of baseline problems.These ca...This paper investigates the capabilities of large language models(LLMs)to leverage,learn and create knowledge in solving computational fluid dynamics(CFD)problems through three categories of baseline problems.These categories include(1)conventional CFD problems that can be solved using existing numerical methods in LLMs,such as lid-driven cavity flow and the Sod shock tube problem;(2)problems that require new numerical methods beyond those available in LLMs,such as the recently developed Chien-physics-informed neural networks for singularly perturbed convection-diffusion equations;and(3)problems that cannot be solved using existing numerical methods in LLMs,such as the ill-conditioned Hilbert linear algebraic systems.The evaluations indicate that reasoning LLMs overall outperform non-reasoning models in four test cases.Reasoning LLMs show excellent performance for CFD problems according to the tailored prompts,but their current capability in autonomous knowledge exploration and creation needs to be enhanced.展开更多
Based on the synchronous joint gravity and magnetic inversion of single interface by Pilkington and the need of revealing Cenozoic and crystalline basement thickness in the new round of oil-gas exploration, we propose...Based on the synchronous joint gravity and magnetic inversion of single interface by Pilkington and the need of revealing Cenozoic and crystalline basement thickness in the new round of oil-gas exploration, we propose a joint gravity and magnetic inversion methodfor two-layer models by concentrating on the relationship between the change of thicknessI and position of the middle layer and anomaly and discuss the effects of the key parameters. Model tests and application to field data show the validity of this method.展开更多
A 1∶8 physical water model was constructed to investigate the fluid flow and mixing phenomena in the basic oxygen furnace(BOF)converter.The particle image velocimetry was employed to measure the velocity distribution...A 1∶8 physical water model was constructed to investigate the fluid flow and mixing phenomena in the basic oxygen furnace(BOF)converter.The particle image velocimetry was employed to measure the velocity distribution of the bath and the high-speed camera was applied to capture the cavity shape in the combined blowing BOF converter.The mixing time for varied operating conditions was measured by the stimulus-response approach.The cavity depth increased with the decrease in the lance height and the increase in the top gas flow rate while the bottom blowing gas had little influence on the cavity depth.The minimum cavity depth was obtained under the condition of a 69.8 m^(3)/h top gas flow rate,a 287.5 mm lance height and a 0.93 m^(3)/h bottom blowing gas flow rate,which was 161.2 mm.The mixing time decreased as the lance height decreased and the top blowing gas flow rate increased.The mixing time was first decreased and then increased with the increase in the bottom gas flow rate.With the condition of 69.8 m^(3)/h gas flow rate of top blowing,the 287.5 mm lance height and the 0.93 m^(3)/h gas flow rate of bottom blowing,the mixing time in the converter was 48.65 s.The empirical formula between the stirring power and the mixing time in the converter was calculated.展开更多
The hydroelastic response of a circular, very large floating structure(VLFS), idealized as a floating circular elastic thin plate, is investigated for the case of time-harmonic incident waves of the surface and interf...The hydroelastic response of a circular, very large floating structure(VLFS), idealized as a floating circular elastic thin plate, is investigated for the case of time-harmonic incident waves of the surface and interfacial wave modes, of a given wave frequency, on a two-layer fluid of finite and constant depth. In linear potential-flow theory, with the aid of angular eigenfunction expansions, the diffraction potentials can be expressed by the Bessel functions. A system of simultaneous equations is derived by matching the velocity and the pressure between the open-water and the platecovered regions, while incorporating the edge conditions of the plate. Then the complex nested series are simplified by utilizing the orthogonality of the vertical eigenfunctions in the open-water region. Numerical computations are presented to investigate the effects of different physical quantities, such as the thickness of the plate, Young’s modulus, the ratios of the densities and of the layer depths, on the dispersion relations of the flexural-gravity waves for the two-layer fluid. Rapid convergence of the method is observed, but is slower at higher wave frequency. At high frequency, it is found that there is some energy transferred from the interfacial mode to the surface mode.展开更多
A weak nonlinear model of a two-layer barotropic ocean with Rayleigh dissipation is built.The analytic asymptotic solution is derived in the mid-latitude stationary wind field,and the physical meaning of the correspon...A weak nonlinear model of a two-layer barotropic ocean with Rayleigh dissipation is built.The analytic asymptotic solution is derived in the mid-latitude stationary wind field,and the physical meaning of the corresponding problem is discussed.展开更多
In the present paper, we endeavor to accomplish a diagram, which demarcates the validity ranges for interfacial wave theories in a two-layer system, to meet the needs of design in ocean engineering. On the basis of th...In the present paper, we endeavor to accomplish a diagram, which demarcates the validity ranges for interfacial wave theories in a two-layer system, to meet the needs of design in ocean engineering. On the basis of the available solutions of periodic and solitary waves, we propose a guideline as principle to identify the validity regions of the interfacial wave theories in terms of wave period T, wave height H, upper layer thickness dl, and lower layer thick-ness d2, instead of only one parameter-water depth d as in the water surface wave circumstance. The diagram proposed here happens to be Le Mehaute＇s plot for free surface waves if water depth ratio r= d1/d2 approaches to infinity and the upper layer water density p1 to zero. On the contrary, the diagram for water surface waves can be used for two-layer interfacial waves if gravity acceleration g in it is replaced by the reduced gravity defined in this study under the condition of σ=（P2 - Pl）/P2 → 1.0 and r 〉 1.0. In the end, several figures of the validity ranges for various interfacial wavetheories in the two-layer fluid are given and compared with the results for surface waves.展开更多
In this paper, the diffraction of water waves by a vertically floating cylinder in a two-layer fluid of a finite depth is studied. Analytical expressions for the hydrodynamic loads on the vertically floating cylinder ...In this paper, the diffraction of water waves by a vertically floating cylinder in a two-layer fluid of a finite depth is studied. Analytical expressions for the hydrodynamic loads on the vertically floating cylinder are obtained by use of the method of eigenfunction expansions. The hydrodynamic loads on the vertically floating cylinder in a two-layer fluid inelude not only the surge, heave and pitch exciting forces due to the incident wave of the surface-wave mode, but also those due to the incident wave of the internal-wave mode. This is different from the case of a homogenous fluid. Some given examples show that, for a two-layer fluid system with a small density difference, the hydrodynamic loads for the surface-wave mode do not differ significantly from those due to surface waves in a single-layer fluid, but the hydrodynamic loads for the internal-wave mode are important over a wide range of frequencies. Moreover, also considered are the free surface and interface elevations generated by the diffraction wave due to the incident wave of the surface-wave and interhal-wave modes, and transfer of energy between modes.展开更多
Many new forms of Boussinesq-type equations have been developed to extend the range of applicability of the classical Boussinesq equations to deeper water in the Study of the surface waves. One approach was used by Nw...Many new forms of Boussinesq-type equations have been developed to extend the range of applicability of the classical Boussinesq equations to deeper water in the Study of the surface waves. One approach was used by Nwogu （1993. J. Wtrw. Port Coastal and Oc. Eng. 119, 618-638） to improve the linear dispersion characteristics of the classical Boussinesq equations by using the velocity at an arbitrary level as the velocity variable in derived equations and obtain a new form of Boussinesq-type equations, in which the dispersion property can be optimized by choosing the velocity variable at an adequate level. In this paper, a set of Boussinesq-type equations describing the motions of the interracial waves propagating alone the interface between two homogeneous incompressible and inviscid fluids of different densities with a free surface and a variable water depth were derived using a method similar to that used by Nwogu （1993. J. Wtrw. Port Coastal and Oc. Eng. 119, 618-638） for surface waves. The equations were expressed in terms of the displacements of free surface and density-interface, and the velocity vectors at arbitrary vertical locations in the upper layer and the lower layer （or depth-averaged velocity vector across each layer） of a two-layer fluid. As expected, the equations derived in the present work include as special cases those obtained by Nwogu （1993, J. Wtrw. Port Coastal and Oc. Eng. 119, 618-638） and Peregrine （1967, J. Fluid Mech. 27, 815-827） for surface waves when the density of the upper fluid is taken as zero.展开更多
A quite general coupled variable coefficient modified KdV (VCmKdV) equation in a two-layer fluid systemis derived by means of the reductive perturbation method.Making use of the CK's direct method,some similarityr...A quite general coupled variable coefficient modified KdV (VCmKdV) equation in a two-layer fluid systemis derived by means of the reductive perturbation method.Making use of the CK's direct method,some similarityreductions of the coupled VCmKdV equation are obtained and their corresponding group explanations are discussed.Some exact solutions of the coupled equations are also presented.展开更多
The problem of oblique wave （internal wave） propagation over a small deformation in a channel flow consisting of two layers was considered. The upper fluid was assumed to be bounded above by a rigid lid, which is an...The problem of oblique wave （internal wave） propagation over a small deformation in a channel flow consisting of two layers was considered. The upper fluid was assumed to be bounded above by a rigid lid, which is an approximation for the free surface, and the lower one was bounded below by an impermeable bottom surface having a small deformation; the channel was unbounded in the horizontal directions. Assuming irrotational motion, the perturbation technique was employed to calculate the first-order corrections of the velocity potential in the two fluids by using Green＇s integral theorem suitably with the introduction of appropriate Green＇s functions. Those functions help in calculating the reflection and transmission coefficients in terms of integrals involving the shape ftmction c（x） representing the bottom deformation. Three-dimensional linear water wave theory was utilized for formulating the relevant boundary value problem. Two special examples of bottom deformation were considered to validate the results. Consideration of a patch of sinusoidal ripples （having the same wave number） shows that the reflection coefficient is an oscillatory function of the ratio of twice the x-component of the wave number to the ripple wave number. When this ratio approaches one, the theory predicts a resonant interaction between the bed and the interface, and the reflection coefficient becomes a multiple of the number of ripples. High reflection of incident wave energy occurs if this number is large. Similar results were observed for a patch of sinusoidal ripples having different wave numbers. It was also observed that for small angles of incidence, the reflected energy is greater compared to other angles of incidence up to π/ 4. These theoretical observations are supported by graphical results.展开更多
A previous study （Song. 2004. Geophys Res Lett, 31 （15）：L15302） of the second-order solutions for random interracial waves is extended in a constant depth, two-layer fluid system with a rigid lid is extended into...A previous study （Song. 2004. Geophys Res Lett, 31 （15）：L15302） of the second-order solutions for random interracial waves is extended in a constant depth, two-layer fluid system with a rigid lid is extended into a more general case of two-layer fluid with a top free surface. The rigid boundary condition on the upper surface is replaced by the kinematical and dynamical boundary conditions of a free surface, and the equations describing the random displacements of free surface, density-interface and the associated velocity potentials in the two-layer fluid are solved to the second order using the same expansion technology as that of Song （2004. Geophys Res Lett, 31 （15）：L15302）. The results show that the interface and the surface will oscillate synchronously, and the wave fields to the first-order both at the free surface and at the density-interface are made up of a linear superposition of many waves with different amplitudes, wave numbers and frequencies. The second-order solutions describe the second-order wave-wave interactions of the surface wave components, the interface wave components and among the surface and the interface wave components. The extended solutions also include special cases obtained by Thorpe for progressive interracial waves （Thorpe. 1968a.Trans R Soc London, 263A：563～614） and standing interracial waves （Thorpe. 1968b. J Fluid Mech, 32：489-528） for the two-layer fluid with a top free surface. Moreover, the solutions reduce to those derived for random surface waves by Sharma and Dean （1979.Ocean Engineering Rep 20） if the density of the upper layer is much smaller than that of the lower layer.展开更多
The derivation of Green function in a two-layer fluid model has been treated in different ways. In a two-layer fluid with the upper layer having a free surface, there exist two modes of waves propagating due to the fr...The derivation of Green function in a two-layer fluid model has been treated in different ways. In a two-layer fluid with the upper layer having a free surface, there exist two modes of waves propagating due to the free surface and the interface. This paper is concerned with the derivation of Green functions in the three dimensional case of a stationary source oscillating. The source point is located either in the upper or lower part of a two-layer fluid of finite depth. The derivation is carried out by the method of singularities. This method has an advantage in that it involves representing the potential as a sum of singularities or multipoles placed within any structures being present. Furthermore, experience shows that the systems of equations resulted from using a singularity method possess excellent convergence characteristics and only a few equations are needed to obtain accurate numerical results. Validation is done by showing that the derived two-layer Green function can be reduced to that of a single layer of finite depth or that the upper Green function coincides with that of the lower, for each case. The effect of the density on the internal waves is demonstrated. Also, it is shown how the surface and internal wave amplitudes are compared for both the wave modes. The fluid in this case is considered to be inviscid and incompressible and the flow is irrotational.展开更多
The modeling of turbulence,especially the high-speed compressible turbulence encountered in aerospace engineering,has always being a significant challenge in terms of balancing efficiency and accuracy.Most traditional...The modeling of turbulence,especially the high-speed compressible turbulence encountered in aerospace engineering,has always being a significant challenge in terms of balancing efficiency and accuracy.Most traditional models typically show limitations in universality,accuracy,and reliance on past experience.The stochastic multi-scale models show great potential in addressing these issues by representing turbulence across all characteristic scales in a reduced-dimensional space,maintaining sufficient accuracy while reducing computational cost.This review systematically summarizes advances in methods related to a widely used and refined stochastic multi-scale model,the One-Dimensional Turbulence(ODT).The advancements in formulations are emphasized for stand-alone incompressible ODT models,stand-alone compressible ODT models,and coupling methods.Some diagrams are also provided to facilitate more readers to understand the ODT methods.Subsequently,the significant developments and applications of stand-alone ODT models and coupling methods are introduced and critically evaluated.Despite the extensively recognized effectiveness of ODT models in low-speed turbulent flows,it is crucial to emphasize that there is still a research gap in the field of ODT coupling methods that are capable of accurately and efficiently simulating complex,three-dimensional,high-speed compressible turbulent flows up to now.Based on an analysis of the advantages and limitations of existing ODT methods,the recent advancement in the conservative compressible ODT model is considered to have provided a promising approach to tackle the modeling challenges of high-speed compressible turbulence.Therefore,this review outlines several recommended new research subjects and challenging issues to inspire further research in simulating complex,three-dimensional,high-speed compressible turbulent flows using ODT models.展开更多
Due to the novel applications of flexible pipes conveying fluid in the field of soft robotics and biomedicine,the investigations on the mechanical responses of the pipes have attracted considerable attention.The fluid...Due to the novel applications of flexible pipes conveying fluid in the field of soft robotics and biomedicine,the investigations on the mechanical responses of the pipes have attracted considerable attention.The fluid-structure interaction(FSI)between the pipe with a curved shape and the time-varying internal fluid flow brings a great challenge to the revelation of the dynamical behaviors of flexible pipes,especially when the pipe is highly flexible and usually undergoes large deformations.In this work,the geometrically exact model(GEM)for a curved cantilevered pipe conveying pulsating fluid is developed based on the extended Hamilton's principle.The stability of the curved pipe with three different subtended angles is examined with the consideration of steady fluid flow.Specific attention is concentrated on the large-deformation resonance of circular pipes conveying pulsating fluid,which is often encountered in practical engineering.By constructing bifurcation diagrams,oscillating shapes,phase portraits,time traces,and Poincarémaps,the dynamic responses of the curved pipe under various system parameters are revealed.The mean flow velocity of the pulsating fluid is chosen to be either subcritical or supercritical.The numerical results show that the curved pipe conveying pulsating fluid can exhibit rich dynamical behaviors,including periodic and quasi-periodic motions.It is also found that the preferred instability type of a cantilevered curved pipe conveying steady fluid is mainly in the flutter of the second mode.For a moderate value of the mass ratio,however,a third-mode flutter may occur,which is quite different from that of a straight pipe system.展开更多
The coefficients embodied in a Boussinesq-type model are very important since they are determined to optimize the linear and nonlinear properties.In most conventional Boussinesq-type models,these coefficients are assi...The coefficients embodied in a Boussinesq-type model are very important since they are determined to optimize the linear and nonlinear properties.In most conventional Boussinesq-type models,these coefficients are assigned the specific values.As for the multi-layer Boussinesq-type models with the inclusion of the vertical velocity,however,the effect of the different values of these coefficients on linear and nonlinear performances has never been investigated yet.The present study focuses on a two-layer Boussinesq-type model with the highest spatial derivatives being 2 and theoretically and numerically examines the effect of the coefficient on model performance.Theoretical analysis show that different values for(0.13≤α≤0.25)do not have great effects on the high accuracy of the linear shoaling,linear phase celerity and even third-order nonlinearity for water depth range of 0<kh≤10(k is wave number and h is water depth).The corresponding errors using different values are restricted within 0.1%,0.1%and 1%for the linear shoaling amplitude,dispersion and nonlinear harmonics,respectively.Numerical tests including regular wave shoaling over mildly varying slope from deep to shallow water,regular wave propagation over submerged bar,bichromatic wave group and focusing wave propagation over deep water are conducted.The comparison between numerical results using different values of,experimental data and analytical solutions confirm the theoretical analysis.The flexibility and consistency of the two-layer Boussinesq-type model is therefore demonstrated theoretically and numerically.展开更多
Internal waves are one of the various phenomena that occur at sea,and they affect acoustic equipment and sea density measurement equipment.In this study,internal waves are simulated using computational fluid dynamics ...Internal waves are one of the various phenomena that occur at sea,and they affect acoustic equipment and sea density measurement equipment.In this study,internal waves are simulated using computational fluid dynamics method in the presence of a submarine in a pre-stratified fluid.Several scenarios were implemented by Froude number changes and submersible velocity by using the Navier-Stokes k-εturbulence model.Results indicate that the realizable k-εturbulence model gives better results than the RNG k-εmodel and the internal waves flow in this model are well represented,which increases the wavelength of the internal waves by increasing the Froude number and floating velocity,while the internal angle of the Kelvin waves is decreased.We also observe that increasing the floating velocity causes the turbulent velocity contours to increase due to the drag coefficient and its relationship with the Reynolds number.The Reynolds number increases with the increasing velocity of the float motion.The results indicate the efficiency of this method in the discovery of subsurface objects.展开更多
Recent industrial explosions globally have intensified the focus in mechanical engineering on designing infras-tructure systems and networks capable of withstanding blast loading.Initially centered on high-profile fac...Recent industrial explosions globally have intensified the focus in mechanical engineering on designing infras-tructure systems and networks capable of withstanding blast loading.Initially centered on high-profile facilities such as embassies and petrochemical plants,this concern now extends to a wider array of infrastructures and facilities.Engineers and scholars increasingly prioritize structural safety against explosions,particularly to prevent disproportionate collapse and damage to nearby structures.Urbanization has further amplified the reliance on oil and gas pipelines,making them vital for urban life and prime targets for terrorist activities.Consequently,there is a growing imperative for computational engineering solutions to tackle blast loading on pipelines and mitigate associated risks to avert disasters.In this study,an empty pipe model was successfully validated under contact blast conditions using Abaqus software,a powerful tool in mechanical engineering for simulating blast effects on buried pipelines.Employing a Eulerian-Lagrangian computational fluid dynamics approach,the investigation extended to above-surface and below-surface blasts at standoff distances of 25 and 50 mm.Material descriptions in the numerical model relied on Abaqus’default mechanical models.Comparative analysis revealed varying pipe performance,with deformation decreasing as explosion-to-pipe distance increased.The explosion’s location relative to the pipe surface notably influenced deformation levels,a key finding highlighted in the study.Moreover,quantitative findings indicated varying ratios of plastic dissipation energy(PDE)for different blast scenarios compared to the contact blast(P0).Specifically,P1(25 mm subsurface blast)and P2(50 mm subsurface blast)showed approximately 24.07%and 14.77%of P0’s PDE,respectively,while P3(25 mm above-surface blast)and P4(50 mm above-surface blast)exhibited lower PDE values,accounting for about 18.08%and 9.67%of P0’s PDE,respectively.Utilising energy-absorbing materials such as thin coatings of ultra-high-strength concrete,metallic foams,carbon fiber-reinforced polymer wraps,and others on the pipeline to effectively mitigate blast damage is recommended.This research contributes to the advancement of mechanical engineering by providing insights and solutions crucial for enhancing the resilience and safety of underground pipelines in the face of blast events.展开更多
基金Fundamental Research Funds for the Central Universities(M23JBZY00050)National Natural Science Foundation of China(22278032)。
文摘The kinetic characteristics of plasma-assisted oxidative pyrolysis of ammonia are studied by using the global/fluid models hybrid solution method.Firstly,the stable products of plasma-assisted oxidative pyrolysis of ammonia are measured.The results show that the consumption of NH_(3)/O_(2)and the production of N_(2)/H_(2)change linearly with the increase of voltage,which indicates the decoupling of nonequilibrium molecular excitation and oxidative pyrolysis of ammonia at low temperatures.Secondly,the detailed reaction kinetics mechanism of ammonia oxidative pyrolysis stimulated by a nanosecond pulse voltage at low pressure and room temperature is established.Based on the reaction path analysis,the simplified mechanism is obtained.The detailed and simplified mechanism simulation results are compared with experimental data to verify the accuracy of the simplified mechanism.Finally,based on the simplified mechanism,the fluid model of ammonia oxidative pyrolysis stimulated by the nanosecond pulse plasma is established to study the pre-sheath/sheath behavior and the resultant consumption and formation of key species.The results show that the generation,development,and propagation of the pre-sheath have a great influence on the formation and consumption of species.The consumption of NH_(3)by the cathode pre-sheath is greater than that by the anode pre-sheath,but the opposite is true for OH and O(1S).However,within the sheath,almost all reactions do not occur.Further,by changing the parameters of nanosecond pulse power supply voltage,it is found that the electron number density,electron current density,and applied peak voltages are not the direct reasons for the structural changes of the sheath and pre-sheath.Furthermore,the discharge interval has little effect on the sheath structure and gas mixture breakdown.The research results of this paper not only help to understand the kinetic promotion of non-equilibrium excitation in the process of oxidative pyrolysis but also help to explore the influence of transport and chemical reaction kinetics on the oxidative pyrolysis of ammonia.
基金supported by the National Natural Science Foundation of China(Grant Nos.52306126,22350710788,12432010,11988102,92270203)the Xplore Prize.
文摘Configuring computational fluid dynamics(CFD)simulations typically demands extensive domain expertise,limiting broader access.Although large language models(LLMs)have advanced scientific computing,their use in automating CFD workflows is underdeveloped.We introduce a novel approach centered on domain-specific LLM adaptation.By fine-tuning Qwen2.5-7B-Instruct on NL2FOAM,our custom dataset of 28,716 natural language-to-OpenFOAM configuration pairs with chain-of-thought(CoT)annotations enables direct translation from natural language descriptions to executable CFD setups.A multi-agent system orchestrates the process,autonomously verifying inputs,generating configurations,running simulations,and correcting errors.Evaluation on a benchmark of 21 diverse flow cases demonstrates state-of-the-art performance,achieving 88.7%solution accuracy and 82.6%first-attempt success rate.This significantly outperforms larger general-purpose models such as Qwen2.5-72B-Instruct,DeepSeek-R1,and Llama3.3-70B-Instruct,while also requiring fewer correction iterations and maintaining high computational efficiency.The results highlight the critical role of domain-specific adaptation in deploying LLM assistants for complex engineering workflows.Our code and fine-tuned model have been deposited at https://github.com/YYgroup/AutoCFD.
基金the Deanship of Research and Graduate Studies at King Khalid University for funding this work through large Research Group Project (Grant No. RGP2/198/45)Project supported by Prince Sattam bin Abdulaziz University (Grant No. PSAU/2025/R/1446)。
文摘This work investigates water-based micropolar hybrid nanofluid(MHNF) flow on an elongating variable porous sheet.Nanoparticles of diamond and copper have been used in the water to boost its thermal conductivity. The motion of the fluid is taken as two-dimensional with the impact of a magnetic field in the normal direction. The variable, permeable, and stretching nature of sheet's surface sets the fluid into motion. Thermal and mass diffusions are controlled through the use of the Cattaneo–Christov flux model. A dataset is generated using MATLAB bvp4c package solver and employed to train an artificial neural network(ANN) based on the Levenberg–Marquardt back-propagation(LMBP) algorithm. It has been observed as an outcome of this study that the modeled problem achieves peak performance at epochs 637, 112, 4848, and 344 using ANN-LMBP. The linear velocity of the fluid weakens with progression in variable porous and magnetic factors.With an augmentation in magnetic factor, the micro-rotational velocity profiles are augmented on the domain 0 ≤ η < 1.5 due to the support of micro-rotations by Lorentz forces close to the sheet's surface, while they are suppressed on the domain 1.5 ≤ η < 6.0 due to opposing micro-rotations away from the sheet's surface. Thermal distributions are augmented with an upsurge in thermophoresis, Brownian motion, magnetic, and radiation factors, while they are suppressed with an upsurge in thermal relaxation parameter. Concentration profiles increase with an expansion in thermophoresis factor and are suppressed with an intensification of Brownian motion factor and solute relaxation factor. The absolute errors(AEs) are evaluated for all the four scenarios that fall within the range 10^(-3)–10^(-8) and are associated with the corresponding ANN configuration that demonstrates a fine degree of accuracy.
基金supported by the National Natural Science Foundation of China Basic Science Center Program for“Multiscale Problems in Nonlinear Mechanics”(Grant No.11988102)the National Natural Science Foundation of China(Grant No.12202451).
文摘This paper investigates the capabilities of large language models(LLMs)to leverage,learn and create knowledge in solving computational fluid dynamics(CFD)problems through three categories of baseline problems.These categories include(1)conventional CFD problems that can be solved using existing numerical methods in LLMs,such as lid-driven cavity flow and the Sod shock tube problem;(2)problems that require new numerical methods beyond those available in LLMs,such as the recently developed Chien-physics-informed neural networks for singularly perturbed convection-diffusion equations;and(3)problems that cannot be solved using existing numerical methods in LLMs,such as the ill-conditioned Hilbert linear algebraic systems.The evaluations indicate that reasoning LLMs overall outperform non-reasoning models in four test cases.Reasoning LLMs show excellent performance for CFD problems according to the tailored prompts,but their current capability in autonomous knowledge exploration and creation needs to be enhanced.
基金Supported by the National Natural Science Foundation of China(Grant No.40674063)National Hi-tech Research and Development Program of China(863Program)(Grant No.2006AA09Z311)
文摘Based on the synchronous joint gravity and magnetic inversion of single interface by Pilkington and the need of revealing Cenozoic and crystalline basement thickness in the new round of oil-gas exploration, we propose a joint gravity and magnetic inversion methodfor two-layer models by concentrating on the relationship between the change of thicknessI and position of the middle layer and anomaly and discuss the effects of the key parameters. Model tests and application to field data show the validity of this method.
基金support from the National Natural Science Foundation of China(U22A20171)the High Steel Center(HSC)at North China University of Technology and University of Science and Technology Beijing,China.
文摘A 1∶8 physical water model was constructed to investigate the fluid flow and mixing phenomena in the basic oxygen furnace(BOF)converter.The particle image velocimetry was employed to measure the velocity distribution of the bath and the high-speed camera was applied to capture the cavity shape in the combined blowing BOF converter.The mixing time for varied operating conditions was measured by the stimulus-response approach.The cavity depth increased with the decrease in the lance height and the increase in the top gas flow rate while the bottom blowing gas had little influence on the cavity depth.The minimum cavity depth was obtained under the condition of a 69.8 m^(3)/h top gas flow rate,a 287.5 mm lance height and a 0.93 m^(3)/h bottom blowing gas flow rate,which was 161.2 mm.The mixing time decreased as the lance height decreased and the top blowing gas flow rate increased.The mixing time was first decreased and then increased with the increase in the bottom gas flow rate.With the condition of 69.8 m^(3)/h gas flow rate of top blowing,the 287.5 mm lance height and the 0.93 m^(3)/h gas flow rate of bottom blowing,the mixing time in the converter was 48.65 s.The empirical formula between the stirring power and the mixing time in the converter was calculated.
基金sponsored by the National Basic Research Program of China(973 Program,Grant No.2014CB046203)the National Natural Science Foundation of China(Grant No.11072140)
文摘The hydroelastic response of a circular, very large floating structure(VLFS), idealized as a floating circular elastic thin plate, is investigated for the case of time-harmonic incident waves of the surface and interfacial wave modes, of a given wave frequency, on a two-layer fluid of finite and constant depth. In linear potential-flow theory, with the aid of angular eigenfunction expansions, the diffraction potentials can be expressed by the Bessel functions. A system of simultaneous equations is derived by matching the velocity and the pressure between the open-water and the platecovered regions, while incorporating the edge conditions of the plate. Then the complex nested series are simplified by utilizing the orthogonality of the vertical eigenfunctions in the open-water region. Numerical computations are presented to investigate the effects of different physical quantities, such as the thickness of the plate, Young’s modulus, the ratios of the densities and of the layer depths, on the dispersion relations of the flexural-gravity waves for the two-layer fluid. Rapid convergence of the method is observed, but is slower at higher wave frequency. At high frequency, it is found that there is some energy transferred from the interfacial mode to the surface mode.
基金Project supported by the National Basic Research Program of China (Grant No. 2011CB403501)the National Natural Science Foundation of China (GrantNos. 41175058,41275062,and 11202106)
文摘A weak nonlinear model of a two-layer barotropic ocean with Rayleigh dissipation is built.The analytic asymptotic solution is derived in the mid-latitude stationary wind field,and the physical meaning of the corresponding problem is discussed.
基金the Knowledge Innovation Project of CAS(KJCX-YW-L02)the National 863 Project of China(2006AAO9A103-4)+1 种基金China National Oil Corporation in Beijing(CNOOC)the National Natural Science Foundation of China(10672056).
文摘In the present paper, we endeavor to accomplish a diagram, which demarcates the validity ranges for interfacial wave theories in a two-layer system, to meet the needs of design in ocean engineering. On the basis of the available solutions of periodic and solitary waves, we propose a guideline as principle to identify the validity regions of the interfacial wave theories in terms of wave period T, wave height H, upper layer thickness dl, and lower layer thick-ness d2, instead of only one parameter-water depth d as in the water surface wave circumstance. The diagram proposed here happens to be Le Mehaute＇s plot for free surface waves if water depth ratio r= d1/d2 approaches to infinity and the upper layer water density p1 to zero. On the contrary, the diagram for water surface waves can be used for two-layer interfacial waves if gravity acceleration g in it is replaced by the reduced gravity defined in this study under the condition of σ=（P2 - Pl）/P2 → 1.0 and r 〉 1.0. In the end, several figures of the validity ranges for various interfacial wavetheories in the two-layer fluid are given and compared with the results for surface waves.
基金This project was financially supported by the National Natural Science Foundation of China(Grant No.10572092)the High Technology Research and Development Programof China(863Program,Grant Nos.2006AA09Z352 and 2006AA09196-6)
文摘In this paper, the diffraction of water waves by a vertically floating cylinder in a two-layer fluid of a finite depth is studied. Analytical expressions for the hydrodynamic loads on the vertically floating cylinder are obtained by use of the method of eigenfunction expansions. The hydrodynamic loads on the vertically floating cylinder in a two-layer fluid inelude not only the surge, heave and pitch exciting forces due to the incident wave of the surface-wave mode, but also those due to the incident wave of the internal-wave mode. This is different from the case of a homogenous fluid. Some given examples show that, for a two-layer fluid system with a small density difference, the hydrodynamic loads for the surface-wave mode do not differ significantly from those due to surface waves in a single-layer fluid, but the hydrodynamic loads for the internal-wave mode are important over a wide range of frequencies. Moreover, also considered are the free surface and interface elevations generated by the diffraction wave due to the incident wave of the surface-wave and interhal-wave modes, and transfer of energy between modes.
文摘Many new forms of Boussinesq-type equations have been developed to extend the range of applicability of the classical Boussinesq equations to deeper water in the Study of the surface waves. One approach was used by Nwogu （1993. J. Wtrw. Port Coastal and Oc. Eng. 119, 618-638） to improve the linear dispersion characteristics of the classical Boussinesq equations by using the velocity at an arbitrary level as the velocity variable in derived equations and obtain a new form of Boussinesq-type equations, in which the dispersion property can be optimized by choosing the velocity variable at an adequate level. In this paper, a set of Boussinesq-type equations describing the motions of the interracial waves propagating alone the interface between two homogeneous incompressible and inviscid fluids of different densities with a free surface and a variable water depth were derived using a method similar to that used by Nwogu （1993. J. Wtrw. Port Coastal and Oc. Eng. 119, 618-638） for surface waves. The equations were expressed in terms of the displacements of free surface and density-interface, and the velocity vectors at arbitrary vertical locations in the upper layer and the lower layer （or depth-averaged velocity vector across each layer） of a two-layer fluid. As expected, the equations derived in the present work include as special cases those obtained by Nwogu （1993, J. Wtrw. Port Coastal and Oc. Eng. 119, 618-638） and Peregrine （1967, J. Fluid Mech. 27, 815-827） for surface waves when the density of the upper fluid is taken as zero.
基金The project supported by National Natural Science Foundation of China under Grant Nos.10547124,10475055,and 90503006the Youth Foundation of Shanghai Jiao Tong University
文摘A quite general coupled variable coefficient modified KdV (VCmKdV) equation in a two-layer fluid systemis derived by means of the reductive perturbation method.Making use of the CK's direct method,some similarityreductions of the coupled VCmKdV equation are obtained and their corresponding group explanations are discussed.Some exact solutions of the coupled equations are also presented.
文摘The problem of oblique wave （internal wave） propagation over a small deformation in a channel flow consisting of two layers was considered. The upper fluid was assumed to be bounded above by a rigid lid, which is an approximation for the free surface, and the lower one was bounded below by an impermeable bottom surface having a small deformation; the channel was unbounded in the horizontal directions. Assuming irrotational motion, the perturbation technique was employed to calculate the first-order corrections of the velocity potential in the two fluids by using Green＇s integral theorem suitably with the introduction of appropriate Green＇s functions. Those functions help in calculating the reflection and transmission coefficients in terms of integrals involving the shape ftmction c（x） representing the bottom deformation. Three-dimensional linear water wave theory was utilized for formulating the relevant boundary value problem. Two special examples of bottom deformation were considered to validate the results. Consideration of a patch of sinusoidal ripples （having the same wave number） shows that the reflection coefficient is an oscillatory function of the ratio of twice the x-component of the wave number to the ripple wave number. When this ratio approaches one, the theory predicts a resonant interaction between the bed and the interface, and the reflection coefficient becomes a multiple of the number of ripples. High reflection of incident wave energy occurs if this number is large. Similar results were observed for a patch of sinusoidal ripples having different wave numbers. It was also observed that for small angles of incidence, the reflected energy is greater compared to other angles of incidence up to π/ 4. These theoretical observations are supported by graphical results.
基金supported by the National Science Foundation for Distinguished Young Scholars of China under contract No.40425015the Cooperative Project of Chinese Academy Sciences and the China National 0ffshore 0il Corporation("Behaviours of internal waves and their roles on the marine stuctures").
文摘A previous study （Song. 2004. Geophys Res Lett, 31 （15）：L15302） of the second-order solutions for random interracial waves is extended in a constant depth, two-layer fluid system with a rigid lid is extended into a more general case of two-layer fluid with a top free surface. The rigid boundary condition on the upper surface is replaced by the kinematical and dynamical boundary conditions of a free surface, and the equations describing the random displacements of free surface, density-interface and the associated velocity potentials in the two-layer fluid are solved to the second order using the same expansion technology as that of Song （2004. Geophys Res Lett, 31 （15）：L15302）. The results show that the interface and the surface will oscillate synchronously, and the wave fields to the first-order both at the free surface and at the density-interface are made up of a linear superposition of many waves with different amplitudes, wave numbers and frequencies. The second-order solutions describe the second-order wave-wave interactions of the surface wave components, the interface wave components and among the surface and the interface wave components. The extended solutions also include special cases obtained by Thorpe for progressive interracial waves （Thorpe. 1968a.Trans R Soc London, 263A：563～614） and standing interracial waves （Thorpe. 1968b. J Fluid Mech, 32：489-528） for the two-layer fluid with a top free surface. Moreover, the solutions reduce to those derived for random surface waves by Sharma and Dean （1979.Ocean Engineering Rep 20） if the density of the upper layer is much smaller than that of the lower layer.
基金supported by the National Natural Science Foundation of China (Grant No. 50779008)
文摘The derivation of Green function in a two-layer fluid model has been treated in different ways. In a two-layer fluid with the upper layer having a free surface, there exist two modes of waves propagating due to the free surface and the interface. This paper is concerned with the derivation of Green functions in the three dimensional case of a stationary source oscillating. The source point is located either in the upper or lower part of a two-layer fluid of finite depth. The derivation is carried out by the method of singularities. This method has an advantage in that it involves representing the potential as a sum of singularities or multipoles placed within any structures being present. Furthermore, experience shows that the systems of equations resulted from using a singularity method possess excellent convergence characteristics and only a few equations are needed to obtain accurate numerical results. Validation is done by showing that the derived two-layer Green function can be reduced to that of a single layer of finite depth or that the upper Green function coincides with that of the lower, for each case. The effect of the density on the internal waves is demonstrated. Also, it is shown how the surface and internal wave amplitudes are compared for both the wave modes. The fluid in this case is considered to be inviscid and incompressible and the flow is irrotational.
基金cosupported by the National Natural Science Foundation of China(No.12202487)。
文摘The modeling of turbulence,especially the high-speed compressible turbulence encountered in aerospace engineering,has always being a significant challenge in terms of balancing efficiency and accuracy.Most traditional models typically show limitations in universality,accuracy,and reliance on past experience.The stochastic multi-scale models show great potential in addressing these issues by representing turbulence across all characteristic scales in a reduced-dimensional space,maintaining sufficient accuracy while reducing computational cost.This review systematically summarizes advances in methods related to a widely used and refined stochastic multi-scale model,the One-Dimensional Turbulence(ODT).The advancements in formulations are emphasized for stand-alone incompressible ODT models,stand-alone compressible ODT models,and coupling methods.Some diagrams are also provided to facilitate more readers to understand the ODT methods.Subsequently,the significant developments and applications of stand-alone ODT models and coupling methods are introduced and critically evaluated.Despite the extensively recognized effectiveness of ODT models in low-speed turbulent flows,it is crucial to emphasize that there is still a research gap in the field of ODT coupling methods that are capable of accurately and efficiently simulating complex,three-dimensional,high-speed compressible turbulent flows up to now.Based on an analysis of the advantages and limitations of existing ODT methods,the recent advancement in the conservative compressible ODT model is considered to have provided a promising approach to tackle the modeling challenges of high-speed compressible turbulence.Therefore,this review outlines several recommended new research subjects and challenging issues to inspire further research in simulating complex,three-dimensional,high-speed compressible turbulent flows using ODT models.
基金Project supported by the National Natural Science Foundation of China (Nos.12072119,12325201,and 52205594)the China National Postdoctoral Program for Innovative Talents (No.BX20220118)。
文摘Due to the novel applications of flexible pipes conveying fluid in the field of soft robotics and biomedicine,the investigations on the mechanical responses of the pipes have attracted considerable attention.The fluid-structure interaction(FSI)between the pipe with a curved shape and the time-varying internal fluid flow brings a great challenge to the revelation of the dynamical behaviors of flexible pipes,especially when the pipe is highly flexible and usually undergoes large deformations.In this work,the geometrically exact model(GEM)for a curved cantilevered pipe conveying pulsating fluid is developed based on the extended Hamilton's principle.The stability of the curved pipe with three different subtended angles is examined with the consideration of steady fluid flow.Specific attention is concentrated on the large-deformation resonance of circular pipes conveying pulsating fluid,which is often encountered in practical engineering.By constructing bifurcation diagrams,oscillating shapes,phase portraits,time traces,and Poincarémaps,the dynamic responses of the curved pipe under various system parameters are revealed.The mean flow velocity of the pulsating fluid is chosen to be either subcritical or supercritical.The numerical results show that the curved pipe conveying pulsating fluid can exhibit rich dynamical behaviors,including periodic and quasi-periodic motions.It is also found that the preferred instability type of a cantilevered curved pipe conveying steady fluid is mainly in the flutter of the second mode.For a moderate value of the mass ratio,however,a third-mode flutter may occur,which is quite different from that of a straight pipe system.
基金supported by the National Natural Science Foundation of China(Grant Nos.51779022,51809053,and 51579034)the Innovation Team Project of Estuary and Coast Protection and Management(Grant No.Y220013)the Open Project Fund of State Key Laboratory of Coastal and Offshore Engineering,Dalian University of Technology(Grant No.LP19015).
文摘The coefficients embodied in a Boussinesq-type model are very important since they are determined to optimize the linear and nonlinear properties.In most conventional Boussinesq-type models,these coefficients are assigned the specific values.As for the multi-layer Boussinesq-type models with the inclusion of the vertical velocity,however,the effect of the different values of these coefficients on linear and nonlinear performances has never been investigated yet.The present study focuses on a two-layer Boussinesq-type model with the highest spatial derivatives being 2 and theoretically and numerically examines the effect of the coefficient on model performance.Theoretical analysis show that different values for(0.13≤α≤0.25)do not have great effects on the high accuracy of the linear shoaling,linear phase celerity and even third-order nonlinearity for water depth range of 0<kh≤10(k is wave number and h is water depth).The corresponding errors using different values are restricted within 0.1%,0.1%and 1%for the linear shoaling amplitude,dispersion and nonlinear harmonics,respectively.Numerical tests including regular wave shoaling over mildly varying slope from deep to shallow water,regular wave propagation over submerged bar,bichromatic wave group and focusing wave propagation over deep water are conducted.The comparison between numerical results using different values of,experimental data and analytical solutions confirm the theoretical analysis.The flexibility and consistency of the two-layer Boussinesq-type model is therefore demonstrated theoretically and numerically.
文摘Internal waves are one of the various phenomena that occur at sea,and they affect acoustic equipment and sea density measurement equipment.In this study,internal waves are simulated using computational fluid dynamics method in the presence of a submarine in a pre-stratified fluid.Several scenarios were implemented by Froude number changes and submersible velocity by using the Navier-Stokes k-εturbulence model.Results indicate that the realizable k-εturbulence model gives better results than the RNG k-εmodel and the internal waves flow in this model are well represented,which increases the wavelength of the internal waves by increasing the Froude number and floating velocity,while the internal angle of the Kelvin waves is decreased.We also observe that increasing the floating velocity causes the turbulent velocity contours to increase due to the drag coefficient and its relationship with the Reynolds number.The Reynolds number increases with the increasing velocity of the float motion.The results indicate the efficiency of this method in the discovery of subsurface objects.
文摘Recent industrial explosions globally have intensified the focus in mechanical engineering on designing infras-tructure systems and networks capable of withstanding blast loading.Initially centered on high-profile facilities such as embassies and petrochemical plants,this concern now extends to a wider array of infrastructures and facilities.Engineers and scholars increasingly prioritize structural safety against explosions,particularly to prevent disproportionate collapse and damage to nearby structures.Urbanization has further amplified the reliance on oil and gas pipelines,making them vital for urban life and prime targets for terrorist activities.Consequently,there is a growing imperative for computational engineering solutions to tackle blast loading on pipelines and mitigate associated risks to avert disasters.In this study,an empty pipe model was successfully validated under contact blast conditions using Abaqus software,a powerful tool in mechanical engineering for simulating blast effects on buried pipelines.Employing a Eulerian-Lagrangian computational fluid dynamics approach,the investigation extended to above-surface and below-surface blasts at standoff distances of 25 and 50 mm.Material descriptions in the numerical model relied on Abaqus’default mechanical models.Comparative analysis revealed varying pipe performance,with deformation decreasing as explosion-to-pipe distance increased.The explosion’s location relative to the pipe surface notably influenced deformation levels,a key finding highlighted in the study.Moreover,quantitative findings indicated varying ratios of plastic dissipation energy(PDE)for different blast scenarios compared to the contact blast(P0).Specifically,P1(25 mm subsurface blast)and P2(50 mm subsurface blast)showed approximately 24.07%and 14.77%of P0’s PDE,respectively,while P3(25 mm above-surface blast)and P4(50 mm above-surface blast)exhibited lower PDE values,accounting for about 18.08%and 9.67%of P0’s PDE,respectively.Utilising energy-absorbing materials such as thin coatings of ultra-high-strength concrete,metallic foams,carbon fiber-reinforced polymer wraps,and others on the pipeline to effectively mitigate blast damage is recommended.This research contributes to the advancement of mechanical engineering by providing insights and solutions crucial for enhancing the resilience and safety of underground pipelines in the face of blast events.
