SEARCH · Engineering Papers

Results for “atmospheric 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 19 records

All sources NASA NTRS DOE OSTI

NASA NTRS1997

Toward Understanding Wake Vortices and Atmospheric Turbulence Interactions Using Large-Eddy Simulation

The vortices produced by an aircraft in flight are a complex phenomena created from a 'sheet of vorticity' leaving the trailing edge of the aircraft surfaces. This sheet tends to roll-up into two counter-rotating vortices. After a few spans downstream of the aircraft, the roll-up process is complete and the vortex pair may be characterized in a simple manner for modeling purposes. Our research will focus on what happens to these post roll-up vortices in the vicinity of an airport terminal. As the aircraft wake vortices descend, they are transported by the air mass which they are embedded and are decayed by both internal and external processes. In the vicinity of the airport, these external influences are usually due to planetary boundary layer (PBL) turbulence. Using large-eddy simulation (LES), one may simulate a variety of PBL conditions. In the LES method, turbulence is generated in the PBL as a response to surface heat flux, horizontal pressure gradient, wind shear, and/or stratification, and may produce convective or unstably stratified, neutral, or stably stratified PBL's. Each of these PBL types can occur during a typical diurnal cycle of the PBL. Thus it is important to be able to characterize these conditions with the LES method. Once this turbulent environment has been generated, a vortex pair will be introduced and the interactions are observed. The objective is to be able to quantify the PBL turbulence vortex interaction and be able to draw some conclusions of vortex behavior from the various scale interactions. This research is ongoing, and we will focus on what has been accomplished to date and the future direction of this research. We will discuss the model being used, show results that validate its use in the PBL, and present a nested-grid method proposed to analyze the entire PBL and vortex pair simultaneously.

DeCroix, David↗

NASA NTRS2003

Dynamic Turbulence Modelling in Large-eddy Simulations of the Cloud-topped Atmospheric Boundary Layer

The use of large eddy simulation, or LES, to study the atmospheric boundary layer dates back to the early 1970s when Deardor (1972) used a three-dimensional simulation to determine velocity and temperature scales in the convective boundary layer. In 1974 he applied LES to the problem of mixing layer entrainment (Deardor 1974) and in 1980 to the cloud-topped boundary layer (Deardor 1980b). Since that time the LES approach has been applied to atmospheric boundary layer problems by numerous authors. While LES has been shown to be relatively robust for simple cases such as a clear, convective boundary layer (Mason 1989), simulation of the cloud-topped boundary layer has proved more of a challenge. The combination of small length scales and anisotropic turbulence coupled with cloud microphysics and radiation effects places a heavy burden on the turbulence model, especially in the cloud-top region. Consequently, over the past few decades considerable effort has been devoted to developing turbulence models that are better able to parameterize these processes. Much of this work has involved taking parameterizations developed for neutral boundary layers and deriving corrections to account for buoyancy effects associated with the background stratification and local buoyancy sources due to radiative and latent heat transfer within the cloud (see Lilly 1962; Deardor 1980a; Mason 1989; MacVean & Mason 1990, for example). In this paper we hope to contribute to this effort by presenting a number of turbulence models in which the model coefficients are calculated dynamically during the simulation rather than being prescribed a priori.

Kirkpatrick, M. P.↗

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 NTRS1996

The Sensitivity of Large-Eddy Simulation to Local and Nonlocal Drag Coefficients at the Lower Boundary

It was found that the homogeneity of the surface drag coefficient plays an important role in the large scale structure of turbulence in large-eddy simulation of the convective atmospheric boundary layer. Particularly when a ground surface temperature was specified, large horizontal anisotropies occurred when the drag coefficient depended upon local velocities and heat fluxes. This was due to the formation of streamwise roll structures in the boundary layer. In reality, these structures have been found to form when shear is approximately balanced by buoyancy. The present cases, however, were highly convective. The formation was caused by particularly low values of the drag coefficient at the entrance to thermal plume structures.

Schowalter, D. G.↗

NASA NTRS1995

Lead-induced atmospheric circulations

