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 199 records · Page 11

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.↗

The Effects of Aircraft Wake Dynamics on Contrail Development

Results of large-eddy simulations of the development of young persistent ice contrails are presented, concentrating on the interactions between the aircraft wake dynamics and the ice cloud evolution over ages front a few seconds to approx. 30 min. The 3D unsteady evolution of the dispersing engine exhausts, trailing vortex pair interaction and breakup, and subsequent Brunt-Vaisala oscillations of the older wake plume are modeled in detail in high-resolution simulations, coupled with it bulk microphysics model for the contrail ice development. The simulations confirm that the early wake dynamics can have a strong influence on the properties of persistent contrails even at late times. The vortex dynamics are the primary determinant of the vertical extent of the contrail (until precipitate ton becomes significant): and this together with the local wind shear largely determines the horizontal extent. The ice density, ice crystal number density, and a conserved exhaust tracer all develop and disperse in different fashions from each other. The total ice crystal number can be significantly reduced due to adiabatic compression resulting from the downward motion of the vortex system, even for ambient conditions that are substantially supersaturated with respect to ice. The fraction of the initial ice crystals surviving, their spatial distribution and the ice mass distribution are all sensitive to the aircraft type, ambient humidity, assumed initial ice crystal number, and ambient turbulence conditions. There is a significant range of conditions for which a smaller transport such as a B737 produces as significant a persistent contrail as a larger transport such as a B747, even though the latter consumes almost five times as much fuel. The difficulties involved in trying to minimize persistent contrail production are discussed.

Lewellen, D. C.↗

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↗

Theoretical and Observational Determination of Global and Regional Radiation Budget, Forcing and Feedbacks at the Top-of-Atmosphere and Surface

Report consists of: 1. List of accomplishments 2. List of publications 3. Abstracts of published or submitted papers and 4. Subject invention disclosure. The accomplishments of the grant listed are: 1. Improved the third-order turbulence closure in cloud resolving models to remove the liquid water oscillation. 2. Used the University of California-Los Angeles (UCLA) large-eddy simulation (LES) model to provide data for radiation transfer testing. 3. Revised shortwave k-distribution models based on HITRAN 2000. 4. Developed a gamma-weighted two-stream radiative transfer model for radiation budget estimate applications. 5. Estimated the effect of spherical geometry to the earth radiation budget. 6. Estimated top-of-atmosphere irradiance over snow and sea ice surfaces. 7. Estimated the aerosol direct radiative effect at the top of the atmosphere. 8. Estimated the top-of-atmosphere reflectance of the clear-sky molecular atmosphere over ocean. 9. Developed and validated new set of Angular Distribution Models for the CERES TRMM satellite instrument (tropical) 10. Developed and validated new set of Angular Distribution Models for the CERES Terra satellite instrument (global) 11. Quantified the top-of-atmosphere direct radiative effect of aerosols over global oceans from merged CERES and MODIS observations 12 Clarified the definition of TOA flux reference level for radiation budget studies 13. Developed new algorithm for unfaltering CERES measured radiances 14. Used multiangle POLDER measurements to produce narrowband angular distribution models and examine the effect of scene identification errors on TOA albedo estimates 15. Developed and validated a novel algorithm called the Multidirectional Reflectance Matching (MRM) model for inferring TOA albedos from ice clouds using multi-angle satellite measurements. 16. Developed and validated a novel algorithm called the Multidirectional Polarized Reflectance Matching (MPRM) model for inferring particle shapes from ice clouds using multi-angle polarized satellite measurements. 17. Developed 4 advanced light scattering models including the three-dimensional (3D) uniaxial perfectly matched layer (UPML) finite-difference time-domain (FDTD) model. 18. Develop sunglint in situ measurement and study reflectance distribution in the sunglint area. 19. Lead a balloon-borne radiometer TOA albedo validation effort. 20. Developed a CERES surface UVB, UVA, and UV index product.

Loeb, Norman G.↗

Limits to the Indirect Aerosol Forcing in Stratocumulus

The indirect radiative forcing of aerosols is poorly constrained by the observational data underlying the simple cloud parameterizations in GCMs. signal of cloud response to increased aerosol concentrations from meteorological noise. Recent satellite observations indicate a significant decrease of cloud water in ship tracks, in contrast to an ensemble of in situ measurements showing no average change in cloud water relative to the surrounding clouds. Both results contradict the expectation of cloud water increasing in polluted 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. The simulations also show that increases in cloud water from suppressing drizzle by increased droplet concentrations are favored at night or at extremely low droplet concentrations. At typical droplet concentrations we find that the Twomey effect on cloud albedo is amplified very little by the secondary indirect effect of drizzle suppression, largely because the absorption of solar radiation by cloud water reduces boundary-layer mixing in the daytime and thereby restricts any possible increase in cloud water from drizzle suppression. The cloud and boundary layer respond to radiative heating variations on a time scale of hours, and on longer time scales respond to imbalances between large-scale horizontal advection and the entrainment of inversion air. We analyze the co-varying response of cloud water, cloud thickness, width of droplet size distributions, and dispersion of the optical depth, as well as the overall response of cloud albedo, to changes in droplet concentrations. We also dissect the underlying physical mechanisms through sensitivity studies. Ship tracks represent an ideal natural laboratory to extricate the

