Engineering Papers⌕ Search

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.

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

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Evaluating the Collision‐Coalescence Process in Idealized Cloud Convection Using Large‐Eddy Simulations With Lagrangian Microphysics

Drizzle initiation through the collision and coalescence of cloud droplets plays a crucial role in warm cloud precipitation. Recent theoretical studies suggest that the influence of collisional growth on the droplet size distribution can be quantified by a non-dimensional drizzle number (Dz). Here, large-eddy simulations with Lagrangian microphysics are employed to evaluate the theory by simulating a tall convection-cloud chamber under various conditions. Results show that the smaller the Dz, the larger the impact of collisions on the right tail of the droplet size distribution, consistent with the theory. The simulations confirm that the collision rate can be estimated from the droplet size distribution interacting only with cloud droplets of the same size at the mode radius. This suggests that the idealized theory can be a useful tool to design a cloud chamber for drizzle investigation, as well as to represent drizzle formation in models of real atmospheric clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Impact of Forest Canopy Structure on Buoyant Plume Dynamics During Wildland Fires

Heterogeneous forest canopies can generate complex turbulent structures, but in the presence of a fire plume, these interactions are not fully understood. This study investigates the influence of forest canopy heterogeneity on buoyant plume dynamics resulting from surface thermal anomalies representing wildland fires, utilizing Large Eddy Simulation (LES). The Parallelized Large-Eddy Simulation Model (PALM) was employed to simulate six canopy configurations: no canopy, homogeneous canopy, external plume-edge canopy, internal plume-edge canopy, 100 m gap canopy, and 200 m gap canopy. Each configuration was analyzed with and without a static surface heat flux patch of 5000 W ∙ m -2 , resulting in a resting buoyant plume. Simulations were conducted under three crosswind speeds: 0, 5, and 10 m ∙ s -1 . Results show that canopy structure significantly modifies plume behavior, mean flow, and turbulent kinetic energy (TKE) budgets. Plume updraft speed and tilt varied with canopy configuration and crosswind speed. Horizontal pressure gradients associated with plume-atmosphere interaction were modified based on the canopy configuration, resulting in varying crosswind speed reductions at the plume region. Strong momentum absorption was observed above the canopy for the crosswind cases, with the greatest enhancement in the gap canopies. Momentum injection from below the canopy due to the heat source was also observed, resulting in plume structure modulation based on canopy configuration. TKE was found to be the largest in the gap canopy configurations. TKE budget analysis revealed that buoyant production dominated over shear production. At the center of the heat patch, the gap canopy configurations showed enhanced buoyancy within the gap. These results improve our knowledge of fire-canopy-atmosphere interactions that can inform fire models on the impacts of canopy heterogeneity on plume dynamics and ember ejections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Simplified Approximations of Direct Cumulus Entrainment and Detrainment

Abstract In recent years, direct calculations of simulated cumulus entrainment and detrainment have facilitated new physical insights into these highly elusive but critically important processes. However, these calculations require substantial computational resources that may limit their widespread usage. To facilitate such calculations, two simplified approximations of direct cumulus entrainment and detrainment are examined herein. The first approximation, termed the “semidirect” method, follows a standard bulk approach but makes more realistic assumptions about the sources of entrained and detrained air near the cloud edges. In contrast, the second approximation (the “projection” method) uses the governing equations of motion to project whether grid points near the cloud edge will entrain or detrain as the mean cloud ascends by one grid point. Verification exercises using large-eddy simulations reveal that both methods generally agree better with corresponding direct entrainment/detrainment estimates than the traditional bulk formulation, with the projection method outperforming the semidirect method. The two methods can be used in a synergistic fashion, with the semidirect method helping to optimize the projection method, to suit a wide range of applications. Because the latter incorporates the essential dynamics of entrainment and detrainment at the local scale, it can be used to gain physical insight into the causal mechanisms regulating these complex processes.

Meteorology & Atmospheric Sciences↗
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.↗
DOE OSTI2025

Exploring the impact of surface topography on Rayleigh-Bénard dry convection in the Pi cloud chamber using OpenFOAM: In cylindrical and rectangular geometries

