Search NASA⌕ Search

SEARCH · Search NASA

Results for “Convective mixing layer”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A two-fluid single-column model of turbulent shallow convection. Part II: Single-column model formulation and numerics

The two-fluid single-column model of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) is extended to include moisture and horizontal wind shear. Turbulent kinetic energy is introduced as a prognostic variable, dependence on a diagnosed boundary-layer height is removed, and subfilter fluxes are approximated using a two-fluid version of a Mellor–Yamada scheme. Three mechanisms for entrainment and detrainment processes are introduced, which represent entrainment of unstable air at the surface, forced detrainment of air at the top of the boundary/cloud layers, and turbulent mixing that relaxes the convective fluid to a reference profile. A semi-implicit Eulerian discretization replaces the semi-implicit semi-Lagrangian implementation of Thuburn et al. (Quart. J. R. Meteorol. Soc., 2019, 145, 1535–1550) to improve numerical stability and conservation. The equations for the implicit time step are solved using a quasi-Newton method, which is shown to perform well in numerical tests for conservation and convergence. The two-fluid single-column model presented in this article will be applied to simulations of shallow cumulus convection in Part III.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A Methodology for Estimating the Energy and Moisture Budget of the Convective Boundary Layer Using Continuous Ground-Based Infrared Spectrometer Observations

Abstract Land–atmosphere interactions play a critical role in both the atmospheric water and energy cycles. Changes in soil moisture and vegetation alter the partitioning of surface water and energy fluxes, influencing diurnal evolution of the planetary boundary layer (PBL). The mixing-diagram framework has proven useful in understanding the evolution of the heat and moisture budget within the convective boundary layer (CBL). We demonstrate that observations from the Department of Energy Atmospheric Radiation Measurement (ARM) Southern Great Plains site provide all of the needed inputs needed for the mixing-diagram framework, allowing us to quantify the impact from the surface fluxes, advection, radiative heating, encroachment, and entrainment on the evolution of the CBL. Profiles of temperature and humidity retrieved from the ground-based infrared spectrometer [Atmospheric Emitted Radiance Interferometer (AERI)] are a critical component in this analysis. Large-eddy simulation results demonstrate that mean mixed-layer values derived are shown to be critical to close the energy and moisture budgets. A novel approach demonstrated here is the use of network of AERIs and Doppler lidars to quantify the advective fluxes of heat and moisture. The framework enables the estimation of the entrainment fluxes as a residual, providing a way to observe the entrainment fluxes without using multiple lidar systems. The high temporal resolution of the AERI observations enables the morning, midday, and afternoon evolution of the CBL to be quantified. This work provides a new way to use observations in this framework to evaluate weather and climate models. Significance Statement The energy and moisture budget of the planetary boundary layer (PBL) is influenced by multiple sources, and accurately representing this evolution in numerical models is critical for weather forecasts and climate predictions. The mixing-diagram approach, driven by profiling observations as illustrated here, provides a powerful way to quantify the contributions from each of these sources. In particular, the energy and moisture mixed into the PBL from above the PBL can be determined accurately from ground-based remote sensors using this approach.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Improved method for temporally interpolating radiosonde profiles in the convective boundary layer

A significantly improved technique for temporally interpolating radiosonde (RS) profiles of potential temperature and water vapor mixing ratio in the planetary boundary layer during daytime is introduced. The key innovation of this technique is its operation on a height grid normalized with the planetary boundary layer height. This study utilized a three-month dataset of three-hourly soundings from the Atmospheric Radiation Measurement Facility's Southern Great Plains site. The technique was evaluated for convective boundary layer cases, with the necessary boundary layer height data obtained from a ground-based infrared spectrometer. A total of 79 comparisons were conducted between reference soundings and interpolated profiles that did and did not employ height normalization. The results demonstrated a substantial improvement in the representation of interpolated profiles using the new technique, characterized by enhanced correlation, improved amplitude representation, and reduced bias for potential temperature, as well as improved correlation and reduced bias for water vapor mixing ratio.

convective boundary layer↗
DOE OSTI2025

Identifying synoptic controls on boundary layer thermodynamic and cloud properties in a regional forecast model

