Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vertical velocities”

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 73 records · Page 4

Collaborative Proposal: Improving Understanding of the Internal Structure and Dynamics of Deep Convection Using ARM Observations and Large Eddy Simulations (Final progress report)

This is the final technical report for the DOE Atmospheric System Research funded project entitled "Collaborative Proposal: Improving Understanding of the Internal Structure and Dynamics of Deep Convection Using ARM Observations and Large Eddy Simulations". This project addressed several science questions related to the structure and behavior of deep convective clouds in the atmosphere. It applied knowledge gained from this work to improve the deep cumulus convection parameterization in the Community Earth System Model (CESM). Specific outcomes related to the project include an improved understanding of: (1) how convective cloud environments regulate the size of thermals reaching the upper troposphere, and hence the shallow-to-deep convective transition; (2) the role of dynamic pressure forcing in thermal and convective cloud evolution; and (3) the role of vertical wind shear on the behavior of moist thermals and the transition from shallow to deep convection. Although the primary motivation of this project was to improve basic understanding of thermal dynamics in atmospheric convection, an important component of this work was using this improved understanding to modify the representation of entrainment and vertical velocity profiles in the Zhang-McFarlane cumulus parameterization used in CESM. This work also had a strong observational component, with analyses of Atmospheric Radiation Measurement (ARM) observations to characterize thermal characteristics and behavior. Observations collected during the ARM-supported CACTI field project were used with high-resolution simulations to explore environmental controls on deep convection initiation. Data from the ARM-supported MC3E field project were also used to evaluate the representation of updraft vertical velocities in the Zhang-McFarlane cumulus parameterization. This project directly supported the training of two post-doctoral researchers and led to 16 peer-reviewed published papers.

54 ENVIRONMENTAL SCIENCES↗

Diffusional growth of cloud droplets in homogeneous isotropic turbulence: DNS, scaled-up DNS, and stochastic model

This paper presents a novel methodology to use direct numerical simulation (DNS) to study the impact of isotropic homogeneous turbulence on the condensational growth of cloud droplets. As shown by previous DNS studies, the impact of turbulence increases with the computational domain size, that is, with the Reynolds number, because larger eddies generate higher and longer-lasting supersaturation fluctuations that affect growth of individual cloud droplets. The traditional DNS can only simulate a limited range of scales because of the excessive computational cost that comes from resolving all scales involved, that is, from large scales at which the turbulent kinetic energy (TKE) is introduced down to the Kolmogorov microscale, and from following every single droplet. The novel approach is referred to as the “scaled-up DNS”. The scaling up is done in two parts, first by increasing both the computational domain and the Kolmogorov microscale and second by using super-droplets instead of real droplets. To ensure proper dissipation of TKE and scalar variance at small scales, molecular transport coefficients are appropriately scaled up with the grid length. For the scaled-up domains, say, meters and tens of meters, one needs to follow billions of real droplets. This is not computationally feasible, and so-called super-droplets are applied in scaled-up DNS simulations. Each super-droplet represents an ensemble of identical real droplets, and the number of real droplets represented by a super-droplet is referred to as the multiplicity attribute. After simple tests showing the validity of the methodology, scaled-up DNS simulations are conducted for five domains, the largest of 64 3 m 3 volume using a DNS of 256 3 grid points and various multiplicities. All simulations are carried out with vanishing mean vertical velocity and with no mean supersaturation, similarly to past DNS studies. As expected, the supersaturation fluctuations as well as the spread in droplet size distribution increase with the domain size, with the droplet radius variance increasing in time t as t 1∕2 as identified in previous DNS studies. Scaled-up simulations with different multiplicities document numerical convergence of the scaled-up solutions. Finally, we compare the scaled-up DNS results with a simple stochastic model that calculates supersaturation fluctuations based on the vertical velocity fluctuations updated using the Langevin equation. Overall, the results document similar scaling to previous small-domain DNS simulations and support the notion that the stochastic subgrid-scale model is a valuable tool for the multi-scale simulation of droplet spectral evolution applying a large-eddy simulation model.

