Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “entrainment rates”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

The diurnal susceptibility of subtropical clouds to aerosols

The diurnal susceptibility of clouds and their radiative properties to aerosols is examined during their Lagrangian transition from subtropical stratocumulus to shallow cumulus regimes. Using large-eddy simulations, we analyze the six-day evolution of an air mass along a 3800 km observed trajectory from the coast of Peru toward the equator. Pristine and polluted scenarios are simulated with forcing imposed from weather reanalysis. The polluted scenario exhibits stronger diurnal variations in cloud water, cloud fraction, and albedo, with enhanced nighttime entrainment and suppressed precipitation. The overall response of cloud properties and outgoing shortwave radiation to droplet number concentration follows a distinct diurnal pattern: strong positive cloud adjustments dominate at night and in the morning, while weak negative adjustments prevail in the afternoon. This cycle is driven by the competition between precipitation suppression, which enhances cloud water and coverage, and entrainment drying, which depletes them. In polluted conditions, enhanced entrainment leads to a deeper and more decoupled boundary layer that cannot be sustained by surface fluxes in the afternoon, resulting in negative cloud adjustments. The enhanced entrainment rate under polluted conditions is caused by the reduced sedimentation of cloud and precipitation water from the entrainment zone. While the Twomey effect dominates the diurnal average albedo response, the diurnal variation in the competing cloud adjustments lead to a near-neutral net susceptibility in the afternoon, highlighting the critical role of diurnally varying processes in aerosol-cloud interactions.

Geosciences↗
DOE OSTI2024

Response of the Current Climate to Land‐Ocean Contrasts in Parameterized Cumulus Entrainment

Abstract Cumulus entrainment substantially regulates the earth's climate but remains poorly constrained in global climate models. Recent studies have shown that cumulus bulk entrainment (or dilution) is particularly sensitive to continentality, with the entrainment rate in simulated maritime cumuli nearly double that of continental cumuli. This study examines the impacts of such land–ocean entrainment contrasts on the current climate using 21‐year simulations with the Geophysical Fluid Dynamics Laboratory's High‐Resolution Atmospheric Model (HIRAM). In response to a 25% reduction in the HIRAM entrainment parameter c 0 over land, precipitation over tropical land regions increases by up to 40%. Along with directly facilitating enhanced convective precipitation, this c 0 reduction induces an increase in soil moisture, which may contribute to a further enhancement of convective precipitation over land. A 25% c 0 reduction over the oceans leads to more widespread modifications of convection patterns, with the strongest signal in the tropical Pacific. Deep convection shifts upstream (eastward) there, inducing enhanced large‐scale ascent over the central Pacific with compensating subsidence and reduced humidity and precipitation over the western Pacific (WP). Land–ocean variations in c 0 project onto the Pacific Walker circulation, with the 25% land reduction strengthening it by 4% and the 25% ocean reduction weakening it by 14%. These changes are driven by variations in convective and large‐scale stratiform heating over the Pacific. While reduced c 0 over land enhances diabatic heating in the Maritime Continent to strengthen the Walker circulation, reduced c 0 over the oceans decreases diabatic heating in the WP to weaken the Walker circulation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Shallow cumulus entrainment: observational retrieval, physical interpretation, and climate impacts (Final grant report)

This document constitutes a final technical report for DE-SC0020083, a research grant awarded by the US Department of Energy Office of Science, to Prof. Daniel Kirshbaum at McGill University. The original research proposal, entitled "Shallow cumulus entrainment: observational retrieval, physical interpretation, and climate impacts", defined four scientific objectives that synthesized ARM observations with turbulence-permitting large-eddy simulations (LES) to quantify and understand the process of cumulus entrainment, whereby a cumulus cloud ingests surrounding air through its lateral edges. These objectives included developing observational climatologies of cumulus bulk entrainment rate at two ARM observatories, physical interpretation of these findings using LES, quantification of the impacts of subcloud circulations on entrainment, and the climate impacts of flow-dependent entrainment variations. In the six years that have elapsed since the grant was awarded (three funded years and three no-cost annual extensions), all four of the objectives were successfully addressed. Four personnel were trained (one B.Sc., one M.Sc. to completion, five years of Ph.D. training, and 2.5 years of postdoctoral training), eight peer-reviewed articles were published in high-impact journals, and 20 presentations were delivered (three invited).

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1998

