Large Eddy Simulation Turbulence Model Research Articles

Comparative Analysis of Hydrodynamic Loads on Bridge Piers: Assessing Standards through Numerical Modeling

Currently, there is an alarming difference in the estimated actions on piers resulting from flash-floods and tsunamis as prescribed by international Structural Design Codes. Previous studies have explored the interaction between surge waves and piers from a fluid-dynamics perspective, but there is no assessment about the relevance of the methods used in Structural Design to determine the hydrodynamic pressure resultants. To determine the degree of conservatism embedded in design, results from the different approaches are compared against validated numerical simulations. Three-dimensional numerical analyses were conducted using the OpenFOAM platform, employing the Large Eddy Simulation turbulence modeling approach to simulate the dam-break wave impacting the pier. The findings revealed that the maximum force was not necessarily caused by the stream’s first strike and could be sustained over time. The hydrodynamic pressure center consistently occurred at two-thirds of the stream depth, significantly impacting the overturning moments experienced by the pier. Identified inconsistencies between Codes underscore the need for revision of non-conservative approaches in predicting the design forces.

Large eddy simulation of micro vortex generator-controlled cavitation across multiple stages

Micro vortex generators (mVGs) control cavitation by altering the boundary layer flow structure. This study employs the wall-adapting local eddy-viscosity large eddy simulation (WALE-LES) turbulence model combined with the Zwart–Gerber–Belamri cavitation model to conduct transient numerical simulations on the National Advisory Committee for Aeronautics 0015 baseline hydrofoil and the hydrofoil equipped with mVGs under various cavitation numbers. The proper orthogonal decomposition method and experiments verify the accuracy and consistency of these simulations regarding cavity scale. The study elucidates mechanisms by which mVGs suppress cloud cavitation at low cavitation numbers and induce vortex cavitation at high cavitation numbers. Results indicate that mVGs maintain sheet cavitation characteristics at low cavitation numbers, reducing wall pressure fluctuations and enhancing flow stability. During cavitation inception, mVG-induced vortex cavitation leads to early cavitation formation. In the sheet cavitation phase, modal energy distribution is more dispersed, while in the inception phase, energy is concentrated with significant dominant modes. Moreover, the counter-rotating vortices generated by mVGs mitigate flow separation, enhance leading-edge flow attachment stability, and reduce high-frequency vibrations caused by bubble shedding. This study significantly advances the understanding of cavitation control by accurately simulating and revealing the cavitation control mechanisms of mVGs across different stages using the WALE-LES model. The findings demonstrate that mVGs can effectively stabilize cavity structures at low cavitation numbers, reducing flow instabilities and enhancing overall hydrofoil performance. These insights will have a significant impact on the design of hydrofoils and the development of cavitation control strategies.

Numerical investigation on the ventilated supercavitating model with propeller

