Search NASA⌕ Search

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 1,711 records · Page 95

Toward Simplification of Dynamic Subgrid-Scale Models

We examine the relationship between the filter and the subgrid-scale (SGS) model for large-eddy simulations, in general, and for those with dynamic SGS models, in particular. From a review of the literature, it would appear that many practitioners of LES consider the link between the filter and the model more or less as a formality of little practical effect. In contrast, we will show that the filter and the model are intimately linked, that the Smagorinsky SGS model is appropriate only for filters of first- or second-order, and that the Smagorinsky model is inconsistent with spectral filters. Moreover, the Germano identity is shown to be both problematic and unnecessary for the development of dynamic SGS models. Its use obscures the following fundamental realization: For a suitably chosen filter, the computible resolved turbulent stresses, property scaled, closely approximate the SGS stresses.

Pruett, C. David↗

On Taylor-Series Approximations of Residual Stress

Although subgrid-scale models of similarity type are insufficiently dissipative for practical applications to large-eddy simulation, in recently published a priori analyses, they perform remarkably well in the sense of correlating highly against exact residual stresses. Here, Taylor-series expansions of residual stress are exploited to explain the observed behavior and "success" of similarity models. Until very recently, little attention has been given to issues related to the convergence of such expansions. Here, we re-express the convergence criterion of Vasilyev [J. Comput. Phys., 146 (1998)] in terms of the transfer function and the wavenumber cutoff of the grid filter.

Pruett, C. David↗

Ocean Turbulence I: One-Point Closure Model Momentum and Heat Vertical Diffusivities

Since the early forties, one-point turbulence closure models have been the canonical tools used to describe turbulent flows in many fields. In geophysics, Mellor and Yamada applied such models using the 1980 state-of-the art. Since then, no improvements were introduced to alleviate two major difficulties: 1) closure of the pressure correlations, which affects the correct determination of the critical Richardson number Ri(sub cr) above which turbulent mixing is no longer possible and 2) the need to express the non-local third-order moments (TOM) in terms of lower order moments rather than via the down-gradient approximation as done thus far, since the latter seriously underestimates the TOMs. Since 1) and 2) are still being dealt with adjustable parameters which weaken the credibility of the models, alternative models, not based on turbulence modeling, have been suggested. The aim of this paper is to show that new information, partly derived from the newest 2-point closure model discussed, can be used to solve these shortcomings. The new one-point closure model, which in its simplest form is algebraic and thus simple to implement, is first shown to reproduce a variety of data. Then, it is used in a Ocean-General Circulation Model (O-GCM) where it reproduces well a large variety of ocean data. While phenomenological models are specifically tuned to ocean turbulence, the present model is not. It is first tested against laboratory data on stably stratified flows and then used in an O-GCM. It is more general, more predictive and more resilient, e.g., it can incorporate phenomena like wave-breaking at the surface, salinity diffusivity, non-locality, etc. One important feature that naturally comes out of the new model is that the predicted Richardson critical value Ri(sub cr) is Ri (sub cr approx. = 1) in agreement with both Large Eddy Simulations (LES) and empirical evidence while all previous models predicted Ri (sub cr approx. = 0.2) which led to a considerable underestimate of the extent of turbulent mixing and thus to an incorrect mixed layer depth. The predicted temperature and salinity profiles (vs. depth) are presented and compared with those of the Kolmogorov-Petruvsky-Piskunuv (KPP) model and Levitus data.

Canuto, V. M.↗

Wake Vortex Transport in Proximity to the Ground

A sensitivity study for aircraft wake vortex transport has been conducted using a validated large eddy simulation (LES) model. The study assumes neutrally stratified and nonturbulent environments and includes the consequences of the ground. The numerical results show that the nondimensional lateral transport is primarily influenced by the magnitude of the ambient crosswind and is insensitive to aircraft type. In most of the simulations, the ground effect extends the lateral position of the downwind vortex about one initial vortex spacing (b(sub o)) in the downstream direction. Further extension by as much as one b(sub o) occurs when the downwind vortex remains 'in ground effect' (IGE) for relatively long periods of time. Results also show that a layer-averaged ambient wind velocity can be used to bound the time for lateral transport of wake vortices to insure safe operations on a parallel runway.

Hamilton, David W.↗

A New Method for Accurate Treatment of Flow Equations in Cylindrical Coordinates Using Series Expansions