Leads provide a significant source of heat and moisture to the Arctic winter atmosphere, and plumes from wide leads have been observed to penetrate the Arctic inversion. We have developed a two-dimensional, high-resolution, deep anelastic numerical model to investigate the atmospheric convection from leads with widths ranging from 100 m to 10 km. A second-order turbulence closure scheme is used to parameterize the atmospheric turbulence in the horizontally inhomogeneous system. This study describes how the lead-induced circulations can enhance the vertical transport of heat into the atmospheric boundary layer. This model is compared with large-eddy simulation results and with lidar observations of a lead-induced ice crystal plume. The model is used to study the effect of varying lead widths and ambient atmospheric conditions on the resultant convection from leads, and some preliminary results are described.

Alam, Afshan↗

NASA NTRS1999

Large Eddy Simulation of Cirrus Clouds

The Regional Atmospheric Modeling System (RAMS) with mesoscale interactive nested-grids and a Large-Eddy Simulation (LES) version of RAMS, coupled to two-moment microphysics and a new two-stream radiative code were used to investigate the dynamic, microphysical, and radiative aspects of the November 26, 1991 cirrus event. Wu (1998) describes the results of that research in full detail and is enclosed as Appendix 1. The mesoscale nested grid simulation successfully reproduced the large scale circulation as compared to the Mesoscale Analysis and Prediction System's (MAPS) analyses and other observations. Three cloud bands which match nicely to the three cloud lines identified in an observational study (Mace et al., 1995) are predicted on Grid #2 of the nested grids, even though the mesoscale simulation predicts a larger west-east cloud width than what was observed. Large-eddy simulations (LES) were performed to study the dynamical, microphysical, and radiative processes in the 26 November 1991 FIRE 11 cirrus event. The LES model is based on the RAMS version 3b developed at Colorado State University. It includes a new radiation scheme developed by Harrington (1997) and a new subgrid scale model developed by Kosovic (1996). The LES model simulated a single cloud layer for Case 1 and a two-layer cloud structure for Case 2. The simulations demonstrated that latent heat release can play a significant role in the formation and development of cirrus clouds. For the thin cirrus in Case 1, the latent heat release was insufficient for the cirrus clouds to become positively buoyant. However, in some special cases such as Case 2, positively buoyant cells can be embedded within the cirrus layers. These cells were so active that the rising updraft induced its own pressure perturbations that affected the cloud evolution. Vertical profiles of the total radiative and latent heating rates indicated that for well developed, deep, and active cirrus clouds, radiative cooling and latent heating could be comparable in magnitude in the cloudy layer. This implies that latent heating cannot be neglected in the construction of a cirrus cloud model. The probability density function (PDF) of w was analyzed to assist in the parameterization of cloud-scale velocities in large-scale models. For the more radiatively-driven, thin cirrus case, the PDFs are approximately Gaussian. However, in the interior of the deep, convectively unstable case, the PDFs of w are multi-modal and very broad, indicating that parameterizing cloud-scale motions for such clouds can be very challenging. The results of this research are described in detail in a paper submitted to the Journal of Atmospheric Science (Wu and Cotton, 1999), which is enclosed as Appendix 2. Using soundings extracted from a mesoscale simulation of the November 26, 1991 cirrus event, the radiative effects on vapor deposition/sublimation of ice crystals was studied using a two-dimensional cloud-resolving model (CRM) version of RAMS, coupled to an explicit bin-resolving microphysics. The CRM simulations of the November 26, 1991 cirrus event demonstrate that the radiative impact on the diffusional growth (or sublimation) of ice crystals is significant. In this case, the ice particles experienced radiative warming. Model results show that radiative feedbacks in the diffusional growth of ice particles can be very complex. Radiative warming of an ice particle will restrict the particle's diffusional growth. In the case of radiative warming, ice particles larger than a certain size will experience so much radiative warming that surface ice saturation vapor pressures become large enough to cause sublimation of the larger crystals, while smaller crystals are growing by vapor deposition. However, ice mass production can be enhanced in the case of radiative cooling of an ice particle. For the November 26, 1991 cirrus event, radiative feedback results in significant reduction in the total ice mass, especially in the production of large ice crystals, and consequently, both radiative and dynamic properties of the cirrus cloud are significantly affected. A complete description of this research has been submitted as a paper to the Journal of Atmospheric Science (Wu et al., 1999), and included as Appendix 3.