To achieve stable and continuous supercavitating flow over underwater vehicles, artificial ventilation is implemented, particularly effective at lower speeds. Previous research on supercavitation primarily focused on analyzing ventilated supercavitating flow with various cavitator types and/or ventilation rates. In this investigation, we examine the behavior of ventilated supercavitating flow over an axisymmetric model featuring both a disk cavitator and a Postdam propeller placed at the bow. Utilizing the Large Eddy Simulation turbulence model, Volume of Fluid method, and Kunz cavitation model, our simulation aims to capture the cavitating flow around the propeller and the ventilated supercavitating flow over the model. Validation of our numerical methods is achieved by comparing our results with experimental data of a ventilated model by Chung and Cho [“Ventilated supercavitation around a moving body in a still fluid: Observation and drag measurement,” J. Fluid Mech. 854, 367–419 (2018)] and the cavitating Postdam propeller by SVA Postdam [S. Potsdam, “PPTC smp'11 Workshop,” in Proceedings of the Workshop on Cavitation and Propeller Performance (2011)]. The results show that with a rotating propeller at the bow of the supercavitating model, the cavitating flow extends and stabilizes compared to configurations utilizing a traditional disk cavitator. The presence of the propeller accelerates the formation of supercavitating flow at a consistent incoming flow speed. Additionally, coupling the propeller with the disk cavitator results in significant increases in propeller thrust, torque, and efficiency. While there is an observed rise in model drag, the impact is not substantial.

A numerical investigation on the dynamics of particle deposition and fouling on a vertical falling film pipe

The zero liquid discharge process is designed to eliminate the discharge of industrial wastewater into the environment, addressing significant environmental concerns that have emerged in recent decades. Falling film evaporator is one of the primary technologies employed in these systems. As observed, there is a noticeable lack of research numerically simulating these phenomena in a vertical pipe including falling film flow. In the present study, an effort has been made to address the existing research gap by conducting numerical simulation of particles deposition and fouling phenomena in a falling film flow in a three-dimensional vertical pipe using computational fluid dynamics. Continuous phases which contain two non-penetrating fluids, i.e., air and water, are simulated using volume of fluid multiphase model and the large eddy simulation turbulence model. Meanwhile, particles motion has been tracked utilizing discrete phase model. Water film enters the pipe with a thickness of 3.7 mm and an inlet velocity of 1 m/s, and particles with a density of 3110 kg/m3 are injected through the water inlet surface in various cases, each with different diameters. As observed, with an increase in particles diameter from 5 to 15 μm, the asymptotic total fouling mass has increased over 14 times, reaching from 0.016 to 0.18 kg/m2, and the total number of deposited particles has also increased. By increasing the fluids inlet velocity from 1 to 1.5 m/s, the total number of deposited particles has decreased throughout the pipe wall, and the total fouling mass has experienced a reduction of 55.8 %.

Aerodynamic performance and flow field characteristics of wind turbines under shear incoming flow conditions

Abstract With the large size of the units, wind farms are evaluated in order to improve the accuracy of the incoming flow measurement and to improve the unit control strategy. In this paper, the turbulent wind field generated by autoregressive linear filtering using a large eddy simulation turbulence model is used to study the wind speed-induced zone under shear incoming flow conditions. At the same time, the wake and aerodynamic performances of the wind turbine under different incoming flow conditions are investigated, and conclusions are obtained. 1) The larger the shear index is, the larger the amplitude is at the first frequency component of the power and the thrust, and the larger the fatigue load is for the wind turbine. 2) For the wind turbine wake flow in the case of smaller wind turbine size, the wake velocity profiles at different shear indices cross and overlap in the transition and far wake regions. At position 0.5 D, the turbulence profiles with different shear indices basically overlap; the larger the shear index is, the more intense the wake flow exchange is and the greater the turbulence intensity in the wake region will be. 3) For the wind speed-induced region, at position −0.1 D, the incoming flow velocity distribution is subject to the combined effect of induction and shear, and the induction effect of the wind turbine wheel is a little larger. The larger the shear index, the greater the turbulence intensity.

Numerical and experimental investigations of swirl-stabilized partially premixed flames using natural gas-hydrogen-air mixtures

This study uses experimental and numerical techniques to examine natural gas and hydrogen mixtures as fuel with a partially premixed swirl-stabilized burner. Variations in the fuel mixture composition were investigated throughout the course of the study. In the experiments, hydrogen was gradually added to the fuel combination, initially made only of natural gas. Chemiluminescence imaging was used to examine how the addition of hydrogen affected the stability and overall structure of the flame. Throughout the experiments, the overall volumetric fuel flow rate remained constant. The change to the flame stabilization mode from a detached flame to a burner attached flame only happened once a certain volumetric hydrogen ratio (H2VOL%) in the fuel was attained. Numerical simulations also evaluated how the lifted flame dimensions and stabilization method changed. The simulations are performed using a commercial compressible, finite volume Large Eddy Simulation software and include the Thickened Flame Model (TFM) for turbulence-chemistry interactions, the Flamelet Generated Manifold (FGM) combustion model, and the Large Eddy Simulation (LES) turbulence model. Experimental and numerical results show that the flame lowers and becomes more compact as the H2VOL% increases. Additionally, experimental and numerical data showed that the flame stabilization mode changes for H2VOL% approaching the value of 70%.

