SEARCH · Search NASA

Results for “Large eddy simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 883 records · Page 49

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Reconciling Chord Length Distributions and Area Distributions for Fields of Fractal Cumulus Clouds

While the total cover of broken cloud fields can in principle be obtained from one-dimensional measurements, the cloud size distribution normally differs between two-dimensional (area) and one-dimensional retrieval (chord length) methods. In this study, we use output from high-resolution Large Eddy Simulations to generate a transfer function between the two. We retrieve chord lengths and areas for many clouds, and plot the one as a function of the other, and vice versa. We find that the cloud area distribution conditional on the chord length behaves like a gamma distribution with well-behaved parameters, with a mean μ=1.1L and a shape parameter β=L−0.645. Using this information, we are able to generate a transfer function that can adjust the chord length distribution so that it comes much closer to the cloud area distribution. Our transfer function improves the error in predicting the mean cloud size, and is performs without strong biases for smaller sample sizes. However, we find that the method is still has difficulties in accurately predicting the frequency of occurrence of the largest cloud sizes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Process-oriented stochastic perturbations applied to the parametrization of turbulence and shallow convection for ensemble prediction

Weather forecasting nowadays often requires some estimation of uncertainties associated with the output of meteorological models, in order to better inform decision making, especially in the context of intense weather events. Ensemble prediction systems provide such information through sets of possible scenarios which are designed to represent various uncertainty sources, including model uncertainties. A wide variety of methods have been proposed to estimate model uncertainties, among which perturbation methods targeting uncertain processes are a promising research field. Here, we focus on the representation of small-scale variability by process-oriented perturbation schemes applied to two key physical processes, namely turbulence and shallow convection. The perturbations are applied to a single-column version of the convection-permitting AROME model, in three idealized boundary-layer cases. Large-eddy simulations (LESs) of the same cases serve as a reference for the subgrid variability that has to be represented, and the results are also compared to those given by the Stochastically Perturbed Parametrization Tendencies (SPPT) method, which is a method commonly used by weather forecast centres to represent model uncertainty. The spread produced by our process-oriented perturbations of turbulence and shallow convection does not represent all the small-scale variability implied by the LESs for temperature and humidity. However, it is of a similar order of magnitude for the wind, thanks to perturbations generated by the stochastic turbulence scheme. The dispersion is structurally different from what is obtained with SPPT. It is non-negligible in the lower levels, where SPPT perturbations are usually suppressed because of numerical instabilities, indicating a possible complementarity between the schemes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2019

Description of the ARM Large-Scale Forcing Data from the Constrained Variational Analysis (VARANAL) Version 2

This technical report represents an update of the previous technical report written by Zhang et al., which described the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility constrained variational analysis (VARANAL) that is used to develop the large-scale forcing data for driving single-column models (SCMs), cloud-resolving models (CRMs), and large-eddy simulation models (LES). The VARANAL algorithm was originally developed by Zhang and Lin and Zhang et al. at Stony Brook University and was migrated to the Lawrence Livermore National Laboratory (LLNL) as the ARM operational objective analysis system in May 1999. Since then, the algorithm has been evolved with time along with the availability of new observations and techniques to meet various modeling needs. Major updates include: 1. The method used to develop multi-year continuous forcing data, 2. The incorporation of eddy correlation flux measurement system (ECOR) turbulent fluxes into the analysis, and 3. Improvements to the workflow (e.g., implementing part of the code into the ARM Data Integrator [ADI]) to increase efficiency. The ARM large-scale forcing data have been widely used for SCM/CRM/LES to understand and improve physical processes in models. Zhang et al. has provided a comprehensive review of the SCM concept, early efforts to derive forcing data for SCM studies, efforts of the ARM constrained variational analysis, and previous SCM studies using ARM cases. This technical report focuses on the constrained variational analysis algorithm and the introduction of the ARM VARANAL products. We also extended the VARANAL algorithm into a three-dimensional constrained variational analysis (3DCVA) and designed an ensemble framework to address the forcing uncertainty. The 3D large-scale forcing data are released as another datastream named “varanal3d”. Please refer to the ARM varanal3d technical report for more information.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Analysis of Compressibility Effects and Nonlinear Property Variations in a Supercritical CO2 Mixing Layer