The Pi convection-cloud chamber can generate steady-state turbulence in both rectangular and cylindrical shapes via Rayleigh-Bénard convection (RBC) by maintaining warm bottom and cold top surfaces. Although most experiments in the Pi chamber were conducted in cylindrical shapes, all previous Pi chamber simulations were conducted in a rectangular shape due to the limitations of those models to discretize a cylindrical domain when using the finite difference method therein. Here, we use OpenFOAM, an open-source finite-volume-based Computational Fluid Dynamics (CFD) software package, to conduct Large-Eddy Simulation (LES) of dry RBC in the Pi chamber at high Rayleigh numbers (10 8 to 10 9 ). Results show that large-scale circulation (LSC) direction varies in the chamber with a constant side wall temperature. Imposing a slight temperature imbalance at the side wall ranging from 0.1 to 0.7 degrees can lock the LSC, aligning better with Pi chamber observations, particularly at higher Rayleigh numbers. In addition, we examine the impact of surface topography on LSC and heat transfer in RBC systems within cylindrical and rectangular shapes under varying conditions. Results show that roughing top/bottom surfaces by adding bars of a few tens millimeters height can strengthen thermal plumes and enhance temperature fluctuations in the chamber. Furthermore, we observe that different bar height configurations lead to notable changes in LSC orientation and thermal stratification, highlighting the complex interactions between surface features and convection patterns. This finding highlights how surface topography and chamber geometry affect Rayleigh-Bénard convection, improving understanding of turbulent heat transfer and atmospheric boundary-layer processes. Direct Numerical Simulations (DNS) are also conducted to validate LES results. In conclusion, while LES effectively captures qualitative behaviors seen in DNS, it tends to underestimate velocity variances near walls, illustrating a trade-off between computational efficiency and accuracy.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2021

Nighttime Convection in Water-ice Clouds at High Northern Latitudes on Mars

We investigate water-ice clouds and their influence on the temperature structure of the Mar-tian atmosphere at high northern latitudes in early summer. New results are obtained throughcoordinated analysis of two types of data from Mars Global Surveyor: atmospheric profiles re-trieved from radio occultation (RO) measurements and wide-angle images from the Mars OrbiterCamera (MOC). Some RO profiles contain a layer of neutral static stability, which indicates thepresence of convective mixing at a local time (about 5 h) when it does not usually occur. Thesenocturnal mixed layers (NMLs) were observed frequently in early summer of Mars year 27 atlatitudes of 53–72◦N and longitudes of 210–330◦E. The base of a typical NML is 3 km abovethe surface, about the same height as the nighttime cloud layer detected by the Phoenix LIDARin early summer of Mars year 29 at 234◦E, 68◦N. The depth of the NMLs ranges from less than1 km to more than 5 km. Comparisons with nearly simultaneous MOC images demonstrate thatNMLs are closely associated with water-ice clouds. There is a dense cluster of NMLs within theannular cloud that appears every year in early summer between Alba Mons and the north polarresidual ice cap. The lighting conditions at this location and season allowed MOC to observe theannular cloud on most orbits, at 118-min intervals. Its appearance varies dramatically with localtime, becoming more symmetrical and better organized at night and dissipating to a crescentshape during the day. According to high-resolution numerical simulations (Spiga et al., 2017),including a large-eddy simulation at the Phoenix landing site, NMLs form when radiative cool-ing by water-ice aerosols causes convective instability; the mixed layer is forced from above bynegative buoyancy. Our results strongly support this conclusion. In addition, MOC images frommidsummer contain eastward-moving frontal clouds. Temperature profiles within these cloudsshow signs of near-surface advection of warm air, which reduces the static stability of the loweratmosphere and contributes, along with cloud radiation, to the formation of an NML.

David Hinson↗
DOE OSTI2026

A Model Intercomparison Study to Investigate Mixing Characteristics in Non‐Precipitating Stratocumulus Clouds