Influence of different turbulence models on tip vortex prediction

Abstract Tip vortex cavitation usually occurs in marine propellers and axial-flow turbines, the cavitation in the tip vortex could potentially result in blade erosion, at the same time, it will cause the pump vibration and noise emission. The tip vortex mainly occurs in rotating machinery with complex three-dimensional flow, but the prediction of the hydrofoil tip vortex can provide a basis for the study of the blade tip vortex in rotating machinery by simplifying the blade model. Hence, in this paper, we use the commercial software CFX to investigate the impact of distinct turbulence models on the configuration of tip vortex distribution in the NACA 16-020 hydrofoil, including SST (Shear Stress Transport), DES (Detached Eddy Simulation), LES (Large Eddy Simulation) turbulence models which are widely used in the field of rotating machinery. At the same time, the axial and radial velocity distributions at different locations downstream of the tip of the hydrofoil are compared to analyse the effects on the velocity distribution near the vortex core line stemming from various turbulence models. By comparing with the experimental data, the most accurate turbulence model for tip vortex prediction is obtained, which will provide guidance for the next step tip vortex prediction and the control of tip vortex flow.

Parametric Study of Flow and Combustion Characteristic in a Cavitied Scramjet with Multi-Position Injection

This study focuses on the three-dimensional flow and combustion characteristics of a cavitied scramjet engine with multi-position injection. A single-equation large eddy simulation (LES) turbulence model is employed, with a detailed reaction mechanism for hydrogen combustion, as described by Jachimowski. The combustion characteristics of hydrogen in the scramjet combustion chamber are analyzed. Based on the combustion chamber model, the influence of different equivalence ratios, injection timing, injection positions, and injection pressures on the flame formation and propagation process are compared. The results indicate that within a certain range, an increase in the equivalence ratio enhances the combustion intensity and chamber pressure. In the case of multi-position injection, the order of injection from different nozzles has little effect on the final flame stabilization mode and pressure distribution. The opposite-side distribution of nozzles can effectively improve the fuel efficiency and the internal pressure. Furthermore, when the nozzles are closely placed in the opposite-side distribution, the combustion efficiency increases, although this leads to a higher total pressure loss. In scenarios where the fuel injection duration is short, an increase in the injection pressure at the upstream nozzles of the cavity results in a higher local equivalence ratio, as well as reduced fuel mixing and ignition time.

Numerical simulation of jet boiling characteristics of high-pressure water injected into high-temperature liquid lead–bismuth in a confined space

The steam generator heat transfer tube rupture (SGTR) accident is potentially devastating and is one of the most serious design basis accidents for lead–bismuth cooled fast reactor (LFR). It is important to study the mechanism of water jet injecting into liquid lead–bismuth eutectic (LBE) for the safety analysis of SGTR accidents. In this paper, based on the volume of fluid (VOF) model, the large eddy simulation (LES) turbulence model and the Lee phase transition model, a three-dimensional numerical model is developed to investigate the typical characteristics of a high-pressure sub-cooled water jet injecting into high-temperature LBE in a confined space. The study focuses on examining the impact of water temperature, inlet pressure, and nozzle diameter on jet development and its surroundings, including the cover gas layer and LBE phase. The results indicate that the central region of the jet is primarily subject to flash boiling, particularly at the nozzle position, whereas the jet interface undergoes heat-transfer boiling. The residual liquid phase migrates upward along the interface and gradually boils. Based on the phase distribution characteristics, a typical jet can be categorized into four zones: the water-steam transition zone, the multiphase flow zone, the end water phase zone, and the steam block zone. Under the conditions analyzed in this paper (water temperature 140–260 °C, injection pressure 2–10 MPa, nozzle diameter 2–10 mm, LBE temperature 400 °C), the jet volume expands approximately exponentially, correlating positively with the three factors studied. It is also found that at high injection pressures and large nozzle diameters, the steam at the top of the jet is prone to destabilization and fragmentation, with steam volume and pressure oscillations being the dominant factors. Steam migration triggers pressure oscillations in the LBE region, and the intensity of these oscillations is closely linked to the migration rate. Under the conditions studied, the maximum pressure oscillation reached 5.5 times the initial pressure. The expansion of the jet steam drives the pressurization of the cover gas layer, which correlates positively with all three factors studied. The nozzle diameter has the greatest effect, increasing the maximum pressure by 168 % from the initial pressure.