Calculations are performed using the large eddy simulation technique to qualitatively assess the effects of fluid compressibility and thermodynamic nonlinearities on the dynamics of a supercritical flow field. Resulting data are also used to interrogate instantaneous subfilter-scale fields, as modeled by the mixed-dynamic Smagorinsky closure and gradient diffusion model. A three-dimensional wall-resolved calculation of a spatially evolving mixing layer composed solely of carbon dioxide is performed. This is followed by a coarse grid calculation to examine model performance as a function of resolution. Theoretical analysis indicates that the partial derivatives of density with respect to pressure and temperature are important modulators of the temperature and pressure field evolution, respectively. Results pertaining to subfilter velocity-velocity and velocity-temperature fields indicate that the models are most stressed in the regions where compressibility and thermo-fluid nonlinearities dominate the physics.

42 ENGINEERING↗

DOE OSTI2023

Turbulent Kinetic Energy Budgets over Gentle Topography Covered by Forests

Large-eddy simulations of flow over a “horizontally” uniform model forest are used to investigate the effects of gentle topography on the turbulent kinetic energy (TKE) budget within the canopy roughness sublayer. Despite significant differences between simulations using idealized sinusoidal topography and real topography of the Amazon forest, results indicate that the effects of topography are located predominantly in the upper canopy and above, and are mostly caused by mean advection of TKE. The “horizontally” averaged TKE budget from idealized and real gentle topographies are almost identical to that for flat terrain, including a clear inertial layer above the roughness sublayer in which shear production is balanced by local dissipation. At topography crests, where observational towers are usually located, mean vertical advection of TKE can be as important as horizontal advection. Here, we propose the use on an approximate TKE balance equation to estimate mean advection from single tower measurements, and introduce a new advection index that can be used as a proxy to quantify the importance of the topography on the TKE budget.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Methods to Evaluate Subcolumn Profiles Based on Two-Point Diagnostics

In atmospheric models, stochastic generation of subgrid-scale profiles or “subcolumns” has been used for a variety of purposes. Such subcolumns can be generated from subgrid probability density functions (PDFs) at different vertical levels, when such PDFs are available. To do so, the generator needs to decide how strongly points should be correlated in the vertical, that is, how much the values should be overlapped. This is sometimes called “PDF overlap.” To assess vertical correlation in a simplified, observable setting, here the vertical correlation of vertical velocity in subcloud layers is examined. Doppler lidar is used to evaluate the vertical profiles of vertical velocity produced by a large-eddy simulation (LES) model and the Subgrid Importance Latin Hypercube Sampler (SILHS) subcolumn generator. In order to diagnose unrealistic features in subcolumn profiles, various statistical diagnostics are examined here, including the bivariate PDF of vertical velocity at two separated points (i.e., altitudes), the two-point velocity correlation, the integral correlation length, the PDF of two-point velocity differences, and the skewness and kurtosis of two-point velocity differences. The profiles produced by LES match lidar well, except that they are too smooth at small scales. The profiles produced by SILHS exhibit sharp jumps from updraft to downdraft that are not observed in the lidar data. To reduce the generation of these unrealistically sharp jumps, the SILHS sampling method is revised. The diagnostics confirm that the revised sampling method reduces the overprediction of sharp jumps.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Turbulence-induced bias in time-averaged laser absorption tomography of correlated concentration and temperature fields with a first-order correction

The influence of correlated scalar fluctuations on time-averaged laser absorption tomography measurements of temperature and species in a piloted turbulent premixed flame was examined using a coupled spectroscopic and fluid-dynamic analysis. To understand bias associated with turbulence, spatio-temporally resolved temperature and species mole fraction profiles predicted by large eddy simulations (LES) were used to synthetically generate time-resolved line-of-sight absorption measurements at short time scales (microsecond) to reflect the unsteady nature of a canonical jet burner across various transverse measurement planes. Inversion methods were employed on the time-averaged line-of-sight data to produce radially-resolved temperature and mole fraction profiles, analogous to those produced by laser absorption tomography performed on a time-averaged axisymmetric flowfield. It is shown that bias in the measurements compared to true time-averaged scalar fields is a function primarily of temperature dependence in absorptivity and non-zero correlation between temperature and species concentration scalars. Finally, a first-order correction to tomography measurements is proposed to account for the bias based on estimated correlations and the known spectroscopic parameters of the probed absorption transitions.