Although most of our understanding of boundary layer cloudiness is based on idealized, subtropical, barotropic marine environments, boundary layer clouds exist across a range of conditions. In this study, we use the Naval Research Laboratory's Coupled Ocean/Atmosphere Mesoscale Prediction System (COAMPS) and an automated cold-front-relative analysis framework to explore the boundary layer structure associated with low clouds across a transect through the cold front of a midlatitude synoptic cyclone. The model credibly captures boundary layer structure in line with conceptual models of midlatitude cyclones and ground-based observations at Graciosa Island in the Azores. The warm sector is conditionally unstable, with clouds that are too shallow and with too little liquid water, compared to cloud property retrievals from satellite and surface-based instruments. The cold-frontal region exhibits convection associated with weak stability and ascent. Northwest of the cold front, the boundary layer is well mixed, deeper, and capped by a strong inversion maintained by large-scale subsidence. Simulated clouds in frontal and post-frontal regions are mostly too thick, with too much liquid water and too little cloud base drizzle. The post-frontal clouds are associated with grid-scale updrafts, which appear to be the model's attempt to represent mesoscale organization of cellular convection typically observed in the cold sector of midlatitude cyclones. The deep, well-mixed post-frontal boundary layers and cloudiness are maintained by strong surface fluxes, as in cold air outbreaks. Our analysis framework serves as a unique approach to model verification, and our results offer insights into differences in boundary layer cloud behavior between subtropical and synoptic cold-sector regimes.