Mixing Enhancement by Tabs in Round Supersonic Jets

The objective of this study was to analyze jet plume mass flow entrainment rates associated with the introduction of counter-rotating streamwise vorticity by prism shaped devices (tabs) located at the lip of the nozzle. We have examined the resulting mixing process through coordinated experimental tests and numerical simulations of the supersonic flow from a model axisymmetric nozzle. In the numerical simulations, the total induced vorticity was held constant while varying the distribution of counter-rotating vorticity around the nozzle lip training edge. In the experiment, the number of tabs applied was varied while holding the total projected area constant. Evaluations were also conducted on initial vortex strength. The results of this work show that the initial growth rate of the jet shear layer is increasingly enhanced as more tabs are added, but that the lowest tab count results in the largest entrained mass flow. The numerical simulations confirm these results.

Seiner, John M.↗
DOE OSTI2026

Examination of entrainment, combustibility, and heat transfer due to cell venting inside simplified battery energy storage enclosures

Battery energy storage systems (ESS) assembled of modules containing lithium-ion batteries can pose fire and explosion hazards during thermal runaway because vented gases may accumulate, combust, or escape enclosures. Here, we propose a prototypical configuration composed of the main components of an ESS rack, where modules are represented as rectangles within an enclosure and vent gas is injected steadily through the top of a module. This idealized geometry allows for exploration of how rack-level characteristics can influence smaller single-cell or module-level scales. Through this geometry, we conduct numerical simulations to understand how variations in the vent velocity, temperature, and vertical location affect entrainment rates, fraction of uncombusted fuel, rich-gas coverage, and heat flux to other modules. At lower vent velocities and vertical positions, buoyancy drives enough entrainment and mixing to react with most of the vent gases on the converging plenum. As jet momentum increases and/or buoyancy decreases via higher vent position in the rack, vent gases can propagate into and burn in the divergence plenum. Heat flux to the impinged module scales strongly with vent momentum and temperature, while heat transfer to upper modules peaks when flames in the diverging plenum reach the surface and declines when combustion is suppressed.

Computational fluid dynamics↗
DOE OSTI2021

A Lagrangian Study of the Transition from Shallow to Deep Convection Using ASR Observations and Les Simulations (Final Report)

With the support from this grant, we have made the following accomplishments. We have found that a formula for fractional entrainment rate, a key aspect for parameterization moist convection, that works for both shallow and deep convection, which provide a promising direction for unified treatment of entrainment for both shallow and deep convection. We also examined the formulation of vertical velocity and showed that the effect of non-hydrostatic pressure and entrainment on vertical velocity are not constants, as assumed in past studies, but can be computed using a simple model for the cloudy updrafts. These results also apply to both deep and shallow convection, which opens the door for a more physically based parameterization of the updraft vertical velocity. We have also conducted extensive analysis of the RACORO aircraft data and have published the results in Wei et al. (2020). We have made extensive effort scrutinizing the data in detail and examining the various hypotheses with the aid of the LES simulations. This was the first application of the modified Paluch diagram on modern aircraft data and the first comprehensive analysis of mixing in clouds using real data (instead of only LES output) in a very long time. These analyses enhanced the science return of the ASR aircraft measurements. This grant also supported a study of the collisions among cold pools and an effort to build parameterizations of moist convection using tools developed in control theory.

42 ENGINEERING↗
NASA NTRS1984

Can cloud-top entrainment promote cloud growth?

