Computational Fluid Dynamics Approach for Predicting Pipeline Response to Various Blast Scenarios: A Numerical Modeling Study

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.

A simplified approach to modelling blasts in computational fluid dynamics (CFD)

This paper presents a time-efficient numerical approach to modelling high explosive(HE)blastwave propagation using Computational Fluid Dynamics(CFD).One of the main issues of using conventional CFD modelling in high explosive simulations is the ability to accurately define the initial blastwave properties that arise from the ignition and consequent explosion.Specialised codes often employ Jones-Wilkins-Lee(JWL)or similar equation of state(EOS)to simulate blasts.However,most available CFD codes are limited in terms of EOS modelling.They are restrictive to the Ideal Gas Law(IGL)for compressible flows,which is generally unsuitable for blast simulations.To this end,this paper presents a numerical approach to simulate blastwave propagation for any generic CFD code using the IGL EOS.A new method known as the Input Cavity Method(ICM)is defined where input conditions of the high explosives are given in the form of pressure,velocity and temperature time-history curves.These time history curves are input at a certain distance from the centre of the charge.It is shown that the ICM numerical method can accurately predict over-pressure and impulse time history at measured locations for the incident,reflective and complex multiple reflection scenarios with high numerical accuracy compared to experimental measurements.The ICM is compared to the Pressure Bubble Method(PBM),a common approach to replicating initial conditions for a high explosive in Finite Volume modelling.It is shown that the ICM outperforms the PBM on multiple fronts,such as peak values and overall overpressure curve shape.Finally,the paper also presents the importance of choosing an appropriate solver between the Pressure Based Solver(PBS)and Density-Based Solver(DBS)and provides the advantages and disadvantages of either choice.In general,it is shown that the PBS can resolve and capture the interactions of blastwaves to a higher degree of resolution than the DBS.This is achieved at a much higher computational cost,showing that the DBS is much preferred for quick turnarounds.

Computational Fluid Dynamics (CFD) Modelling to Estimate Fluvial Bank Erosion—A Case Study

River bank erosion models are an important prerequisite for understanding the development of river meanders and for estimating likely land-loss and potential danger to floodplain infrastructure. Although bank erosion models have been developed that consider large-scale mass failure, the contribution of fluvial erosion (the process of particle-by-particle erosion due to the shearing action of the river flow) to bank retreat has not received as much consideration. In principle, such fluvial bank erosion rates can be quantified using excess shear stress formulations, but in practice, it has proven difficult to estimate the parameters involved. In this study, a series of three-dimensional Computational Fluid Dynamics (CFD) simulations for a meander loop on the River Asker (200 m long) at Bridport in southern England were undertaken to elucidate the overall flow structures and in particular to provide estimates of the applied fluid shear stress exerted on the riverbanks. The CFD models, which simulated relatively low and relatively high flow conditions, were established using Fluent 6.2 software. The modelling outcomes show that the key qualitative features of the flow endure even as flow discharge varies. At bank full, the degrees of velocity and simulated shear stresses within the inner bank separation zones are shown to be higher than those observed under low flow conditions, and that these elevated shear stresses may be sufficient to result in the removal of accumulated sediments into the main downstream flow.

Numerical Study on the Aerodynamic and Fluid−Structure Interaction of An NREL-5MW Wind Turbine