54 ENVIRONMENTAL SCIENCES↗

An Investigation of LES Wall Modeling for Rayleigh–Bénard Convection via Interpretable and Physics-Aware Feedforward Neural Networks with DNS

Abstract The traditional approach of using the Monin–Obukhov similarity theory (MOST) to model near-surface processes in large-eddy simulations (LESs) can lead to significant errors in natural convection. In this study, we propose an alternative approach based on feedforward neural networks (FNNs) trained on output from direct numerical simulation (DNS). To evaluate the performance, we conduct both a priori and a posteriori tests. In the a priori (offline) tests, we compare the statistics of the surface shear stress and heat flux, computed from filtered DNS input variables, to the stress and flux obtained from the filtered DNS. Additionally, we investigate the importance of various input features using the Shapley additive explanations value and the conditional average of the filter grid cells. In the a posteriori (online) tests, we implement the trained models in the System for Atmospheric Modeling (SAM) LES and compare the LES-generated surface shear stress and heat flux with those in the DNS. Our findings reveal that vertical velocity, a traditionally overlooked flow quantity, is one of the most important input features for determining the wall fluxes. Increasing the number of input features improves the a priori test results but does not always improve the model performance in the a posteriori tests because of the differences in input variables between the LES and DNS. Last, we show that physics-aware FNN models trained with logarithmic and scaled parameters can well extrapolate to more intense convection scenarios than in the training dataset, whereas those trained with primitive flow quantities cannot. Significance Statement The traditional near-surface turbulence model, based on a shear-dominated boundary layer flow, does not represent near-surface turbulence in natural convection. Using a feedforward neural network (FNN), we can construct a more accurate model that better represents the near-surface turbulence in various flows and reveals previously overlooked controlling factors and process interactions. Our study shows that the FNN-generated models outperform the traditional model and highlight the importance of the near-surface vertical velocity. Furthermore, the physics-aware FNN models exhibit the potential to extrapolate to convective flows of various intensities beyond the range of the training dataset, suggesting their broader applicability for more accurate modeling of near-surface turbulence.

54 ENVIRONMENTAL SCIENCES↗

An Analytic Formula for Entraining CAPE in Midlatitude Storm Environments

Abstract This article introduces an analytic formula for entraining convective available potential energy (ECAPE) with an entrainment rate that is determined directly from an environmental sounding, rather than prescribed by the formula user. Entrainment is connected to the background environment using an eddy diffusivity approximation for lateral mixing, updraft geometry assumptions, and mass continuity. These approximations result in a direct correspondence between the storm-relative flow and the updraft radius and an inverse scaling between the updraft radius squared and entrainment rate. The aforementioned concepts, combined with the assumption of adiabatic conservation of moist static energy, yield an explicit analytic equation for ECAPE that depends entirely on state variables in an atmospheric profile and a few constant parameters with values that are established in past literature. Using a simplified Bernoulli-like equation, the ECAPE formula is modified to account for updraft enhancement via kinetic energy extracted from the cloud’s background environment. CAPE and ECAPE can be viewed as predictors of the maximum vertical velocity w max in an updraft. Hence, these formulas are evaluated using w max from past numerical modeling studies. Both of the new formulas improve predictions of w max substantially over commonly used diagnostic parameters, including undiluted CAPE and ECAPE with a constant prescribed entrainment rate. The formula that incorporates environmental kinetic energy contribution to the updraft correctly predicts instances of exceedance of by w max , and provides a conceptual explanation for why such exceedance is rare among past simulations. These formulas are potentially useful in nowcasting and forecasting thunderstorms and as thunderstorm proxies in climate change studies. Significance Statement Substantial mixing occurs between the upward-moving air currents in thunderstorms (updrafts) and the surrounding comparatively dry environmental air, through a process called entrainment. Entrainment controls thunderstorm intensity via its diluting effect on the buoyancy of air within updrafts. A challenge to representing entrainment in forecasting and predictions of the intensity of updrafts in future climates is to determine how much entrainment will occur in a given thunderstorm environment without a computationally expensive high-resolution simulation. To address this gap, this article derives a new formula that computes entrainment from the properties of a single environmental profile. This formula is shown to predict updraft vertical velocity more accurately than past diagnostics, and can be used in forecasting and climate prediction to improve predictions of thunderstorm behavior and impacts.