Validation Assessment of Turbulent Reacting Flow Model Using the Area-Validation Metric on Medium-Scale Methanol Pool Fire Results

Abstract Accident analysis and ensuring power plant safety are pivotal in the nuclear energy sector. Significant strides have been achieved over the past few decades regarding fire protection and safety, primarily centered on design and regulatory compliance. Yet, after the Fukushima accident a decade ago, the imperative to enhance measures against fire, internal flooding, and power loss has intensified. Hence, a comprehensive, multilayered protection strategy against severe accidents is needed. Consequently, gaining a deeper insight into pool fires and their behavior through extensive validated data can greatly aid in improving these measures using advanced validation techniques. A model validation study was performed at Sandia National Laboratories (SNL) in which a 30-cm diameter methanol pool fire was modeled using the SIERRA/Fuego turbulent reacting flow code. This validation study used a standard validation experiment to compare model results against, and conclusions have been published. The fire was modeled with a large eddy simulation (LES) turbulence model with subgrid turbulent kinetic energy closure. Combustion was modeled using a strained laminar flamelet library approach. Radiative heat transfer was accounted for with a model utilizing the gray-gas approximation. In this study, additional validation analysis is performed using the area validation metric (AVM). These activities are done on multiple datasets involving different variables and temporal/spatial ranges and intervals. The results provide insight into the use of the area validation metric on such temporally varying datasets and the importance of physics-aware use of the metric for proper analysis.

Evolution of cavitation clouds under cavitation impinging jets based on three-view high-speed visualization

Cavitation cloud is an important index for evaluating the unsteady behaviour of cavitation water jets. However, when cavitation jets are applied in a restricted space, the optimal standoff distance is not applied. To study the evolution and collapse mechanisms of cavitation clouds at small standoff distances, this study conducted visualization research on the behaviour of submerged cavitation water jet impacting the target surface at three standoff distances from three-view using a high-speed visualization system. The effect of the optimal standoff distance and the small standoff distance on the morphology of cavitation clouds was explored by observing their shedding behaviour and evolution on the target surface, obtaining a pattern for the morphological evolution of cavitation clouds. The evolution of cavitation clouds under impact flow conditions can be divided into stages: incipient cavitation, ring vortex cavitation, spindle-shaped cavitation cloud, mushroom-type cavitation cloud, cavitation ring, ring vortex cavitation, and collapse. The periodicity of cavitation clouds on the target surface was confirmed by monitoring the grayscale change at the center of the target. Furthermore, the changing of cavitation jet structures is illustrated through Proper Order Decomposition (POD) analysis, showing that the frequency of cavitation clouds increased with the decrease of standoff distance within a certain range. The behaviour of cavitation flow impinging on the target surface at optimal standoff distance and small standoff distance was simulated using the Large Eddy Simulation (LES) turbulence model combined with the Zwart-Gerber-Belamri (ZGB) cavitation model. The results show that the numerical simulation can capture the high-pressure zone generated when the cavitation cloud collapses, and the area of the high-pressure zone increases with the standoff distance within a certain range, the concentration of cavitation bubbles on the target surface did not represent the intensity of erosion. The transient vortex structures of the flow field at different standoff distances were investigated based on the Q-criterion. The findings of this study not only provide critical technical support for the application of cavitation jets in confined spaces, but also contribute to a broader understanding of the erosion mechanism of cavitation jets.

Modeling of the influence of local heat sources on a light gas stratification formation and erosion in a large-scale experimental facility using eddy resolving numerical approach