Ackerman, Andrew↗

Evidence for Tropopause Layer Moistening by Convection During CRYSTAL-FACE

Measurements and analysis of the impact of deep convection on tropopause layer moisture are easily confounded by difficulties making precise observations with sufficient spatial coverage before and after convective events and difficulties distinguishing between changes due to local convection versus large-scale advection. The interactions between cloud microphysics and dynamics in the convective transport of moisture into the tropopause layer also result in a sufficiently complex and poorly characterized system to allow for considerable freedom in theoretical models of stratosphere-troposphere exchange. In this work we perform detailed large-eddy simulations with an explicit cloud microphysics model to study the impact of deep convection on tropopause layer moisture profiles observed over southern Florida during CRYSTALFACE. For four days during the campaign (July 11, 16, 28, and 29) we initialize a 100-km square domain with temperature and moisture profiles measured prior to convection at the PARSL ground site, and initiate convection with a warm bubble that produces an anvil at peak elevations in agreement with lidar and radar observations on that day. Comparing the moisture field after the anvils decay with the initial state, we find that convection predominantly moistens the tropopause layer (as defined by minimum temperature and minimum potential temperature lapse rate), although some drying is also predicted in localized layers. We will also present results of sensitivity tests designed to separate the roles of cloud microphysics and dynamics.

Ackerman, A.↗

The NASA Turbulent Heat Flux Experiments: Summary and Lessons Learned

The Turbulent Heat Flux (THX) experiments were conducted at NASA Glenn Research Center (GRC) in order to collect measurements of velocities and temperatures for computational fluid dynamics (CFD) validation of heated flows, with a focus on propulsion system components. The experiments spanned 5 phases; four of which were conducted in the GRC AeroAcoustic Propulsion Laboratory (AAPL) using the Small Hot Jet Flow Rig (SHJAR). In addition to making velocity measurements with Particle Image Velocimetry (PIV), the THX experiments introduced a new Raman-scattering based capability to measure temperatures. Computational studies were also conducted for each of the experimental configurations, in order to provide a baseline of expected CFD results and conduct an assessment of the capability of various CFD approaches for calculating flows where the turbulent transport of heat was important. Two of the collected sets of data were used for American Institute of Aeronautics and Astronautics (AIAA) Propulsion Aerodynamic Workshops (PAWs). The data set from the 5th phase, collected for heated supersonic jets, was also used to construct new validation cases for the NASA Turbulence Model Resource (TMR). This paper provides an overview of the experiments and associated computations for each of the 5 test phases. Key experimental findings are presented. Lessons learned are provided concerning the effect of computational modeling choice on accuracy of predicting turbulent flows where thermal transport is important. Emphasis is placed on comparing Reynolds- averaged Navier-Stokes approaches with large-eddy simulation approaches. The benefits of utilizing a conjugate heat transfer method in conjunction with CFD solver for film cooling is demonstrated.

RANS↗

The NASA Turbulent Heat Flux (THX) Experiments: Summary and Lessons Learned

The Turbulent Heat Flux (THX) experiments were conducted at NASA Glenn Research Center (GRC) in order to collect measurements of velocities and temperatures for computational fluid dynamics (CFD) validation of heated flows, with a focus on propulsion system components. The experiments spanned 5 phases; four of which were conducted in the GRC AeroAcoustic Propulsion Laboratory (AAPL) using the Small Hot Jet Flow Rig (SHJAR). In addition to making velocity measurements with Particle Image Velocimetry (PIV), the THX experiments introduced a new Raman-scattering based capability to measure temperatures. Computational studies were also conducted for each of the experimental configurations, in order to provide a baseline of expected CFD results and conduct an assessment of the capability of various CFD approaches for calculating flows where the turbulent transport of heat was important. Two of the collected sets of data were used for American Institute of Aeronautics and Astronautics (AIAA) Propulsion Aerodynamic Workshops (PAWs). The data set from the 5th phase, collected for heated supersonic jets, was also used to construct new validation cases for the NASA Turbulence Model Resource (TMR). This paper provides an overview of the experiments and associated computations for each of the 5 test phases. Key experimental findings are presented. Lessons learned are provided concerning the effect of computational modeling choice on accuracy of predicting turbulent flows where thermal transport is important. Emphasis is placed on comparing Reynolds-averaged Navier-Stokes approaches with large-eddy simulation approaches. The benefits of utilizing a conjugate heat transfer method in conjunction with CFD solver for film cooling is demonstrated.

