Large Eddy Simulation Model Research Articles

Numerical Simulation of Flow and Flame Dynamics of a Pool Fire Under Combined Effects of Wind and Slope

Wind has a significant effect on pool fire behavior, which is relevant to many fire conditions, such as wildfires, building fires, and oil transportation fires. Although fire behavior and morphology changes have received considerable attention and been widely researched, there are few works concerning the flow and flam dynamics of pool fire. A large eddy simulation model is adopted to investigate the flow and flame dynamics of a rectangular pool fire considering the combined effects of wind and slope. The results show that, with a wind speed of 0.5 m/s, a flame develops immediately downstream of the fire source and sustains two flanks of plume. Further downstream, the plume starts to rise due to buoyant force. Temperature, velocity, and vorticity distributions show significantly different shapes at different streamwise locations. Near the fire source, the flame is confined to a small region around the fire source. The air circulation downstream shows a cylindrical spiring pattern. When the wind speed increases, the temperature and velocity become more parallel to the surface and their maximum values increase. On the contrary, the temperature fluctuations and turbulent kinetic energy decrease with the wind speed, and they are more frequent near the flame tails.

Improving the estimate of higher-order moments from lidar observations near the top of the convective boundary layer

Abstract. Ground-based lidar data have proven extremely useful for profiling the convective boundary layer (CBL). Many groups have derived higher-order moments (e.g., variance, skewness, fluxes) from high-temporal-resolution lidar data using an autocovariance approach. However, these analyses are highly uncertain near the CBL top when the depth of the CBL (zi) is changing during the analysis period. This is because the autocovariance approach is usually applied to constant height levels and the character of the eddies is changing on either side of the changing CBL top. Here, a new approach is presented wherein the autocovariance analysis is performed on a normalized height grid, with a temporally smoothed zi. Output from a large eddy simulation model demonstrates that deriving higher-order moments from time series on a normalized height grid has better agreement with the slab-averaged quantities than the moments derived from the original height grid.

Experimental and numerical investigations into cold flow characteristics of multiple micro-mixing jets for hydrogen-rich gas turbines

To explore the flow field of micro-mixing jets, cold flow characteristics of a model Micromix burner were investigated by particle image velocimetry (PIV) system and Large-eddy simulation (LES) model. Results show that LES results are in good agreement with experimental results. In the flow field of multiple micro-mixing jets, the jet velocities of nozzles farther away from the burner center have a high increase and decay rate. When the outlet Reynolds number increases, the Reynolds stress increases first and then decreases in the merging region indicating that the velocity fluctuation disappears in the second jet half, but it has little effect on the flow field structure. Comparing the flow fields of round multiple micro-mixing jets, the merging point and combined point in the elliptical jets flow field move backward. Moreover, the maximum velocity for elliptical jets is also faster than the round jets, which is caused by the high turbulent kinetic energy in the elliptical jet flow field. When the tube spacing increases from 2 to 3 times the tube diameter, positions of the two feature points change linearly. Further, the surrounding jets can decrease the velocity attenuation of the center nozzle and elongate the axial length of the two feature regions.

A biparametric analysis on the steady performance of bidirectional pumping hydrodynamic-static hybrid dry gas seal

PurposeThe purpose of this paper is to obtain the matching relationship between spiral groove, equalizing groove and operation parameters through the biparametric analysis for the bidirectional pumping hydrodynamic-static hybrid dry gas seal (BP-HHDGS).Design/methodology/approachThe large eddy simulation (LES) model in Fluent is used to simulate the flow field of BP-HHDGS, and the biparameter variables method is chosen to analyze the effects of different parameters on the performance of BP-HHDGS.FindingsBP-HHDGS has a greater opening force than hydrostatic dry gas seal (HDGS); the vortex is formed after lubricating gas is exhausted from the throttle. Increasing the depth of the equalizing groove and spiral groove has a synergistic enhancement effect on the opening force and leakage of BP-HHDGS. There is a matching relationship between spiral angle and rotational speed. The preferred parameter ranges in current conditions are found as follows: spiral angle αa = 15°–24°; groove-dam ratio λ = 0.4–0.7; equalizing groove depth hj > 35 µm; spiral groove depth hg = 5-10 µm.Originality/valueThe high starting capacity of HDGS is given to the hydrodynamic type seal, and thus the application promotion of HDGS in high-speed working condition is realized at the same time. This work also provides precise and quick theoretical guidance for the selection and design of hydrodynamic-static dry gas seal and further promotion.

