Results for “Large eddy simulation” — Page 21

NASA NTRSDate not supplied

Wall-modeled LES of Flow Over a Gaussian Bump

We perform wall-modeled large eddy simulations (WMLES) of turbulent flow over a Gaussian-shaped bump geometry, to assess its performance in the accelerating and separation regions of the flow. The flow conditions are based on the ongoing CFD validation experiments of Slotnick [NATO STO-MP-AVT-307, 2019]. The oncoming flow Mach number is 0.176, and two Reynolds numbers are simulated that are about 10000 and 36000 based on boundary layer thickness upstream of the bump. Preliminary Reynolds-averaged Navier-Stokes simulations are first performed to assess the effects of Mach number, Reynolds number, tunnel top and sidewall effects. Finally, WMLES results with an equilibrium wall model will be presented at two Reynolds numbers to assess their performance for this flow by making detailed comparisons with available experimental and higher-fidelity numerical data.

turbulence, wall model, large eddy simulation↗

NASA NTRS1996

Planetary Boundary Layer Simulation Using TASS

Boundary conditions to an existing large-eddy simulation model have been changed in order to simulate turbulence in the atmospheric boundary layer. Several options are now available, including the use of a surface energy balance. In addition, we compare convective boundary layer simulations with the Wangara and Minnesota field experiments as well as with other model results. We find excellent agreement of modelled mean profiles of wind and temperature with observations and good agreement for velocity variances. Neutral boundary simulation results are compared with theory and with previously used models. Agreement with theory is reasonable, while agreement with previous models is excellent.

Schowalter, David G.↗

NASA NTRS2017

Prediction of Turbulent Temperature Fluctuations in Hot Jets

Large-eddy simulations were used to investigate turbulent temperature fluctuations and turbulent heat flux in hot jets. A high-resolution finite-difference Navier-Stokes solver, WRLES, was used to compute the flow from a 2-inch round nozzle. Several different flow conditions, consisting of different jet Mach numbers and temperature ratios, were examined. Predictions of mean and fluctuating velocities were compared to previously obtained particle image velocimetry data. Predictions of mean and fluctuating temperature were compared to new data obtained using Raman spectroscopy. Based on the good agreement with experimental data for the individual quantities, the combined quantity turbulent heat flux was examined.

large-eddy simulations↗

NASA NTRS2017

Prediction of Turbulent Temperature Fluctuations in Hot Jets

Large-eddy simulations were used to investigate turbulent temperature fluctuations and turbulent heat flux in hot jets. A high-resolution finite-difference Navier-Stokes solver, WRLES, was used to compute the flow from a 2-inch round nozzle. Several different flow conditions, consisting of different jet Mach numbers and temperature ratios, were examined. Predictions of mean and fluctuating velocities were compared to previously obtained particle image velocimetry data. Predictions of mean and fluctuating temperature were compared to new data obtained using Raman spectroscopy. Based on the good agreement with experimental data for the individual quantities, the combined quantity turbulent heat flux was examined.

hot jets↗

NASA NTRSDate not supplied

Wall-modeled LES of Turbulent Flow Over a Two-dimensional Gaussian Bump

We perform wall-modeled large eddy simulations (WMLES) of turbulent flow over a nominally two-dimensional Gaussian-shaped bump geometry to assess its performance in the accelerating and separation regions of the flow. The flow conditions are based on high-fidelity numerical simulation of Uzun & Malik (AIAA Journal 2022). The oncoming flow Mach number is 0.2, with the bump length-based Reynolds number (ReL) of 2 million. In our previous study, while WMLES with the constant coefficient Vreman subgrid scale model and equilibrium wall model performed satisfactorily at lower ReL = 1 million, it failed to predict flow separation at the higher Re, contrary to the experimental observations. We investigate the sensitivity of WMLES to different subgrid scale models, wall models, and grid resolution and topology on flow separation by comparing with available data.

Computational Fluid Dynamics↗

NASA NTRS2017

Prediction of Turbulent Temperature Fluctuations in Hot Jets