Recent aircraft observations of marine stratocumulus clouds consistently showed that cloud microphysical relationships vary with altitude, indicating inhomogeneous mixing characteristics near cloud top and homogeneous mixing characteristics in mid-levels of clouds. Here, we conduct model intercomparison of an idealized, non-precipitating stratocumulus cloud to evaluate model consistency and examine whether simulations can reproduce the observed mixing characteristics. The results show that eleven large-eddy simulations with various dynamics and microphysics schemes show good agreement on the thermodynamical, microphysical, and dynamical properties of the stratocumulus-topped boundary layer in a steady state. The inter-model spread in steady-state liquid water path is significantly reduced compared to previous model intercomparison studies. This improvement might be due to better models and more consistent initial conditions than those used decades ago. In addition, most simulations, including a low-dimensional simulation, capture inhomogeneous mixing characteristics near the cloud top and homogeneous mixing characteristics inside the cloud. Moreover, simulations using Lagrangian microphysics schemes agree better with the observed mixing characteristics compared with those using the bin microphysics schemes. Since most simulations do not fully resolve the entrainment process, the apparent mixing characteristics arise from the variations in the resolved cloud properties. Our results support the vertical circulation mixing hypothesis, which suggests that homogeneous mixing characteristics in mid-levels of clouds are due to the vertical circulation of entrainment-affected and diluted parcels from the cloud top moved to lower levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Toward Understanding the Differences between Mesoscale and Large-Eddy Simulations of Tropical Cyclones

In this work, we investigate the ability of mesoscale and large-eddy simulation (LES) model configurations to predict the mean wind speed profile within the boundary layer of tropical cyclones (TCs). To this end, we perform idealized simulations of five hypothetical intense storms ranging from categories 1 to 5 on the Saffir–Simpson scale and extract time-averaged quantities near the eyewall region. We compare the model-generated data against mean wind speed profiles compiled from dropsondes launched from reconnaissance aircraft operating in the North Atlantic basin. Our analysis shows that mesoscale- and LES-generated mean wind fields display important differences in the boundary layer, including the magnitude of shear as well as the height where their low-level wind speed maxima are located. In addition, a comparison between the two model configurations with the dropsonde data shows that both modeling approaches are unable to capture the typical structure of mean winds in the lower part of the TC boundary layer (10–500 m), calling into question the use of simulations of near-axisymmetric storms for investigating the wind structure of past events. To better understand these differences, we conduct a momentum-budget analysis and show that modeled turbulent fluxes are underestimated in the mesoscale boundary layer parameterization compared to the LES model. Based on the analysis of the horizontal turbulent fluxes and their potential impact on mean flow quantities, a TC-specific boundary layer parameterization may be needed.

17 WIND ENERGY↗
NASA NTRS2021

Large-Eddy Simulations of Marine Boundary-Layer Clouds Associated With Cold Air Outbreaks During the Activate Campaign – Part 1: Case Setup and Sensitivities to Large-Scale Forcings

Large-eddy simulation (LES) is able to capture key boundary-layer (BL) turbulence and cloud processes. Yet, large-scale forcing and surface turbulent fluxes of sensible and latent heat are often poorly prescribed for LES simulations. We derive these quantities from measurements and reanalysis obtained for two cold air outbreak (CAO) events during Phase I of the Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment (ACTIVATE) in February-March 2020. We study the two contrasting CAO cases by performing LES and test the sensitivity of BL structure and clouds to large-scale forcings and turbulent heat fluxes. Profiles of atmospheric state and large-scale divergence and surface turbulent heat fluxes obtained from the reanalysis data ERA5 agree reasonablywell with those derived fromACTIVATE field measurements for both cases at the sampling time and location. Therefore, we adopt the time evolving heat fluxes, wind and advective tendencies profiles from ERA5 reanalysis data to drive the LES.We find that large-scale thermodynamic advective tendencies and wind relaxations are important for the LES to capture the evolving observed BL meteorological states characterized by the hourly ERA5 reanalysis data and validated by the observations. We show that the divergence (or vertical velocity) is important in regulating the BL growth driven by surface heat fluxes in LES simulations. The evolution of liquid water path is largely affected by the evolution of surface heat fluxes. The liquid water path imulated in LES agrees reasonably well with the ACTIVATE measurements. This study paves the path to investigate aerosol-cloud-meteorology interactions using LES informed and evaluated by ACTIVATE field measurements.