42 ENGINEERING↗

DOE OSTI2024

On the Departure from Monin–Obukhov Surface Similarity and Transition to the Convective Mixed Layer

Large-eddy simulations are used to evaluate mean profile similarity in the convective boundary layer (CBL). Particular care is taken regarding the grid sensitivity of the profiles and the mitigation of inertial oscillations in the simulation spin-up. The nondimensional gradients Φ for wind speed and air temperature generally align with Monin–Obukhov similarity across cases but have a steeper slope than predicted within each profile. The same trend has been noted in several other recent studies. The Businger-Dyer relations are modified here with an exponential cutoff term to account for the decay in Φ to first-order approximation, yielding improved similarity from approximately 0.05z i to above 0.3z i , where z i is the CBL depth. The necessity for the exponential correction is attributed to an extended transition from surface scaling to zero gradient in the mixed layer, where the departure from Monin–Obukhov similarity may be negligible at the surface but becomes substantial well below the conventional surface layer height of 0.1 z i .

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

The Effect of Surface Heating Heterogeneity on Boundary Layer Height and Its Dependence on Background Wind Speed

The planetary boundary layer height (PBLH) is a fundamental variable of the planetary boundary layer (PBL). The effect of surface heterogeneity is challenging in the PBL parameterization for numerical weather and climate models. Here we use large eddy simulation data to investigate the impact of surface heating heterogeneity on the PBLH spatial variations and mean values over a domain of typical mesoscale model grid size. It is found that this impact depends on background wind condition. Variance decomposition and Fourier spectra both reveal that for the fully developed sheared convective PBL, surface heterogeneity contributes little to the spatial variation of PBLH, except for the case with heterogeneity scale of 14.4 km. This suggests that surface heterogeneity with length scales less than ~10 km, which is typical for mesoscale modeling, does not induce significant subgrid spatial variation of PBLH. The domain-averaged PBLH is found to generally increase as surface heterogeneity scale decreases, until reaching the scale on the order of PBLH. The findings can be used to inform the parameterization development of PBL processes over heterogeneous surface.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

ASSESSING FUEL PROPERTY EFFECTS ON CAVITATION AND EROSION PROPENSITY USING A COMPUTATIONAL FUEL SCREENING TOOL

To advance compression ignition combustion strategies, researchers have evaluated fuel property effects and their impact on achieving higher efficiencies and lower emissions levels relative to current capabilities. Within the Department of Energy's Co-Optima initiative, there has been a recent focus on understanding the influence of fuel properties on fuel injection performance. To help identify candidate fuels that can meet desired injector performance metrics, a computational fuel screening tool is under development that can link fuel properties with the tendency of a given fuel to cavitate and lead to cavitation-induced erosion. In the initial development of this tool, five liquid fuel properties were selected to represent candidate fuels, namely density, viscosity, vapor pressure, surface tension, and heat of vaporization. A design of experiments methodology was employed to generate a set of pseudo fuel cases, which can represent the main effects and interactions among the fuel properties and be related to cavitation erosion predictions. Large eddy simulations were performed using a mixture modeling approach to predict the cavitation and erosion propensity of these pseudo-fuels in terms of the mean fuel vapor mole fraction and stored impact energy from repeated cloud collapse events. The low order regression models generated from this study revealed the importance of liquid fuel density on cavitation formation, whereas liquid viscosity was found to have a strong negative correlation with erosion severity. The surrogate models were then used in the fuel screening tool to rank the cavitation and erosion tendency of four candidate single-component fuels: methyl decanoate, isooctane, ethanol and n-dodecane. While the fuel screening tool was not able to quantitatively predict the cavitation and erosion response metrics, the tool was able to accurately rank the relative cavitation and erosion propensity of the four fuels. Overall, ethanol and iso-octane were observed to produce the highest levels of cavitation and erosion, respectively.