The motivation of this work is the ongoing effort at the Center for Turbulence Research (CTR) to use large eddy simulation (LES) techniques to calculate the noise radiated by jet engines. The focus on engine exhaust noise reduction is motivated by the fact that a significant reduction has been achieved over the last decade on the other main sources of acoustic emissions of jet engines, such as the fan and turbomachinery noise, which gives increased priority to jet noise. To be able to propose methods to reduce the jet noise based on results of numerical simulations, one first has to be able to accurately predict the spatio-temporal distribution of the noise sources in the jet. Though a great deal of understanding of the fundamental turbulence mechanisms in high-speed jets was obtained from direct numerical simulations (DNS) at low Reynolds numbers, LES seems to be the only realistic available tool to obtain the necessary near-field information that is required to estimate the acoustic radiation of the turbulent compressible engine exhaust jets. The quality of jet-noise predictions is determined by the accuracy of the numerical method that has to capture the wide range of pressure fluctuations associated with the turbulence in the jet and with the resulting radiated noise, and by the boundary condition treatment and the quality of the mesh. Higher Reynolds numbers and coarser grids put in turn a higher burden on the robustness and accuracy of the numerical method used in this kind of jet LES simulations. As these calculations are often done in cylindrical coordinates, one of the most important requirements for the numerical method is to provide a flow solution that is not contaminated by numerical artifacts. The coordinate singularity is known to be a source of such artifacts. In the present work we use 6th order Pade schemes in the non-periodic directions to discretize the full compressible flow equations. It turns out that the quality of jet-noise predictions using these schemes is especially sensitive to the type of equation treatment at the singularity axis. The objective of this work is to develop a generally applicable numerical method for treating the singularities present at the polar axis, which is particularly suitable for highly accurate finite-differences schemes (e.g., Pade schemes) on non-staggered grids. The main idea is to reinterpret the regularity conditions developed in the context of pseudo-spectral methods. A set of exact equations at the singularity axis is derived using the appropriate series expansions for the variables in the original set of equations. The present treatment of the equations preserves the same level of accuracy as for the interior scheme. We also want to point out the wider utility of the method, proposed here in the context of compressible flow equations, as its extension for incompressible flows or for any other set of equations that are solved on a non-staggered mesh in cylindrical coordinates with finite-differences schemes of various level of accuracy is straightforward. The robustness and accuracy of the proposed technique is assessed by comparing results from simulations of laminar forced-jets and turbulent compressible jets using LES with similar calculations in which the equations are solved in Cartesian coordinates at the polar axis, or in which the singularity is removed by employing a staggered mesh in the radial direction without a mesh point at r = 0.

Constantinescu, G.S.↗

GAP Noise Computation By The CE/SE Method

A typical gap noise problem is considered in this paper using the new space-time conservation element and solution element (CE/SE) method. Implementation of the computation is straightforward. No turbulence model, LES (large eddy simulation) or a preset boundary layer profile is used, yet the computed frequency agrees well with the experimental one.

Loh, Ching Y.↗

Icehouse Effect: A Selective Arctic Cooling Trend Current Models are Missing

The icehouse effect is a hypothesized climate feedback mechanism which could result in human-caused surface cooling trends in polar regions. Once understood in detail, it becomes apparent that these trends, which are discernable in the literature, but have been largely dismissed, do not conflict with the consensus assessment of the evidence, which infers century-scale Arctic warming. In fact, confirmation of the hypothesis would substantially strengthen the argument that there is a detectable human influence on today's climate. This apparent enigma is resolved only through careful attention to the detail of the hypothesis and the data supporting it. The posited surface cooling is entirely dependent on the existence of climate warming in layers capping the stable boundary layer. Also, the cooling is not pandemic, but is selective. It is readily revealed in properly sorted data by making use of the principles of micrometeorological similarity. Specifically, the cooling is manifest under a range of favorable turbulence conditions which can develop and disappear locally on time scales of minutes to hours because of the intrinsically intermittent nature of stable boundary layer turbulence. Because of the fine-scale nature of the processes which produce the cooling, modeling it is a difficult proposition. Vertical resolution on the order of 1 meter is required. Adequate models of intermittent surface fluxes coupled with radiation exchange do not currently exist, not as parameterizations for aggregated systems, nor in large eddy simulation (LES) models. This presentation will introduce the theory. An important testable null hypothesis emerges: the icehouse effect produces a unique signature or "fingerprint" which could not be produced by any other known process. The presence of this signature will be demonstrated using nearly all available Arctic temperature observations. Its aggregate effect is clearly found in Arctic monthly surface temperature trends when sorted by climatological stability. Using all available Arctic rawinsonde ascents - about 1.1 million profiles, "frozen moments" of the icehouse processes are captured in various states. Because turbulent time scales are so short in the stable boundary layer. each of these snapshots can be treated as independent -- their chronology is irrelevant. Micrometeorological similarity is invoked to reassemble the soundings into bins of similar stability and it is in a wide, coherent range of these stability bins where the cooling effect is revealed.

Wetzel, Peter J.↗

On the Two-Dimensionalization of Homogeneous Rotating Turbulence