Large-eddy simulation↗
DOE OSTI2025

Ensemble‐Based, Large‐Eddy Reconstruction of Wind Turbine Inflow in a Near‐Stationary Atmospheric Boundary Layer Through Generative Artificial Intelligence

ABSTRACT To validate the second‐by‐second dynamics of turbines in field experiments, it is necessary to accurately reconstruct the winds going into the turbine. Current time‐resolved inflow reconstruction techniques estimate wind behavior in unobserved regions using relatively simple spectral‐based models of the atmosphere. Here, we develop a technique for time‐resolved inflow reconstruction that is rooted in a large‐eddy simulation model of the atmosphere. Our “large‐eddy reconstruction” technique blends observations and atmospheric model information through a diffusion model machine learning algorithm, allowing us to generate probabilistic ensembles of reconstructions for a single 10‐min observational period. Our generated inflows can be used directly by aeroelastic codes or as inflow boundary conditions in a large‐eddy simulation. We verify the second‐by‐second reconstruction capability of our technique in three synthetic field campaigns, finding positive Pearson correlation coefficient values () between ground‐truth and reconstructed streamwise velocity, as well as smaller positive correlation coefficient values for unobserved fields (spanwise velocity, vertical velocity, and temperature). We validate our technique in three real‐world case studies by driving large‐eddy simulations with reconstructed inflows and comparing to independent inflow measurements. The reconstructions are visually similar to measurements, follow desired power spectra properties, and track second‐by‐second behavior ().

17 WIND ENERGY↗
DOE OSTI2024

Development of the ARM Lagrangian Large-Scale Forcing Data (ARMLAGTRAJ) Value-Added Product Based on the lagtraj Framework

The Atmospheric Radiation Measurement (ARM) large-scale forcing data developed based on the constrained variational analysis (VARANAL) value-added product (VAP) (Zhang and Lin 1997, Zhang et al. 2001, Xie et al. 2004, Tang et al. 2019) has been widely used for single-column models (SCMs), cloud-resolving models (CRMs), and large-eddy simulation models (LESs) to understand and improve physical processes in models. Recently, the U.S. Department of Energy (DOE) ARM user facility conducted several major field campaigns using ship-based moving observational platforms. For example, the Marine ARM GPCI Investigation of Clouds (MAGIC) field campaign focused on the role of subtropical marine-boundary layer (MBL) clouds, and the Multidisciplinary Drifting Observatory for the Study of Arctic Climate (MOSAiC) field campaign aimed to improve understanding of the coupled climate systems in the Arctic. Observations from moving platforms are critical to provide a comprehensive characterization of coupled-system processes associated with all stages of the cloud and/or sea-ice life cycle. Traditional ARM large-scale forcing data have been developed at fixed locations. They need to be extended to include these moving platforms to address data needs for ship-based field campaigns or to support LES modeling in a Lagrangian framework. With these considerations in mind, we develop ARM-type Lagrangian large-scale forcing data sets based on the lagtraj framework (Boeing et al. 2020) with notable enhancements in generating forcings that are more suitable for ARM field campaigns. The lagtraj is a novel tool that generates forcings for LES and SCM simulation in both Lagrangian and Eulerian perspective. This technical report focuses on the major changes we performed on the lagtraj algorithm and provides an overview of the ARM Lagrangian Large-Scale Forcing Data (ARMLAGTRAJ) value-added products.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Nonsteady Load Responses of Wind Turbines to Atmospheric and Mountain-Generated Turbulence Eddies, With Impacts on the Main Bearing: A Validation Study