Assessment of Thermal Boundary Models for Large Eddy Simulations of Natural Convection

Assessment of Thermal Boundary Models for Large Eddy Simulations of Natural Convection

LES and RANS Modeling of an 84-Pin Hexagonal Rod Bundle with Spacer Grid and Experimental Validation

As part of an ongoing integrated research project (IRP), detailed investigations of the flow characteristics of an 84-pin hexagonal rod bundle representing a potential modular gas-cooled fast reactor design have been performed. The rod bundle features a pitch-to-diameter ratio of 1.5 and a large central guide tube, along with simple spacer grids and an outer enclosure. The primary focus of the current work is modeling and simulation of this geometry using multiple computational techniques. These include high-fidelity large eddy simulation (LES) and advanced Reynolds-averaged Navier-Stokes (RANS) approaches. The flow predictions are validated at a Reynolds number of 12 000 using matched index of refraction particle image velocimetry experimental data that were also generated by Texas A&M University as part of the IRP. The comparison data include velocity magnitude and turbulent kinetic energy (TKE) at different lateral locations and at multiple distances downstream of the spacer grid. Many characteristics of the experimental flow are replicated in the LES and RANS runs, and a general agreement between simulation and experiment is shown. Except for the immediate vicinity of the grid, LES is able to capture the velocity magnitude within a 10% error and the TKE to within the experimental uncertainty. The LES shows notably better agreement with the experimental data than RANS, particularly regarding the TKE decay behavior downstream of the grid. Both methods show significant departures from the experiment in their predictions close to the central guide tube. The experimental and high-fidelity data will be used to inform closures for novel multiscale methodologies. The long-term goal of the work is to use the high-fidelity data to yield improved multiscale thermal analysis techniques for solving fuel performance problems of direct relevance to industry.

Experimental characterization and 3D simulations of turbulent flames assisted by nanosecond plasma discharges

This paper presents quantitative experimental data generated for the validation of plasma-assisted combustion (PAC) simulations. These data are then used to validate the phenomenological model of Castela et al. They are also useful to test other PAC models. In the experiment presented here, Nanosecond Repetitively Pulsed (NRP) discharges are applied to a lean-premixed turbulent methane-air flame initially near the lean blow-off limit. The discharges significantly enhance the combustion and stabilize the flame after a few pulses. Electrical and optical diagnostics are employed to extensively quantify the transient and steady state of the plasma-stabilization process. The flame shape is characterized by OH* chemiluminescence imaging. In the discharge region, OH density profiles are obtained by 1D laser-induced fluorescence, and the gas temperature is measured by optical emission spectroscopy measurements. The local gas temperature increases by 1250 K, and the OH number density rises sevenfold when NRP discharges are applied. These results evidence the cumulative thermal and chemical effects of NRP discharges, which are especially challenging to replicate numerically. A Large Eddy Simulation (LES) of the experiment is performed. Combustion chemistry is modeled by an analytically reduced mechanism, while the plasma discharge is described by the low-CPU cost phenomenological model of Castela et al., which aims to capture the main thermal and chemical effects induced by the discharges. The model of Castela et al. is validated in the burnt gases by the remarkable agreement between the simulations and the experiments regarding the flame shape, the local gas temperature, and the OH number density. More generally, this work demonstrates the relevance of simplified plasma models in LES solvers to simulate complex plasma-assisted burners.

Turbulence model adaptability at critical Reynolds numbers and applications in wake control via fairings