RANS↗

Scalable High Performance Computing: Direct and Large-Eddy Turbulent Flow Simulations Using Massively Parallel Computers

This final report contains reports of research related to the tasks "Scalable High Performance Computing: Direct and Lark-Eddy Turbulent FLow Simulations Using Massively Parallel Computers" and "Devleop High-Performance Time-Domain Computational Electromagnetics Capability for RCS Prediction, Wave Propagation in Dispersive Media, and Dual-Use Applications. The discussion of Scalable High Performance Computing reports on three objectives: validate, access scalability, and apply two parallel flow solvers for three-dimensional Navier-Stokes flows; develop and validate a high-order parallel solver for Direct Numerical Simulations (DNS) and Large Eddy Simulation (LES) problems; and Investigate and develop a high-order Reynolds averaged Navier-Stokes turbulence model. The discussion of High-Performance Time-Domain Computational Electromagnetics reports on five objectives: enhancement of an electromagnetics code (CHARGE) to be able to effectively model antenna problems; utilize lessons learned in high-order/spectral solution of swirling 3D jets to apply to solving electromagnetics project; transition a high-order fluids code, FDL3DI, to be able to solve Maxwell's Equations using compact-differencing; develop and demonstrate improved radiation absorbing boundary conditions for high-order CEM; and extend high-order CEM solver to address variable material properties. The report also contains a review of work done by the systems engineer.

Morgan, Philip E.↗

Limits to Cloud Susceptibility

1-kilometer AVHRR observations of ship tracks in low-level clouds off the west coast of the U S. were used to determine limits for the degree to which clouds might be altered by increases in anthropogenic aerosols. Hundreds of tracks were analyzed to determine whether the changes in droplet radii, visible optical depths, and cloud top altitudes that result from the influx of particles from underlying ships were consistent with expectations based on simple models for the indirect effect of aerosols. The models predict substantial increases in sunlight reflected by polluted clouds due to the increases in droplet numbers and cloud liquid water that result from the elevated particle concentrations. Contrary to the model predictions, the analysis of ship tracks revealed a 15-20% reduction in liquid water for the polluted clouds. Studies performed with a large-eddy cloud simulation model suggested that the shortfall in cloud liquid water found in the satellite observations might be attributed to the restriction that the 1-kilometer pixels be completely covered by either polluted or unpolluted cloud. The simulation model revealed that a substantial fraction of the indirect effect is caused by a horizontal redistribution of cloud water in the polluted clouds. Cloud-free gaps in polluted clouds fill in with cloud water while the cloud-free gaps in the surrounding unpolluted clouds remain cloud-free. By limiting the analysis to only overcast pixels, the current study failed to account for the gap-filling predicted by the simulation model. This finding and an analysis of the spatial variability of marine stratus suggest new ways to analyze ship tracks to determine the limit to which particle pollution will alter the amount of sunlight reflected by clouds.

Coakley, James A., Jr.↗

Effect of rotation on isotropic turbulence - Computation and modelling

This paper uses numerical simulation to analyse the effects of uniform rotation on homogeneous turbulence. Both large-eddy and full simulations were made. The results indicate that the predominant effect of rotation is to decrease the rate of dissipation of the turbulence and increase the lengthscales, especially those along the axis of rotation. These effects are a consequence of the reduction, due to the generation of inertial waves, of the net energy transfer from large eddies to small ones. Experiments are also influenced by a more complicated interaction between the rotation and the wakes of the turbulence-generating grid which modifies the nominal initial conditions in the experiment. The latter effect is accounted for in simulations by modifying the initial conditions. Finally, a two-equation model is proposed that accounts for the effects of rotation and is able to reproduce the experimental decay of the turbulent kinetic energy.

Bardina, J.↗

Cost-effective accurate coarse-grid method for highly convective multidimensional unsteady flows

