Search NASASearch

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.

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

At least 19 records

Turbulent flow simulation - A large eddy simulator's viewpoint

The state of the art of the data base on turbulent structures and numerical turbulence modelling is assessed. The failure of statistical methods to produce accurate results revealed that turbulence is more than simply random noise. Both randomness and coherence may be characteristics of turbulence, which may also feature noise that is a remnant of previous coherence. Five methods are presently used for numerically predicting turbulent flow behavior: correlations, integral methods, the k-epsilon model, and phenomenological and exact solutions to the Navier-Stokes equations. The higher level models include large eddy simulation and full simulation with numerical solutions to the Navier-Stokes equations; the two methods are applicable to experimental, rather than design work.

Ferziger, J. H.

A Nonlinear Dynamic Subscale Model for Partially Resolved Numerical Simulation (PRNS)/Very Large Eddy Simulation (VLES) of Internal Non-Reacting Flows

A brief introduction of the temporal filter based partially resolved numerical simulation/very large eddy simulation approach (PRNS/VLES) and its distinct features are presented. A nonlinear dynamic subscale model and its advantages over the linear subscale eddy viscosity model are described. In addition, a guideline for conducting a PRNS/VLES simulation is provided. Results are presented for three turbulent internal flows. The first one is the turbulent pipe flow at low and high Reynolds numbers to illustrate the basic features of PRNS/VLES; the second one is the swirling turbulent flow in a LM6000 single injector to further demonstrate the differences in the calculated flow fields resulting from the nonlinear model versus the pure eddy viscosity model; the third one is a more complex turbulent flow generated in a single-element lean direct injection (LDI) combustor, the calculated result has demonstrated that the current PRNS/VLES approach is capable of capturing the dynamically important, unsteady turbulent structures while using a relatively coarse grid.

Shih, Tsan-Hsing

Large-eddy simulation of an atmospheric bore and associated gravity wave effects on wind farm performance in the southern Great Plains

Gravity waves are a common occurrence in the atmosphere, with a variety of generation mechanisms. Their impact on wind farms has only recently gained attention, with most studies focused on wind farm-induced gravity waves. In this study, the interaction between a wind farm and gravity waves generated by an atmospheric bore event is assessed using multiscale large-eddy simulations. The atmospheric bore is created by a thunderstorm downdraft from a nocturnal mesoscale convective system (MCS). The associated gravity waves impact the wind resource and power production at a nearby wind farm during the American Wake Experiment (AWAKEN) in the US southern Great Plains. A two-domain nested setup (Δx=300 and 20 m) is used in the Weather Research and Forecasting (WRF) model, forced with data from the High-Resolution Rapid Refresh model, to capture both the formation of the bore and its interaction with individual wind turbines. The MCS is resolved on the large outer domain, where the structure of the bore and the associated gravity waves are found to be especially sensitive to parameterized microphysics processes. On the finer inner domain, gravity wave interactions with individual wind turbines are resolved; wake dynamics are captured using a generalized actuator disk parameterization in WRF. The gravity waves are found to have a strong effect on the atmosphere above the wind farm; however, the effect of the waves is more nuanced closer to the surface where there is additional turbulence, both ambient and wake-generated. Notably, the gravity waves modulate the mesoscale environment by weakening and dissipating the preexisting low-level jet, which reduces hub-height wind speed and hence the simulated power output, which is confirmed by the observed supervisory control and data acquisition (SCADA) power data. Additionally, the gravity waves induce local wind direction variations correlated with fluctuations in pressure, which lead to fluctuations in the simulated power output as various turbines within the farm are subjected to waking from nearby turbines.

17 WIND ENERGY

Collaborative Proposal: Improving understanding of the internal structure and dynamics of deep convection using ARM observations and large eddy simulations