The water-drop shaped fairings with varying shape angles are attached to a circular cylinder to achieve wake control and vortex suppression at critical Reynolds numbers. To ensure the capability of Reynolds averaged Navier–Stokes (RANS), detached eddy simulation (DES) and large eddy simulation (LES) models at the critical Reynolds number region, three representative turbulence models are employed: LES with σ subgrid-scale (SGS) model, delayed DES model with improved wall-modeling capability (IDDES) and shear stress transport (SST) k−ω RANS model. These models are utilized to simulate flow around a circular cylinder at Reynolds number Re=2.5×105. The solver used in this paper is further developed based on the high-resolution algorithm platform for incompressible flow (HRAPIF). The comparative analysis of the results from the three turbulence models has been rigorously validated and investigated. An exhaustive examination of the mean flow field, Reynolds stresses, characteristic lengths, and instantaneous flow fields among the models reveals instructive insights. The IDDES and σ-LES models predict the hydrodynamic forces, the so-called ‘drag crisis’, alongside the pressure distribution and skin friction coefficient with high precision. The σ-LES model stands out for its superior accuracy, while the IDDES model is also a viable alternative, offering commendable accuracy with a reduced demand for mesh density. Subsequently, the IDDES model is selected for wake control calculations using fairings with five distinct shape angles (30∘,45∘,60∘,75∘ and 90∘) at Re=2.5×105. In-depth comparisons to the bare cylinder and subcritical results reveal that the wake control effect varies at the critical Reynolds region. The fairing with a 30° shape angle substantially suppresses hydrodynamic forces. The lift coefficient experiences a remarkable decrease of approximately 96%, while the drag coefficient diminishes by about 90%. Concurrently, fairings with angles from 45° to 90° lead to reductions in drag coefficient of 11.6%, 10%, 3% and 4%, respectively. A 75% lowering in the lift force coefficient is found even with a short fairing with shape angle α=90∘, which may be a meaningful instruction to industrial design.

Subgrid-scale model for large eddy simulations of incompressible turbulent flows within the lattice Boltzmann framework.

Large eddy simulations are a popular method for turbulent simulations because of their accuracy and efficiency. In this paper, a coupling algorithm is proposed that combines nonequilibrium moments (NM) and the volumetric strain-stretching (VSS) model within the framework of the lattice Boltzmann method (LBM). This algorithm establishes a relation between the NM and the eddy viscosity by using a special calculation form of the VSS model and Chapman-Enskog analysis. The coupling algorithm is validated in three typical flow cases: freely decaying homogeneous isotropic turbulence, homogeneous isotropic turbulence with body forces, and incompressible turbulent channel flow at Re_{τ}=180. The results show that the coupling algorithm is accurate and efficient when compared with the results of direct numerical simulations. Using calculation format of the eddy viscosity, a uniform calculation format is used for each grid point of the flow field during the modeling process. The modeling process uses only the local distribution function to obtain the local eddy viscosity coefficients without any additional processing on the boundary, while optimizing the memory access process to fit the inherent parallelism of the LBM. The efficiency of the calculation is improved by about 20% compared to the central difference method within the lattice Boltzmann framework for calculating the eddy viscosity.

Study of surface layer characteristics in the presence of suspended snow particles using observational data and Large-Eddy Simulation

The snowdrift is a two-phase flow consisting of air and suspended particles. In the presence of snow particles in the air, additional stability appears in the surface layer due to the density gradient. The density gradient reduces turbulence and affects the properties of the surface layer. Therefore, to describe the properties of the flow with included snow particles, additional clarifications are required. A description of the surface layer parameterization with the presence of suspended snow particles is presented in this paper. The formulation of the effect of snow particles consists in reformulation of the Obukhov turbulent length scale. The novel surface layer parameterization allows to take into account the effect of snow particles on turbulent flow and may improve the estimates of friction velocity and boundary-layer height.The parameterization was successfully tested on the observational data. Description of snow particles influence was included in the Large-Eddy Simulation (LES) model. The numerical experiments confirmed an increase in the stability of the surface layer. Mechanism of suspended particles influence on the surface layer is analogous to a thermal stabilization of the turbulent flow, in which negative buoyancy acts to reduce the turbulent kinetic energy.

Development of a Multi‐Scale Meteorological Large‐Eddy Simulation Model for Urban Thermal Environmental Studies: The “City‐LES” Model Version 2.0