The main objective of this work is to evaluate the results of application of Implicit Large Eddy Simulation (ILES) turbulence modeling approach to simulate the complex hydrogen safety problem. The paper presents the results of numerical simulation of the OECD/NEA HYMERES (HYdrogen Mitigation Experiments for Reactor Safety) HP2 series tests in the PANDA facility, aimed at studying the thermal effect of passive autocatalytic recombiner (PAR) operation on mixing the containment atmosphere, using CABARET-SC1 code. The code is based on the eddy resolving CABARET technique, which allows implicit modeling of the subgrid turbulence scales without using tuning parameters (ILES approximation). The natural convection flow inside the PAR model was calculated explicitly using a simplified model of the heating section and the housing. Good agreement of calculated transient with experimental data was obtained, which indicates the effectiveness of such approach.

Application of passive vortex generators to enhance vertical mixing in an open raceway pond

A novel use of vortex generators in a raceway pond is described here, with the flow field quantitatively simulated using computational fluid dynamics using the large eddy simulation turbulence model. Persistence lengths of the swirling motion generated by the vortex generators indicate that significant vertical mixing can be achieved by placing vortex generators in the straight section opposite the paddle wheel, downstream of the first hairpin bend. Relatively simple vortex generators are capable of creating stronger swirling motions that persist for a longer distance than those created by the paddle wheel. For optimal performance, vortex generators are positioned side by side but in opposite directions, and their diameters should be equal to or slightly less than the desired liquid depth. The optimal length of a 0.18 m diameter vortex generator in a 0.2 m deep pond was determined to be 0.3 m. Furthermore, it has been demonstrated that a longer persistence length is achieved by inducing a swirling motion with its rotational axis parallel to the primary flow direction.

Thermal and fluid aspects of dual piezoelectric jet pump of double chambers

The rapid development of electronics has led to a significant reduction in the size of electronic devices, while increasing their power consumption. This has resulted in the need for small-scale thermal management solutions. Piezoelectric pump is one of the efficient devices for local thermal management. This paper presents a novel configuration of a piezoelectric pump, which separates the pump chamber with a rigid wall to create two independent pulsating flows. This configuration aims to reduce eddies within each chamber. Furthermore, the proposed design incorporates two pairs of open-sided channels to guide the pulsating flow towards the hot target. The electrostatics physic of the piezoelectric patches along with energy and momentum equations in three dimensions are manipulating mathematically and solved numerically using the finite element method (FEM) and the Arbitrary Lagrangian-Eulerian (ALE) technique. The implicit Large Eddy Simulation (iLES) turbulence model is adopted for this purpose. Several parameters, including frequency, voltage, and the radius of the piezoelectric patch are investigated to better understanding their influence on both the zero-net volumetric flow rate and the heat transfer process. The findings indicate that by separating the pump chamber, it becomes feasible to obtain a 3D solution with the iLES using a reasonable mesh size. Nusselt number increases with all the examined parameters, where it enhances by 42.48 % at V = 100 Volt when the pump is switched on compared to natural convection mode. The peak-peak volumetric flowrates rises by 50 % and 350 % when the applied voltage is raised from 50 to 100 Volts and the frequency is raised from 20 to 50.50 kHz, respectively.

Large-Eddy Simulation Study of Injector Geometry and Parcel Injection Location on Spray Simulation of the Engine Combustion Network Spray G Injector

Abstract Recent improvements in computing power and numerical techniques have enabled us resolve minute details of spray plume behavior and its wall boundary interactions, and made detailed spray simulations using large-eddy simulation (LES) turbulence models available to many more engineers. However, guidelines for parcel-based spray simulation boundary and initial conditions are still based on results from lower-resolution and Reynolds averaged Navier–Stokes (RANS) simulations. Hence, it is necessary to critically examine those assumptions and compare their impact with results using a RANS turbulence model. Three different parameters, including whether a simulation includes a detailed injector tip geometry or a flat surface, whether parcels are initialized at the counterbore exit, which is more common, or at the nozzle exit, and whether to use an experimentally derived rate of injection or one-way coupling with a separate internal nozzle volume of fluid simulation, were examined with an LES turbulence model. Both local data close to the injector and global penetration results were used to compare simulations. Local data, such as the local liquid volume fraction, showed greater variation between the conditions, which may have an impact on mixing and combustion predictions in engine applications. Spray penetration and other global measures demonstrated limited sensitivity to the boundary conditions/initialization procedure. Results were also compared with prior results that used a RANS turbulence model. RANS simulations had overall smoother responses to the changes, as would be expected, but LES simulations showed similar trends in the effects for the measured variables.

