Large Eddy Simulation Model Research Articles

Dry powder inhaler deposition in the larynx and the risk of steroid inhaler laryngitis: a computational fluid dynamics study

Dry Powder Inhalers (DPIs) are a mainstay in the treatment of obstructive respiratory diseases, including asthma and chronic obstructive pulmonary disease (COPD). Deposition of inhaled corticosteroids in the larynx elicits local side effects, potentially leading to steroid inhaler laryngitis. The objective of this study was to estimate the dose of DPIs that are deposited in the larynx relative to other regions of the respiratory tract using computational fluid dynamics (CFD). An anatomically accurate model of the airways (mouth to main bronchi) was constructed based on medical imaging of a healthy adult. Respiratory airflow and particle transport were simulated for constant inhalation rates of 30, 45, and 60 L/min. Two turbulence models were compared, namely the large eddy simulation (LES) and the models. DPIs were assumed to generate an aerosol cloud with a log-normal particle size distribution characterized by the mass median aerodynamic diameter () and geometric standard deviation (). We compared two commercial DPIs, namely DPI 1 had a large particle size ( = 50 μm, ) and DPI 2 had a small particle size ( = 2 μm, ). The laryngeal dose was 1.6-to-3.8-fold higher than the bronchial dose for DPI 1, while the laryngeal and bronchial doses (units of mass per unit surface area) were similar for DPI 2 for both turbulence models and all inhalation rates. The choice of turbulence model had little impact on the total extrathoracic deposition, but a significant impact on regional doses, with the LES model predicting higher larynx-to-bronchi relative doses than the model. Our prediction that the larynx is a hotspot for DPI deposition is consistent with the observation of laryngeal side effects in DPI users. Importantly, our simulations suggest that DPIs with larger particles ( = 50 μm) may increase the risk of steroid inhaler laryngitis.

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

Purpose The 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/approach The 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. Findings BP-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/value The 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

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

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.

Well-planned greenery improves air urban quality - Modelling the effect of altered airflow and pollutant deposition

Urban air quality is influenced by vegetation through alterations in airflow and pollutant deposition processes. We investigated these interactions by integrating the Vegetation Impact Dynamic Assessment model (VIDA) with the Large-Eddy Simulation model PALM. Our analysis focus on nitrogen dioxide (NO₂) and particulate matter (PM) concentrations at the local scale, considering three tree genera. Our findings reveal the necessity of accounting for both gaseous pollutants and particles separately due to their differing mechanisms of deposition onto leaves. The coupled PALM-VIDA model demonstrates a significant reduction in PM levels across the modelling domain and within street canyons when deposition to vegetation is incorporated. Reduction in NO₂ through deposition to vegetation is lower but human NO2 exposure can still be decreased if tree species selection and placement leads to desirable effects on air flow. Sparse tree arrangements or species with sparse crowns facilitate ventilation and are often better at reducing NO₂ concentrations in street canyons compared to denser vegetation with higher deposition but negative effects on ventilation. Our study informs urban planning and green infrastructure design, underscoring the multifaceted role of urban greenery in air pollution mitigation strategies. Its main conclusion is that both deposition processes and the influence of air mixing and ventilation need to be considered to accurately assess the effects of urban trees on local air quality. Ill-considered placement and species selection may cause a net increase in pollutants underneath the trees. However, careful planning can address this risk and instead improve overall air quality.

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.

Building-block-flow computational model for large-eddy simulation of external aerodynamic applications

Computational fluid dynamics is an essential tool for accelerating the discovery and adoption of transformative designs across multiple engineering disciplines. Despite its many successes, no single approach consistently achieves high accuracy for all flow phenomena of interest, primarily due to limitations in the modeling assumptions. Here, we introduce a closure model for wall-modeled large-eddy simulation to address this challenge. The model, referred to as the Building-block Flow Model (BFM), rests on the premise that a finite collection of simple flows encapsulates the essential missing physics necessary to predict more complex scenarios. The BFM is designed to: (1) predict multiple flow regimes, (2) unify the closure model at solid boundaries and the rest of the flow, (3) ensure consistency with numerical schemes and gridding strategies by accounting for numerical errors, (4) be directly applicable to arbitrary complex geometries, and (5) be scalable to model additional flow physics in the future. The BFM is utilized to predict key quantities in five cases, including an aircraft in landing configuration, demonstrating similar or superior capabilities compared to previous state-of-the-art models. The design of BFM opens up new opportunities for developing closure models that can accurately represent various flow physics across different scenarios.