AbstractTo bridge the gaps between meteorological large‐eddy simulation (LES) models and computational fluid dynamics (CFD) models for microscale urban climate simulations, the present study has developed a meteorological LES model for urban areas. This model simulates urban climates across both mesoscale (city scale) and microscale (city‐block scale). The paper offers an overview of this LES model, which distinguishes itself from standard numerical weather prediction models by resolving buildings and trees at the microscale simulations. It also differs from standard CFD models by accounting for atmospheric stratification and physical processes. Noteworthy features of this model include: (a) the calculation of long‐ and short‐wave radiations in three dimensions, incorporating multiple reflections within urban canopy layers using the radiosity method, and accounting for building and tree shadows in the simulations; (b) the provision of various heat stress indices (Universal Thermal Climate Index, Wet Bulb Globe Temperature, MRT, THI); (c) the assessment of the efficacy of heat stress mitigation measures such as dry‐mist spraying, roadside trees, cool pavements, and green/cool roofs strategies; (d) the capability to run on supercomputers, with the code parallelized in a three‐dimensional manner, and the model can also run on a graphics processing unit cluster. Following the introduction of this model, the study confirms its basic performance through various numerical experiments, including simulations of thermals in the convective boundary layer, coherent structure of turbulence over urban canopy, and thermal environment and heat stress indices in urban districts. The model developed in this study is intended to serve as a community tool for addressing both fundamental and applied studies in urban climatology.

Study on the characteristics of fluid exciting force on impeller of high-speed centrifugal pump under the condition of poiseuille inflow

Abstract The paper adopted the turbulence model of Wall-Model Large Eddy Simulation with improved algebraic form (WMLES S-Omega) to carry out the study on the fluid exciting force on impeller of high-speed centrifugal pump under the conditions of Poiseuille inflow and uniform inflow. The external characteristics and the vibration displacement of the shaft of high-speed centrifugal pump are tested. The results indicated that it is the high matching that the test and numerical external data of high-speed centrifugal pump under the condition of Poiseuille inflow. The test results of vibration displacement of shaft match the fluid exciting force better in terms of waveform and phase difference under the condition of Poiseuille inflow. The contribution degree of force and torque from the shroud inner surface in the X and Y direction is 76.3%, 80.7%, 79.5% and 85% of the resultant force and torque, and its force from the inner and outer surfaces of shroud and hub in the Z direction is 94.3%, and its torque from blades’ surfaces in the Z direction is 98.3%. The study results can provide some theoretical support for the research of fluid exciting force and hydraulic optimization design of centrifugal pump.

Eddy viscosity enhanced temporal direct deconvolution model for temporal large-eddy simulation of turbulent channel flow

In this work, a novel eddy viscosity enhanced temporal direct deconvolution model (TDDM) is proposed for temporal large-eddy simulation (TLES) of turbulent channel flow at large filter widths. To improve the accuracy of the constant eddy viscosity (CEV) model, particularly in the near-wall region, a damping function is incorporated to refine its performance. Moreover, a spatial filtering strategy is introduced to reduce the aliasing errors associated with the computation of subfilter-scale (SFS) stress, thereby enhancing numerical stability. In the a posteriori study, the accuracy of the CEV model is assessed comprehensively by comparing the TLES results with corresponding temporally filtered direct numerical simulation data. The results demonstrate that the CEV-enhanced TDDM provides accurate predictions across various statistical properties of velocity, instantaneous flow structures, kinetic energy spectra, and SFS energy fluxes. The coefficient sensitivity analysis of the CEV model reveals that the model coefficient significantly influences low Reynolds number flows, while its impact on high Reynolds number flows is relatively small. TLES on coarse grids demonstrate that the CEV-enhanced TDDM exhibits strong robustness and accuracy at different grid resolutions. Additionally, the CEV-enhanced TDDM in high Reynolds number flows is stable and accurate at remarkably large filter widths.

Quantifying the Wind‐Induced Bias of Rainfall Measurements for the Thies CLIMA Optical Disdrometer