A 5-MW wind turbine has been modeled and analyzed for fluid-structure interaction and aerodynamic performance.In this study, a full-scale model of a 5-MW wind turbine is first developed based on a computational fluid dynamics(CFD) approach, in which the unsteady, noncompressible Reynolds Averaged Navier-Stokes(RANS) method is used. The main focus of the study is to analyze the tower shadow effect on the aerodynamic performance of the wind turbine under different inlet flow conditions. Subsequently, the finite element model is established by considering fluid/structure interactions to study the structural stress, displacement, strain distributions and flow field information of the structure under the uniform wind speed. Finally, the fluid-structure interaction model is established by considering turbulent wind and the tower shadow effect. The variation rules of the dynamic response of the one-way and two-way fluid-structure interaction(FSI) models under different wind speeds are analyzed, and the numerical calculation results are compared with those of the centralized mass model. The results show that the tower shadow effect and structural deformation are the main factors affecting the aerodynamic load fluctuation of the wind turbine, which in turn affects the aerodynamic performance and structural stability of the blades. The structural dynamic response of the coupled model shows significant similarity, while the structural displacement response of the former exhibits less fluctuation compared with the conventional centralized mass model. The one-way fluid-structure interaction(FSI)model shows a higher frequency of stress-strain and displacement oscillations on the blade compared with the two-way FSI model.

Numerical Solutions of the Classical and Modified Buckley-Leverett Equations Applied to Two-Phase Fluid Flow

Climate change is a reality. The burning of fossil fuels from oil, natural gas and coal is responsible for much of the pollution and the increase in the planet’s average temperature, which has raised discussions on the subject, given the emergencies related to climate. An energy transition to clean and renewable sources is necessary and urgent, but it will not be quick. In this sense, increasing the efficiency of oil extraction from existing sources is crucial, to avoid waste and the drilling of new wells. The purpose of this work was to add diffusive and dispersive terms to the Buckley-Leverett equation in order to incorporate extra phenomena in the temporal evolution between the water-oil and oil-water transitions in the pipeline. For this, the modified Buckley-Leverett equation was discretized via essentially weighted non-oscillatory schemes, coupled with a three-stage Runge-Kutta and a fourth-order centered finite difference methods. Then, computational simulations were performed and the results showed that new features emerge in the transitions, when compared to classical simulations. For instance, the dispersive term inhibits the diffusive term, adding oscillations, which indicates that the absorption of the fluid by the porous medium occurs in a non-homogeneous manner. Therefore, based on research such as this, decisions can be made regarding the replacement of the porous medium or the insertion of new components to delay the replacement.

Numerical Modelling of Ore-forming Dynamics of Fractal Dispersive Fluid Systems

Based on an analysis of the fractal structures and mass transport mechanism of typical shear-fluid-ore formation system, the fractal dispersion theory of the fluid system was used in the dynamic study of the ore formation system. The model of point-source diffusive illuviation of the shear-fluid-ore formation system was constructed, and the numerical simulation of dynamics of the ore formation system was finished. The result shows that: (1) The metallogenic system have nested fractal structure. Different fractal dimension values in different systems show unbalance and inhomogeneity of ore-forming processes in the geohistory. It is an important parameter to symbolize the process of remobilization and accumulation of ore-forming materials. Also it can indicate the dynamics of the metallogenic system quantitatively to some extent. (2) In essence, the fractal dispersive ore-forming dynamics is a combination of multi-processes dominated by fluid dynamics and supplemented by molecule dispersion in fluids and fluid-rock interaction. It changes components and physico-chemical properties of primary rocks and fluids, favouring deposition and mineralization of ore-forming materials. (3) Gold ore-forming processes in different types of shear zones are quite different. (1) In a metallogenic system with inhomogeneous volumetric change and inhomogeneous shear, mineralization occurs in structural barriers in the centre of a shear zone and in geochemical barriers in the shear zone near its boundaries. But there is little possibility of mineralization out of the shear zone. (2) As to a metallogenic system with inhomogeneous volumetric change and simple shear, mineralization may occur only in structural barriers near the centre of the shear zone. (3) In a metallogenic system with homogeneous volumetric change and inhomogeneous shear, mineralization may occur in geochemical barriers both within and out of the shear zone.

A Computational Fluid Dynamics (CFD) Analysis of an Undulatory Mechanical Fin Driven by Shape Memory Alloy