Previous computational and field experiments identify three characteristic time scales in the aerodynamic responses of utility-scale wind turbine loads to atmospheric boundary layer (ABL) turbulence: a 30-90 second time scale for the passage of high/low speed "streaks" through the rotor plane, the blade and rotor rotation time scales (approximately 1 to 5 seconds), and a sub-second time scale created by blade rotation through gradients within eddy coherent structure. In the current study we compare aerodynamic load responses from daytime ABL turbulence quantified with large-eddy simulation and a actuator line model of the NREL 5 MW wind turbine with analysis of field data from the NREL/GE 1.5 MW wind turbine 5 kilometers east of the Rocky Mountain Front Range in Colorado. In addition, we contrast the responses to the passage of the mountain-generated eddies embedded within the westerly winds with the ABL eddies embedded within northerly/southerly winds. These analyses are in context with the nonsteady forcing of the main bearing by the aerodynamic generation of nontorque bending moments on the main shaft. Potentially relevant to main bearing failure mechanisms, both computational and field data show that the magnitudes of turbulence-generated nontorque bending moments, that we show generate nonsteady force on the main bearing, are of order, and often larger than, torque (which underlies power). However, the temporal variations in these two responses are uncorrelated, implying that the aerodynamic mechanisms that drive power and main bearing response are fundamentally different. We find this to be the case in the field with both mountain-generated eddies (westerly winds) and ABL-generated eddies (northerly/southerly winds). Whereas the time and length scales are comparable, the mountain eddies were somewhat more energetic than the northerly/southerly ABL eddies. Interestingly, however, the fluctuations in nontorque bending moment that force the main bearing were found to be stronger when forced by the ABL eddies than the mountain eddies. The field studies validate the key results from the computational study and show even stronger response in the nontorque bending moment than in the computer simulations. In all cases, the torque and nontorque bending moments are temporally uncorrelated, torque and power are driven by time variations in rotor-averaged horizontal wind velocity and nontorque bending moments are driven by time changes in the degree of nonuniformity in the distribution of velocity over the rotor plane. Thus the results generalize the mechanisms underlying nonsteady aerodynamic forcing to classes of turbulence eddy types with strength of order or stronger than ABL eddies with transverse scale of order the wind turbine rotor. These include atmospheric turbulence eddies, topography-generated turbulence eddies and, by extension, impacts of turbine-wake-scale turbulence eddies on downstream wind turbine rotors.

17 WIND ENERGY↗
DOE OSTI2026

The remarkable inefficiency of stratocumulus

Marine stratocumulus clouds play a central role in Earth's climate system by reflecting incoming solar radiation and exerting a strong cooling effect. Their organization into open and closed mesoscale cellular morphologies can be thought of as an example of bistable dynamics driven by aerosol–cloud interactions and mesoscale processes. From the perspective of non-equilibrium thermodynamics, these structures are an example of a far-from-equilibrium open system that continuously produces and exports entropy. While entropy production has been studied in idealized deep convective systems, it has not yet been quantified for shallow clouds. Here, we compute and decompose the internal entropy production of open- and closed-cell stratocumulus using an ensemble of large-eddy simulations. We show that the overall entropy production of stratocumulus is low, reflecting the limited vertical extent and corresponding reduced ability to utilize the energy fluxes at the system's boundaries. Moist processes dominate the overall irreversibility, which, combined with their low entropy production, leads to a mechanical efficiency about an order of magnitude smaller than in deep convective systems. Although the dominant irreversible processes differ between open- and closed-cell regimes, the distributions of total entropy production largely overlap across the ensemble, limiting the ability to distinguish the dynamics of individual cases based solely on total entropy production.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Influence of simple terrain on the spatial variability of a low-level jet and wind farm performance in the AWAKEN field campaign