AbstractThe wind‐induced bias of rainfall measurements obtained from non‐catching instruments is addressed in this work with reference to the Laser Precipitation Monitor (LPM) optical disdrometer manufactured by Thies CLIMA. A numerical simulation approach is adopted to quantify the expected bias, involving three different models with increasing complexity. Computational Fluid‐Dynamics simulation of the airflow field around the instrument with an embedded Lagrangian particle‐tracking module to obtain raindrop trajectories are performed by solving the Unsteady Reynolds Averaged Navier‐Stokes (URANS) equations and a Large Eddy Simulation (LES) model. URANS‐uncoupled, LES‐uncoupled, and LES‐coupled approaches are tested to assess the impact of modeling the airflow turbulent fluctuations in detail. Due to the non‐radially symmetric external shape of the instrument, various combinations of the wind speed and direction are considered. Catch ratios for monodisperse rain are obtained as a function of the particle Reynolds number and the wind direction and fitted to obtain site‐independent curves to support application of the simulation results. Based on literature expressions to link the drop size distribution of real rainfall events with the rainfall intensity (which instead depend on the local rainfall climatology at the measurement site), sample collection efficiency curves are obtained from the catch ratios of monodisperse rain. The resulting adjustment curves allow rainfall measurements to be corrected using either a real‐time or post‐processing approach. However, at high wind speed and assuming that the wind blows parallel to the instrument sensing area, the instrument may fail to report precipitation altogether.

Impacts of building modifications on the turbulent flow and heat transfer in urban surface boundary layers

This study examines turbulent airflow and upward heat transport in real urban environments using a building-resolving large-eddy simulation model to understand the characteristics of turbulent airflow and upward heat transport when geometrical distributions of buildings are modified. The target areas were two real urban districts within Osaka City, Japan, having different morphological features. In the numerical experiments, the initial condition was set to a neutral condition in which temperature is uniformly distributed vertically, and buildings emitted heat at a constant rate. The results in the two districts indicated that the features of turbulence and heat transport distinctly differed with different building arrangement. Specifically, taller buildings significantly decelerated airflows and induced warming behind buildings. More high-rise buildings (which resulted in a larger building variability) in a district with a larger building density caused a large heat flux and warming at higher levels. The sensitivity experiments in which a density and height variability of buildings were modified showed that a building density at higher levels and a building height variability significantly influenced warming at upper levels. An increased building height variability weakened wind speed and disturbed horizontal heat advection, whereas a large building density caused numerous heat sources.

An improved detached eddy simulation method for cavitation multiphase flow

The evolution of cavitation flow involves intense transient characteristics and complex multi-scale motions, significantly increasing the difficulty and computational cost of solving in separation zones. This study proposes an Improved Detached Eddy Simulation (IDES) method, based on a single-equation eddy-viscosity model, aimed at enhancing the resolution of small-scale cavitation flows. By incorporating the broad-spectrum von Karman scale concept, the method effectively improves the resolution capability for small-scale flows in separated zones. The performance of the IDES model was validated through three computational cases, and the main conclusions are as follows: The results of the triangular cylinder flows show that the IDES model demonstrates a resolution capability for large-scale turbulence comparable to the Large Eddy Simulation (LES) model. The orifice flow results indicate that activating the LES method with insufficient grid resolution is highly discouraged, whereas the IDES model performs more robustly under such conditions. For hydrofoil flows, the IDES model more accurately captures the temporal evolution of cavitation, avoiding the premature cavity shedding predicted by LES under coarse grid conditions. The decomposition results of the vorticity transport equation show that both the dilatation term and viscous term are significant contributors in the vorticity transport process, with their average contributions to vorticity reaching 49.25% and 39.99%, respectively. Vortices can be considered the primary controllers of cavitation dynamics, while the dilatation term and viscous term can be regarded as the main controllers of vorticity. Overall, the energy compensation mechanism of the IDES model, based on local flow information, gives it unique advantages in handling small-scale turbulence and cavitation phenomena, providing both high computational efficiency and accuracy.

Study on the hydrogen-air premixed flame propagation characteristics in semi-open space with obstacle

In semi-open space, obstacles have the potential to accelerate flame propagation and increase hydrogen-air deflagration pressure. Therefore, this paper is dedicated to exploring the influence of obstacles on the premixed hydrogen deflagration, which is crucial for enhancing the safety of industrial production and energy utilization. By using Large Eddy Simulation (LES) model in OpenFOAM, this study investigates the deflagration characteristics of premixed hydrogen in the presence of three different shaped obstacles. The analysis results reveal that under obstacle conditions, the flame shape can be categorized into four phases: the hemispherical phase, finger-shaped phase, jet phase, and vortex phase. The velocity of the flame front is nearly same for elliptical and rectangular obstacle condition, but it is 36% higher compared to triangular obstacle condition. The impact of triangular and rectangular obstacles on explosion overpressure is less than that of triangular obstacles on explosion overpressure, but it is 16% higher than that of elliptical obstacle. Analyzing the vorticity generated by different obstacle reveals that the vorticity produced by rectangular obstacle is twice as much as that produced by elliptical obstacle, whereas the vorticity produced by triangular obstacle is 2.4 times greater than that produced by elliptical obstacle. The acceleration of hydrogen-air explosion process occurs due to the narrow space created by obstacle and pipeline walls, and the shape of obstacle significantly influences this acceleration.