Wu, Ting↗

NASA NTRS1999

Large Eddy Simulation of Aircraft Wake Vortices in a Homogeneous Atmospheric Turbulence: Vortex Decay and Descent

The effects of ambient turbulence on decay and descent of aircraft wake vortices are studied using a validated, three-dimensional: large-eddy simulation model. Numerical simulations are performed in order to isolate the effect of ambient turbulence on the wake vortex decay rate within a neutrally-stratified atmosphere. Simulations are conducted for a range of turbulence intensities, by injecting wake vortex pairs into an approximately homogeneous and isotropic turbulence field. The decay rate of the vortex circulation increases clearly with increasing ambient turbulence level, which is consistent with field observations. Based on the results from the numerical simulations, simple decay models are proposed as functions of dimensionless ambient turbulence intensity (eta) and dimensionless time (T) for the circulation averaged over a range of radial distances. With good agreement with the numerical results, a Gaussian type of vortex decay model is proposed for weak turbulence: while an exponential type of Tortex decay model can be applied for strong turbulence. A relationship for the vortex descent based on above vortex decay model is also proposed. Although the proposed models are based on simulations assuming neutral stratification, the model predictions are compared to Lidar vortex measurements observed during stable, neutral, and unstable atmospheric conditions. In the neutral and unstable atmosphere, the model predictions appear to be in reasonable agreement with the observational data, while in the stably-stratified atmosphere, they largely underestimate the observed circulation decay with consistent overestimation of the observed vortex descent. The underestimation of vortex decay during stably-stratified conditions suggests that stratification has an important influence on vortex decay when ambient levels of turbulence are weak.

Han, Jongil↗

NASA NTRS2002

The Role of Turbulence in Chemical and Dynamical Processes in the Near-Field Wake of Subsonic Aircraft

During this grant, covering the period from September 1998 to December 2001, we continued the investigation of the role of turbulent mixing in the wake of subsonic aircraft initiated in 1994 for NASA's Atmospheric Effects of Aviation Project. The goal of the research has been to provide sufficient understanding and quantitative analytical capability to assess the dynamical, chemical, and microphysical interactions in the near-field wake that have the greatest potential to influence the global atmospheric impact of the projected fleet of subsonic aircraft. Through large-eddy simulations we have shown that turbulence in the early wake dynamics can have a strong effect on both the ice microphysics of contrail evolution and on wake chemistry. The wake vortex dynamics are the primary determinant of the vertical extent of the contrail; this together with the local wind shear largely determines the horizontal extent. The fraction of the initial ice crystals surviving the wake vortex dynamics, their spatial distribution, and the ice mass distribution are all sensitive to the aircraft type, assumed initial ice crystal number, and ambient humidity and turbulence conditions. Our model indicates that there is a significant range of conditions for which a smaller aircraft such as a B737 produces as significant a persistent contrail as a larger aircraft such as a B747, even though the latter consumes almost five times as much fuel. Large-eddy simulations of the near wake of a B757 provided a fine-grained chemical-dynamical representation of simplified NOx - HOx chemistry in wakes of ages from a few seconds to several minutes. By sampling the simulated data in a manner similar to that of in situ aircraft measurements it was possible to provide a likely explanation for a puzzle uncovered in the 1996 SUCCESS flight measurements of OH and HO2 The results illustrate the importance of considering fluid dynamics effects in interpreting chemistry results when mixing rates and species fluctuations are large, and demonstrate the feasibility of using 3D unsteady LES with coupled chemistry to study such phenomena.

Lewellen, D. C.↗

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 NTRS1999

Numerical Modeling Studies of Wake Vortices: Real Case Simulations

A three-dimensional large-eddy simulation model, TASS, is used to simulate the behavior of aircraft wake vortices in a real atmosphere. The purpose for this study is to validate the use of TASS for simulating the decay and transport of wake vortices. Three simulations are performed and the results are compared with the observed data from the 1994-1995 Memphis field experiments. The selected cases have an atmospheric environment of weak turbulence and stable stratification. The model simulations are initialized with appropriate meteorological conditions and a post roll-up vortex system. The behavior of wake vortices as they descend within the atmospheric boundary layer and interact with the ground is discussed.