Large-eddy simulation of the incompressible Navier-Stokes equations has been used to examine the long-time development of initially isotropic turbulence subjected to solid-body rotation. The simulations were carried out using a pseudo-spectral method with 128 x 128 x 512 collocation points in a computational domain that is four times larger along the rotation axis than in the other directions; subgrid-scale motions were parameterized using a spectral eddy viscosity model modified for system rotation. Simulation results show that the correlation length along the rotation am's of velocities orthogonal to the rotation vector exhibits rapid growth while the integral length-scale of velocities aligned with the rotation axis is relatively unaffected by rotation. Examination of the energy spectrum of two-dimensional, two-component motions indicates the presence of an inverse cascade of energy. System rotation also causes an alignment of vorticity along the rotation axis with relatively stronger cyclonic vorticity than anticyclonic. The onset of anisotropic effects are well characterized by Rossby numbers defined in terms of both macroscopic and microscopic quantities.

Squires, K. D.↗

Cloud-Scale Numerical Modeling of the Arctic Boundary Layer

The research objective of this NASA grant-funded project was to determine in detail how large-scale processes. in combination with cloud-scale radiative, microphysical, and dynamical processes, govern the formation and multi-layered structure of Arctic stratus clouds. This information will be useful for developing and improving 1D (one dimensional) boundary layer models for the Arctic. Also, to quantitatively determine the effects of leads on the large-scale budgets of sensible heat, water vapor, and condensate in a variety of Arctic winter conditions. This information will be used to identify the most important lead-flux processes that require parameterization in climate models. Our approach was to use a high-resolution numerical model, the 2D (two dimensional) University of Utah Cloud Resolving Model (UU CRM), and its 1D version, the University of Utah Turbulence Closure Model (UU TCM), a boundary layer model based on third-moment turbulence closure, as well as a large-eddy simulation (LES) model originally developed by C.H. Moeng.

Kruegen, Steven K.↗

Anisotropic Structure of Rotating Homogeneous Turbulence at High Reynolds Numbers

Large eddy simulation is used to investigate the development of anisotropies and the evolution towards a quasi two-dimensional state in rotating homogeneous turbulence at high Reynolds number. The present study demonstrates the existence of two transitions in the development of anisotropy. The first transition marks the onset of anisotropy and occurs when a macro-Rossby number (based on a longitudinal integral lengthscale) has decreased to near unity while the second transition occurs when a micro-Rossby number (defined in this work as the ratio of the rms fluctuating vorticity to background vorticity) has decreased to unity. The anisotropy marked by the first transition corresponds to a reduction in dimensionality while the second transition corresponds to a polarization of the flow, i.e., relative dominance of the velocity components in the plane normal to the rotation axis. Polarization is reflected by emergence of anisotropy measures based on the two-dimensional component of the turbulence. Investigation of the vorticity structure shows that the second transition is also characterized by an increasing tendency for alignment between the fluctuating vorticity vector and the background angular velocity vector with a preference for corrotative vorticity.

Cambon, Claude↗

A Unified LES/RANS Approach Using the CE/SE Method

The goal of the proposed research is to develop advanced large eddy simulations (LES) and Reynolds average Navier-Stokes (RANS) capabilities for complex problems using the conservation element and solution element (CE/SE) method. This proposed use of the CE/SE method is justified by the many demonstrated advantages of the CE/SE method over the traditional computational fluid dynamics (CFD) methods. The proposed research will produce improved engineering modeling for near-term applications for focus programs as well as long-term impacts on the fundamental understanding of turbulence.

Liou, William W.↗

Enhancement of Cloud Cover and Suppression of Nocturnal Drizzle in Stratocumulus Polluted by Haze

Recent satellite observations indicate a significant decrease of cloud water in ship tracks, in contrast to an ensemble of in situ ship-track measurements that show no average change in cloud water relative to the surrounding clouds. We find through large-eddy simulations of stratocumulus that the trend in the satellite data is likely an artifact of sampling only overcast clouds. The simulations instead show cloud cover increasing with droplet concentrations. Our simulations also show that increases in cloud water from drizzle suppression (by increasing droplet concentrations) are favored at night or at extremely low droplet concentrations.

Ackerman, Andrew S.↗

Influence of Idealized Heterogeneity on Wet and Dry Planetary Boundary Layers Coupled to the Land Surface