Recent observational and large eddy simulation (LES) modeling studies have nearly unanimously supported the view of deep cumulus convection being composed of a series of quasi-spherical bubbles of buoyant air, known as moist thermals. Despite the prevalence of moist thermals in deep convection, a comprehensive theory for the dynamics of these structures is lacking. Most current conceptual models for cumulus convection are based on canonical scaling theories for dry thermals or plumes; however, there is considerable evidence that the behavior of moist thermals differs markedly from these theories. Furthermore, the theoretical basis for most cumulus parameterizations originates from the plume conceptual model, and therefore these parameterizations are inconsistent with the real structure of moist convection. Motivated by the aforementioned knowledge gaps, this “end-to-end” research effort use theory, observations, numerical simulations, and direct improvements to the Zhang-McFarlane (ZM) convection scheme in the global climate Community Atmosphere Model (CAM) to address the following research questions: What key environmental parameters determine whether or not shallow convection will transition into deep convection, in the context of thermal-like updrafts? What factors regulate the size of thermals within cumulus updrafts? How does vertical wind shear influence thermal behavior, and as a consequence, vertical velocity and mass flux profiles and the shallow-to-deep convective transition? What are the critical processes that determine updraft vertical velocities and their connection to the vertical mass flux profile for thermal-like updrafts? Idealized LES modeling will be used in conjunction with theoretical models for the core properties of thermal-like updrafts to better understand key processes that regulate thermal ascent rates and entrainment properties. Thermal-tracking procedures will be used to characterize the behavior of thermals within the LES, and recently developed direct measures of entrainment and detrainment will be used to quantify entrainment/detrainment rates. Building from these results, we will analyze the structure of moist thermals from hemispheric range-height indicator scans taken during the Atmospheric Radiation Measurement Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign, and from “real case” LES of CACTI events. This combined modeling and observational analysis will provide essential validation for the existing body of research on moist thermal dynamics, which is based primarily on modeling studies. With the insight gained from the aforementioned activities, we will modify the Zhang-McFarlane convection scheme to improve its representation of updraft vertical velocity and entrainment rate profiles. These process-level changes will be tested in the Community Atmosphere Model to assess the impact on global climate simulations.

54 ENVIRONMENTAL SCIENCES

Transonic Lift and Drag Predictions using Wall Modelled Large Eddy Simulations

Wall-modelled Large Eddy Simulations (WMLES) of the NASA Common Research Model (CRM) at transonic conditions and various angles of attacks leading up to and including shock-induced flow separation are performed using the LAVA computational framework. The simulations are shown to accurately predict the lift curve slope and the onset of separation characterized by the break in the pitching moment. Furthermore, careful assessment of skin-friction drag at cruise condition and its subsequent decrease with increasing angles of attack are shown to be in agreement with viscous sublayer resolving Reynolds Averaged Navier Stokes(RANS) simulations. The small differences between WMLES and RANS appear to be of the same order as differences seen between two RANS models at the cruise-point. Some sensitivity is observed to the coefficient used in the subgrid scale model, although this can be reconciled by noting the low chord- and shock-incidence Reynolds numbers (Rex≈106) seen in the outboard regions of the wing, along with the uncertainties associated with tripping and the numerical transition that occurs near the leading edge. Among the primary drawbacks seen in the WMLES predictions is the subdued increase in wave-drag with increasing angle-of-attack and the shock intensity when compared with experimental data. Preliminary assessment suggests that additional span- and stream-wise grid resolution is likely needed in mid- and out-board portions of the wing to better resolve the shock-induced separated flow and to further investigate prediction accuracy of unsteady temporal characteristics of the problem.

ARMD

Analysis of subgrid models using direct and large-eddy simulations of isotropic turbulence

Direct and large eddy simulations of forced and decaying isotropic turbulence have been performed using a pseudospectral and a finite-difference code. Subgrid models that include a one-equation subgrid kinetic energy model with and without a stochastic backscatter forcing term and a new scale similarity model have been analyzed in both Fourier space and physical space. The Fourier space analysis showed that the energy transfer across the cutoff wavenumber k(sub c) is dominated by local interaction. The correlation between the exact and the modeled (by a spectral eddy viscosity) nonlinear terms and the subgrid energy transfer in physical space was found to be quite low. In physical space, a similar correlation analysis was carried out using top hat filtering. Results show that the subgrid stress and the energy flux predicted by the subgrid models correlates very well with the exact data. The scale similarity model showed very high correlation for reasonable grid resolution. However, with decrease in grid resolution, the scale similarity model became more uncorrelated, when compared to the kinetic energy subgrid model. The subgrid models were then used for large-eddy simulations for a range of Reynolds number. It was determined that the dissipation was modeled poorly and that the correlation with the exact results was quite low for all the models. In general, for coarse grid resolution, the scale similarity model consistently showed very low correlation while the kinetic energy model showed a relatively higher correlation. These results suggest that to use the scale similarity model relatively fine grid resolution may be required, whereas, the kinetic energy model could be used even in coarse grid.