Shen, Shao-Hua↗

NASA NTRS2002

Sensitivity of Stratocumulus Optical Depths to Droplet Concentrations: Satellite Observations and Large-Eddy Simulations

A number of observations and simulations have shown that increased droplet concentrations in ship tracks increase their total cross-sectional area, thereby enhancing cloud albedo and providing a negative (cooling) radiative forcing at the surface and the top of the atmosphere. In some cases cloud water has been found to be enhanced in ship tracks, which has been attributed to suppression of drizzle and implies an enhanced susceptibility of cloud albedo to droplet concentrations. However, observations from aircraft and satellite indicate that on average cloud water is instead reduced in daytime ship tracks. Such a reduction in liquid water may be attributable to cloud-burning caused by solar heating by soot within the ship exhaust, or by increased precipitation resulting from giant nuclei in the ship exhaust. We will summarize the observational evidence and present results from large-eddy simulations that evaluate these mechanisms. Along the way we will present our insights into the interpretation of satellite retrievals of cloud microphysical properties.

Ackerman, A. S.↗

NASA NTRS2001

Reduction of Cloud Water in Ship Tracks: Observations and Large-Eddy Simulations

Ship tracks represent a natural laboratory to study the effects of aerosols on clouds. A number of observations and simulations have shown that increased droplet concentrations in ship tracks increase their total cross-sectional area, thereby enhancing cloud albedo and providing a negative radiative forcing at the surface and the top of the atmosphere. In some cases, cloud water has been found to be enhanced in ship tracks, which has been attributed to suppression of drizzle and implies an enhanced susceptibility of cloud albedo to droplet concentrations. However, more recently compiled observations indicate that cloud water is instead reduced in daytime ship tracks on average. Such a response is consistent with cloud-burning due to solar absorption by soot (the semi-direct radiative forcing of aerosols), recently suggested to be suppressing trade cumulus cloud coverage over the Indian Ocean. We will summarize observational evidence and present large-eddy simulations that consider these competing mechanisms in the effects of aerosols on cloud albedo.

Ackerman, A. S.↗

NASA NTRS1993

Local isotropy in distorted turbulent boundary layers at high Reynolds number

This is a report on the continuation of our experimental investigations of the hypothesis of local isotropy in shear flows. This hypothesis, which states that at sufficiently high Reynolds numbers the small-scale structures of turbulent motions are independent of large-scale structures and mean deformations, has been used in theoretical studies of turbulence and computational methods such as large-eddy simulation. Since Kolmogorov proposed his theory, there have been many experiments, conducted in wakes, jets, mixing layers, a tidal channel, and atmospheric and laboratory boundary layers, in which attempts have been made to verify - or refute - the local-isotropy hypothesis. However, a review of the literature over the last five decades indicated that, despite all these experiments in shear flows, there was no consensus in the scientific community regarding this hypothesis, and, therefore, it seemed worthwhile to undertake a fresh experimental investigation into this question.

Saddoughi, Seyed G.↗

NASA NTRS1998

The NASA-Langley Wake Vortex Modelling Effort in Support of an Operational Aircraft Spacing System

Two numerical modelling efforts, one using a large eddy simulation model and the other a numerical weather prediction model, are underway in support of NASA's Terminal Area Productivity program. The large-eddy simulation model (LES) has a meteorological framework and permits the interaction of wake vortices with environments characterized by crosswind shear, stratification, humidity, and atmospheric turbulence. Results from the numerical simulations are being used to assist in the development of algorithms for an operational wake-vortex aircraft spacing system. A mesoscale weather forecast model is being adapted for providing operational forecast of winds, temperature, and turbulence parameters to be used in the terminal area. This paper describes the goals and modelling approach, as well as achievements obtained to date. Simulation results will be presented from the LES model for both two and three dimensions. The 2-D model is found to be generally valid for studying wake vortex transport, while the 3-D approach is necessary for realistic treatment of decay via interaction of wake vortices and atmospheric boundary layer turbulence. Meteorology is shown to have an important affect on vortex transport and decay. Presented are results showing that wake vortex transport is unaffected by uniform fog or rain, but wake vortex transport can be strongly affected by nonlinear vertical change in the ambient crosswind. Both simulation and observations show that atmospheric vortices decay from the outside with minimal expansion of the core. Vortex decay and the onset three-dimensional instabilities are found to be enhanced by the presence of ambient turbulence.