The primary significance of Cloud Deepening through Entrainment (CDE) is that it can prevent the cloud top entrainment instability from destroying a cloud deck. Without suppressing the instability, CDE transforms it from a cloud destroyer to a cloud builder. The analysis does not depend on an entrainment hypothesis. Moreover, it is not restricted to PBL stratocumulus sheets. Stratiform clouds in the free atmosphere can be subject to CDE we need only reinterpret Ps as the pressure at the base of an elevated turbulent mixed layer. Modest departures from well mixedness will alter the results quantitatively but not qualitatively. Processes other than entrainment, such as surface evaporation, radiative cooling, and advection will often work with CDE to build a cloud layer; but of course they can also oppose CDE by reducing the relative humidity. If we make the weak assumption that the deepening of a cloud layer favors an increase in the cloud top entrainment rate (without specifying any particular functional relationship) we are led to speculate that CDE can cause runaway cloud growth, even in the absence of cloud top entrainment instability. through CDE entrainment leads to a deeper cloud, which leads to stronger entrainment.

Randall, D. A.↗
NASA NTRS1985

Lidar observations of vertically organized convection in the planetary boundary layer over the ocean

Observations of a convective planetary boundary layer (PBL) were made with an airborne, downward-looking lidar system over the Atlantic Ocean during a cold air outbreak. The lidar data revealed well-organized, regularly spaced cellular convection with dominant spacial scales between two and four times the height of the boundary layer. It is demonstrated that the lidar can accurately measure the structure of the PBL with high vertical and horizontal resolution. Parameters important for PBL modeling such as entrainment zone thickness, entrainment rate, PBL height and relative heat flux can be inferred from the lidar data. It is suggested that wind shear at the PBL top may influence both entrainment and convective cell size.

Melfi, S. H.↗
DOE OSTI2020

Dynamics of Subsiding Shells in Actively Growing Clouds with Vertical Updrafts

The dynamics of a subsiding shell at the edges of actively growing shallow cumulus clouds with updrafts is analyzed using direct numerical simulation. The actively growing clouds have a fixed in-cloud buoyancy and velocity. Turbulent mixing and evaporative cooling at the cloud edges generate a subsiding shell that grows with time. A self-similar regime is observed for first- and second-order moments when normalized with respective maximum values. Internal scales derived from integral properties of the flow problem are identified. A self-similarity analysis using these scales reveals that contrary to classical self-similar flows, the turbulent kinetic energy budget terms and velocity moments scale according to the buoyancy and not with the mean velocity. The shell thickness is observed to increase linearly with time. The buoyancy scale remains time invariant and is set by the initial cloud–environment thermodynamics. The shell accelerates ballistically with a magnitude set by the saturation value of the buoyancy of the cloud–environment mixture. In this regime, the shell is buoyancy driven and independent of the in-cloud velocity. Relations are obtained for predicting the shell thickness and minimum velocities by linking the internal scales with external flow parameters. The values thus calculated are consistent with the thickness and velocities observed in typical shallow cumulus clouds. The entrainment coefficient is a function of the initial state of the cloud and the environment, and is shown to be on the same order of magnitude as fractional entrainment rates calculated for large-scale models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Cumulonimbus Clouds Convert a Smaller Fraction of CAPE into Kinetic Energy in a Warmer Atmosphere

Abstract This study investigates how entrainment’s diluting effect on cumulonimbus updraft buoyancy is affected by the temperature of the troposphere, which is expected to increase by the end of the century. A parcel model framework is constructed that allows for independent variations in the temperature ( T ), the entrainment rate ε , the free-tropospheric relative humidity (RH), and the convective available potential energy (CAPE). Using this framework, dilution of buoyancy is evaluated with T and RH independently varied and with CAPE either held constant or increased with temperature. When CAPE is held constant, buoyancy decreases as T increases, with parcels in warmer environments realizing substantially smaller fractions of their CAPE as kinetic energy (KE). This occurs because the increased moisture difference between an updraft and its surroundings at warmer temperatures drives greater updraft dilution. Similar results are found in midlatitude and tropical conditions when CAPE is increased with temperature. With the expected 6%–7% increase in CAPE per kelvin of warming, KE only increases at 2%–4% K −1 in narrow updrafts but tracks more closely with CAPE at 4%–6% in wider updrafts. Interestingly, the rate of increase in the KE with T becomes larger than that of CAPE when the later quantity increases at more than 10% K −1 . These findings emphasize the importance of considering entrainment in studies of moist convection’s response to climate change, as the entrainment-driven dilution of buoyancy may partially counteract the influence of increases in CAPE on updraft intensity. Significance Statement Cumulonimbus clouds mix air with their surrounding environment through a process called entrainment, which controls how efficiently environmental energy is converted into upward speed in thunderstorm updrafts. Our research shows that warmer temperatures will exacerbate the moisture difference between cumulonimbus updrafts and their surroundings, leading to greater mixing and less efficient conversion of environmental energy into updraft speeds. This effect should be considered in future research that investigates how climate change will affect cumulonimbus clouds.