In wind energy research, scientific challenges are often associated with complex terrain sites, where orography, vegetation, and buildings disrupt flow uniformity. However, even sites characterized as simple terrain can exhibit significant spatial variability in wind speed, particularly during stable boundary layers (SBLs) and low-level jets (LLJs). This study investigates these terrain interactions using both simulations and observations from the American WAKe ExperimeNt (AWAKEN). We employ a multiscale Weather Research and Forecasting (WRF) model simulation, integrating mesoscale forcing in the coarse domains and representing three rows of turbines from the King Plains wind farm as generalized actuator disks (GAD) in the large-eddy simulation (LES) domains. During a nocturnal LLJ event on 3 April 2023, the downstream, wake-affected turbine rows outperformed the upstream, unwaked row by 25 %–51 %. This counterintuitive result arises from terrain-induced streamwise variations in hub-height wind speed of approximately 4 m s −1 over 5 km – equivalent to ∼50 % of the upstream reference speed. This enhancement outweighs the wake-induced reduction in mean wind speed (∼12 %) and global blockage effects reported in the literature (∼1 %–3.4 %). The multiscale simulations capture the intra-farm spatial variability in power performance observed in SCADA data. Terrain-induced vertical displacement of the LLJ, coupled with large wind shear below the jet maximum, drives the substantial streamwise acceleration within the wind farm. These findings underscore the importance of accounting for spatial variability related to terrain, even in simple landscapes, particularly during LLJ conditions. Incorporating such effects into reduced-order modeling frameworks for wind farm design and control could significantly enhance their effectiveness.

17 WIND ENERGY↗
DOE OSTI2025

JHTDB-wind: a web-accessible large-eddy simulation database of a wind farm with virtual sensor querying

This paper introduces JHTDB-wind (https://turbulence.idies.jhu.edu/datasets/windfarms, last access: 11 November 2025), a publicly accessible database containing large-eddy simulation (LES) data from wind farms. Building on the framework of the Johns Hopkins Turbulence Database (JHTDB), which hosts direct numerical simulation (DNS) and some LES datasets of canonical turbulent flows, JHTDB-wind stores the 4D space–time history of the flow and provides users the ability to access and query the data via a web-based virtual sensor interface. The initial dataset comprises LES results from a large wind farm with 10×6 turbines, modeled using a filtered actuator line method, under conventionally neutral atmospheric conditions. These data comprise 1 h (hour) of flow field data (velocity, pressure, potential temperature deviation, subgrid-scale (SGS) eddy viscosity, and turbine forces, approximately 15 TB (terabytes) and wind turbine data – including both turbine-level operational quantities and blade-level aerodynamic quantities (approximately 1.3 TB) – stored in Zarr and Parquet formats, respectively. Data retrieval is facilitated by the giverny Python package, allowing remote users to query the database in Python or MATLAB (C and Fortran support are available for flow field data). This paper details the simulation setup and demonstrates data access through examples that analyze wind farm flow structures and turbine performance. The framework is extensible to future datasets, including the JHTDB-wind diurnal cycle simulation analyzed in Xiao et al. (2025).

17 WIND ENERGY↗
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.↗
NASA NTRS1998

Large-Eddy Simulations and Lidar Measurements of Vortex-Pair Breakup in Aircraft Wakes

Results of large-eddy simulations of an aircraft wake are compared with results from ground-based lidar measurements made at NASA Langley Research Center during the Subsonic Assessment Near-Field Interaction Flight Experiment field tests. Brief reviews of the design of the field test for obtaining the evolution of wake dispersion behind a Boeing 737 and of the model developed for simulating such wakes are given. Both the measurements and the simulations concentrate on the period from a few seconds to a few minutes after the wake is generated, during which the essentially two-dimensional vortex pair is broken up into a variety of three-dimensional eddies. The model and experiment show similar distinctive breakup eddies induced by the mutual interactions of the vortices, after perturbation by the atmospheric motions.

Lewellen, D. C.↗
NASA NTRS2004

Single-column Model Intercomparison for a Stably Stratified Atmospheric Boundary Layer

The parameterization of the stably stratified atmospheric boundary layer is a difficult issue, which has a large impact on the medium-range weather forecasts and on climate integrations. A non-strongly stratified arctic case is simulated by nineteen single-column turbulence schemes. The statistics from the Large-eddy simulation (LES) intercomparison made for the same case by eight different models are used as a guiding reference. The single-column parameterizations include research schemes and operational schemes from major forecast and climate research centres. First order schemes, a large number of turbulence kinetic energy closures, and other proposals have submitted results. There is a large spread in the results; in general, the operational schemes mix more efficiently than the research ones, and the TKE and other higher order closures give results closer to the LES statistics. The sensitivities of the schemes to the parameters of their turbulence closures are partially explored.

Cuxart, J.↗