Proctor, Fred H.↗

NASA NTRS1992

Direct simulation of the stably stratified turbulent Ekman layer

The Navier-Stokes equations and the Boussinesq approximation were used to compute a 3D time-dependent turbulent flow in the stably stratified Ekman layer over a smooth surface. The simulation data are found to be in very good agreement with atmospheric measurements when nondimensionalized according to Nieuwstadt's local scaling scheme. Results suggest that, when Reynolds number effects are taken into account, the 'constant Froud number' stable layer model (Brost and Wyngaard, 1978) and the 'shearing length' stable layer model (Hunt, 1985) for the dissipitation rate of turbulent kinetic energy are both valid. It is concluded that there is good agreement between the direct numerical simulation results and large-eddy simulation results obtained by Mason and Derbyshire (1990).

Coleman, G. N.↗

NASA NTRS2003

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↗

NASA NTRS2000

Wake Vortex Transport and Decay in Ground Effect: Vortex Linking with the Ground

Numerical simulations are carried out with a three-dimensional Large-Eddy Simulation (LES) model to explore the sensitivity of vortex decay and transport in ground effect (IGE). The vortex decay rates are found to be strongly enhanced following maximum descent into ground effect. The nondimensional decay rate is found to be insensitive to the initial values of circulation, height, and vortex separation. The information gained from these simulations is used to construct a simple decay relationship. This relationship compares well with observed data from an IGE case study. Similarly, a relationship for lateral drift due to ground effect is constructed from the LES data. In the second part of this paper, vortex linking with the ground is investigated. Our numerical simulations of wake vortices for IGE show that a vortex may link with its image beneath the ground, if the intensity of the ambient turbulence is moderate to high. This linking with the ground (which is observed in real cases)gives the appearance of a vortex tube that bends to become vertically oriented and which terminates at the ground. From the simulations conducted, the linking time for vortices in the free atmosphere; i.e., a function of ambient turbulence intensity.

Proctor, Fred H.↗

NASA NTRS1991

Numerical simulations of turbulent thermal convection with differential rotation

The solar nebula, from which the planets in our solar system formed, featured a disk of gas and dust grains in rapid, differential rotation, and at some stage was likely to have been unstable to thermal convection. This situation is suspected by many to lead to significant turbulent Reynolds stress production and angular momentum transport in such systems, and estimates of transport rates have been attempted from unsubstantiated phenomenological models. In order to determine the circumstances and physical conditions under which our own planetary system formed and to explain recent observations of young stellar systems, it is necessary to develop realistic models of heat and angular momentum transport for such flows. Developing an understanding of complicated flows featuring thermal convection, rotation, and shear is also of wide interest in stellar astrophysics and in planetary and terrestrial atmospheric studies. The ultimate objective is to develop workable models based on the numerical simulations for constructing global solar nebula models; viz., relatively simple prescriptions for heat and angular momentum fluxes from given system parameters (e.g., ratios of rotation, shear, and convective lapse rates) are characterized, quantified, and developed. Toward this end, our program has been to attempt to understand the behavior of the direct numerical simulations of Boussinesq convection, which, despite the complexity of the results, is still an overly simplified approximation to the real system and should be more amenable to analysis. These results are also intended to be tested against turbulence models, especially those designed for atmospheric boundary layers, and may provide a basis for subgrid-scale models. In order to make the numerical simulations more realistic with regard to the solar nebula problem, a fully compressible code that will allow incorporation of large density stratifications and realistic thermodynamic and radiative properties is developed. In order to explore the properties of these flows at the very high values of Re found in natural systems and the very low values of Pr found in most astrophysical contexts, we will need to employ large-eddy simulations for which we want to determine the most appropriate subgrid-scale model to incorporate.

Cabot, William H.↗