Peters, John M.↗
DOE OSTI2024

Fluid-Dynamic Mechanisms Underlying Wind Turbine Wake Control with Strouhal-Timed Actuation

A reduction in wake effects in large wind farms through wake-aware control has considerable potential to improve farm efficiency. This work examines the success of several emerging, empirically derived control methods that modify wind turbine wakes (i.e., the pulse method, helix method, and related methods) based on Strouhal numbers on the O(0.3). Drawing on previous work in the literature for jet and bluff-body flows, the analyses leverage the normal-mode representation of wake instabilities to characterize the large-scale wake meandering observed in actuated wakes. Idealized large-eddy simulations (LES) using an actuator-line representation of the turbine blades indicate that the n=0 and ±1 modes, which correspond to the pulse and helix forcing strategies, respectively, have faster initial growth rates than higher-order modes, suggesting these lower-order modes are more appropriate for wake control. Exciting these lower-order modes with periodic pitching of the blades produces increased modal growth, higher entrainment into the wake, and faster wake recovery. Modal energy gain and the entrainment rate both increase with streamwise distance from the rotor until the intermediate wake. This suggests that the wake meandering dynamics, which share close ties with the relatively well-characterized meandering dynamics in jet and bluff-body flows, are an essential component of the success of wind turbine wake control methods. A spatial linear stability analysis is also performed on the wake flows and yields insights on the modal evolution. In the context of the normal-mode representation of wake instabilities, these findings represent the first literature examining the characteristics of the wake meandering stemming from intentional Strouhal-timed wake actuation, and they help guide the ongoing work to understand the fluid-dynamic origins of the success of the pulse, helix, and related methods.

17 WIND ENERGY↗
NASA NTRS1986

Fire control method and analytical model for large liquid hydrocarbon pool fires

The dominate parameter governing the behavior of a liquid hydrocarbon (JP-5) pool fire is wind speed. The most effective method of controlling wind speed in the vicinity of a large circular (10 m dia.) pool fire is a set of concentric screens located outside the perimeter. Because detailed behavior of the pool fire structure within one pool fire diameter is unknown, an analytical model supported by careful experiments is under development. As a first step toward this development, a regional pool fire model was constructed for the no-wind condition consisting of three zones -- liquid fuel, combustion, and plume -- where the predicted variables are mass burning rate and characteristic temperatures of the combustion and plume zones. This zone pool fire model can be modified to incorporate plume bending by wind, radiation absorption by soot particles, and a different ambient air flow entrainment rate. Results from the zone model are given for a pool diameter of 1.3 m and are found to reproduce values in the literature.

Fenton, D. L.↗
NASA NTRS2000

The Turbulent/Non-Turbulent Interface Bounding a Far-Wake