Large-eddy simulations (LES) were used to investigate turbulent temperature fluctuations and turbulent heat flux in hot jets. A high-resolution finite-difference Navier-Stokes solver was used to compute the flow from a 2-inch round nozzle. Three different flow conditions of varying jet Mach numbers and temperature ratios were examined. The LES results showed that the temperature field behaves similar to the velocity field, but with a more rapidly spreading mixing layer. Predictions of mean, mu-bar(sub i), and fluctuating, mu'(sub i), velocities were compared to particle image velocimetry data. Predictions of mean, T-bar, and fluctuating, T', temperature were compared to data obtained using Rayleigh scattering and Raman spectroscopy. Very good agreement with experimental data was demonstrated for the mean and fluctuating velocities. The LES correctly predicts the behavior of the turbulent temperature field, but over-predicts the levels of the fluctuations. The turbulent heat flux was examined and compared to Reynolds-averaged Navier-Stokes (RANS) results. The LES and RANS simulations produced very similar results for the radial heat flux. However, the axial heat flux obtained from the LES differed significantly from the RANS result in both structure and magnitude, indicating that the gradient diffusion type model in RANS is inadequate. Finally, the LES data was used to compute the turbulent Prandtl number and verify that a constant value of 0.7 used in the RANS models is a reasonable assumption.

hot jets↗

NASA NTRS2023

A Methodology for Estimating the Energy and Moisture Budget of the Convective Boundary Layer Using Continuous Ground-Based Infrared Spectrometer Observations

Land-atmosphere interactions play a critical role in both the atmospheric water and energy cycles. Changes in soil moisture and vegetation alter the partitioning of surface water and energy fluxes, influencing diurnal evolution of the planetary boundary layer (PBL). The mixing diagram framework has proven useful in understanding the evolution of the heat and moisture budget within the convective boundary layer (CBL). We demonstrate that observations from the Department of Energy Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site provide all of the needed inputs needed for the mixing diagram framework, allowing us to quantify the impact from the surface fluxes, advection, radiative heating, encroachment, and entrainment on the evolution of the CBL. Profiles of temperature and humidity retrieved from the ground-based infrared spectrometer (called the Atmospheric Emitted Radiance Interferometer, or AERI) are a critical component in this analysis. Large eddy simulation results demonstrate that mean mixed-layer values derived are shown to be critical to close the energy and moisture budgets. A novel approach demonstrated here is the use of network of AERIs and Doppler lidars to quantify the advective fluxes of heat and moisture. The framework enables the estimation of the entrainment fluxes as a residual, providing a way to observe the entrainment fluxes without using multiple lidar systems. Finally, the high temporal resolution of the AERI observations enable the morning, midday, and afternoon evolution of the CBL to be quantified. This work provides a new way to use observations in this framework to evaluate weather and climate models.

Boundary layer↗

NASA NTRSDate not supplied

WMLES of Boundary-Layer Transition

The potential for the accurate prediction of laminar-turbulent transition by Wall-Modeled Large-Eddy Simulation (WMLES) is investigated through analysis of Direct Numerical Simulation (DNS) and WMLES results. The analysis indicates satisfactory prediction of transition by WMLES is possible, provided a) a sufficiently accurate wall model for the initial stage of transition is employed, b) a satisfactory transition is made between this initial-stage wall model and the wall model for the downstream fully-developed turbulence, and c) limitations on temporal and spatial resolution, and on the thickness of the wall-modeled region, imposed by the need to accurately resolve the disturbance waves and the transition region are observed. This is illustrated with a WMLES computation that yields reasonable agreement with DNS results; however, the wall-modeled region thickness, and thus the cost savings, is less than is typical for WMLES of fully developed turbulence.

wall-modeled large-eddy simulation, boundary-layer↗

NASA NTRSDate not supplied

Wall-modeled LES of Flow Over a Gaussian Bump

We perform wall-modeled large eddy simulations (WMLES) of turbulent flow over a Gaussian-shaped bump geometry, to assess its performance in the acceleration and separation regions of the flow. An unstructured finite-volume solver is used along with an equilibrium wall model. The Mach number of the oncoming flow selected for the simulation is 0.176, and two Reynolds numbers (Re) are simulated that are about 10000 and 36000 based on boundary layer thickness upstream of the bump. Spanwise periodic results for the lower Re case are compared with the available DNS data, while the full three-dimensional simulation results for the higher Re case are compared to available experimental data. Sensitivity of WMLES results are assessed for a number of factors including grid resolution, wall model exchange location, type of wall model, and unstructured grid topology. The WMLES results for the lower Re case agree well with available DNS data in terms of the wall pressure variation, velocity and turbulent stress profile comparisons. The skin-friction predictions show a reduced tendency to separate. This appears to be accentuated for the higher Re case, for the medium resolution grid used in this study, where WMLES does not predict any flow separation in the mean, which is in contrast to the large separation bubble observed in the experiments.