A fundamentally multidimensional convection scheme is described based on vector transient interpolation modeling rewritten in conservative control-volume form. Vector third-order upwinding is used as the basis of the algorithm; this automatically introduces important cross-difference terms that are absent from schemes using component-wise one-dimensional formulas. Third-order phase accuracy is good; this is important for coarse-grid large-eddy or full simulation. Potential overshoots or undershoots are avoided by using a recently developed universal limiter. Higher order accuracy is obtained locally, where needed, by the cost-effective strategy of adaptive stencil expansion in a direction normal to each control-volume face; this is controlled by monitoring the absolute normal gradient and curvature across the face. Higher (than third) order cross-terms do not appear to be needed. Since the wider stencil is used only in isolated narrow regions (near discontinuities), extremely high (in this case, seventh) order accuracy can be achieved for little more than the cost of a globally third-order scheme.

Leonard, B. P.↗

Boundary layer transition in hypersonic flows

This paper summarizes some of the recent progress made at NASA Langley Research Center in the understanding, prediction and modeling of high speed boundary-layer transition. Linear and nonlinear theories together with large-eddy and direct numerical simulations have been used to understand various aspects of the transition problem while low disturbance 'quiet' tunnels provide means for validating the theoretical results.

Malik, M. R.↗

Aerosol-Cloud Interactions during Tropical Deep Convection: Evidence for the Importance of Free Tropospheric Aerosols

NASA's 2002 CRYSTAL-FACE field experiment focused on the formation and evolution of tropical cirrus cloud systems in southern Florida. Multiple aircraft extensively sampled cumulonimbus dynamical and microphysical properties, as well as characterizing ambient aerosol populations both inside and outside the full depth of the convective column. On July 18, unique measurements were taken when a powerful updraft was traversed directly by aircraft, providing a window into the primary source region of cumulonimbus anvil crystals. Observations of the updraft, entered at approximately l0 km altitude and -34 C, indicated more than 200 cloud particles per mL at vertical velocities exceeding 20 m/s and the presence of significant condensation nuclei and liquid water within the core. In this work, aerosol and cloud phase observations are integrated by simulating the updraft conditions using a large-eddy resolving model with 3 explicit multiphase microphysics, including treatment of size-resolved aerosol fields, aerosol activation and freezing, and evaporation of cloud particles back to the aerosol phase. Simulations were initialized with observed thermodynamic and aerosol size distributions profiles and convection was driven by surface fluxes assimilated from the ARPS forecast model. Model results are consistent with the conclusions that most crystals are homogeneously frozen droplets and that entrained free tropospheric aerosols may contribute a significant fraction of the crystals. Thus most anvil crystals appear to be formed aloft in updraft cores, well above cloud base. These conclusions are supported by observations of hydrometeor size distribution made while traversing the dore, as well as aerosol and cloud particle size distributions generally observed by aircraft below 4km and crystal properties generally observed by aircraft above 12km.

Ackerman, A.↗

Streamwise Vorticity Generation in Laminar and Turbulent Jets

Complex streamwise vorticity fields are observed in the evolution of non-circular jets. Generation mechanisms are investigated via Reynolds-averaged (RANS), large-eddy (LES) and direct numerical (DNS) simulations of laminar and turbulent rectangular jets. Complex vortex interactions are found in DNS of laminar jets, but axis-switching is observed only when a single instability mode is present in the incoming mixing layer. With several modes present, the structures are not coherent and no axis-switching occurs, RANS computations also produce no axis-switching. On the other hand, LES of high Reynolds number turbulent jets produce axis-switching even for cases with several instability modes in the mixing layer. Analysis of the source terms of the mean streamwise vorticity equation through post-processing of the instantaneous results shows that, complex interactions of gradients of the normal and shear Reynolds stresses are responsible for the generation of streamwise vorticity which leads to axis-switching. RANS computations confirm these results. k - epsilon turbulence model computations fail to reproduce the phenomenon, whereas algebraic Reynolds stress model (ASM) computations, in which the secondary normal and shear stresses are computed explicitly, succeeded in reproducing the phenomenon accurately.

Demuren, Aodeji O.↗

Lumley's PODT definition of large eddies and a trio of numerical procedures

Lumley's 1967 Moscow paper provided, for the first time, a completely rational definition of the physically-useful term 'large eddy', popular for a half-century. The numerical procedures based upon his results are: (1) PODT (Proper Orthogonal Decomposition Theorem), which extracts the Large Eddy structure of stochastic processes from physical or computer simulation two-point covariances, and 2) LEIM (Large-Eddy Interaction Model), a predictive scheme for the dynamical large eddies based upon higher order turbulence modeling. Earlier Lumley's work (1964) forms the basis for the final member of the triad of numerical procedures: this predicts the global neutral modes of turbulence which have surprising agreement with both structural eigenmodes and those obtained from the dynamical equations. The ultimate goal of improved engineering design tools for turbulence may be near at hand, partly due to the power and storage of 'supermicrocomputer' workstations finally becoming adequate for the demanding numerics of these procedures.

Payne, Fred R.↗