Menon, S.

Statistical Ensemble of Large Eddy Simulations

A statistical ensemble of large eddy simulations (LES) is run simultaneously for the same flow. The information provided by the different large scale velocity fields is used to propose an ensemble averaged version of the dynamic model. This produces local model parameters that only depend on the statistical properties of the flow. An important property of the ensemble averaged dynamic procedure is that it does not require any spatial averaging and can thus be used in fully inhomogeneous flows. Also, the ensemble of LES's provides statistics of the large scale velocity that can be used for building new models for the subgrid-scale stress tensor. The ensemble averaged dynamic procedure has been implemented with various models for three flows: decaying isotropic turbulence, forced isotropic turbulence, and the time developing plane wake. It is found that the results are almost independent of the number of LES's in the statistical ensemble provided that the ensemble contains at least 16 realizations.

Carati, Daniele

Investigating the Impact of Surface Heterogeneity on the Convective Boundary Layer Over Urban Areas Through Coupled Large-Eddy Simulation and Remote Sensing

Large-eddy Simulation (LES) was used to study convective boundary layer (CBL) flow through suburban regions with both large and small scale heterogeneities in surface temperature. Constant remotely sensed surface temperatures were applied at the surface boundary at resolutions of 10 m, 90 m, 200 m, and 1 km. Increasing the surface resolution from 1 km to 200 m had the most significant impact on the mean and turbulent flow characteristics as the larger scale heterogeneities became resolved. While previous studies concluded that scales of heterogeneity much smaller than the CBL inversion height have little impact on the CBL characteristics, we found that further increasing the surface resolution (resolving smaller scale heterogeneities) results in an increase in mean surface heat flux, thermal blending height, and potential temperature profile. The results of this study will help to better inform sub-grid parameterization for meso-scale meteorological models. The simulation tool developed through this study (combining LES and high resolution remotely sensed surface conditions) is a significant step towards future studies on the micro-scale meteorology in urban areas.

Dominguez, Anthony

Time-Domain Filtering for Spatial Large-Eddy Simulation

An approach to large-eddy simulation (LES) is developed whose subgrid-scale model incorporates filtering in the time domain, in contrast to conventional approaches, which exploit spatial filtering. The method is demonstrated in the simulation of a heated, compressible, axisymmetric jet, and results are compared with those obtained from fully resolved direct numerical simulation. The present approach was, in fact, motivated by the jet-flow problem and the desire to manipulate the flow by localized (point) sources for the purposes of noise suppression. Time-domain filtering appears to be more consistent with the modeling of point sources; moreover, time-domain filtering may resolve some fundamental inconsistencies associated with conventional space-filtered LES approaches.

Pruett, C. David

Predictions of LAGOON Nose Landing Gear Flow and Noise Using Wall-Modeled Large-Eddy Simulations

Wall-modeled large-eddy simulations (WMLESs) of the LAGOON nose landing gear are conducted with compressible Navier–Stokes equations and immersed boundary technique using the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework. The simulations are conducted using six different Cartesian octree meshes for the grid sensitivity analysis of the near-field and far-field numerical predictions, where the far-field noise results are computed with the Ffowcs Williams–Hawkings acoustic analogy. The effects of numerical tripping induced at the exact locations of the tripping devices in the experiments are also examined. In general, better comparison with the experimental results are shown for the the near-field results obtained with the simulations under the effects of numerical tripping. The effects of tripping are not significant on the far-field noise calculations and the results have reasonable comparison with the experimental data in the low and medium frequency ranges when an impermeable formulation of the acoustic analogy is used.

CST

Large eddy simulation - A predictive approach to turbulent flow computation

The filtering approach (Leonard, 1974) to large eddy simulation is described. Large eddy simulation is then surveyed with reference to subgrid scale models, numerical methods, channel flow, and free shear flows. It is concluded that the premise of LES (i.e., that the small eddies are much easier to model than the large ones) seems to be borne out by computational experience to date. The greatest obstacle in LES at present is the lack of computer power which effectively denies the ability to simulate all the significant eddies in flows of interest.

Ferziger, J. H.

