SEARCH · Engineering Papers
Results for “BOUNDARY LAYER TRANSITION - EFFECT OF TURBULENCE”
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.
Assessment and Improvement of the SST-Gamma Transition Model in Nalu-Wind
We conduct laminar–turbulent boundary-layer transition simulations using a local correlation-based transition model for two-dimensional incompressible flow and present enhancements to improve the accuracy of transition predictions. Menter’s Galilean-invariant 𝛾 transition model is implemented in the incompressible, unstructured-grid flow solver Nalu-Wind and is validated against experimental data and results from NASA’s flow solvers. The test cases of the AIAA Transition Prediction and Modeling Workshop are investigated, namely, the T3A/T3B flat plates and the NLF(1)-0416 and S809 airfoils. Based on the results, best practices for transition simulations, particularly for an unstructured-grid flow solver, are identified. Additional airfoil simulations are conducted for two wind turbine airfoils, S822 at Reynolds numbers of 𝑂(10 5 ) and DU00-W-212 at Reynolds numbers of 𝑂(10 7 ), to assess the model at low and high Reynolds numbers. Furthermore, through this work, we propose several approaches to enhance transition simulations, including 1) enforcing positivity of the implicit operator for the source terms of the transition model, 2) employing a constant turbulence intensity in stationary external flow simulations, and 3) recommending meshing for unstructured-grid flow solvers. Finally, we provide detailed documentation of the validation and data for the canonical cases to the transition modeling community.
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.
A Model for Turbulence Spectra in the Equilibrium Range of the Stable Atmospheric Boundary Layer
Stratification can cause turbulence spectra to deviate from Kolmogorov's isotropic -5/3 power law scaling in the universal equilibrium range at high Reynolds numbers. However, a consensus has not been reached with regard to the exact shape of the spectra. Here we propose a shape of the turbulent kinetic energy and temperature spectra in horizontal wavenumber for the equilibrium range that consists of three regimes at small Froude number: the buoyancy subrange, a transition region, and the isotropic inertial subrange through dimensional analysis and substantial revision of previous theoretical approximation. These spectral regimes are confirmed by various observations in the atmospheric boundary layer. The representation of the transition region in direct numerical simulations will require large-scale separation between the Dougherty-Ozmidov scale and the Kolmogorov scale for strongly stratified turbulence at high Reynolds numbers, which is still challenging computationally. In addition, we suggest that the failure of Monin-Obukhov similarity theory in the very stable atmospheric boundary layer is due to the fact that it does not consider the buoyancy scale that characterizes the transition region.
Zooming in: SCREAM at 100 m using regional refinement over the San Francisco Bay Area
Pushing global climate models to large-eddy simulation (LES) scales over complex terrain has remained a major challenge. This study presents the first known implementation of a global model – SCREAM (Simple Cloud-Resolving E3SM Atmosphere Model) – at 100 m horizontal resolution using a regionally refined mesh (RRM) over the San Francisco Bay Area. Two hindcast simulations were conducted to test performance under both strong synoptic forcing and weak, boundary-layer-driven conditions. We demonstrate that SCREAM can stably run at LES scales while realistically capturing topography, surface heterogeneity, and coastal processes. The 100 m SCREAM-RRM substantially improves near-surface wind speed, temperature, humidity, and pressure biases compared to the baseline 3.25 km simulation, and better reproduces fine-scale wind oscillations and boundary-layer structures. These advances leverage SCREAM's scale-aware SHOC turbulence parameterization, which transitions smoothly across scales without tuning. Performance tests show that while CPU-only simulations remain costly, GPU acceleration with SCREAMv1 on NERSC's Perlmutter system enables two-day hindcasts to complete in under two wall-clock days. Our results open the door to LES-scale studies of orographic flows, boundary-layer turbulence, and coastal clouds within a fully comprehensive global modeling framework.
Classical 1/3 scaling of convection holds up to Ra = 10 15
Significance The heat transport law in turbulent convection remains central to current research in the field. Our present knowledge of the heat transport law for R a > 1 0 12 is inconclusive, where the Rayleigh number R a is a measure of the strength of convection. Massively parallel simulations of the three-dimensional convection have progressed to R a = 1 0 15 in slender cells. We resolve velocity gradients inside thin boundary layers and show that the turbulent heat transport continues to follow the classical 1 / 3 scaling law with no transition to the so-called “ultimate” state that is variously argued to have 1 / 2 scaling. Our work suggests that the boundary layers remain marginally stable and continue to act as the bottleneck for global heat transport.
Coupling of Laminar-Turbulent Transition with RANS Computational Fluid Dynamics
This project combines several new concepts to create a boundary layer transition prediction capability that is suitable for analyzing modern hypersonic flight vehicles. The first new concept is the use of ''optimization'' methods to detect the hydrodynamic instabilities that cause boundary layer transition; the use of this method removes the need for many limiting assumptions of other methods and enables quantification of the interactions between boundary layer instabilities and the flow field imperfections that generate them. The second new concept is the execution of transition analysis within a conventional hypersonics CFD code, using the same mesh and numerical schemes for the transition analysis and the laminar flow simulation. This feature enables rapid execution of transition analysis with less user oversight required and no interpolation steps needed.
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.
Assessing CLUBB PDF Closure Assumptions for a Continental Shallow-to-Deep Convective Transition Case Over Multiple Spatial Scales
Assumed-PDF (probability density function) higher-order turbulence closures (APHOCs) are now widely used for parameterizing boundary layer turbulence and shallow convection in Earth system models (ESMs). A better understanding of the resolution-dependent behavior of APHOCs is essential for improving the performance of next-generation ESMs with intended horizontal resolutions finer than 10 km. In this study, we evaluate the PDF family of Analytic double-Gaussian 1 implemented in Cloud Layers Unified By Binormals (CLUBB) over a range of spatial scales (Dx) from 2 to 100 km. A 120-km-wide large eddy simulation (LES) for a continental convection case during 2016 Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field campaign serves as benchmark to evaluate the PDF closure using an off-line approach. We find during the shallow convection period, the CLUBB PDF closure tends to produce positive biases of cloud properties and liquid water flux near cloud base for all scales of analysis. It produces negative biases for these variables near cloud top that are more severe for Dx larger than 25 km. Results show that replacing the CLUBB-parameterized moisture and temperature skewnesses with LES-derived ones can fix most of the biases if clipping of input moments is allowed to prevent the occurrence of unrealizable solutions. Overall, the performance of the PDF closure is better for smaller Dx = 2–5 km than for larger Dx = 50–100 km; for a given grid spacing, it is better when the convective clouds become deeper in the late afternoon. Likely causes for the resolution dependence and implications for improving the PDF closure are discussed.
Scrape-off layer stability, turbulence and transport in MAST-U (Final Report)
U.S. Department of Energy (DOE) supported work at Lodestar Research Corporation under grant DE-SC0019270 is summarized in this final report. The research conducted under this grant explored the properties of the edge plasma in spherical torus magnetic fusion devices. Harnessing the energy that is released in fusion reactions would provide a safe and abundant source of power to meet the growing energy needs of the world population. The edge and scrape-off layer (SOL) region of fusion plasmas is the region where the hot core plasma transitions to a cooler surrounding boundary plasma. The properties of this boundary plasma determine how interactions with material surfaces will impact both the overall device performance and the lifetime of crucial components such as the walls of the device chamber and other hardware. The research carried out studied filamentary plasma instabilities and turbulent transport in the edge and SOL of the MAST U spherical torus.
Pressure-Gradient-Based RANS Model for Predicting Separation in Transitional and Turbulent Flows
Predicting flow separation poses a significant challenge for RANS models, particularly in transitional flows over airfoils. We propose a novel improvement to RANS models to predict incipient separation in both transitional and fully turbulent flows. Our approach modifies the eddy viscosity model in regions indicated by a pressure-gradient criterion that accounts for intermittency - determining whether the boundary layer is laminar or turbulent. This model demonstrates robust generalization across various airfoil shapes and Reynolds numbers. Applied to the NREL Phase VI wind turbine rotor, our model shows improved aerodynamic performance predictions compared to the baseline RANS model.
Determination of confinement regime boundaries via separatrix parameters on Alcator C-Mod based on a model for interchange-drift-Alfvén turbulence
The separatrix operational space (SepOS) model (Eich and Manz 2021 Nucl. Fusion 61 086017) is shown to predict the L–H transition, the L-mode density limit, and the ideal magnetohydrodynamic ballooning limit in terms of separatrix parameters for a wide range of Alcator C-Mod plasmas. The model is tested using Thomson scattering measurements across a wide range of operating conditions on C-Mod, spanning $\overline{n}_{\mathrm{e}}$ = 0.3–5.5 $\times 10^{20}\,$m −3 , $B_{\mathrm{t}} = 2.5$–8.0 T, and $B_{\mathrm{p}}$= 0.1–1.2 T. An empirical regression for the electron pressure gradient scale length, $\lambda_{{p}_{\mathrm{e}}}$, against a turbulence control parameter, $\alpha_{\mathrm{t}}$, and the poloidal fluid gyroradius, $\rho_{\mathrm{s,p}}$, is constructed for H-modes and found to require positive exponents for both regression parameters, indicating turbulence widening of near-scrape-off layer widths at high $\alpha_{\mathrm{t}}$ and an inverse scaling with $B_{\mathrm{p}}$, consistent with results on ASDEX Upgrade. The SepOS model is also tested in the unfavorable drift direction and found to apply well to all three boundaries, including the L–H transition as long as a correction to the Reynolds energy transfer term, $\alpha_\mathrm{RS} \lt 1$ is applied. I-modes typically exist in the unfavorable drift direction for values of $\alpha_{\mathrm{t}} \lesssim 0.3$. Finally, an experiment studying the transition between the Type-I ELMy and EDA H-mode is analyzed using the same framework. It is found that a recently identified boundary $\alpha_{\mathrm{t}} = 0.55$ at excludes most EDA H-modes but that the balance of wavenumbers responsible for the L-mode density limit, namely $k_\mathrm{EM} = k_\mathrm{RBM}$, may better describe the transition on C-Mod. The ensemble of boundaries validated and explored is then applied to project regime access and limit avoidance for the SPARC primary reference discharge parameters.
Investigating the Impacts of Daytime Boundary Layer Clouds on Surface Energy Fluxes and Boundary Layer Structure During CHEESEHEAD19
Studies of land-atmosphere interactions under a clear sky and low cumulus cloud conditions are common from long-term observatories like at the southern great plains. How well the relationships and responses of surface radiative and turbulent heat fluxes determined from these investigations hold for more heterogeneous surfaces in other climate regimes, however, is uncertain. In this study, detailed observations of the surface energy budget and daytime boundary layer properties are analyzed using measurements from the Chequamegon Heterogenous Ecosystem Energy-Balance Study Enabled by a High-Density Extensive Array of Detectors 2019 (CHEESEHEAD19) field campaign, July-October 2019, across a heterogeneous forested landscape of northern Wisconsin. A cloud regime framework is employed to classify consecutive periods of clear skies from lower atmosphere stratiform and cumulus clouds. A seasonal transition from low cumulus to low stratiform periods occurred, together with a diurnal pattern in cloudy or clear sky period dominance. Radiative forcing was highly dependent on sky conditions, leading to changes in the redistribution efficiency of radiative energy by the surface turbulent heat fluxes. During CHEESEHEAD19, small Bowen ratios dominated with daytime latent heat fluxes three times as large as sensible heat fluxes for all sky conditions studied; the forested region, therefore, falls within an energy-limited regime. The depth of the daytime mixed layer depended upon the sky condition and thermodynamic setting; deeper mixed layers occurred during periods of low cumulus and not clear skies. Profiles of vertical velocity were found to have enhanced variance under low cumulus compared to clear sky periods, suggesting potential for cloud feedbacks on boundary layer structure and surface energy fluxes.
Numerical Simulations of Orographic Convection across Multiple Gray Zones
Idealized simulations are used to determine the sensitivity of moist orographic convection to horizontal grid spacing Δ h . In simulated mechanically (MECH) and thermally (THERM) forced convection over an isolated ridge, Δ h is varied systematically over both the deep-convection (Δ h ~ 10 –1 km) and turbulence (Δ h ~ 1 km–100 m) gray zones. To aid physical interpretation, a new parcel-based bulk entrainment/detrainment diagnosis for horizontally heterogeneous flows is developed. Within the deep-convection gray zone, the Δ h sensitivity is dominated by differences in parameterized versus explicit convection; the former initiates convection too far upstream of the ridge (MECH) and too early in the diurnal heating cycle (THERM). These errors stem in part from a large underprediction of parameterized entrainment and detrainment. Within the turbulence gray zone, sensitivities to Δh arise from the representation of both subcloud- and cloud-layer turbulence. As Δ h is decreased, MECH exhibits stronger cloud-layer entrainment to enhance the convective mass flux M co , while THERM exhibits stronger detrainment to suppress Mco and delay convection initiation. The latter is reinforced by increased subcloud turbulence at smaller Δ h , which leads to drying and diffusion of the central updraft responsible for initiating moist convection. Numerical convergence to a robust solution occurs only in THERM, which develops a fully turbulent flow with a resolved inertial subrange (for Δ h ≤ 250 m). In MECH, by contrast, turbulent transition occurs within the orographic cloud, the details of which depend on both physical location and Δ h .
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”.
University of Nebraska unmanned aerial system (UAS) profiling during the LAPSE-RATE field campaign
Abstract. This paper describes the data collected by the University of Nebraska-Lincoln (UNL) as part of the field deployments during the Lower Atmospheric Process Studies at Elevation – a Remotely-piloted Aircraft Team Experiment (LAPSE-RATE) flight campaign in July 2018. The UNL deployed two multirotor unmanned aerial systems (UASs) at multiple sites in the San Luis Valley (Colorado, USA) for data collection to support three science missions: convection initiation, boundary layer transition, and cold air drainage flow. We conducted 172 flights resulting in over 21 h of cumulative flight time. Our novel design for the sensor housing onboard the UAS was employed in these flights to meet the aspiration and shielding requirements of the temperature and humidity sensors and to separate them from the mixed turbulent airflow from the propellers. Data presented in this paper include timestamped temperature and humidity data collected from the sensors, along with the three-dimensional position and velocity of the UAS. Data are quality-controlled and time-synchronized using a zero-order-hold interpolation without additional post-processing. The full dataset is also made available for download at https://doi.org/10.5281/zenodo.4306086 (Islam et al., 2020).
Turbulence Structures in the Very Stable Boundary Layer Under the Influence of Wind Profile Distortion
Abstract In very stable boundary layers (VSBL), a “cocktail” of submeso motions routinely result in elevated mean wind speed maxima above the ground, acting as a new source of turbulence generation. This new source of turbulent kinetic energy enhances turbulent mixing and causes mean wind profile distortion (WPD). As a results, this transient distortion in the wind profile adjusts the classical log‐law. Addressing how WPD‐induced turbulence regulates flow structures, turbulent fluxes, and transitions in stability regimes across layers remains a challenge. Eddy covariance data measured at four levels on a 62‐m tower are employed to address these questions. It is shown that the WPD initiates large turbulent eddies that penetrate downward, leading to enhanced vertical mixing and comparable turbulent transport efficiencies across layers. As a consequence, turbulence intensity and fluxes are increased. As the WPD is intensified, turbulent fluxes and turbulent flux transport caused by large eddies are also enhanced, leading to a transition from very stable to weakly stable regimes. Due to the influence of WPD‐induced large eddies, the large‐eddy turbulent Prandtl number does not deviate appreciably from unity and the partitioning between turbulent kinetic and potential energies is linearly related to the gradient Richardson number.
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.