54 ENVIRONMENTAL SCIENCES↗

Remote Influence of Andean Convection on Amazonian Rainfall and Its Mechanisms

Models from Coupled Model Intercomparison Project Phase 6 produce too much precipitation over the Andes but too little over the Amazon or the Wet Andes-Dry Amazon (WADA) bias pattern. Unlike the conventional view that convection parameterization and land model deficiencies can contribute to Amazonian rainfall biases, we approach this long-standing biased model behavior through the lens of Andean convection. Using Community Earth System Model v1.1 and focusing on the wet season, our mechanism-denial experiments demonstrate that Andean convection notably reduces precipitation over the Amazon during austral summer. The Andean forced Amazonian response operates on weather timescale. Furthermore, the reduction of Amazonian rainfall is detectable within a few hours after initial Andean forcing. The precipitation response is primarily driven by variations in the moisture budget and is moderated by changes in convective available potential energy over the Amazon. Changes in the total advection of moisture over the Amazon are dominated by the vertical advection term and can be attributed to discrepancies in the dynamic omega field. In the experiments, the Andean east flank region is scrutinized where the vertical velocity and moisture fields play an intermediary role for the Andean driven WADA connection. The Andean forcing induces descending anomalies on the Andean east flank. The disturbances of wind and geopotential fields over the Andean east flank propagate eastward via Kelvin waves. Over the Amazon, descending anomalies and advective drying lead to reduction of mid-to-high level cloud, increase of shortwave cloud forcing and surface net radiation, and enhancement of themodynamic stability and rainfall reduction.

58 GEOSCIENCES↗

Aerosol-induced closure of marine cloud cells: enhanced effects in the presence of precipitation

Abstract. The Weather Research Forecasting (WRF) version 4.3 model is configured within a Lagrangian framework to quantify the impact of aerosols on evolving cloud fields. Kilometer-scale simulations utilizing meteorological boundary conditions are based on 10 case study days offering diverse meteorology during the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA). Measurements from aircraft, the ground-based Atmosphere Radiation Measurement (ARM) site at Graciosa Island in the Azores, and A-Train and geostationary satellites are utilized for validation, demonstrating good agreement with the WRF-simulated cloud and aerosol properties. Higher aerosol concentration leads to suppressed drizzle and increased cloud water content in all case study days. These changes lead to larger radiative cooling rates at cloud top, enhanced vertical velocity variance, and increased vertical and horizontal wind speed near the base of the lower-tropospheric inversion. As a result, marine cloud cell area expands, narrowing the gap between shallow clouds and increasing cloud optical thickness, liquid water content, and the top-of-atmosphere outgoing shortwave flux. While similar aerosol effects are observed in lightly to non-raining clouds, they tend to be smaller by comparison. These simulations show a relationship between cloud cell area expansion and the radiative adjustments caused by liquid water path and cloud fraction changes. The adjustments are positive and scale as 74 % and 51 %, respectively, relative to the Twomey effect. While higher-resolution large-eddy simulations may provide improved representation of cloud-top mixing processes, these results emphasize the importance of addressing mesoscale cloud-state transitions in the quantification of aerosol radiative forcing that cannot be attained from traditional climate models.

54 ENVIRONMENTAL SCIENCES↗

MLBS Halo scanning Lidar / Reviewed Data / Fixed-point vertical scans