CFD Modeling and Analysis of a Heat-Generating Melt Pool Under External Vessel Cooling and Top Flooding

In some nuclear reactors, under accidental conditions, core debris forms a molten pool, which is later located in a core catcher. The core catcher proposed by the authors uses special refractory material to absorb enthalpy of corium so that temperatures are within 1500 K, which is possible to cool with side cooling and top flooding. Since performing a full-scale prototypic experiment is extremely challenging and complex because of the involvement of very high temperatures and the presence of radioactive materials, it is important to develop a Computational Fluid Dynamics (CFD) model capable of simulating coolability of the melt pool with the above cooling strategy. In the present work, a CFD model was developed for the above purpose and was benchmarked with experiments conducted under simulated conditions by the authors. The experiment involved the melting of about 25 L of sodium borosilicate glass at about 1473 K and cooling it in a scaled-down core catcher model. In the presence of decay heat inside the melt pool, turbulent natural convection plays an important role in the temperature distribution inside the melt pool and on the vessel walls. For this, we used different turbulence models. Comparisons among the Standard k-ε, Shear Stress Transport (SST) k-ω, and two-dimensional (2D) Large Eddy Simulation (LES) turbulence models show that SST k-ω and 2D LES turbulences are found to be in good agreement with the experimental results for the temperature distribution in the melt pool, and SST k-ω is found to be computationally less expensive than 2D LES. In general, the CFD model is capable of simulating heat transfer with phase changes inside the heat-generating melt pool. In view of this, the model can be further extended to include cooling of the melt pool in the prototype core catcher. The evolution of crust formation has been investigated in detail using a CFD model.

Effect of pitch ratio on the cavitation of controllable pitch propeller

This paper numerically investigated the effects of the pitch ratio on the cavitation of a controllable pitch propeller. The geometric model of marine propellers is obtained by changing the pitch of the propeller VP1304. The computational fluid dynamics (CFD) solver STAR-CCM+ and large eddy simulation (LES) turbulent model are used to analyze marine propellers. The open water characteristics and cavitation behavior obtained by numerical simulation are in good agreement with available experimental data. Because the controllable pitch propeller always operates under high-efficiency conditions, the open water characteristic curves of four propellers with different pitch ratios are first calculated to obtain the corresponding advance coefficients at the highest efficiency point. Then the cavitation behavior at the highest efficiency points of four propellers with different pitch ratios is analyzed. As the pitch increases, the main type of cavitation gradually changes from leading edge cavitation to hub vortex cavitation. In addition, the causes of different cavitation types are analyzed from the aspects of pressure and velocity distribution.

Numerical study of motorbike aerodynamic wing kit

Abstract This study aims to design the configuration of an aerodynamic wing kit (AWK) on a racing motorbike to achieve the highest downforce-to-drag ratio. The numerical study involves a motorbike traveling in a straight line, where the AWK improves performance and safety by generating downforce to prevent lift. The geometry of the AWK uses a NACA 4412 airfoil with a span of 0.6 m. The computational mesh is generated using SnappyHexMesh and installed on a simplified motorbike to minimize the mesh skewness. The Navier–Stokes equations are solved with OpenFOAM computational fluid dynamics using the Reynolds-averaged Navier–Stokes k-ω SST and large eddy simulation (LES) turbulence models. Case 1 compares a motorbike with and without a dummy, both equipped with the AWK, varying the angle of attack (AoA) from 0 to −41°. Case 2 studies the single wing at different wind speeds (i.e. 20, 60 and 100 m/s) to determine the highest downforce-to-drag ratio at an AoA of −37°. These results serve as the basis for Case 3, which investigates non-parallel wing configurations with a fixed upper wing and a rotating lower wing. In Case 4, where both upper and lower wings rotate simultaneously as parallel wings, the peak downforce-to-drag ratio occurs at an AoA of −41°. Finally, Case 5 modifies the AoA of −41° of the parallel wing of Case 4 to a closed-wing version to comply with Fédération Internationale de Motocyclisme safety regulations. With the LES turbulence model, unsteady and complex turbulence structures can be visualized using the Q-criterion. A comparison of the time-averaged lift coefficient between the wingless and closed-wing configurations shows an increase in downforce of approximately 360%. Subsequently, the popularity of AWK will contribute to the safety of racing motorbike driving.