The velocity fields of a turbulent wake behind a flat plate obtained from the direct numerical simulations of Moser et al. are used to study the structure of the flow in the intermittent zone where there are, alternately, regions of fully turbulent flow and non-turbulent velocity fluctuations either side of a thin randomly moving interface. Comparisons are made with a wake that is 'forced' by amplifying initial velocity fluctuations. There is also a random temperature field T in the flow; T varies between constant values of 0.0 and 1.0 on the sides of the wake. The value of the Reynolds number based on the centreplane mean velocity defect and halfwidth b of the wake is Re approx. = 2000. It is found that the thickness of the continuous interface is about equal to 0.07b, whereas the amplitude of fluctuations of the instantaneous interface displacement y(sub I)(t) is an order of magnitude larger, being about 0.5b. This explains why the mean statistics of vorticity in the intermittent zone can be calculated in terms of the probability distribution of y(sub I) and the instantaneous discontinuity in vorticity across the interface. When plotted as functions of y - y(sub I), the conditional mean velocity (U) and temperature (T) profiles show sharp jumps Delta(U) and Delta(T) at the interface adjacent to a thick zone where (U) and (T) vary much more slowly. Statistics for the vorticity and velocity variances, available in such detail only from DNS data, show how streamwise and spanwise components of vorticity are generated by vortex stretching in the bulges of the interface. Flow fields around the interface, analyzed in terms of the local streamline pattern, confirm previous results that the advancement of the vortical interface into the irrotational flow is driven by large-scale eddy motion. It is argued that because this is an inviscid mechanism the entrainment process is not sensitive to the value of Re, and that small-scale nibbling only plays a subsidiary role. While mean Reynolds stresses decrease gradually in the intermittent zone, conditional stresses are found to decrease sharply towards zero at the interface. Using one-point turbulence models applied to either unconditional or conditional statistics for the turbulent region and then averaged, the entrainment rate E(sub b) would, if calculated exactly, be zero. But if computed with standard computational methods, E(sub b) would be non-zero because of numerical diffusion. It is concluded that the current practice in statistical models of approximating entrainment by a diffusion process is computationally arbitrary and physically incorrect. An analysis shows how E(sub b) is related to Delta(U) and the jump in shear stress at the interface, and correspondingly to Delta(T) and the heat flux.

Bisset, David K.↗
DOE OSTI2025

The Effect of Updraft Entrainment on Convective Cell Deepening in Realistic Large-Eddy Simulations

Entrainment of surrounding cooler and drier air into convective updrafts is one of the key processes that influence deep convection initiation and growth. Numerous studies have investigated the effect of entrainment on isolated convective cloud growth in idealized simulations, but the importance of this effect in realistic conditions with many interacting convective clouds remains uncertain. We examine the impact of entrainment on the depth reached by convective clouds in realistic large-eddy simulations (LES) over central Argentina during the Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign. Cloudy updrafts and their associated properties are assigned to convective cells tracked with radar reflectivity signatures. Several thousand convective cells are tracked over two high convective available potential energy (CAPE) and two low CAPE cases that support cells of varying depths and intensities. Entrainment is calculated explicitly as the fluxes of air into the outer surface of each cloudy updraft. Single-predictor logistic regression models are used to determine the relative importance of updraft, near-updraft, and preconvective initiation atmospheric conditions in predicting whether convective cells become deep. We then build a multiple-predictor regression model pairing important updraft and meteorological metrics with fractional entrainment rate. The probability of cells transitioning to deep convection is most sensitive to ambient 600-hPa relative humidity (42% of total metric contribution to cloud depth predictability), followed by low-level CAPE (28%), cloud-base updraft width (19%), and fractional entrainment (11%). Thus, the initial width of the updraft along with potential buoyancy and its dilution through the midtroposphere collectively determine whether deep convection will result from shallower clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Influences of Environmental Relative Humidity and Horizontal Scale of Subcloud Ascent on Deep Convective Initiation

This study examines two factors impacting initiation of moist deep convection: free-tropospheric environmental relative humidity (φ E ) and horizontal scale of subcloud ascent (R sub ), the latter exerting a dominant control on cumulus cloud width. A simple theoretical model is used to formulate a “scale selection” hypothesis: that a minimum R sub is required for moist convection to go deep, and that this minimum scale decreases with increasing φ E . Specifically, the ratio of R 2 sub to saturation deficit (1 - φ E ) must exceed a certain threshold value that depends on cloud-layer environmental lapse rate. Idealized, large-eddy simulations of moist convection forced by horizontally varying surface fluxes show strong sensitivity of maximum cumulus height to both φ E and R sub consistent with the hypothesis. Increasing R sub by only 300–400 m can lead to a large increase (>5 km) in cloud height. A passive tracer analysis shows that the bulk fractional entrainment rate decreases rapidly with R sub but depends little on φ E . However, buoyancy dilution increases as either R sub or φ E decreases; buoyancy above the level of free convection is rapidly depleted in dry environments when R sub is small. While deep convective initiation occurs with an increase in relative humidity of the near environment from moistening by earlier convection, the importance of this moisture preconditioning is inconclusive as it is accompanied by an increase in R sub . Overall, it is concluded that small changes to R sub driven by external forcing or by convection itself could be a dominant regulator of deep convective initiation.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1980