This dataset contains high-frequency vertical velocity recorded by fixed-point vertical scans done by the University of Virginia Halo Streamline scanning lidar at the MLBS station. The quality control is performed according to the algorithm of Goring & Nikora (2002).

17 WIND ENERGY↗

TRACER SPARC Doppler Lidar Vertical Stare Data

These files contain 24 hour periods of data collected from the SPARC Halo Streamline XR+ Doppler lidar. While not conducting other scans, the lidar directs the beam to zenith, allowing for the measurement of vertical velocity. These data were collected during the TRACER project.

54 ENVIRONMENTAL SCIENCES↗

Characterizing warm atmospheric boundary layer over land by combining Raman and Doppler lidar measurements

PBL plays a critical role in the atmosphere by transferring heat, moisture, and momentum. The warm PBL has a distinct diurnal cycle including daytime convective mixing layer (ML) and nighttime residual layer developments. Thus, for PBL characterization and process study, simultaneous determinations of PBL height (PBLH) and ML height (MLH) are necessary. Here, new approaches are developed to provide reliable PBLH and MLH to characterize warm PBL evolution. The approaches use Raman lidar (RL) water vapor mixing ratio (WVMR) and Doppler lidar (DL) vertical velocity measurements at the Southern Great Plains (SGP) atmospheric observatory, which was established by the Atmospheric Radiation Measurement (ARM) user facility. Compared with widely used lidar aerosol measurements for PBLH, WVMR is a better trace for PBL vertical mixing. For PBLH, the approach classifies PBL water vapor structures into a few general patterns, then uses a slope method and dynamic threshold method to determine PBLH. For MLH, wavelet analysis is used to re-construct 2-D variance from DL vertical wind velocity measurements according to the turbulence eddy size to minimize the impacts of gravity wave and eddy size on variance calculations; then, a dynamic threshold method is used to determine MLH. Remotely-sensed PBLHs and MLHs are compared with radiosonde measurements based on the Richardson number method. Good agreements between them confirm that the proposed new algorithms are reliable for PBLH and MLH characterization. The algorithms are applied to warm seasons’ RL and ML measurements at the SGP site for five years to study warm season PBL structure and processes. The weekly composited diurnal evolutions of PBLHs and MLHs in warm climate were provided to illustrate diurnal and seasonal PBL evolutions. This reliable PBLH and MLH dataset will be valuable for PBL process study, model evolution, and PBL parameterization improvement.

54 ENVIRONMENTAL SCIENCES↗

ARM SGP PBLH and MLH datasets from Raman lidar and Doppler lidar

The planetary boundary layer (PBL) plays a critical role in the atmosphere by transferring heat, moisture, and momentum. The warm PBL has a distinct diurnal cycle including the daytime convective mixing layer (ML) and nighttime residual layer developments. Thus, simultaneous determinations of PBL height (PBLH) and ML height (MLH) are necessary for studying PBL characterization and processes. Here, new approaches are developed to provide reliable PBLH and MLH estimates to characterize warm PBL evolution. The approaches use Raman lidar (RL) water vapor mixing ratio (WVMR) and Doppler lidar (DL) vertical velocity measurements at the Southern Great Plains (SGP) atmospheric observatory, which was established by the Atmospheric Radiation Measurement (ARM) User Facility. Compared to widely used lidar aerosol measurements for PBLH, WVMR is a better tracer for PBL vertical mixing. For PBLH, the approach classifies PBL water vapor structures into a few general patterns, then uses a slope method and dynamic threshold method to determine PBLH. For MLH, wavelet analysis is used to reconstruct 2D variance from DL vertical wind velocity measurements according to the turbulence eddy size to minimize the impacts of gravity wave and eddy size on variance calculations; then, a dynamic threshold method is used to determine MLH. Remotely-sensed PBLHs and MLHs are compared with radiosonde measurements based on the Richardson number method. Good agreements between them confirm that the proposed new algorithms are reliable for PBLH and MLH characterization. The algorithms are applied to warm-season RL and ML measurements at the SGP site for five years to study warm-season PBL structure and processes. The weekly composited diurnal evolutions of PBLHs and MLHs in a warm climate were provided to illustrate diurnal and seasonal PBL evolutions. This reliable data set of PBLH and MLH values will be valuable for studying PBL processes, model evolution, and PBL parameterization improvements. The MLH dataset includes the MLH in values of km above ground level. The PBLH dataset includes the PBLH in values of km above ground level, along with a flag ("situation_PBLH") to determine the state of the PBL (1 = Cloudy Condition, 2 = Stable Layer, 3 = Multi-layer WVMR structure, 4 = Well-Mixed PBL, 5 = A de-coupled layer, 6 = Other).