turbulence↗

NASA NTRSDate not supplied

WMLES of Boundary-Layer Transition

An analysis of the potential for the accurate prediction of laminar-turbulent transition by Wall-Modeled Large-Eddy Simulation (WMLES) indicates that satisfactory prediction of transition by WMLES is possible, provided a) a sufficiently accurate wall model for the initial stage of transition is employed, b) a satisfactory transition is made between this initial-stage wall model and the wall model for the downstream fully-developed turbulence, and c) limitations on temporal and spatial resolution, and on the thickness of the wall-modeled region, imposed by the need to accurately resolve the disturbance waves and the transition region are observed. This is illustrated with WMLES computations of subharmonic transition in a flat-plate boundary layer that yield reasonable agreement with DNS results. However, the wall-modeled region thickness, and thus the cost savings, is less than is typical for WMLES of fully developed turbulence.

wall-modeled large-eddy simulation↗

NASA NTRSDate not supplied

WMLES of Boundary-Layer Transition

An analysis of the potential for the accurate prediction of laminar-turbulent transition by Wall-Modeled Large-Eddy Simulation (WMLES) indicates that satisfactory prediction of transition by WMLES is possible, provided a) a sufficiently accurate wall model for the initial stage of transition is employed, b) a satisfactory transition is made between this initial-stage wall model and the wall model for the downstream fully-developed turbulence, and c) limitations on temporal and spatial resolution, and on the thickness of the wall-modeled region, imposed by the need to accurately resolve the disturbance waves and the transition region are observed. This is illustrated with WMLES computations of subharmonic transition in a flat-plate boundary layer that yield reasonable agreement with DNS results. However, the wall-modeled region thickness, and thus the cost savings, is less than is typical for WMLES of fully developed turbulence.

Wall-modeled large-eddy simulation↗

NASA NTRS2023

Numerical Study on the Aerodynamics of an Iced Airfoil With Scale-Resolving Simulations

Scale-resolving simulations of the NACA 23012 airfoil with horn ice accretion on the leading edge are conducted using the hybrid Reynolds-averaged Navier–Stokes/large-eddy simulation (hybrid RANS/LES) and wall-modeled large-eddy simulation (WMLES) approaches implemented in the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework. Aerodynamic results at the Reynolds number of 1.8 million show good comparison with the experimental measurements at different angles of attack from pre-stall to post-stall regimes. The pressure plateaus caused by the flow separation and the recovery of pressure inside the separation bubble around the iced leading edge are well predicted with the scale-resolving simulations when sufficient grid resolution is used around the accreted ice. The unsteadiness of the turbulent flows around the iced airfoil is also examined through the turbulent kinetic energy with the Reynolds normal stress anisotropy. Kelvin–Helmholtz instability (KHI) arises at the shear layer triggered by the upper ice horn and leads to rapid laminar-to-turbulent transition over a large range of angle of attack. With the increase of the angle of attack, the region with high turbulence intensity induced by the unstable shear layer spreads quickly over the entire upper surface of the airfoil. The coherent KHI modes from the upper and lower ice horns are extracted using the spectral proper orthogonal decomposition (SPOD) technique. The SPOD modes extracted from the upper shear layer have large-scale variations in the spanwise direction and low-rank behavior where the energy of the leading SPOD mode at each Strouhal number of the KHI largely represents the total energy when the mode number in the spanwise direction is small.

ARMD↗

NASA NTRS2023

Numerical Study on the Aerodynamics of an Iced Airfoil With Scale-Resolving Simulations