Multi-dimensional modeling of solitary wave–structure interaction problems by using a [formula omitted]-LES-SPH model

The interaction mechanisms between waves and marine structures are a popular research topic. This paper applies the weakly compressible smoothed particle hydrodynamics (WCSPH) method to study the dynamics of green water overtopping. To enhance the accuracy of the simulations, the SPH method coupled with the large eddy simulation (LES) model is employed for numerical investigations. Initially, we validate the effectiveness of the model by simulating the generation of solitary waves and irregular waves, as well as numerically reproducing the water surface morphology during the interaction between solitary waves and the deck. Subsequently, the validated model is used to study the dynamic characteristics of different types of waves overtopping, revealing significant variations in their motion. Furthermore, we investigate the effect of deck roughness during the entire green water overtopping process in terms of both protrusions extent and distribution, confirming that a reasonable setting of the protrusions can greatly reduce the wave impact loads on the deck, thereby protecting the structure. Additionally, a three-dimensional model is developed to study the green water problem, and we find that the turbulence phenomenon is more pronounced in the three-dimensional scenario.

Aerosol composition, air quality, and boundary layer dynamics in the urban background of Stuttgart in winter

Abstract. Aerosol distributions are of great relevance for air quality, especially for cities like Stuttgart, which has limited air exchange due to its location in a basin. We collected a comprehensive set of data from remote sensing and in situ methods including radiosondes for the urban background of downtown Stuttgart to determine the impact of boundary layer mixing processes on local air quality and to evaluate the simulation results of the high-resolution large eddy simulation (LES) model PALM-4U at 10 m grid spacing. Stagnant meteorological conditions caused accumulation of aerosols, and chemical composition analysis shows that ammonium nitrate (37 ± 9 %) and organic aerosol (OA; 34 ± 9 %) dominated during this winter study. Case studies show that clouds during previous nights can weaken temperature inversion and accelerate boundary layer mixing after sunrise by up to 3 h. This is important for ground-level aerosol dilution during the morning rush hour. Furthermore, our observations validate results of the LES model PALM-4U in terms of boundary layer heights and aerosol mixing for 48 h. The simulated aerosol concentrations follow the trend of our observations but are still underestimated by a factor of 4.5 ± 2.1 due to missing secondary aerosol formation processes and uncertainties of emissions and boundary conditions in the model. This paper firstly evaluates the PALM-4U model performance in simulating aerosol spatio-temporal distributions, which can help to improve the LES model and to better understand sources and sinks for air pollution as well as the role of horizontal and vertical transport.

Numerical Study of High-g Combustion Characteristics in a Channel with Backward-Facing Steps

High gravity (high-g) combustion can significantly increase flame propagation speed, thereby potentially shortening the axial length of aero-engines and increasing their thrust-to-weight ratio. In this study, we utilized the large eddy simulation model to investigate the combustion characteristics and flame morphology evolution of premixed propane–air flames in a channel with a backward-facing step. The study reveals that both the increase in centrifugal force and flow velocity can enhance pressure fluctuations during combustion and increase the turbulence intensity. The presence of centrifugal force promotes the occurrence of Rayleigh–Taylor instability (RTI) between hot and cold fluids. The combined effects of RTI and Kelvin–Helmholtz instability (KHI) enhance the disturbance between hot and cold fluids, shorten the fuel combustion time, and intensify the dissipation of large-scale vortices. The increase in fluid flow velocity can raise the flame front’s hydrodynamic stretch rate, thereby enhancing the turbulence level during combustion to a certain extent and increasing the fuel consumption rate. When a strong centrifugal force is applied, the global flame propagation speed can be more than doubled. Within a certain range, the increase in high-g field strength can enhance the intensity of RTI and accelerate the transition of RTI to the nonlinear stage.