Many fishes use undulatory fin to propel themselves in the underwater environment. These locomotor mechanisms have a popular interest to many researchers. In the present study, we perform a three-dimensional unsteady computation of an undulatory mechanical fin that is driven by Shape Memory Alloy （SMA）. The objective of the computation is to investigate the fluid dynamics of force production associated with the undulatory mechanical fin. An unstructured, grid-based, unsteady Navier-Stokes solver with automatic adaptive remeshing is used to compute the unsteady flow around the fin through five complete cycles. The pressure distribution on fin surface is computed and integrated to provide fin forces which are decomposed into lift and thrust. The velocity field is also computed throughout the swimming cycle. Finally, a comparison is conducted to reveal the dynamics of force generation according to the kinematic parameters of the undulatory fin （amplitude, frequency and wavelength）.

Modeling of gas-solid flow in a CFB riser based on computational particle fluid dynamics

A three-dimensional model for gas-solid flow in a circulating fluidized bed（CFB） riser was developed based on computational particle fluid dynamics（CPFD）.The model was used to simulate the gas-solid flow behavior inside a circulating fluidized bed riser operating at various superficial gas velocities and solids mass fluxes in two fluidization regimes,a dilute phase transport（DPT） regime and a fast fluidization（FF） regime.The simulation results were evaluated based on comparison with experimental data of solids velocity and holdup,obtained from non-invasive automated radioactive particle tracking and gamma-ray tomography techniques,respectively.The agreement of the predicted solids velocity and holdup with experimental data validated the CPFD model for the CFB riser.The model predicted the main features of the gas-solid flows in the two regimes;the uniform dilute phase in the DPT regime,and the coexistence of the dilute phase in the upper region and the dense phase in the lower region in the FF regime.The clustering and solids back mixing in the FF regime were stronger than those in the DPT regime.

COMPUTATIONAL FLUID DYNAMICS(CFD) SIMULATIONS OF DRAG REDUCTION WITH PERIODIC MICRO-STRUCTURED WALL

Computational fluid dynamics（CFD） simulations are adopted to investigate rectangular microchannel flows with various periodic micro-structured wall by introducing velocity slip boundary condition at low Reynolds number. The purpose of the current study is to numerically find out the effects of periodic micro-structured wall on the flow resistance in rectangular microchannel with the different spacings between microridges ranging from 15 to 60 pm. The simulative results indicate that pressure drop with different spacing between microridges increases linearly with flow velocity and decreases monotonically with slip velocity; Pressure drop reduction also increases with the spacing between microridges at the same condition of slip velocity and flow velocity. The results of numerical simulation are compared with theoretical predictions and experimental results in the literatures. It is found that there is qualitative agreement between them.

Computational fluid dynamics simulation of gas-liquid two phases flow in 320 m^3 air-blowing mechanical flotation cell using different turbulence models

According to the recently developed single-trough floating machine with the world's largest volume(inflatable mechanical agitation flotation machine with volume of 320 m3) in China, the gas-fluid two-phase flow in flotation cell was simulated using computational fluid dynamics method. It is shown that hexahedral mesh scheme is more suitable for the complex structure of the flotation cell than tetrahedral mesh scheme, and a mesh quality ranging from 0.7 to 1.0 is obtained. Comparative studies of the standard k-ε, k-ω and realizable k-ε turbulence models were carried out. It is indicated that the standard k-ε turbulence model could give a result relatively close to the practice and the liquid phase flow field is well characterized. In addition, two obvious recirculation zones are formed in the mixing zones, and the pressure on the rotor and stator is well characterized. Furthermore, the simulation results using improved standard k-ε turbulence model show that surface tension coefficient of 0.072, drag model of Grace and coefficient of 4, and lift coefficient of 0.001 can be achieved. The research results suggest that gas-fluid two-phase flow in large flotation cell can be well simulated using computational fluid dynamics method.

Optimization study of a PEM fuel cell performance using 3D multi-phase computational fluid dynamics model

