SEARCH · Engineering Papers
Results for “Boundary Layer Transition”
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.
Using Satellite and ARM Observations to Evaluate Cold Air Outbreak Cloud Transitions in E3SM Global Storm‐Resolving Simulations
Abstract This study examines marine boundary layer cloud regime transition during a cold air outbreak (CAO) over the Norwegian Sea, simulated by a global storm‐resolving model (GSRM) known as the Simple Cloud‐Resolving Energy Exascale Earth System Model Atmosphere Model (SCREAM). By selecting observational references based on a combination of large‐scale conditions rather than strict time‐matched comparisons, this study finds that SCREAM qualitatively captures the CAO cloud transition, including boundary layer growth, cloud mesoscale structure, and phase partitioning. SCREAM also accurately locates the greatest ice and liquid in the mesoscale updrafts, however, underestimates supercooled liquid water in cumulus clouds. The model evaluation approach adopted by this study takes advantages of the existing computational‐expensive global simulations of GSRM and the available observations to understand model performance and can be applied to assessments of other cloud regimes in different regions. Such practice provides valuable guidance on the future effort to correct and improve biased model behaviors.
On the enhancement of boundary layer skin friction by turbulence: an angular momentum approach
Turbulence enhances the wall shear stress in boundary layers, significantly increasing the drag on streamlined bodies. Other flow features such as free stream pressure gradients and streamwise boundary layer growth also strongly influence the local skin friction. In this paper, an angular momentum integral (AMI) equation is introduced to quantify these effects by representing them as torques that alter the shape of the mean velocity profile. This approach uniquely isolates the skin friction of a Blasius boundary layer in a single term that depends only on the Reynolds number most relevant to the flow's engineering context, so that other torques are interpreted as augmentations relative to the laminar case having the same Reynolds number. The AMI equation for external flows shares this key property with the so-called FIK relation for internal flows (Fukagata et al. , Phys. Fluids , vol. 14, 2002, pp. L73–L76). Without a geometrically imposed boundary layer thickness, the length scale in the Reynolds number for the AMI equation may be chosen freely. After a brief demonstration using Falkner–Skan boundary layers, the AMI equation is applied as a diagnostic tool on four transitional and turbulent boundary layer direct numerical simulation datasets. Regions of negative wall-normal velocity are shown to play a key role in limiting the peak skin friction during the late stages of transition, and the relative strengths of terms in the AMI equation become independent of the transition mechanism a very short distance into the fully turbulent regime. The AMI equation establishes an intuitive, extensible framework for interpreting the impact of turbulence and flow control strategies on boundary layer skin friction.
Observations Pertaining to Precipitation within the Northeast Pacific Stratocumulus-to-Cumulus Transition
Three genuine stratocumulus-to-cumulus transitions sampled during the Cloud System Evolution over the Trades (CSET) campaign are documented. The focus is on Lagrangian evolution of in situ precipitation, thought to exceed radar/lidar retrieved values because of Mie scattering. Two of the three initial stratocumulus cases are pristine [cloud droplet number concentrations ($N_d$) of ~22 cm –3 ] but occupied boundary layers of different depths, while the third is polluted ($N_d$ ~ 225 cm –3 ). Hourly satellite-derived cloud fraction along Lagrangian trajectories indicate that more quickly deepening boundary layers tend to transition faster, into more intense but more occasional precipitation. These transitions begin either in the morning or late afternoon, suggesting that preceding night processes can precondition or delay the inevitable transition. The precipitation shifts toward larger drop sizes throughout the transition as the boundary layers deepen, with aerosol concentrations only diminishing in two of the three cases. Ultraclean ($N_d$ < 1 cm –3 ) cumulus clouds evolved from pristine stratocumulus cloud with unusually high precipitation rates occupying a shallow, well-mixed boundary layer. Results from a simple one-dimensional evaporation model and from radar/lidar retrievals suggest subcloud evaporation likely increases throughout the transition. This, coupled with larger drop sizes capable of lowering the latent cooling profile, facilitates the transition to more surface-driven convection. The coassociation between boundary layer depth and precipitation does not provide definitive conclusions on the isolated effect of precipitation on the pace of the transition. Differences between the initial conditions of the three examples provide opportunities for further modeling studies.
Low-Reynolds-number oscillating boundary layers on adiabatic slopes
We investigate the instabilities and transition mechanisms of Boussinesq stratified boundary layers on sloping boundaries when subjected to oscillatory body forcing parallel to the slope. We examine idealized forms of boundary layers on hydraulically smooth abyssal slopes in tranquil mid- to low-latitude regions, where low-wavenumber internal tides gently heave isopycnals up and down adiabatic slopes in the absence of mean flows, high-wavenumber internal tides, shelf breaks, resonant tide–bathymetry interactions (critical slopes) and other phenomena associated with turbulence ‘hot spots’. In non-rotating low-Reynolds-number flow, increased stratification on the downslope phase has a relaminarizing effect, while on the upslope phase we find transition-to-turbulence pathways arise from shear production triggered by gravitational instabilities. When rotation is significant (low slope Burger numbers) we find that boundary layer turbulence is sustained throughout the oscillation period, resembling stratified Stokes–Ekman layer turbulence. Simulation results suggest that oscillating boundary layers on smooth slopes at low Reynolds number ($Re\leqslant 840$), unity Prandtl number and slope Burger numbers greater than unity do not cause significant irreversible turbulent buoyancy flux (mixing), and that flat-bottom dissipation rate models derived from the tide amplitude are accurate within an order of magnitude.
Diagnosing Predicted Boundary Layer Evolution Errors during the Morning Transition using Observations and the Common Community Physics Package Single Column Model (CCPP SCM)
Explore the source record for details and available documents.
Observed Covariations in Boundary Layer and Cumulus Cloud Layer Processes
In this paper we examine variations in boundary-layer processes spanning the shallow to deep cumulus transition. This is accomplished by differentiating boundary layer properties on the basis of convective outcomes, ranging from shallow to deep, as observed at the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site in Oklahoma, USA. Doppler lidar, radar, and radiosonde data are combined to determine statistical differences in boundary layer and cloud layer properties using a large sample (236) of days with a range of convective outcomes: shallow, congestus, and deep convection. In these analyses, the radar characterizes diurnal cloud depth, the lidar quantifies updraft and downdraft properties in the subcloud layer, and daily radiosonde data provides the convective inhibition (CIN). Combined, these data are used to test the hypothesis that deep convection occurs when the strength of the boundary layer turbulence (i.e., TKE) exceeds the strength of the energy barrier (i.e., CIN) at the top of the CBL. Results show that days with deep convective clouds have significantly lower vertical velocity variance and weaker updrafts within the subcloud layer. However, CIN values are also found to be significantly lower on deep convective days, allowing for these weaker updrafts to penetrate the energy barrier and reach the level of free convection (LFC). In contrast, shallow convective outcomes occur when the updrafts are strong in an absolute sense, but are weak when compared to the strength of the energy barrier. These findings support the use of the CIN/TKE framework in parameterizing convection in coarse resolution models.
Evaluation of methods to determine the surface mixing layer height of the atmospheric boundary layer in the central Arctic during polar night and transition to polar day in cloudless and cloudy conditions
This study evaluates methods to derive the surface mixing layer (SML) height of the Arctic atmospheric boundary layer (ABL) using in situ measurements inside the Arctic ABL during winter and the transition period to spring. An instrumental payload carried by a tethered balloon was used for the measurements between December 2019 and May 2020 during the year-long Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition. Vertically highly resolved (centimeter scale) in situ profile measurements of mean and turbulent parameters were obtained, reaching from the sea ice to several hundred meters above ground. Two typical conditions of the Arctic ABL over sea ice were identified: cloudless situations with a shallow surface-based inversion and cloudy conditions with an elevated inversion. Both conditions are associated with significantly different SML heights whose determination as accurately as possible is of great importance for many applications. We used the measured turbulence profile data to define a reference of the SML height. With this reference, a more precise critical bulk Richardson number of 0.12 was derived, which allows an extension of the SML height determination to regular radiosoundings. Furthermore, we have tested the applicability of the Monin–Obukhov similarity theory to derive SML heights based on measured turbulent surface fluxes. The application of the different approaches and their advantages and disadvantages are discussed.
Assessing Physical Relationships Between Atmospheric State, Fluxes, and Boundary Layer Stability at McMurdo Station, Antarctica
Observations at McMurdo Station, Antarctica from 24 November 2015 through 3 January 2017 were used to characterize the physical relationships between boundary layer stability and atmospheric state and fluxes. The basis of this analysis was self-organizing maps (SOMs), a neural network algorithm, used to identify the range of potential temperature profiles present in the twice-daily radiosonde data during the ARM (Atmospheric Radiation Measurement) West Antarctic Radiation Experiment (AWARE) campaign. The SOM identified profiles ranging from strongly stable to weakly stable regimes over the lowest 500 m of the atmosphere. It was found that in the winter (MJJA), moderate and strongly stable regimes occur most frequently (61%), while weakly stable regimes dominate in the summer (DJ, 83.4%). The mechanisms responsible for the dominance of different stability regimes in each season were analyzed to determine why these regimes occur with varying frequency throughout the year. This analysis found that wind speed variations and radiative cooling are responsible for the stability observed in the winter, radiative warming, as well as weaker wind speeds, are responsible for summer weak stability, and stability variations in the transition seasons (FMA, SON) are characterized by a change in sign of net radiation with increasing stability, as wind speed changes little across stability regimes. Low-level jets were observed to occur about 50% of the time below areas of enhanced stability aloft and were observed most frequently in the transition seasons. The boundary layer depth, as determined by the Bulk Richardson number, was found to decrease with increasing stability.
Characterizing Land-Atmosphere Interactions During the Afternoon-to-Evening Transition Using ARM SGP Observations
The land-atmosphere (LA) system consists of the soil, vegetation, and lowest layer of the atmosphere. The energy and water transport processes in this system contribute to the structure and the daily cycle of the atmospheric boundary layer (ABL), which is the layer between the land surface and the free troposphere that controls the exchange of momentum, energy, and water. The state of the ABL is the result of complex two-way interactions between local processes such as surface fluxes and entrainment at the ABL top as well as large-scale and mesoscale processes such as advection and subsidence.
Boundary Layer Controls on the Shallow-to-Deep Cumulus Transition (Final Technical Report)
This project advanced our understanding of the processes governing cumulus cloud formation and provides an improved observational basis for validating earth system models. To be specific, the project used laser- and radar-remote sensors to examine the physical properties of updrafts that rise from earth’s surface and initiate clouds deeper in the atmosphere. These updrafts comprise “thermals” and “plumes” and occur at small spatial and temporal scales (e.g., 10s of minutes, 100s of meters). These small scales preclude explicit representation in most earth system and climate models, and thus necessitate “sub-grid-scale” parameterization of updraft processes. The innovation of this project was to directly measure the size, shape, strength and water vapor content of these updrafts with Doppler and Raman lidars, respectively, and to link these updraft properties to cloud processes using vertically pointed weather radars. The resulting data sets comprise 100s of thousands of updrafts and thousands of clouds, which far exceeds previous efforts, and thereby provides a robust statistical and physical representation of these processes. From these large datasets the project produced a sequence of scientific analyses that: (1) Elucidate how variations in the turbulent structure of the convective boundary layer control shallow cumulus convection, (2) Quantify the upward transport of water vapor to cloud base via thermals and plumes, (3) Validate large-eddy simulations of updrafts and shallow convective clouds, (4) Demonstrate a size-to-strength relationship between updraft width and updraft speed, and (5) Demonstrate how updrafts interact with the stability at the top of the convective boundary layer to modulate the depth and vigor of convective clouds. These results have been disseminated via several published journal articles, academic theses, and conference presentations. Collectively these results contribute to the Atmospheric System Research (ASR) program’s goal to “improve understanding of the key cloud, aerosol, precipitation, and radiation processes that affect the Earth’s radiative balance and hydrological cycle, particularly processes that limit the predictive ability of regional and global earth system models”.
Imaging gate-tunable Tomonaga–Luttinger liquids in 1H-MoSe 2 mirror twin boundaries
One-dimensional electron systems exhibit fundamentally different properties than higher-dimensional systems. For example, electron-electron interactions in one-dimensional electron systems have been predicted to induce Tomonaga-Luttinger liquid behaviour. Naturally occurring grain boundaries in single-layer transition metal dichalcogenides exhibit one-dimensional conducting channels that have been proposed to host Tomonaga-Luttinger liquids, but charge density wave physics has also been suggested to explain their behaviour. Clear identification of the electronic ground state of this system has been hampered by an inability to electrostatically gate such boundaries and tune their charge carrier concentration. Here we present a scanning tunnelling microscopy and spectroscopy study of gate-tunable mirror twin boundaries in single-layer 1H-MoSe 2 devices. Gating enables scanning tunnelling microscopy and spectroscopy for different mirror twin boundary electron densities, thus allowing precise characterization of electron-electron interaction effects. Visualization of the resulting mirror twin boundary electronic structure allows unambiguous identification of collective density wave excitations having two velocities, in quantitative agreement with the spin-charge separation predicted by finite-length Tomonaga-Luttinger liquid theory.
Modeling Arctic boundary layer cloud streets at grey-zone resolutions
To better understand how model resolution affects the formation of Arctic boundary layer clouds, we investigated the influence of grid spacing on simulating cloud streets that occurred near Utqiagvik (formerly Barrow), Alaska on 2 May 2013 and were observed by the Moderate Resolution Imaging Spectroradiometer (MODIS). The Weather Research and Forecasting model was used to simulate the clouds using nested domains with increasingly fine resolution ranging from a horizontal grid spacing of 27 km in the boundary-layer-parameterized mesoscale domain to a grid spacing of 0.111 km in the large-eddy-permitting domain. We investigated the model-simulated mesoscale environment, horizontal and vertical cloud structures, boundary layer stability, and cloud properties, all of which were subsequently used to interpret the observed roll-cloud case. Increasing model resolution led to a transition from a more buoyant boundary layer to a more shear-driven turbulent boundary layer. The clouds were stratiform-like in the mesoscale domain, but as the model resolution increased, roll-like structures, aligned along the wind field, appeared with ever smaller wavelengths. A stronger vertical water vapor gradient occurred above the cloud layers with decreasing grid spacing. With fixed model grid spacing at 0.333 km, changing the model configuration from a boundary layer parameterization to a large-eddy-permitting scheme produced a more shear-driven and less unstable environment, a stronger vertical water vapor gradient below the cloud layers, and the wavelengths of the rolls decreased slightly. In this study, only the large-eddy-permitting simulation with gird spacing of 0.111 km was sufficient to model the observed roll clouds.
Local correlation-based transition models for high-Reynolds-number wind-turbine airfoils
Abstract. Modern wind-turbine airfoil design requires robust performance predictions for varying thicknesses, shapes, and appropriate Reynolds numbers. The airfoils of current large offshore wind turbines operate with chord-based Reynolds numbers in the range of 3–15 million. Turbulence transition in the airfoil boundary layer is known to play an important role in the aerodynamics of these airfoils near the design operating point. While the lack of prediction of lift stall through Reynolds-averaged Navier–Stokes (RANS) computational fluid dynamics (CFD) is well known, airfoil design using CFD requires the accurate prediction of the glide ratio (L/D) in the linear portion of the lift polar. The prediction of the drag bucket and the glide ratio is greatly affected by the choice of the transition model in RANS CFD of airfoils. We present the performance of two existing local correlation-based transition models – one-equation model (γ− SA) and two-equation model (γ-Reθt‾- SA) coupled with the Spalart–Allmaras (SA) RANS turbulence model – for offshore wind-turbine airfoils operating at a high Reynolds number. We compare the predictions of the two transition models with available experimental and CFD data in the literature in the Reynolds number range of 3–15 million including the AVATAR project measurements of the DU00-W-212 airfoil. Both transition models predict a larger L/D compared to fully turbulent results at all Reynolds numbers. The two models exhibit similar behavior at Reynolds numbers around 3 million. However, at higher Reynolds numbers, the one-equation model fails to predict the natural transition behavior due to early transition onset. The two-equation transition model predicts the aerodynamic coefficients for airfoils of various thickness at higher Reynolds numbers up to 15 million more accurately compared to the one-equation model. As a result, the two-equation model predictions are more comparable to the predictions from eN transition model. However, a limitation of this model is observed at very high Reynolds numbers of around 12–15 million where the predictions are very sensitive to the inflow turbulent intensity. The combination of the two-equation transition model coupled with the Spalart–Allmaras (SA) RANS turbulence model is a good method for performance prediction of modern wind-turbine airfoils using CFD.
High Free-Tropospheric Aitken-Mode Aerosol Concentrations Buffer Cloud Droplet Concentrations in Large-Eddy Simulations of Precipitating Stratocumulus
A new Aitken mode aerosol microphysics scheme is developed for a large eddy simulation model in order to better investigate cloud-aerosol interactions in the marine boundary layer and to study the Aitken buffering hypothesis of McCoy et al. (2021), https://doi.org/10.1029/2020jd033529. This scheme extends the single-mode two-moment prognostic aerosol scheme of Berner et al. (2013), https://doi.org/10.5194/acp-13-12549-2013. Seven prognostic variables represent accumulation and Aitken log-normal aerosol modes in air and droplets as well as 3 gas species. Scavenging of interstitial and other unactivated aerosol by cloud and rain drops are treated using the scheme described in Berner et al. (2013), https://doi.org/10.5194/acp-13-12549-2013. The scheme includes coagulation of unactivated aerosol and a simple chemistry model with gas phase H 2 SO 4 , SO 2 , and DMS as prognostic variables to capture basic influences of sulfur chemistry on the model aerosols. Nucleation of H 2 SO 4 aerosol particles from gas-phase H 2 SO 4 is neglected. A deep, precipitating stratocumulus case (VAMOS Ocean Cloud Atmosphere Land Study RF06) is used to test the new scheme. The presence of the Aitken mode aerosol increases the cloud droplet concentration through activation of the larger Aitken particles and delays the creation of an ultraclean, strongly precipitating cumulus state. Scavenging of unactivated accumulation and Aitken particles by cloud and precipitation droplets accelerates the collapse. Increasing either the above-inversion Aitken concentration or the surface Aitken flux increases the Aitken population in the boundary layer and prevents the transition to an ultraclean state.
Radiative heating rate profiles over the southeast Atlantic Ocean during the 2016 and 2017 biomass burning seasons
Marine boundary layer clouds, including the transition from stratocumulus to cumulus, are poorly represented in numerical weather prediction and general circulation models. Further uncertainties in the cloud structure arise in the presence of biomass burning carbonaceous aerosol, as is the case over the southeast Atlantic Ocean, where biomass burning aerosol is transported from the African continent. As the aerosol plume progresses across the southeast Atlantic Ocean, radiative heating within the aerosol layer has the potential to alter the thermodynamic environment and therefore the cloud structure; however, limited work has been done to quantify this along the trajectory of the aerosol plume in the region. The deployment of the first Atmospheric Radiation Measurement (ARM) Mobile Facility (AMF1) in support of the Layered Atlantic Smoke Interactions with Clouds field campaign provided a unique opportunity to collect observations of cloud and aerosol properties during two consecutive biomass burning seasons during July through October of 2016 and 2017 over Ascension Island (7.96° S, 14.35° W). Using observed profiles of temperature, humidity, and clouds from the field campaign alongside aerosol optical properties from Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), as input for the Rapid Radiation Transfer Model (RRTM), profiles of the radiative heating rate due to aerosols and clouds were computed. Radiative heating is also assessed across the southeast Atlantic Ocean using an ensemble of back trajectories from the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model. Idealized experiments using the RRTM with and without aerosols and a range of values for the single-scattering albedo (SSA) demonstrate that shortwave (SW) heating within the aerosol layer above Ascension Island can locally range between 2 and 8 K d -1 depending on the aerosol optical properties, though impacts of the aerosol can be felt elsewhere in the atmospheric column. When considered under clear conditions, the aerosol has a cooling effect at the TOA, and based on the observed cloud properties at Ascension Island, the cloud albedo is not large enough to overcome this. Shortwave radiative heating due to biomass burning aerosol is not balanced by additional longwave (LW) cooling, and the net radiative impact results in a stabilization of the lower troposphere. However, these results are extremely sensitive to the single-scattering albedo assumptions in models.
Demonstrating Hierarchical System Development With the Common Community Physics Package Single‐Column Model: A Case Study Over the Southern Great Plains
This study demonstrates a specific application of the hierarchical system development (HSD) approach to investigate, analyze, and attribute model issues within the Unified Forecast System (UFS), with a focus on process isolation. By evaluating a non‐precipitating, shallow cumulus case at the Atmospheric Radiation Measurement Southern Great Plains site in the UFS global forecast against the observation, the investigation identifies a warmer and deeper daytime convective planetary boundary layer (PBL) and misrepresented nocturnal PBL transition. Hypothesis testing, which employs the Common Community Physics Package (CCPP) single‐column model (SCM) and uses the same physics as the UFS global model, confirms that these issues are attributed to the model physics and initialization. Specifically, misrepresented PBL processes are linked to problematic surface condition and a lack of cloud formation, which may stem from deficiencies in PBL and cloud microphysics parameterizations and their interactions. The UFS initial condition contributes to an earlier, excessively collapsed daytime convective boundary layer and a lack of decoupling between the stable boundary layer and residual layer late in the afternoon. This work introduces an avenue for the community to engage with the application of HSD, along with the CCPP and CCPP SCM, to understand the interplay of model physics, disentangle the roles of model components, as well as facilitate model and forecast improvement.
Numerical Investigation of Wall-Cooling Effect on Aero-Optical Distortions for Hypersonic Boundary Layer
Compressible wall-modeled large-eddy simulations of Mach 8 turbulent boundary-layer flows over a flat plate were carried out for the conditions of the hypersonic wind tunnel at Sandia National Laboratories. The simulations provide new insight into the effect of wall cooling on the aero-optical path distortions for hypersonic turbulent boundary-layer flows. Four different wall-to-recovery temperature ratios, 0.3, 0.48, 0.71, and 0.89, are considered. Despite the much lower grid resolution, the mean velocity, temperature, and resolved Reynolds stress profiles from the simulation for a temperature ratio of 0.48 are in good agreement with those from a reference direct numerical simulation. The normalized root-mean-square optical path difference obtained from the present simulations is compared with that from reference direct numerical simulations, Sandia experiments, as well as predictions obtained with a semi-analytical model by Notre Dame University. Here the present analysis focuses on the effect of wall cooling on the wall-normal density correlations, on key underlying assumptions of the aforementioned model such as the strong Reynolds analogy, and on the elevation angle effect on the optical path difference. Wall cooling is found to increase the velocity fluctuations and decrease the density fluctuations, resulting in an overall reduction of the normalized optical path distortion. Compared to the simulations, the basic strong Reynolds analogy overpredicts the temperature fluctuations for cooled walls. Also different from the strong Reynolds analogy, the velocity and temperature fluctuations are not perfectly anticorrelated. Finally, as the wall temperature is raised, the density correlation length, away from the wall but inside the boundary layer, increases significantly for beam paths tilted in the downstream direction.