We examine the influence of surface heterogeneity on boundary layers using a large-eddy simulation coupled to a land-surface model. Heterogeneity, imposed in strips varying from 2-30 km (1 less than lambda/z(sub i) less than 18), is found to dramatically alter the structure of the free convective boundary layer by inducing significant organized circulations. A conditional sampling technique, based on the scale of the surface heterogeneity (phase averaging), is used to identify and quantify the organized surface fluxes and motions in the atmospheric boundary layer. The impact of the organized motions on turbulent transport depends critically on the scale of the heterogeneity lambda, the boundary layer height zi and the initial moisture state of the boundary layer. Dynamical and scalar fields respond differently as the scale of the heterogeneity varies. Surface heterogeneity of scale 4 less than lamba/z(sub i) less than 9 induces the strongest organized flow fields (up, wp) while heterogeneity with smaller or larger lambda/z(sub i) induces little organized motion. However, the organized components of the scalar fields (virtual potential temperature and mixing ratio) grow continuously in magnitude and horizontal scale, as lambda/z(sub i) increases. For some cases, the organized motions can contribute nearly 100% of the total vertical moisture flux. Patch-induced fluxes are shown to dramatically impact point measurements that assume the time-average vertical velocity to be zero. The magnitude and sign of this impact depends on the location of the measurement within the region of heterogeneity.

Houser, Paul↗

Transition-Sensitized Turbulence Models for Compressible and Incompressible Flows

During the grant period from January 1,2002 to December 31,2002 work was carried out on three projects to extend the range of applicability of advanced turbulence models. First, a new transition-sensitized turbulence model was tested and refined. Second, the influence of compressibility on the pressure-strain rate correlation was studied. Third, the relationship between time-filtered large eddy simulation (TLES) and Reynolds-averaged Navier Stokes (RANS) modeling was investigated leading to submission of the article. The transition-sensitized turbulence model encompasses the early-stage transition and turbulent flow regimes describing the evolution of the ensemble mean disturbance energy and dissipation rate. It is founded on a consistent mathematical description of the laminar regime with its linear disturbances and the fully turbulent regime with its stochastic fluctuations. The unified description is provided by the ensemble viewpoint.

Thacker, William D.↗

An Aircraft Encounter with Turbulence in the Vicinity of a Thunderstorm

Large eddy simulations of three convective turbulence events are investigated and compared with observational data. Two events were characterized with severe turbulence and the other with moderate turbulence. Two of the events occurred during NASA s turbulence flight experiments during the spring of 2002, and the third was an event identified by the Flight Operational Quality Assurance (FOQA) Program. Each event was associated with developing or ongoing convection and was characterized by regions of low to moderate radar reflectivity. Model comparisons with observations are favorable. The data sets from these simulations can be used to test turbulence detection sensors.

Hamilton, David W.↗

Pressure Oscillations and Structural Vibrations in Space Shuttle RSRM and ETM-3 Motors

The complex interactions between internal motor pressure oscillations resulting from vortex shedding, the motor's internal acoustic modes, and the motor's structural vibration modes were assessed for the Space Shuttle four-segment booster Reusable Solid Rocket Motor and for the five-segment engineering test motor ETM-3. Two approaches were applied 1) a predictive procedure based on numerically solving modal representations of a solid rocket motor s acoustic equations of motion and 2) a computational fluid dynamics two-dimensional axi-symmetric large eddy simulation at discrete motor burn times.

Mason, D. R.↗

Numerical Study of Wake Vortex Behavior in Turbulent Domains with Ambient Stratification

A three-dimensional large eddy simulation model is used to investigate the sensitivity of ambient stratification with turbulence on the behavior of aircraft wake vortices. Modeled ambient turbulence levels range from very weak to moderate, and stratification levels range from strongly stable to unstable. The results of profound significance from this study are: 1) very little sensitivity between vortex linking time and the level of stratification, 2) the mean vortex separation remained nearly constant regardless of stratification and turbulence (at least prior to linking), 3) the wake vortices did not rise regardless of the level of stratification, and 4) for very strong stratification, the vortex stopped descending and quickly dissipated even before vortex linking could occur. These results are supported by experimental data and are contrary to conclusions from other numerical studies that assume laminar flow and/or relatively-low Reynolds numbers.

Switzer, George F.↗

Microphysics of Pyrocumulonimbus Clouds

The intense heat from forest fires can generate explosive deep convective cloud systems that inject pollutants to high altitudes. Both satellite and high-altitude aircraft measurements have documented cases in which these pyrocumulonimbus clouds inject large amounts of smoke well into the stratosphere (Fromm and Servranckx 2003; Jost et al. 2004). This smoke can remain in the stratosphere, be transported large distances, and affect lower stratospheric chemistry. In addition recent in situ measurements in pyrocumulus updrafts have shown that the high concentrations of smoke particles have significant impacts on cloud microphysical properties. Very high droplet number densities result in delayed precipitation and may enhance lightning (Andrew et al. 2004). Presumably, the smoke particles will also lead to changes in the properties of anvil cirrus produces by the deep convection, with resulting influences on cloud radiative forcing. In situ sampling near the tops of mature pyrocumulonimbus is difficult due to the high altitude and violence of the storms. In this study, we use large eddy simulations (LES) with size-resolved microphysics to elucidate physical processes in pyrocumulonimbus clouds.

Jensen, Eric↗