An optimization study using a comprehensive 3D, multi-phase, non-isothermal model of a PEM （proton exchange membrane） fuel cell that incorporates significant physical processes and key parameters affecting fuel cell performance is presented and discussed in detail. The model accounts for both gas and liquid phase in the same computational domain, and thus allows for the implementation of phase change inside the gas diffusion layers. The model includes the transport of gaseous species, liquid water, protons, energy, and water dissolved in the ion-conducting polymer. Water is assumed to be exchanged among three phases： liquid, vapottr, and dissolved, with equilibrium among these phases being assumed. This model also takes into account convection and diffusion of different species in the channels as well as in the porous gas diffusion layer, heat transfer in the solids as well as in the gases, and electrochemical reactions. The results showed that the present multi-phase model is capable of identifying important parameters for the wetting behaviour of the gas diffusion layers and can be used to identify conditions that might lead to the onset of pore plugging, which has a detrimental effect on the fuel cell performance. This model is used to study the effects of several operating, design, and material parameters on fuel cell performance. Detailed analyses of the fuel cell performance under various operating conditions have been conducted and examined.

Flow Field and Temperature Field in GaN-MOCVD Reactor Based on Computational Fluid Dynamics Modeling

Metal organic chenlical vapor deposition （AIOCVD） growth systems arc one of the. main types of equipment used for growing single crystal materials, such as GaN. To obtain fihn epitaxial materials with uniform performanee, the flow field and ternperature field in a GaN-MOCVD reactor are investigated by modeling and simulating. To make the simulation results more consistent with the actual situation, the gases in the reactor are considered to be compressible, making it possible to investigate the distributions of gas density and pressure in the reactor. The computational fluid dynamics method is used to stud,v the effects of inlet gas flow velocity, pressure in the reactor, rotational speed of graphite susceptor, and gases used in the growth, which has great guiding~ significance for the growth of GaN fihn materials.

NUMERICAL SIMULATION BY COMPUTATIONAL FLUID DYNAMICS AND EXPERIMENTAL STUDY ON STIRRED BIOREACTOR WITH PUNCHED IMPELLER

Instantaneous flow field and temperature field of the two-phase fluid are measured by particle image velocimetry （PIV） and steady state method during the state of onflow. A turbulent two-phase fluid model of stirred bioreactor with punched impeller is established by the computational fluid dynamics （CFD）, using a rotating coordinate system and sliding mesh to describe the relative motion between impeller and baffles. The simulation and experiment results of flow and temperature field prove their warps are less than 10% and the mathematic model can well simulate the fields, which will also provide the study on optimized-design and scale-up of bioreactors with reference value.

Computational fluid dynamics modeling in intracranial atherosclerotic disease

Intracranial atherosclerotic disease(ICAD)is an important cause for ischemic stroke and transient ischemic stroke(TIA)throughout the world,especially in Asians,which is not fully appreciated,partly due to its inaccessibility and limitations of current neuroimaging methods.The computational fluid dynamics(CFD)modeling technique provides a novel approach to reveal the hemodynamic characteristics in ICAD,e.g.,the distributions of pressure,wall shear stress and flow velocity.In this review article,we aim to provide an overview of the general methodology of CFD modeling in arterial stenotic diseases,the established application of this technique in coronary artery disease,and more importantly,perspectives and challenges of this technique in the investigation of ICAD.Promising findings of preliminary studies using a CFD model for hemodynamic analysis in ICAD warrant verifications.Further studies in this area will help rectify loopholes in the current secondary prevention strategy,and inform individualized treatment for ICAD patients in the near future.

Computational Fluid Dynamics(CFD) Analysis and Optimization of Reconstructed Intake System of Cylinder Head Based on Slicing Reverse Method

To find out and improve the flow characteristics inside the intake system of cylinder head,the application of computational fluid dynamics(CFD)in the evaluation and optimization of the reconstructed intake system based on slicing reverse method was proposed.The flow characteristics were found out through CFD,and the velocity vector field,pressure field and turbulent kinetic energy field for different valve lifts were discussed,which were in good agreement with experimental data,and the quality of reconstruction was evaluated.In order to improve its flow characteristic,an optimization plan was proposed.The results show that the flow characteristics after optimization are obviously improved.The results can provide a reference for the design and optimization of the intake system of cylinder head.