Eissner, Jordan M. [Univ. of Kansas, Lawrence, KS ↗
DOE OSTI2025

A large-eddy simulation exploration of the assumptions used in retrieving entrainment from a mixing diagram approach with ground-based remote sensors

Entrainment is a crucial component of the atmospheric boundary layer (BL) moisture and heat budget. While usually thought of as only entrainment flux, entrainment within the mixed layer budget equation is really composed of two terms: the flux of a property across the boundary separating the BL from the free troposphere and the change in the concentration of a property as the depth of the BL changes. In a recent study, Wakefield et al. (2023) used ground-based remote-sensing observations to estimate entrainment flux as the residual of a mixing diagram framework that was applied to the daytime convective boundary layer. This present work uses large-eddy simulation (LES) to examine how well this residual assumption for entrainment fluxes alone compares to the actual sum of those two entrainment terms derived from spatial averages of the LES output. We highlight the importance of the second entrainment term in closing the mixed layer budget and show that the residual assumption does not represent entrainment flux only but rather a total entrainment term when the boundary layer depth is changing.

Rosenberger, Tessa E. [Cleveland State Univ., Clev↗
DOE OSTI2024

Aerosol Vertical Turbulent Mass Flux Retrievals Through Novel Remote Sensing Algorithm

Abstract Integrated measurements of aerosol, radiation, cloud, and turbulent transport in the planetary boundary layer (PBL) are essential for understanding and modeling climate and air quality. Here, we developed a new technique for the identification of convective turbulent regions and deriving the vertical distribution of aerosol turbulent mass fluxes within PBL. The algorithm uses retrievals from coherent Doppler lidars and a high spectral resolution lidar. The technique was applied to study particle mass fluxes over 2 months (November–December 2020) during the campaign conducted at the DOE Atmospheric Radiation Measurement Southern Great Plains (SGP) site in Lamont, Oklahoma. The algorithm developed here is capable of continuously deriving vertically resolved (curtains) aerosol mass fluxes. Our data analysis shows that at the site, the 30‐min averaged fluxes at 135 m above the surface were mainly positive (upward) at ∼1 μg m −2 s −1 , suggesting that the surface is the primary source of the particle mass supplied to the boundary layer at the SGP site. Analyses of the individual case studies have revealed that not all the derived fluxes can be linked to surface emissions. Both positive and negative values in a range of ±5 μg m −2 s −1 can be caused by convective thermals interacting between the residual layer and the mixed layer and by rotation of the horizontal wind with the height. Large erroneous negative fluxes can also be caused by drizzling/precipitating clouds. We anticipate that the application of the current technique will lead to a more realistic representation of aerosol mass budgets and bidirectional mixing rates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Influence of Surface Aerosol Injection on Stratocumulus‐to‐Cumulus Transition: Cloud‐Surface Coupling and Background Aerosol Concentrations

The influence of surface aerosol injection on the stratocumulus‐to‐cumulus transition (SCT) is explored using large‐eddy simulations. We examine how cloud‐surface coupling (or the strength of the marine boundary layer (MBL) stratification that limits vertical turbulent mixing and convection) impacts the vertical transport of aerosols, and how injected aerosols influence cloud properties and associated cloud radiative effects during the SCT. By injecting aerosols at different stages of the SCT, noting that cloud is more decoupled from the surface over time due to entrainment warming, we find that cloud‐surface coupling significantly affects aerosol vertical transport. However, injection timing (before drizzle if any) does not notably affect the SCT and the efficiency of marine cloud brightening, because aerosol number concentrations due to injections at different times rapidly converge before the transition onset. By varying the background aerosol concentration, we find that injected aerosols can significantly extend the persistence of stratocumulus decks by suppressing precipitation in clean environments but have little impact on stratocumulus breakup with higher background aerosol concentrations due to saturated aerosol effects. In clean MBLs, the SCT‐delay‐induced increase in cloud fraction dominates the overall cooling effects in response to aerosols, followed by Twomey effects. These cooling effects are slightly offset by decreased liquid water path (LWP) due to entrainment drying. In polluted MBLs, the Twomey effect is more dominant, followed by cloud fraction adjustments, and these coolings are also partly offset by LWP adjustments. All the simulations are made in relatively small domains in which injected aerosols are homogenized over a short time scale.

Zhang, Haipeng [Univ. of Maryland, College Park, M↗
DOE OSTI2024

Hydrodynamical Simulations of Proton Ingestion Flashes in Type I X-Ray Bursts

We perform the first multidimensional fluid simulations of thermonuclear helium ignition underneath a hydrogen-rich shell. This situation is relevant to Type I X-ray bursts on neutron stars that accrete from a hydrogen-rich companion. Using the low-Mach number fluid code MAESTROeX, we investigate the growth of the convection zone due to nuclear burning, and the evolution of the chemical abundances in the atmosphere of the star. We also examine the convective boundary mixing processes that cause the evolution to differ significantly from previous one-dimensional simulations that rely on mixing-length theory. We find that the convection zone grows outward as penetrating fluid elements cool the overlying radiative layer, rather than directly from the increasing entropy of the convection zone itself. Simultaneously, these flows efficiently mix composition, carrying carbon out of and protons into the convection zone even before contact with the hydrogen shell. We discuss the implications of these effects for future modeling of these events and observations.

astrophysical fluid dynamics↗
DOE OSTI2023

Impact of Seasonal Snow-Cover Change on the Observed and Simulated State of the Atmospheric Boundary Layer in a High-Altitude Mountain Valley

The structure and evolution of the atmospheric boundary layer (ABL) under clear-sky fair weather conditions over mountainous terrain is dominated by the diurnal cycle of the surface energy balance and thus strongly depends on surface snow cover. We use data from three passive ground-based infrared spectrometers deployed in the East River Valley in Colorado's Rocky Mountains to investigate the response of the thermal ABL structure to changes in surface energy balance during the seasonal transition from low to high snow cover. Temperature profiles were retrieved from the infrared radiances using the optimal estimation physical retrieval Tropospheric Remotely Observed Profiling via Optimal Estimation. A nocturnal surface inversion formed in the valley during clear-sky days, which was subsequently mixed out during daytime with the development of a convective boundary layer when snow cover was low. Over high snow cover, a very shallow convective boundary layer formed, above which the inversion persisted through the daytime hours. We compare these observations to NOAA's operational High-Resolution-Rapid-Refresh model and find large warm biases on clear-sky days resulting from the model's inability to form strong nocturnal inversions and to maintain the stable stratification in the valley during daytime when there was snow on the ground. We suggest several factors contributing to the large model errors. These are (a) the inability of the model to represent well-developed thermally driven flows likely due to the too coarse horizontal grid spacing (3 km), (b) too much convective mixing during daytime, and (c) too strong vertical coupling between the valley atmosphere and the free troposphere.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Vertically Resolved Convective–Stratiform Echo-Type Identification and Convectivity Retrieval for Vertically Pointing Radars

Using data from the airborne HIAPER Cloud Radar (HCR), a partitioning algorithm (ECCO-V) that provides vertically resolved convectivity and convective versus stratiform radar-echo classification is developed for vertically pointing radars. The algorithm is based on the calculation of reflectivity and radial velocity texture fields that measure the horizontal homogeneity of cloud and precipitation features. The texture fields are translated into convectivity, a numerical measure of the convective or stratiform nature of each data point. The convective–stratiform classification is obtained by thresholding the convectivity field. Subcategories of low, mid-, and high stratiform, shallow, mid-, deep, and elevated convective, and mixed echoes are introduced, which are based on the melting-layer and divergence-level altitudes. As the algorithm provides vertically resolved classifications, it is capable of identifying different types of vertically layered echoes, and convective features that are embedded in stratiform cloud layers. Its robustness was tested on data from four HCR field campaigns that took place in different meteorological and climatological regimes. The algorithm was adapted for use in spaceborne and ground-based radars, proving its versatility, as it is adaptable not only to different radar types and wavelengths, but also different research applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Effect of Formation Heterogeneity on CO 2 Dissolution in Subsurface Porous Media

Dissolution trapping is one of the most dominant mechanisms for the secure storage of CO 2 injected in porous subsurface formations saturated with brine. This trapping mechanism is enhanced by convective mixing, which occurs due to the gravitational instability between the different fluid layers in the aquifer. The reservoir permeability also plays a crucial role in the dissolution rate and overall fluid flow dynamics during the density-driven convection in porous media. Herein this study investigates the role of complex heterogeneity, i.e., irregular permeability distribution in CO 2 dissolution, using a novel experimental approach to create medium permeability heterogeneity in Hele–Shaw cells. Complex subsurface transport phenomena such as a preferential dissolution path, CO 2 sweep efficiency, changes in finger morphology, and CO 2 concentration distribution are visualized by creating heterogeneous media. Experimental results showed that reservoir permeability heterogeneity causes significant channeling effects and poor sweep efficiency. A scaling relationship between average finger growth rate (Gr) and permeability (k) was obtained as Gr [m s –1 ] = 266.8k [m 2 ] + 1.20 × 10 –6 . Furthermore, the mass of CO 2 dissolved is calculated using the spectrophotometric method to characterize the convective instability. The convective flux was analyzed by comparing the experimental dissolution flux with the theoretical diffusion flux, calculating a maximum Sherwood number of 6.8. The study’s findings improve the current understanding of the CO 2 convection morphology in heterogeneous media, allowing better assessment of long-term CO 2 storage.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

TRACER-Sonde Ozonesonde and radiosonde data

The TRACER-Sonde campaign used electrochemical cell ozonesondes to measure vertical profiles of ozone, which can be a tracer for convective mixing, in conjunction with meteorological variables measured by a radiosonde. There were 32 days in which TRACER-Sonde had twice-daily ozonesondes from the TRACER M1 site at La Porte at ~11Z and ~15Z (6 am and 10 am LT). The ozonesonde profiles provide an additional method to validate boundary-layer height determinations and test assumptions about vertical mixing in pre-convective environments. There are six additional days with ~15Z radiosonde only (no ozone data) profiles.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Modification of wind turbine wakes by large-scale, convective atmospheric boundary layer structures

In this study, we consider the impact of large-scale, convective structures in an unstable atmospheric boundary layer on wind turbine wakes. Simulation data from a high-fidelity large-eddy simulation (LES) of the AWAKEN wind farm site matching unstable atmospheric conditions were analyzed, and both turbine performance and wake behavior were affected based on their location relative to the convective structures. Turbines located in updraft regions of the flow experienced lower inflow velocity and generated less power, but their wakes were observed to recover faster and saw greater turbulent kinetic energy mixing higher in the boundary layer. The opposite effect was found for turbines in the downdraft regions of the convective structures. A simplified model of this wake behavior was also developed based on a two-dimensional k–ε Reynolds-Averaged Navier–Stokes formulation. This simplified model included the effects of vertical transport, but could be efficiently solved as a parabolic system, and was found to capture similar wake modifications observed in the high-fidelity LES computations.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Large eddy simulation of a supersonic lifted hydrogen flame: Impacts of Lewis, turbulent Schmidt and Prandtl numbers

Parametric large eddy simulations (LES) of a supersonic lifted hydrogen flame are reported. The emphases are on two aspects: impacts of (1) Lewis number (Lei of the ith species) and (2) turbulent Schmidt and Prandtl numbers (⁠Sc t and Pr t ⁠) on supersonic turbulent flame and flow structures. Five cases are considered: species-specific Le i , Sc t = Sr t = 1.0 (C 0 ); unity Le i , Sc t = Pr t = 1.0 (C 1 ); species-specific Sc t = 0.5, Pr t = 1.0 (C 2 ); species-specific Le i , Sc t = 1.0, Pr t = 0.5 (C 3 ); and species-specific Le i , Sc t = Pr t = 0.5 (C 4 ). Numerical results of instantaneous and/or time-averaged species mole fractions, mixture fraction, heat release rate, flame base location, and mixed modes of premixed and diffusion combustion are compared between cases C 0 and C 1 . Differences in auto-ignition locations and strengths and flame structures and stabilization specify the impacts of Lewis number. They are triggered by different predictions of species mass and thermal diffusions at fuel-coflow and/or coflow-ambient air mixing layers. These differences are rationalized by a scale analysis of mass/thermal diffusion and convection for case C 0 , which suggests the relatively low but non-negligible former against the latter. Cases C 0 and C 2 –C 4 barely see differences in terms of instantaneous and/or time-averaged temperature, velocity, and mixed combustion modes except for further downstream areas where combustion occurs. Both Sc t and Pr t impose less significant influences than Lewis number, as sub-grid scale mass/thermal diffusion is subordinate to its resolved counterpart according to their scale analysis for case C 4 .

42 ENGINEERING↗
DOE OSTI2021

Connection Between Mass Flux Transport and Eddy Diffusivity in Convective Atmospheric Boundary Layers

Turbulence parameterizations for convective boundary layer in coarse-scale atmospheric models usually consider a combination of the eddy-diffusive transport and a non-local transport, typically in the form of a mass flux term, such as the widely adopted eddy-diffusivity mass-flux (EDMF) approach. These two types of turbulent transport are generally considered to be independent of each other. Using results from large-eddy simulations, here, we show that a Taylor series expansion of the updraft and downdraft mass-flux transport can be used to approximate the eddy-diffusivity transport in the atmospheric surface layer and the lower part of the mixed layer, connecting both eddy-diffusivity and mass-flux transport theories in convective conditions, which also quantifies departure from the Monin-Obukhov similarity (MOS) in the surface layer. This study provides a theoretical support for a unified EDMF parameterization applied to both the surface layer and mixed layer and highlights important correction required for surface models relying on MOS

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A Decentralized Approach for Modeling Organized Convection Based on Thermal Populations on Microgrids

Abstract In this study, a spectral model for convective transport is coupled to a thermal population model on a two‐dimensional horizontal “microgrid,” covering the typical gridbox size of general circulation models. The goal is to explore new ways of representing impacts of spatial organization in cumulus cloud fields. The thermals are considered the smallest building block of convection, with thermal life cycle and movement represented through binomial functions. Thermals interact through two simple rules, reflecting pulsating growth and environmental deformation. Long‐lived thermal clusters thus form on the microgrid, exhibiting scale growth and spacing that represent simple forms of spatial organization and memory. Size distributions of cluster number are diagnosed from the microgrid through an online clustering algorithm, and provided as input to a spectral multiplume eddy‐diffusivity mass flux scheme. This yields a decentralized transport system, in that the thermal clusters acting as independent but interacting nodes that carry information about spatial structure. The main objectives of this study are (a) to seek proof of concept of this approach, and (b) to gain insight into impacts of spatial organization on convective transport. Single‐column model experiments demonstrate satisfactory skill in reproducing two observed cases of continental shallow convection. Metrics expressing self‐organization and spatial organization match well with large‐eddy simulation results. We find that in this coupled system, spatial organization impacts convective transport primarily through the scale break in the size distribution of cluster number. The rooting of saturated plumes in the subcloud mixed layer plays a key role in this process.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Large Eddy Simulation of the Diurnal Cycle of Shallow Convection in the Central Amazon

Climate models often face challenges in accurately simulating the daily precipitation cycle over tropical land areas, particularly in the Amazon. One contributing factor may be the incomplete representation of the diurnal evolution of shallow cumulus (ShCu) clouds. This study aimed to enhance the understanding of the diurnal cycles of ShCu clouds—from formation to maturation and dissipation—over the Central Amazon (CAMZ). Using observational data from the Green Ocean Amazon 2014 (GoAmazon) campaign and large eddy simulation (LES) modeling, we analyzed the diurnal cycles of six selected pure ShCu cases and their composite behavior. Our results revealed a well-defined cycle, with cloud formation occurring between 10 and 11 local time (LT), maturity from 13 to 15 LT, and dissipation by 17–18 LT. The vertical extent of the liquid water mixing ratio and the intensity of the updraft mass flux were closely associated with increases in turbulent kinetic energy (TKE), enhanced buoyancy flux within the cloud layer, and reduced large-scale subsidence. We further analyzed the diurnal cycles of the convective available potential energy (CAPE), the convective inhibition (CIN), the Bowen ratio (BR), and the vertically integrated TKE in the mixed layer (ITKE-ML), exploring their relationships with the cloud base mass flux (Mb) and cloud depth across the six ShCu cases. ITKE-ML and Mb exhibited similar diurnal trends, peaking at approximately 14–15 LT. However, no consistent relationships were found between CAPE (or BR) and Mb. Similarly, comparisons of the cloud depth with CAPE, BR, ITKE-ML, CIN, and Mb revealed no clear relationships. Smaller ShCu clouds were sometimes linked to higher CAPE and lower CIN. It is important to emphasize that these findings are preliminary and based on a limited sample of ShCu cases. Further research involving an expanded dataset and more detailed analyses of the TKE budget and synoptic conditions is necessary. Such efforts would yield a more comprehensive understanding of the factors influencing ShCu clouds’ vertical development.

54 ENVIRONMENTAL SCIENCES↗