SIMULATING SPRAY DYNAMICS WITH A FINITE ELEMENT METHOD FOR INTERNAL COMBUSTION ENGINES USING LARGE EDDY SIMULATIONS

Spray dynamics in an internal combustion engine is comprised of complex phenomena while interacting with unsteady turbulence. Physics requires detailed modeling of the dynamics for spray and carrier gases to accurately predict a spray's fate. Large eddy simulation (LES) turbulence modeling approaches are capable predictors of the turbulent processes and of dynamically modeling sub-grid scales, therefore enabling calculation of model coefficients for the smallest resolved scale. Dynamic LES methods are also well suited for unsteady flows associated with spray injection and engine fluid dynamics. In this study, for the first time, a dynamic Verman LES scheme employed in a stabilized a finite element (FE) framework is used to model the spray process with emphasis on injected fuels for simulating internal combustion engines. Spray modeling often comprises a coupled Eulerian-Lagrangian approach to capture the droplet/particle dynamics, where the droplets are modeled in the Lagrangian frame, interacting with the gas. The momentum and heat exchange between the fluid gases and the evaporating and atomizing spray droplets are modeled in a two-way coupling system as described in this paper. Direct injected liquid is modeled in this paper as a spherical ligament of fuel and ligament breakup to atomization use our version of the Kelvin Helmholtz break-up scheme. Discussed are models and methods of the whole system in some detail. The method for simulation of the fluid's momentum, heat transfer, and turbulence are discussed as is the system to evaluate droplet or ligament properties. Validation or results of the modeling are presented on test cases as determined by the Engine Combustion Network, base cases G along with some benchmarking cases. We demonstrate in the paper by using the FE method combined with the Verman LES model, and our version of the Kelvin Helmholtz-Rayleigh Taylor spray break-up model combined with two-way coupling is quite capable of modeling sprays, having distinct and significant advantages over the more traditional finite volume discretization framework using similar spray modeling methods.

Analysis of Dam Break Wave Using Analytical, Computational Fluid Dynamics, and Experimental Approaches

This research aims to examine the capability of the Computational Fluid Dynamics (CFD) method in simulating the behavior of dam break waves. It begins by building a 2D numerical simulation using OpenFOAM. To overcome the influence of turbulence, we employed the Large Eddy Simulation (LES) turbulent model, specifically the k-Equation and Smagorinsky model. The simulation was developed by applying the Navier-Stokes equations using the finite volume method in OpenFOAM. The analysis focuses on the free surface of a dam break. The results are in good accordance with both analytical and experimental results. The simulation has followed the trend of experimental and analytical free surface profiles at the dam break’s early and late conditions. The low mesh number on the computational domain caused significant differences in the wavefront of the dam break. It reduced the accuracy of the calculation between the water and air interface. This study highlights the importance of understanding dam break wave behavior as part of risk mitigation for dam leakage. The behavior of dam break waves can be observed by determining observation positions at different locations, with the water gate of a dam serving as the reference point. These highly accurate numerical results indicate that the CFD approach employing OpenFOAM can be relatively cost-effective yet accurate in analyzing multiphase problems, such as dam breaks. This CFD approach is expected to contribute to developing mitigation and disaster prevention in the future.