mixing layer height↗

dlppi2 (b1)

The Doppler lidar (DL) is an active remote-sensing instrument that provides range- and time-resolved measurements of the line-of-sight component of air velocity (i.e., radial velocity) and attenuated aerosol backscatter. The DL operates in the near-infrared and is sensitive to backscatter from atmospheric aerosol, which are assumed to be ideal tracers of atmospheric wind fields. The DL works by transmitting short pulses of infrared laser light into the atmosphere. Atmospheric aerosols scatter a small fraction of that light energy back to the transceiver, where it is collected and recorded as a time-resolved signal. From the delay between the outgoing pulse and the backscattered signal, the instrument infers the distance to the scattering volume. Coherent detection is used to measure the Doppler frequency shift of the backscatter signal. This is accomplished by mixing the backscatter signal with a reference laser beam (i.e., local oscillator) of known frequency. The onboard signal processor then determines the Doppler frequency shift from the spectrum of the mixed signal. The Doppler frequency shift and thus the radial air velocity is determined from the peak of the Doppler spectrum. The attenuated backscatter is determined from the energy content of the Doppler spectra. The DL provides accurate measurements of radial velocity in regions of the atmosphere where aerosol concentrations are high enough to ensure good signal-to-noise ratio. Thus, valid data are usually limited to the atmospheric boundary layer where aerosol is ubiquitous. Valid measurements can also be obtained in elevated aerosol layers or in optically thin clouds above the boundary layer. Most of the ARM DLs have full upper-hemispheric scanning capability, enabling 3D mapping of turbulent flows in the atmospheric boundary layer. With the scanner pointed vertically, the DL provides height- and time-resolved measurements of vertical velocity.

54 ENVIRONMENTAL SCIENCES↗

Simultaneous velocity and density measurements of fully developed Rayleigh-Taylor mixing

Here, the dynamics of molecular mixing and the energy transfer process in the Rayleigh-Taylor instability (RTI) are studied through the collection of simultaneous density-velocity measurements. These experiments provide simultaneous density-velocity field measurements, in contrast to previous point measurements. Statistically stationary experiments are performed in a “convective-type” gas tunnel facility, with density contrast achieved through the injection of helium into the bottom stream. Three experiments at Atwood number ≈ 0.1 are captured at three outer scale Reynolds numbers Re = 520, 2260, and 4050. Particle image velocimetry and laser induced fluorescence are employed simultaneously. Statistics of the density and velocity show self-similar collapse of RTI profiles at large Reynolds number Re >2000. Flat velocity profiles indicate homogeneous turbulence characteristics in the core of the mixing region. Significant anisotropy develops in the flow, with horizontal velocity fluctuations being only 60% of the vertical velocity fluctuations. The turbulent mass flux is found to be asymmetric about the centerline, with increased peak towards the spike. Measurements of the molecular mixing show that mixing is maximized at the core of the flow and increases with increased Reynolds number. The transport equation of density-specific-volume correlation b shows that it is mostly produced in the core of the mixing region, and that the spatial evolution of its profile is the result of transport by bulk motion of the bubble and spike. Energy transfer from gravitational potential energy to turbulent kinetic energy and viscous dissipation is observed to occur in the experiment with a ratio of dissipated energy to potential energy released of 38%. The analysis of the turbulent kinetic energy transport equation budget reveals that production is the primary mechanism towards the growth of turbulent kinetic energy in the core of the flow, and is asymmetrically slightly skewed towards the spike. However, it is through the transport that the strong advection at the edges of the mixing region is maintained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Analysis of aerosol cloud interactions with a consistent signal of meteorology and other influencing parameters