cavitation-induced erosion↗

DOE OSTI2019

Description of the ARM Large-Scale Forcing Data from the Constrained Variational Analysis (VARANAL) Version 2

This technical report represents an update of the previous technical report written by Zhang et al. (2001a) (available at http://www.arm.gov/publications/tech_reports/arm-tr-005.pdf), which described the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility constrained variational analysis (VARANAL) that is used to develop the large-scale forcing data for driving single-column models (SCMs), cloud-resolving models (CRMs), and large-eddy simulation models (LES). The VARANAL algorithm was originally developed by Zhang and Lin (1997) and Zhang et al. (2001b) at Stony Brook University and was migrated to the Lawrence Livermore National Laboratory (LLNL) as the ARM operational objective analysis system in May 1999. Since then, the algorithm has been evolved with time along with the availability of new observations and techniques to meet various modeling needs. Major updates include: (1) The method used to develop multi-year continuous forcing data (Xie et al. 2004), (2) The incorporation of eddy correlation flux measurement system (ECOR) turbulent fluxes into the analysis (Tang et al. 2019), and (3) Improvements to the workflow (e.g., implementing part of the code into the ARM Data Integrator [ADI]) to increase efficiency. The ARM large-scale forcing data have been widely used for SCM/CRM/LES to understand and improve physical processes in models. Zhang et al. (2016) has provided a comprehensive review of the SCM concept, early efforts to derive forcing data for SCM studies, efforts of the ARM constrained variational analysis, and previous SCM studies using ARM cases. This technical report focuses on the constrained variational analysis algorithm and the introduction of the ARM VARANAL products.

58 GEOSCIENCES↗

DOE OSTI2023

Effect of the integral length scales of turbulent inflows on wind turbine loads

As wind turbines become larger, the fluctuations in the inflow become increasingly influential in the turbine structural loading. These fluctuations are characterized by the integral length scale, which represents the average size of the largest energy-containing turbulent eddies. Current design standards neglect the varying integral length scales that characterize inflows of wind turbines in operation. Using large-eddy simulations, we generate turbulent inflows of varying integral length scales and quantify the loads of the IEA 15-MW reference wind turbine. Results illustrate that the impact of turbulence on rotor and tower loads is up to 10 times greater than the impact of the mean shear profile. Increasing integral length scales from 0.3x to 0.5x the rotor diameter reduces blade root flapwise moments and rotor and tower loads. Increasing integral length scales from 0.5x to 1.7x the rotor diameter increases the rotor aerodynamic thrust force and the blade root flapwise shear loads and decreases the tower base torsional moment and the tilting and yawing rotor aerodynamic moments. Additionally, turbulence intensity has a greater impact on wind turbine loads than integral length scales. In conclusion, findings indicate that design standards should consider varying integral length scales for accurate wind turbine loading characterization in turbulent inflow conditions.

17 WIND ENERGY↗

DOE OSTI2020

The effect of dynamic near-wake modulation on utility-scale wind turbine wake development

High-resolution field-scale experiments using flow visualization with natural snowfall and high-fidelity large eddy simulations are combined to investigate the effect of dynamic turbine operation and atmospheric conditions on wind turbine wake mixing and recovery in the wake of a 2.5 MW wind turbine. Instantaneous near-wake expansion and deflection in response to changes in blade pitch and wind direction, termed dynamic wake modulation, is quantified using both techniques, demonstrating excellent agreement. The simulations are used to extend these results by calculating the energy flux into the wake 7 rotor diameters downstream, showing that dynamic turbine-atmospheric interactions enhance mixing in the far-wake. This finding is exhibited under both uniform and turbulent inflow conditions. Under turbulent flow, a synergistic relationship is also observed between dynamic wake modulation and wake meandering, as wake recovery can be further accelerated when the two phenomena occur together. The results of this study have implications for the development of more realistic far-wake models that include the significant impact of dynamic wake modulation on wake mixing and development. Additionally, the findings from the current study can be used to develop advanced control algorithms to speed up wake breakdown and recovery, further improving wind farm efficiency.

17 WIND ENERGY↗

DOE OSTI2024

Evaluation of a New Approach for Entrainment and Detrainment Rate Estimation

Entrainment and detrainment rates (ε and δ) constitute the most critical free parameters in mass flux schemes commonly employed for cumulus parameterizations. Recently, Zhu et al. (2021) introduced a new approach that utilizes aircraft observations to simultaneously estimate ε and δ for cumulus clouds, overcoming the limitation of other observation-based approaches that solely yield ε without offering insights into δ. This study aims to comprehensively evaluate the reliability of this new approach. First, evaluation using an Explicit Mixing Parcel Model demonstrates the capability of the new approach to back-calculate predetermined ε and δ based on the physical properties before and after the entrainment mixing. Second, evaluation using large-eddy simulations illustrates that the new approach yields consistent ε and δ profiles compared to the traditional approach. Sensitivity tests indicate a weak sensitivity of the estimated δ with the new approach to the entrained air source. A decrease in the proportion of cloudy air in the assumed detrained air leads to a reduction in the estimated δ, while ε remains unaffected. Finally, the most appropriate assumptions for entrained and detrained air are discussed. Estimating ε for cumulus parameterizations involves acquiring ambient air more than 500 m away from the cloud edge as entrained air. Due to implicit mean field approximations in the traditional approach, determining the optimal assumption for detrained air properties proves challenging. Finally, this study confirms the reliability of the new approach in estimating ε and δ, providing confidence in its application to extensive observational data and advancement in parameterization.

54 ENVIRONMENTAL SCIENCES↗

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

ARMD↗

DOE OSTI2020

Changing the rotational direction of a wind turbine under veering inflow: a parameter study

All current-day wind-turbine blades rotate in clockwise direction as seen from an upstream perspective. The choice of the rotational direction impacts the wake if the wind profile changes direction with height. Here, we investigate the respective wakes for veering and backing winds in both hemispheres by means of large-eddy simulations. We quantify the sensitivity of the wake to the strength of the wind veer, the wind speed, and the rotational frequency of the rotor in the Northern Hemisphere. A veering wind in combination with counterclockwise-rotating blades results in a larger streamwise velocity output, a larger spanwise wake width, and a larger wake deflection angle at the same downwind distance in comparison to a clockwise-rotating turbine in the Northern Hemisphere. In the Southern Hemisphere, the same wake characteristics occur if the turbine rotates counterclockwise. These downwind differences in the wake result from the amplification or weakening or reversion of the spanwise wind component due to the effect of the superimposed vortex of the rotor rotation on the inflow's shear. An increase in the directional shear or the rotational frequency of the rotor under veering wind conditions increases the difference in the spanwise wake width and the wake deflection angle between clockwise- and counterclockwise-rotating actuators, whereas the wind speed lacks a significant impact.

17 WIND ENERGY↗

DOE OSTI2026

Inverse Mapping of the Collision Kernel and Wall Flux Scaling in a Tall Convection‐Cloud Chamber Using Local Sensors and Knowledge‐Informed Deep Learning

Droplet collision–coalescence is a crucial process in cloud physics, but accurately representing this process under different dynamical conditions remains challenging. A proposed future convective‐cloud chamber aims to investigate this key process, but the method for observing it remains unclear, even though it is theoretically established that collision‐coalescence will occur. This study serves as a proof‐of‐concept demonstration of how knowledge‐informed deep learning, combined with measurement data from local sensors in the chamber, can be used to estimate the collision kernels, which determine how the droplet size distribution evolves during collision‐coalescence. In addition to estimating the collision kernel, we also address wall fluxes, another uncertain but important process that acts as a source of heat and moisture in the chamber. Ensemble runs of large‐eddy simulations are conducted by scaling the wall fluxes and the collision kernel, while the measured flow and cloud properties are used as inputs for a neural network. Results indicate that this approach successfully maps the scaling of wall fluxes and the collision kernel with biases of approximately 1% or less relative to the range of the target data. This proof‐of‐concept lays the groundwork for future applications; when the real measurements are available, real sensor data combined with the trained model presented in this work will enable estimation of the actual wall fluxes and collision kernel.

cloud chamber↗