Numerical simulation of a direct internal reforming solid oxide fuel cell using computational fluid dynamics method

A detailed mathematical model of a direct internal reforming solid oxide fuel cell(DIR-SOFC) incorporating with simulation of chemical and physical processes in the fuel cell is presented. The model is developed based on the reforming and electrochemical reaction mechanisms,mass and energy conservation,and heat transfer. A computational fluid dynamics(CFD) method is used for solving the complicated multiple partial differential equations(PDEs) to obtain the numerical approximations. The resulting distributions of chemical species concentrations,temperature and current density in a cross-flow DIR-SOFC are given and analyzed in detail. Further,the influence between distributions of chemical species concentrations,temperature and current density during the simulation is illustrated and discussed. The heat and mass transfer,and the kinetics of reforming and electrochemical reactions have significant effects on the parameter distributions within the cell. The results show the particular characteristics of the DIR-SOFC among fuel cells,and can aid in stack design and control.

The State Key Laboratory of Numerical Modelling for Atmospheric Science and Geophysical Fluid Dynamics

The State Key Laboratory of Numerical Modelling for Atmospheric Science and Geophysical Fluid Dynamics (LASG) was established in 1985. It was considered an Excellent Opening Laboratory in the national assessment organized by the former State Science and Technology Committee (now the Ministry of Science and Technology) in 1988. Because of this, LASG was upgraded to a State Key Laboratory in 1989. It won the honor of "Excellent

State Key Lab of Numerical Modeling for Atmospheric Science and Geophysical Fluid Dynamics

The State Key Laboratory of Numerical Modeling for Atmospheric Science and Geophysical Fluid Dynamics (LASG) was set up as early as in 1985. It has been appraised as an Excellent State Key Lab in each of the three national assessments organized by State Economic Planning Commission, National Natural Science Foundation (NSFC) and the Ministry of Science and Technology (MOST) in 1991, 1996 and 2000, respectively.

基于CFD-DEM耦合的梯级溜槽的设计与分析

针对传统物料转运过程中溜槽和输送带磨损严重、出口处粉尘浓度过高的问题,建立含臂架的梯级溜槽几何模型,采用基于计算流体力学与离散单元法(CFD-DEM)耦合的数值模拟方法,分析了臂架对转运溜槽的磨损以及对其出口处粉尘排放浓度的影响。仿真结果表明:含臂架的梯级溜槽可以有效控制物料流的速度和方向,降低对溜槽内表面的冲击磨损,降低出口处的粉尘量。

Particulate flow modelling in a spiral separator by using the Eulerian multi-fluid VOF approach

The Euler-Euler model is less effective in capturing the free surface of flow film in the spiral separator,and thus a Eulerian multi-fluid volume of fluid(VOF)model was first proposed to describe the particulate flow in spiral separators.In order to improve the applicability of the model in the high solid concentration system,the Bagnold effect was incorporated into the modelling framework.The capability of the proposed model in terms of predicting the flow film shape in a LD9 spiral separator was evaluated via comparison with measured flow film thicknesses reported in literature.Results showed that sharp air–water and air-pulp interfaces can be obtained using the proposed model,and the shapes of the predicted flow films before and after particle addition were reasonably consistent with the observations reported in literature.Furthermore,the experimental and numerical simulation of the separation of quartz and hematite were performed in a laboratory-scale spiral separator.When the Bagnold lift force model was considered,predictions of the grade of iron and solid concentration by mass for different trough lengths were more consistent with experimental data.In the initial development stage,the quartz particles at the bottom of the flow layer were more possible to be lifted due to the Bagnold force.Thus,a better predicted vertical stratification between quartz and hematite particles was obtained,which provided favorable conditions for subsequent radial segregation.