Entrainment into a stratocumulus layer with distributed radiative cooling

It is shown that the radiative cooling of a cloud layer strongly influences the turbulent flux profiles and the entrainment rate, and that the radiative cooling should be modeled as acting inside the turbulent layer. Numerical experiments demonstrate that a cloud-topped mixed-layer model, similar to that of Lilly (1968), is quite sensitive to delta (p) sub R, the depth of the radiatively cooled layer near cloud top. As delta (p) sub R increases, the model's sensitivity to the entrainment assumption is markedly heightened; for large delta (p) sub R the cloud top and cloud base rise as the entrainment parameter k is increased, while for small delta (p) sub R an increase in k has almost no effect. The model is most sensitive to delta (p) sub R for the cold-water, strong-divergence regime of greatest interest.

Randall, D. A.↗
DOE OSTI2024

Development of optimal methods for collection, transportation, holding, handling, and tagging of juvenile American shad

American shad (Alosa sapidissima) are an anadromous fish species native to the United States that have an extensive range along the east coast, but their populations are declining worldwide. Acoustic telemetry can play a vital role in better understanding the behavior and survival of this sensitive species. A stepwise, systematic approach was used to assess juvenile shad survival during multiple stressors typically associated with an acoustic telemetry study. Minimizing out-of-water handling and the use of saltwater during collection increased 24-hour survival (78 to 99%) between the first two shad transports. Several experiments were conducted to determine the best transmitter design method for implanting juvenile shad. Saltwater was fundamental in keeping tagged shad alive overnight and during the tagging process. Shad as small as 50 mm were tagged using the pectoral incision method with no suture. In a 60-d holding evaluation, tagged fish survived at a comparable rate as their non-tagged counterparts (81.5% for tagged, 70% for untagged). Also, when tagged and untagged shad were exposed to a tank of predators, there was no difference in their survival. Results from this study are encouraging for improving conservation efforts for small, sensitive species of fish, like American shad. The techniques identified may help broaden the application of telemetry studies by informing researchers of optimal methods to ensure survival of tagged fish in future field studies looking to evaluate survival and entrainment rates of juvenile American shad.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Using Machine Learning to Generate a GISS ModelE Calibrated Physics Ensemble (CPE)

A neural network (NN) surrogate of the NASA GISS ModelE atmosphere (version E3) is trained on a perturbed parameter ensemble (PPE) spanning 45 physics parameters and 36 outputs. The NN is leveraged in a Markov Chain Monte Carlo (MCMC) Bayesian parameter inference framework to generate a second posterior constrained ensemble coined a “calibrated physics ensemble,” or CPE. The CPE members are characterized by diverse parameter combinations and are, by definition, close to top-of-atmosphere radiative balance, and must broadly agree with numerous hydrologic, energy cycle and radiative forcing metrics simultaneously. Global observations of numerous cloud, environment, and radiation properties (provided by global satellite products) are crucial for CPE generation. The inference framework explicitly accounts for discrepancies (or biases) in satellite products during CPE generation. We demonstrate that product discrepancies strongly impact calibration of important model parameter settings (e.g., convective plume entrainment rates; fall speed for cloud ice). Structural improvements new to E3 are retained across CPE members (e.g., stratocumulus simulation). Notably, the framework improved the simulation of shallow cumulus and Amazon rainfall while not degrading radiation fields, an upgrade that neither default parameters nor Latin Hypercube parameter searching achieved. Analyses of the initial PPE suggested several parameters were unimportant for output variation. However, many “unimportant” parameters were needed for CPE generation, a result that brings to the forefront how parameter importance should be determined in PPEs. From the CPE, two diverse 45-dimensional parameter configurations are retained to generate radiatively-balanced, auto-tuned atmospheres that were used in two E3 submissions to CMIP6.

54 ENVIRONMENTAL SCIENCES↗