Quantifying the impact of aerosols on cloud micro/macro physical properties and estimating the signature of Aerosol Cloud Interactions (ACI) is one of the challenging tasks in atmospheric sciences. The Moderate Resolution Imaging Spectroradiometer and the European Centre for Medium-Range Weather Forecasts ERA-5 reanalysis data are employed to systematically study the ACI over the monsoon region in Pakistan. Based on the monsoon occurrence and rainfall intensity, the whole region is divided into three sub-regions labeled as highly intensive (R1), moderately intensive (R2) and weak (R3) monsoon region. The results indicate that the monthly mean Aerosol Optical Depth (AOD) peaks in the summer monsoon months (Jun, Jul, Aug, Sep). Here, the well-known Twomey effect whereby the Cloud Droplet Radius (CDR) decreases with increasing AOD holds only over R3; the opposite effects (Anti-Twomey effect) are found over R1 and R2, all passing the test of statistical significance (p<0.05). The multi-year AOD is found to be positively correlated with Cloud Liquid Water Path (CLWP) and Cloud Optical Depth (COD) over R1 and R2, suggesting that thicker clouds containing more water droplets are formed in polluted atmosphere. Over R3, decreases in CLWP and COD are found with increasing AOD only when AOD is less than~0.325. The analysis of ACI over R1 and R2 during the winter months shows similar but stronger responses of CDR, CLWP and COD to the variation in AOD. The weaker responses during the summer monsoon season may attributed to the occurrence of high level cloud and unstable atmospheric condition. Further investigation of the influences of Relative Humidity and pressure vertical velocity on the CDR-AOD relationships shows that although the magnitude of the CDR-AOD correlations change with meteorological conditions, the sign of correlations remain unchanged with meteorological conditions.

54 ENVIRONMENTAL SCIENCES↗

Interpreting the Diurnal Cycle of Clouds and Precipitation in the ARM GoAmazon Observations: Shallow to Deep Convection Transition

The Green Ocean Amazon (GoAmazon) 2014/5 field campaign data are used to study the diurnal cycle of clouds and precipitation. Through a careful classification of days with shallow cumulus, congestus and deep convection, we investigate the major differences among locally generated convection regimes and the most important environmental factors governing the shallow-to-deep convection transition. On shallow cumulus days, a greater sensible heat flux drives deeper boundary layer growth, which entrains drier free-tropospheric air and lowers the relative humidity, thus leading to a significantly higher cloud base than those on days with deeper convection. Congestus and deep convection regimes exhibit distinct cloud top height distributions with noticeable differences in the vertical wind shear in the mid-troposphere, suggesting an important role of wind shear in limiting the vertical extent of convection. On deep convection days, with preexisting nocturnal convection or cold-pools from external disturbances, the timing of peak surface precipitation (12:00–13:00 LST) tends to be in-phase with the diurnal variation in surface fluxes. However, it takes longer for local deep convection to develop without these disturbances. A plume model with thermodynamic and dynamical constraints is developed to explore the relative importance of various convection-controlling factors. Initial cloud-base vertical velocity and buoyancy are important in helping parcels ascend to the level of free convection (LFC). After parcels reach the LFC, entrainment of environmental air and lower free troposphere humidity become crucial in determining cloud top. Entrainment rate differentiates among convection regimes, which may be tied to the cloud size distribution at cloud base.

54 ENVIRONMENTAL SCIENCES↗