Effects of Eddy Viscosity on Time Correlations in Large Eddy Simulation

Subgrid-scale (SGS) models for large. eddy simulation (LES) have generally been evaluated by their ability to predict single-time statistics of turbulent flows such as kinetic energy and Reynolds stresses. Recent application- of large eddy simulation to the evaluation of sound sources in turbulent flows, a problem in which time, correlations determine the frequency distribution of acoustic radiation, suggest that subgrid models should also be evaluated by their ability to predict time correlations in turbulent flows. This paper compares the two-point, two-time Eulerian velocity correlation evaluated from direct numerical simulation (DNS) with that evaluated from LES, using a spectral eddy viscosity, for isotropic homogeneous turbulence. It is found that the LES fields are too coherent, in the sense that their time correlations decay more slowly than the corresponding time. correlations in the DNS fields. This observation is confirmed by theoretical estimates of time correlations using the Taylor expansion technique. Tile reason for the slower decay is that the eddy viscosity does not include the random backscatter, which decorrelates fluid motion at large scales. An effective eddy viscosity associated with time correlations is formulated, to which the eddy viscosity associated with energy transfer is a leading order approximation.

He, Guowei

High-Order Finite-Difference Large-Eddy Simulations of a Supersonic Nozzle

Large-eddy simulations of the flow emanating from supersonic nozzles were carried out for the 6th Propulsion Aerodynamics Workshop to be held at the AIAA SciTech Conference. The purpose of the work is to detail the ability of a high-order of accuracy finite-difference based numerical scheme to accurately compute this flow. Results for both a perfectly expanded and over-expanded flow case are presented. Several complicating factors are discussed: 1) the limitations on the computational grid imposed by the numerical method, 2) the effect of additional dissipation necessary to resolve shock waves and 3) the effect of the nozzle boundary layer on the flow. Results indicate that the simulations replicate the basic flow physics, but there is room for improvement in accuracy. Areas for future study are outlined.

turbulence

Large Eddy Simulation of a Turbulent Jet

Here we present the results of a Large Eddy Simulation of a non-buoyant jet issuing from a circular orifice in a wall, and developing in neutral surroundings. The effects of the subgrid scales on the large eddies have been modeled with the dynamic large eddy simulation model applied to the fully 3D domain in spherical coordinates. The simulation captures the unsteady motions of the large-scales within the jet as well as the laminar motions in the entrainment region surrounding the jet. The computed time-averaged statistics (mean velocity, concentration, and turbulence parameters) compare well with laboratory data without invoking an empirical entrainment coefficient as employed by line integral models. The use of the large eddy simulation technique allows examination of unsteady and inhomogeneous features such as the evolution of eddies and the details of the entrainment process.

Webb, A. T.

A posterirori study of models for large eddy simulations of drop-laden flows

Large Eddy Simulation (LES) is conducted of a three-dimensional temporal mixing layer whose stream is initially laden with liquid drops which may evaporate during the simulation. The gas-phase equations are written in Eulerian frame for two perfect gas species (carrier gas and vapor emanating from the drops), while the liquid-phase equations are written in a Lagrangian frame.

spray modeling

A High-Resolution Capability for Large-Eddy Simulation of Jet Flows

A large-eddy simulation (LES) code that utilizes high-resolution numerical schemes is described and applied to a compressible jet flow. The code is written in a general manner such that the accuracy/resolution of the simulation can be selected by the user. Time discretization is performed using a family of low-dispersion Runge-Kutta schemes, selectable from first- to fourth-order. Spatial discretization is performed using central differencing schemes. Both standard schemes, second- to twelfth-order (3 to 13 point stencils) and Dispersion Relation Preserving schemes from 7 to 13 point stencils are available. The code is written in Fortran 90 and uses hybrid MPI/OpenMP parallelization. The code is applied to the simulation of a Mach 0.9 jet flow. Four-stage third-order Runge-Kutta time stepping and the 13 point DRP spatial discretization scheme of Bogey and Bailly are used. The high resolution numerics used allows for the use of relatively sparse grids. Three levels of grid resolution are examined, 3.5, 6.5, and 9.2 million points. Mean flow, first-order turbulent statistics and turbulent spectra are reported. Good agreement with experimental data for mean flow and first-order turbulent statistics is shown.

DeBonis, James R.