Scale-resolving simulations of the NACA 23012 airfoil with horn ice accretion on the leading edge are conducted using the hybrid Reynolds-averaged Navier–Stokes/large-eddy simulation (hybrid RANS/LES) and wall-modeled large-eddy simulation (WMLES) approaches implemented in the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework. Aerodynamic results at the Reynolds number of 1.8 million show good comparison with the experimental measurements at different angles of attack from pre-stall to post-stall regimes. The pressure plateaus caused by the flow separation and the recovery of pressure inside the separation bubble around the iced leading edge are well predicted with the scale-resolving simulations when sufficient grid resolution is used around the accreted ice. The unsteadiness of the turbulent flows around the iced airfoil is also examined through the turbulent kinetic energy with the Reynolds normal stress anisotropy. Kelvin–Helmholtz instability (KHI) arises at the shear layer triggered by the upper ice horn and leads to rapid laminar-to-turbulent transition over a large range of angle of attack. With the increase of the angle of attack, the region with high turbulence intensity induced by the unstable shear layer spreads quickly over the entire upper surface of the airfoil. The coherent KHI modes from the upper and lower ice horns are extracted using the spectral proper orthogonal decomposition (SPOD) technique. The SPOD modes extracted from the upper shear layer have large-scale variations in the spanwise direction and low-rank behavior where the energy of the leading SPOD mode at each Strouhal number of the KHI largely represents the total energy when the mode number in the spanwise direction is small.

ARMD↗

NASA NTRSDate not supplied

New Capabilities and Improvements to the High-Order Glenn Flux Reconstruction Code

The Glenn Flux Reconstruction (GFR) code is a computational fluid dynamics (CFD) code under development at NASA Glenn Research Center. GFR is based on the high-order flux reconstruction (FR) method and provides a large-eddy simulation (LES) capability that is both accurate and efficient for complex aeropropulsion flows. Three significant new capabilities have been added to the code that improve its performance and functionality. First, a variety of explicit Runge-Kutta methods, including some with adaptive time stepping, were added to GFR with two methods offering a 33% improvement in time-to-solution. Second, GFR can now utilize fully unstructured, mixed-element meshes to more easily facilitate the grid generation process for complex geometries. Finally, a rotating reference frame capability has been added to GFR for solving rotating turbomachinery problems. A selection of results demonstrating these new capabilities are presented in this work. The Taylor-Green vortex problem is used to verify the new unstructured capability by showing similar accuracy and resolution for all element types. LES of the Turbulent Heat Flux Phase III (THX3) experiment with comparison to another high-order LES code and a popular Reynolds-averaged Navier-Stokes (RANS) code demonstrates the accuracy of the code for complex aeropropulsion flows. Finally, LES of a spacecraft cabin ventilation fan shows the ability of GFR to efficiently establish a fan performance map and identify operating points for further analysis at high orders of accuracy.

High-Order Methods↗

NASA NTRSDate not supplied

Evaluation of Wall-Modeled Les for Flow Over A Multi-Element Airfoil

We perform wall-modeled large-eddy simulations (WMLES) of turbulent flow over a 30P30N multielement airfoil at various angles of attack to assess its ability to accurately predict lift using an unstructured node-based finite-volume methodology. The flow conditions are based on the experiments reported in Klausmeyer & Lin [NASA /TM-112858]. The oncoming flow Mach number is 0.2, and the stowed chord-based Reynolds number (𝑅𝑒𝑐) is 9 million. The angle of attack (𝛼) is varied between 8◦ and 23◦ with the expected stall angle of around 21◦. The simulated geometry is periodic in the span, thus any three-dimensional effects present in the experiment are ignored. We first analyze the predictions with baseline isotropic and anisotropic grids containing over 16 grid points per boundary-layer thickness. The near-wall streamwise/nominal spacing ratio is 0.4 at the leading/trailing edges of the slat, main element and flap, and 1 elsewhere, and the first wall-normal/nominal spacing ratio is 0.375 over most of the airfoil for the anisotropic grid, while these are unity for the isotropic grid. It is found that the predicted lift and wall pressures, near-wall velocity profiles, and the location of onset of resolved turbulence are sensitive to the near-wall grid anisotropy. Further sensitivities were explored at lower angles of attack by independently varying the streamwise and wall-normal grid anisotropies. The effects of both off-wall and near-wall grid refinement were also quantified. For the grid resolutions used here, our results suggest that the first wall-normal spacing, which coincides with the WMLES exchange location, has a dominant effect on the predictions for this flow configuration.

Computational Fluid Dynamics↗

NASA NTRS2021

Preconditioning of overcast-to-broken cloud transitions by riming in marine cold air outbreaks

Marine cold air outbreaks (CAOs) commonly form overcast cloud decks that transition into broken cloud fields downwind, dramatically altering the local radiation budget. In this study, we investigate the impact of frozen hydrometeors on these transitions. We focus on a CAO case in the NW Atlantic, the location of the multi-year flight campaign ACTIVATE (Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment). We use MERRA-2 (Modern-Era Retrospective analysis for Research and Applications, version 2) reanalysis fields to drive large eddy simulations with mixed-phase two-moment microphysics in a Lagrangian framework. We find that transitions are triggered by substantial rain (rainwater paths >25 g m−2), and only simulations that allow for aerosol depletion result in sustained breakups, as observed. Using a range of diagnostic ice nucleating particle concentrations, Ninp, we find that increasing ice progressively accelerates transitions, thus abbreviating the overcast state. Ice particles affect the cloud-topped boundary layer evolution, primarily through riming-related processes prior to substantial rain, leading to (1) a reduction in cloud liquid water, (2) early consumption of cloud condensation nuclei, and (3) early and light precipitation cooling and moistening below cloud. We refer to these three effects collectively as “preconditioning by riming”. Greater boundary layer aerosol concentrations available as cloud condensation nuclei (CCN) delay the onset of substantial rain. However, cloud breakup and low CCN concentration final stages are found to be inevitable in this case, due, primarily, to liquid water path buildup. An ice-modulated cloud transition speed suggests the possibility of a negative cloud–climate feedback. To address prevailing uncertainties in the model representation of mixed-phase processes, the magnitude of ice formation and riming impacts and, thereby, the strength of an associated negative cloud–climate feedback process, requires further observational evaluation by targeting riming hot spots with in situ imaging probes that allow for both the characterization of ice particles and abundance of supercooled droplets.

Marine cold air outbreaks↗

NASA NTRS1998

Two Dimensional Wake Vortex Simulations in the Atmosphere: Preliminary Sensitivity Studies

A numerical large-eddy simulation model is currently being used to quantify aircraft wake vortex behavior with meteorological observables. The model, having a meteorological framework, permits the interaction of wake vortices with environments characterized by crosswind shear, stratification, and humidity. The addition of grid-scale turbulence as an initial condition appeared to have little consequence. Results show that conventional nondimensionalizations work very well for vortex pairs embedded in stably stratified flows. However, this result is based on simple environments with constant Brunt-Vaisala frequency. Results presented here also show that crosswind profiles exert important and complex interactions on the trajectories of wake vortices. Nonlinear crosswind profiles tended to arrest the descent of wake vortex pairs. The member of the vortex pair with vorticity of same sign as the vertical change in the ambient along-track vorticity may be deflected upwards.

Proctor, F. H.↗

NASA NTRS2021

Nighttime Convection in Water-ice Clouds at High Northern Latitudes on Mars

We investigate water-ice clouds and their influence on the temperature structure of the Mar-tian atmosphere at high northern latitudes in early summer. New results are obtained throughcoordinated analysis of two types of data from Mars Global Surveyor: atmospheric profiles re-trieved from radio occultation (RO) measurements and wide-angle images from the Mars OrbiterCamera (MOC). Some RO profiles contain a layer of neutral static stability, which indicates thepresence of convective mixing at a local time (about 5 h) when it does not usually occur. Thesenocturnal mixed layers (NMLs) were observed frequently in early summer of Mars year 27 atlatitudes of 53–72◦N and longitudes of 210–330◦E. The base of a typical NML is 3 km abovethe surface, about the same height as the nighttime cloud layer detected by the Phoenix LIDARin early summer of Mars year 29 at 234◦E, 68◦N. The depth of the NMLs ranges from less than1 km to more than 5 km. Comparisons with nearly simultaneous MOC images demonstrate thatNMLs are closely associated with water-ice clouds. There is a dense cluster of NMLs within theannular cloud that appears every year in early summer between Alba Mons and the north polarresidual ice cap. The lighting conditions at this location and season allowed MOC to observe theannular cloud on most orbits, at 118-min intervals. Its appearance varies dramatically with localtime, becoming more symmetrical and better organized at night and dissipating to a crescentshape during the day. According to high-resolution numerical simulations (Spiga et al., 2017),including a large-eddy simulation at the Phoenix landing site, NMLs form when radiative cool-ing by water-ice aerosols causes convective instability; the mixed layer is forced from above bynegative buoyancy. Our results strongly support this conclusion. In addition, MOC images frommidsummer contain eastward-moving frontal clouds. Temperature profiles within these cloudsshow signs of near-surface advection of warm air, which reduces the static stability of the loweratmosphere and contributes, along with cloud radiation, to the formation of an NML.

David Hinson↗