SEARCH · Search NASA

Results for “large eddy simulation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 703 records · Page 39

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Comparison of wind-farm control strategies under realistic offshore wind conditions: wake quantities of interest

Wind-farm control strategies aim to increase the efficiency, and therefore lower the levelized cost of energy, of wind farms. This is done by using turbine settings such as the yaw angle, blade pitch angles, or generator torque to manipulate the wake that negatively affects downstream turbines in the farm. Two inherently different wind-farm control methods have been identified in the literature: wake steering (WS) and active wake mixing (AWM). As one of two companion papers focused on understanding practical aspects of these two wind-farm control strategies using large-eddy simulation (LES), we below analyze the wake quantities of interest for a single wind turbine performing WS and AWM, while the companion article (Frederik et al., 2025) focuses on turbine quantities of interest including power and structural loads for the same computational setup and also includes two-turbine arrays with full and partial wake overlap. The simulations, which are based in the LES solver AMR-Wind, are tailored to have inflow conditions representative of measurements from a site off the East Coast of the US, including with strong veer and low turbulence. The turbine, which is modeled in OpenFAST and coupled to the LES, is the IEA 15 MW, an open-source offshore design. After presenting an overview of the wake recovery for the different wake-control cases, the analysis probes the fluid-dynamic causes for the different performance of the arrays reported in the companion article by examining control volumes around the wakes and the budget of the mean-flow kinetic energy (MKE) within these volumes. In the high veer environment considered, the MKE recovery is dominated by mean convection, and this is shown to especially benefit the WS strategy when a neighboring turbine is directly downstream: there is ≈65 % more available power for a downstream turbine than in the baseline case, and this power is gained primarily through mean convection on the left-tip and top-tip faces of the control volume. However, the case with imperfect knowledge of the exact wind direction favors the pulse-type AWM strategy, largely because of ≈9 % increased turbulent entrainment from aloft versus the baseline that could be related to an apparent resistance to skewing in the pulsed wake. The general reduced effectiveness of helix-type and other individual-pitch-based AWM strategies for inflow with high veer and low turbulence as reported in the companion paper is due, in part, to low magnitudes of phase-averaged turbulent entrainment. Two main findings of this study are thus that veer has a significant impact on the effectiveness of different wake-control strategies and that pulse-type AWM may be a useful strategy when the objective is power maximization in realistic, offshore flow environments with imperfect knowledge of the exact wake overlap position on the downstream turbine.

17 WIND ENERGY↗

NASA NTRS2014

NASA Tech Briefs, June 2014

Topics include: Real-Time Minimization of Tracking Error for Aircraft Systems; Detecting an Extreme Minority Class in Hyperspectral Data Using Machine Learning; KSC Spaceport Weather Data Archive; Visualizing Acquisition, Processing, and Network Statistics Through Database Queries; Simulating Data Flow via Multiple Secure Connections; Systems and Services for Near-Real-Time Web Access to NPP Data; CCSDS Telemetry Decoder VHDL Core; Thermal Response of a High-Power Switch to Short Pulses; Solar Panel and System Design to Reduce Heating and Optimize Corridors for Lower-Risk Planetary Aerobraking; Low-Cost, Very Large Diamond-Turned Metal Mirror; Very-High-Load-Capacity Air Bearing Spindle for Large Diamond Turning Machines; Elevated-Temperature, Highly Emissive Coating for Energy Dissipation of Large Surfaces; Catalyst for Treatment and Control of Post-Combustion Emissions; Thermally Activated Crack Healing Mechanism for Metallic Materials; Subsurface Imaging of Nanocomposites; Self-Healing Glass Sealants for Solid Oxide Fuel Cells and Electrolyzer Cells; Micromachined Thermopile Arrays with Novel Thermo - electric Materials; Low-Cost, High-Performance MMOD Shielding; Head-Mounted Display Latency Measurement Rig; Workspace-Safe Operation of a Force- or Impedance-Controlled Robot; Cryogenic Mixing Pump with No Moving Parts; Seal Design Feature for Redundancy Verification; Dexterous Humanoid Robot; Tethered Vehicle Control and Tracking System; Lunar Organic Waste Reformer; Digital Laser Frequency Stabilization via Cavity Locking Employing Low-Frequency Direct Modulation; Deep UV Discharge Lamps in Capillary Quartz Tubes with Light Output Coupled to an Optical Fiber; Speech Acquisition and Automatic Speech Recognition for Integrated Spacesuit Audio Systems, Version II; Advanced Sensor Technology for Algal Biotechnology; High-Speed Spectral Mapper; "Ascent - Commemorating Shuttle" - A NASA Film and Multimedia Project DVD; High-Pressure, Reduced-Kinetics Mechanism for N-Hexadecane Oxidation; Method of Error Floor Mitigation in Low-Density Parity-Check Codes; X-Ray Flaw Size Parameter for POD Studies; Large Eddy Simulation Composition Equations for Two-Phase Fully Multicomponent Turbulent Flows; Scheduling Targeted and Mapping Observations with State, Resource, and Timing Constraints;

Source record↗

NASA NTRSDate not supplied

High-Lift Common Research Model: RANS, HRLES and WMLES Perspectives for CLmax Prediction Using LAVA

A unified assessment of three turbulence treatments: Reynolds Averaged Navier-Stokes(RANS), Hybrid RANS/LES (HRLES) and Equilibrium Wall-Modelled Large Eddy Simulation(WMLES) is presented for the High-Lift Common Research Model (CRM-HL). For the free-air configuration, steady-state RANS simulations show very accurate drag polar predictions in the low-𝛼linear regime. However, strong grid sensitivity is reported near the maximum lift-state(𝐶𝐿max), with finer-grids showing larger errors and predicting erroneous flow topologies on the wing. Our RANS simulations show that several corrections for the Spalart-Allmaras (SA)turbulence model widely used in the community lead to more erroneous results compared to the baseline closure, without exception. Both scale-resolving methods (HRLES and WMLES)address these drawbacks and predict an outboard separation pattern on the main element that is in good agreement with the oil flow photographs taken from the QinetiQ wind tunnel experiments, when LES-appropriate grids and numerical discretizations are used. While RANS simulations with the baseline SA closure do not show any wing-root separation post𝐶𝐿max, both HRLES and WMLES show onset of corner flow separation with varying degrees of progression, along with a weak pitch break in the wing-contribution of the overall pitching moment. This post-𝐶𝐿max pitch break seen in the free-air simulations is weaker than the break observed in experiments, with a weaker break reported in WMLES for each iteration of grid-refinement. In-tunnel simulations using both SA-baseline RANS and WMLES show a much stronger post-𝐶𝐿max break with the WMLES predictions showing excellent agreement with the experiment in terms of both the flow-topology observed and the pressure-coefficients at various spanwise stations. Sensitivity to the tunnel wall boundary layer is characterized via comparisons between viscous and inviscid treatments for the tunnel walls. WMLES predictions show moderate sensitivity at the predicted inboard flow-state at 𝐶𝐿max along with the progression towards a post-𝐶𝐿max stall; however, this stalled state at 𝛼≈20◦(inside the tunnel) obtained with both tunnel wall treatments appears to be largely identical.

TTT↗

NASA NTRS2011

Noise of Embedded High Aspect Ratio Nozzles

A family of high aspect ratio nozzles were designed to provide a parametric database of canonical embedded propulsion concepts. Nozzle throat geometries with aspect ratios of 2:1, 4:1, and 8:1 were chosen, all with convergent nozzle areas. The transition from the typical round duct to the rectangular nozzle was designed very carefully to produce a flow at the nozzle exit that was uniform and free from swirl. Once the basic rectangular nozzles were designed, external features common to embedded propulsion systems were added: extended lower lip (a.k.a. bevel, aft deck), differing sidewalls, and chevrons. For the latter detailed Reynolds-averaged Navier-Stokes (RANS) computational fluid dynamics (CFD) simulations were made to predict the thrust performance and to optimize parameters such as bevel length, and chevron penetration and azimuthal curvature. Seventeen of these nozzles were fabricated at a scale providing a 2.13 inch diameter equivalent area throat." ! The seventeen nozzles were tested for far-field noise and a few data were presented here on the effect of aspect ratio, bevel length, and chevron count and penetration. The sound field of the 2:1 aspect ratio rectangular jet was very nearly axisymmetric, but the 4:1 and 8:1 were not, the noise on their minor axes being louder than the major axes. Adding bevel length increased the noise of these nozzles, especially on their minor axes, both toward the long and short sides of the beveled nozzle. Chevrons were only added to the 2:1 rectangular jet. Adding 4 chevrons per wide side produced some decrease at aft angles, but increased the high frequency noise at right angles to the jet flow. This trend increased with increasing chevron penetration. Doubling the number of chevrons while maintaining their penetration decreased these effects. Empirical models of the parametric effect of these nozzles were constructed and quantify the trends stated above." Because it is the objective of the Supersonics Project that future design work be done more by physics-based computations and less by experiments, several codes under development were evaluated against these test cases. Preliminary results show that the RANS-based code JeNo predicts the spectral directivity of the low aspect ratio jets well, but has no capability to predict the non-axisymmetry. An effort to address this limitations, used in the RANS-based code of Leib and Goldstein, overpredicted the impact of aspect ratio. The broadband shock noise code RISN, also limited to axisymmetric assumptions, did a good job of predicting the spectral directivity of underexpanded 2:1 cold jet case but was not as successful on high aspect ratio jets, particularly when they are hot. All results are preliminary because the underlying CFD has not been validated yet. An effort using a Large Eddy Simulation code by Stanford University predicted noise that agreed with experiments to within a few dB.

Bridges, James E.↗

NASA NTRS2011

Small-Scale Dissipation in Binary-Species Transitional Mixing Layers

Motivated by large eddy simulation (LES) modeling of supercritical turbulent flows, transitional states of databases obtained from direct numerical simulations (DNS) of binary-species supercritical temporal mixing layers were examined to understand the subgrid-scale dissipation, and its variation with filter size. Examination of the DSN-scale domain- averaged dissipation confirms previous findings that, out of the three modes of viscous, temperature and species-mass dissipation, the species-mass dissipation is the main contributor to the total dissipation. The results revealed that the percentage of species-mass by total dissipation is nearly invariant across species systems and initial conditions. This dominance of the species-mass dissipation is due to high-density-gradient magnitude (HDGM) regions populating the flow under the supercritical conditions of the simulations; such regions have also been observed in fully turbulent supercritical flows. The domain average being the result of both the local values and the extent of the HDGM regions, the expectations were that the response to filtering would vary with these flow characteristics. All filtering here is performed in the dissipation range of the Kolmogorov spectrum, at filter sizes from 4 to 16 times the DNS grid spacing. The small-scale (subgrid scale, SGS) dissipation was found by subtracting the filtered-field dissipation from the DNS-field dissipation. In contrast to the DNS dissipation, the SGS dissipation is not necessarily positive; negative values indicate backscatter. Backscatter was shown to be spatially widespread in all modes of dissipation and in the total dissipation (25 to 60 percent of the domain). The maximum magnitude of the negative subgrid- scale dissipation was as much as 17 percent of the maximum positive subgrid- scale dissipation, indicating that, not only is backscatter spatially widespread in these flows, but it is considerable in magnitude and cannot be ignored for the purposes of LES modeling. The Smagorinsky model, for example, is unsuited for modeling SGS fluxes in the LES because it cannot render backscatter. With increased filter size, there is only a modest decrease in the spatial extent of backscatter. The implication is that even at large LES grid spacing, the issue of backscatter and related SGS-flux modeling decisions are unavoidable. As a fraction of the total dissipation, the small-scale dissipation is between 10 and 30 percent of the total dissipation for a filter size that is four times the DNS grid spacing, with all OH cases bunched at 10 percent, and the HN cases spanning 24 30 percent. A scale similarity was found in that the domain-average proportion of each small-scale dissipation mode, with respect to the total small-scale dissipation, is very similar to equivalent results at the DNS scale. With increasing filter size, the proportion of the small-scale dissipation in the dissipation increases substantially, although not quite proportionally. When the filter size increases by four-fold, 52 percent for all OH runs, and 70 percent for HN runs, of the dissipation is contained in the subgrid-scale portion with virtually no dependence on the initial conditions of the DNS. The indications from the dissipation analysis are that modeling efforts in LES of thermodynamically supercritical flows should be focused primarily on mass-flux effects, with temperature and viscous effects being secondary. The analysis also reveals a physical justification for scale-similarity type models, although the suitability of these will need to be confirmed in a posteriori studies.

Bellan, Josette↗

NASA NTRSDate not supplied

A Wall-Modeled LES Perspective for the High Lift Common Research Model Using LAVA

A new immersed boundary Wall-Modelled Large Eddy Simulation (WMLES) formulation is developed to study high-lift aerodynamics on the NASA High-Lift Common Research Model (HL-CRM). A sequence of Cartesian Octree grids with sizes ranging from 100 Million through 2.02 Billion grid points is utilized to systematically assess grid-sensitivity and convergence for the in-tunnel (QinetiQ) configuration of the model, and remarkable agreement between the immersed boundary and the curvilinear body-aligned WMLES formulations is reported on grids with comparable resolutions. In the free-air configuration of the model, consistent predictions between the Curvilinear Overset and the Cartesian Octree formulations are reported for angles of attack up to C(L,max) at a=19.57. However, some differences in the onset of stall are seen between the two methods for a>20°: while the curvilinear WMLES experiences wing-root separation with increasing angle of attack (Topology A), the Cartesian Octree formulation shows a different flow topology characterized by boundary layer weakness on the main element, emanating from the pylon-wing attachment (Topology B). In order to obtain further insight into the two-distinct topologies, carefully designed numerical experiments to isolate effects of the model standoff and the tunnel wall-boundary layers are conducted using the immersed boundary WMLES formulation. The increased incidence angle-of-attack on the inboard portion of the wing due to the standoff is shown to be sufficient for triggering a switch from Topology B to Topology A in Cartesian WMLES. The role of the floor boundary layer is further examined in detail by identification of additional corner-flow vorticity generated by the viscous juncture flow interactions between the floor boundary layer and the standoff leading to formation of a strong coherent and persistent vortex on the belly-side of the fuselage. The intensity of this vortex is shown to increase with the thickness of the floor boundary layer. A further increase in the incidence angle of attack near the leading-edge strake caused by the presence of this belly-side vortex is quantified for two-distinct floor boundary layers. Both of the floor boundary layers considered result in the onset of large scale wing-root separation at a=21.47in non-confined (free-air) configurations.

TTT↗

NASA NTRS2015

Investigation of Unsteady Tip Clearance Flow in a Low-Speed One and Half Stage Axial Compressor with LES And PIV

The primary focus of this paper is to investigate how a rotor's unsteady tip clearance flow structure changes in a low speed one and half stage axial compressor when the rotor tip gap size is increased from 0.5 mm (0.49% of rotor tip blade chord, 2% of blade span) to 2.4 mm (2.34% chord, 4% span) at the design condition are investigated. The changes in unsteady tip clearance flow with the 0.62 % tip gap as the flow rate is reduced to near stall condition are also investigated. A Large Eddy Simulation (LES) is applied to calculate the unsteady flow field at these three flow conditions. Detailed Stereoscopic PIV (SPIV) measurements of the current flow fields were also performed at the Johns Hopkins University in a refractive index-matched test facility which renders the compressor blades and casing optically transparent. With this setup, the unsteady velocity field in the entire flow domain, including the flow inside the tip gap, can be measured. Unsteady tip clearance flow fields from LES are compared with the PIV measurements and both LES and PIV results are used to study changes in tip clearance flow structures. The current study shows that the tip clearance vortex is not a single structure as traditionally perceived. The tip clearance vortex is formed by multiple interlaced vorticities. Therefore, the tip clearance vortex is inherently unsteady. The multiple interlaced vortices never roll up to form a single structure. When phased-averaged, the tip clearance vortex appears as a single structure. When flow rate is reduced with the same tip gap, the tip clearance vortex rolls further upstream and the tip clearance vortex moves further radially inward and away from the suction side of the blade. When the tip gap size is increased at the design flow condition, the overall tip clearance vortex becomes stronger and it stays closer to the blade suction side and the vortex core extends all the way to the exit of the blade passage. Measured and calculated unsteady flow fields inside the tip gap agree fairly well. Instantaneous velocity vectors inside the tip gap from both the PIV and LES do show flow separation and reattachment at the entrance of tip gap as some earlier studies suggested. This area at the entrance of tip gap flow (the pressure side of the blade) is confined very close to the rotor tip section. With a small tip gap (0.5mm), the gap flow looks like a simple two-dimensional channel flow with larger velocity near the casing for both flow rates. A small area with a sharp velocity gradient is observed just above the rotor tip. This strong shear layer is turned radially inward when it collides with the incoming flow and forms the core structure of the tip clearance vortex. When tip gap size is increased to 2.4 mm at the design operation, the radial profile of the tip gap flow changes drastically. With the large tip gap, the gap flow looks like a two-dimensional channel flow only near the casing. Near the rotor top section, a bigger region with very large shear and reversed flow is observed.

compressor↗

NASA NTRS2015

Investigation of Unsteady Tip Clearance Flow in a Low-Speed One and Half Stage Axial Compressor with LES And PIV

unsteady↗

DOE OSTI2025

Investigation of Main Bearing Fatigue Estimate Sensitivity to Synthetic Turbulence Models Using a Novel Drivetrain Model Implemented in OpenFAST

ABSTRACT A coupled medium‐fidelity drivetrain model is developed and implemented in OpenFAST for a 10‐MW land‐based reference turbine. The implementation is verified against a fully coupled multibody wind turbine model, including a detailed drivetrain. The new model can simultaneously and accurately estimate main bearing loads and represent elastic bending of the drivetrain. It has low computational cost and is useful for early design phases, sensitivity analyses and complex systems like wind farms (where computational expense must be expended elsewhere). Here, the model is implemented for a monopile offshore wind turbine and used to investigate the sensitivity of main bearing basic rating life to different synthetic turbulence models. Large‐eddy simulations (LES) targeting stable, neutral, and unstable atmospheric conditions at below‐, near‐ and above‐rated wind speeds are used as a reference. The turbulence models recommended by the International Electrotechnical Commission, the Mann spectral tensor model, and the Kaimal spectral model with exponential coherence are fitted to the LES data. Additionally, a constrained turbulence generator, PyConTurb (short for Python Constrained Turbulence ), based on LES data, is applied in the aero‐hydro‐servo‐elastic simulations. Taking PyConTurb as the baseline, the Kaimal model significantly underestimates fatigue of the downwind main bearing, with between 10% and 40% less damage. The Mann model also underestimates the downwind main bearing fatigue by up to 30%. The upwind main bearing damage is driven by mean loads, and differences between models are less significant, although the trends are similar. Reasons for these discrepancies are investigated and attributed to differences in spatial and temporal variations among the turbulence models.

17 WIND ENERGY↗

DOE OSTI2025

Main bearing response in a waked 15-MW floating wind turbine in below-rated conditions

Increased wind turbine size raises unknowns related to structural flexibility. Moreover, moving to deeper waters, component reliability becomes more critical. This work investigates main bearing response dependence on drivetrain flexibility and wake impingement in a two-turbine wind farm. A 15-MW floating direct-drive turbine is considered. Large eddy simulations (LES) are employed to model neutral, stable and unstable atmospheric conditions at below-rated mean wind speed, while the engineering codes OpenFAST and FAST.Farm simulate turbine and wake behavior. Results indicate significant sensitivities in fatigue estimates to lateral distance between the upstream and downstream turbine. The trends are most substantial in stable conditions, where the waked downwind main bearing sees twice the fatigue damage estimates of the upstream turbine for one position and 50% for another. Main bearing fatigue sensitivity to drivetrain flexibility is minor, while properly including generator rotor inertia loads is important for the axial forces of the locating (axially fixed) bearing, especially in stable conditions.

17 WIND ENERGY↗

DOE OSTI2026

Combined effects of horizontal and vertical resolution on reliable turbulence prediction at tidal energy sites: A systematic study in the Salish Sea, WA

Predicting turbulence characteristics with coastal ocean models is essential for tidal energy converter deployment. While large eddy simulation provides a detailed turbulence representation, computational limitations restrict its use to smaller domains. We systematically evaluate whether well-configured coastal models can provide reliable turbulence prediction through progressive refinement of 3D model representation. We implemented four model configurations (Levels 1–4) using terrain-following coordinates, isolating the impacts of horizontal resolution, vertical resolution, and layer distribution. We validated all configurations against field measurements from the Salish Sea, WA. Results show that tidal current velocity predictions remain unchanged regarding model configuration, but turbulence properties are sensitive to resolution refinement. Increasing vertical resolution alone proved insufficient; even with vertical sigma-levels rising from 11 to 41, significant underprediction persisted until finer horizontal resolution better captured bathymetric variations. The Level 4 configuration, incorporating geometric sigma-levels distribution, achieved turbulence prediction skill scores exceeding 0.90. Turbulence closure comparison revealed Mellor- Yamada 2.5 outperformed k-epsilon in TKE prediction (skill scores 0.84–0.94 versus 0.72–0.81) due to better boundary layer parameterization. This study shows that well-configured coastal models effectively bridge the gap between simplified tools and costly high-fidelity modeling, offering the tidal energy industry practical and cost-effective turbulence data at commercially relevant scales.

Marine Energy↗

DOE OSTI2024

Langmuir turbulence in suspended kelp farms

This study investigates the influence of suspended kelp farms on ocean mixed layer hydrodynamics in the presence of currents and waves. We use the large eddy simulation method, where the wave effect is incorporated by solving the wave-averaged equations. Distinct Langmuir circulation patterns are generated within various suspended farm configurations, including horizontally uniform kelp blocks and spaced kelp rows. Intensified turbulence arises from the farm-generated Langmuir circulation, as opposed to the standard Langmuir turbulence observed without a farm. The creation of Langmuir circulation within the farm is attributed to two primary factors depending on farm configuration: (i) enhanced vertical shear due to kelp frond area density variability, and (ii) enhanced lateral shear due to canopy discontinuity at lateral edges of spaced rows. Both enhanced vertical and lateral shear of streamwise velocity, representing the lateral and vertical vorticity components, respectively, can be tilted into downstream vorticity to create Langmuir circulation. This vorticity tilting is driven by the Craik–Leibovich vortex force associated with the Stokes drift of surface gravity waves. In addition to the farm-generated Langmuir turbulence, canopy shear layer turbulence is created at the farm bottom edge due to drag discontinuity. The intensity of different types of turbulence depends on both kelp frond area density and the geometric configuration of the farm. The farm-generated turbulence has substantial consequences for nutrient supply and kelp growth. These findings also underscore the significance of the presence of obstacle structures in modifying ocean mixed layer characteristics.

Bo, Tong (ORCID:0000000260300561)↗

DOE OSTI2025

Unsteady aerodynamic loads on pitching aerofoils represented by Gaussian body force distributions

The actuator line model (ALM) is an approach commonly used to represent lifting and dragging devices like wings and blades in large-eddy simulations (LES). The crux of the ALM is the projection of the actuator point forces onto the LES grid by means of a Gaussian regularisation kernel. The minimum width of the kernel is constrained by the grid size; however, for most practical applications like LES of wind turbines, this value is an order of magnitude larger than the optimal value that maximises accuracy. This discrepancy motivated the development of corrections for the actuator line, which, however, neglect the effect of unsteady spanwise shed vorticity. In this work we develop a model for the impact of spanwise shed vorticity on the unsteady loading of an aerofoil modelled as a Gaussian body force distribution, where the model is applicable within the regime of unsteady attached flow. The model solution is derived both in the time and frequency domain and features an explicit dependence on the Gaussian kernel width. We verify the model with ALM-LES for both pitch steps and periodic pitching. The model solution is compared with Theodorsen theory and validated with both computational fluid dynamics using body fitted grids and experiment. It is concluded that the optimal kernel width for unsteady aerodynamics is approximately 40 % of the chord. The ALM is able to predict the magnitude of the unsteady loading up to a reduced frequency of 𝑘 ≈ 0.2.

17 WIND ENERGY↗

DOE OSTI2026

Revisiting the validity of eddy viscosity models for predicting airflow over water waves

In this study, we revisit the validity of eddy viscosity models for predicting wave-induced airflow disturbances over ocean surface waves. We first derive a turbulence curvilinear model for the phase-averaged Navier–Stokes equations, extending the work of Cao, Deng & Shen (2020 J. Fluid Mech. 901, A27), by incorporating turbulence stress terms previously neglected in the linearised viscous curvilinear model. To verify our formulation, we perform a priori tests by numerically solving the model using mean wind and turbulence stress profiles from large-eddy simulations (LES) of airflow over waves across various wave ages. Results show that including turbulence stress terms improves wave-induced airflow predictions compared with the previous viscous curvilinear model. We further show that using a standard mixing-length eddy viscosity yields inaccurate predictions at certain wave ages, as it fails to capture wave-induced turbulence, which fundamentally differs from mean shear-driven turbulence. The LES data show that accurate representations of wave-induced stresses require a complex-valued eddy viscosity. The maximum magnitude of this eddy viscosity scales as ∼𝑢 𝜏 ⁢𝜁 𝑖𝑛𝑛𝑒𝑟 , where 𝑢 𝜏 is the friction velocity and 𝜁 𝑖𝑛𝑛𝑒𝑟 is the inner-layer thickness, the height at which the eddy-turnover time matches the wave advection time scale. This scaling aligns with the prediction by Belcher & Hunt (1993 J. Fluid Mech. 251, 109–148). Overall, the findings demonstrate that traditional eddy viscosity models are inadequate for capturing wave-induced turbulence. More sophisticated turbulence models are essential for the accurate prediction of airflow disturbances and form drag in wind–wave interaction models.

16 TIDAL AND WAVE POWER↗

DOE OSTI2025

Oceanic mixing and waves in the presence of a suspended canopy

Large-eddy simulations are analysed to determine the influence of suspended canopies, such as those formed in macroalgal farms, on ocean mixed layer (OML) deepening and internal wave generation. In the absence of a canopy, we show that Langmuir turbulence, when compared with wind-driven shear turbulence, results in a deeper OML and more pronounced internal waves beneath the OML. Subsequently, we examine simulations with suspended canopies of varying densities located in the OML, in the presence of a background geostrophic current. Intensified turbulence occurs in the shear layer at the canopy’s bottom edge, arising from the interaction between the geostrophic current and canopy drag. Structures resembling Kelvin–Helmholtz (KH) instability emerge as the canopy shear layer interacts with the underlying stratification, radiating internal waves beneath the OML. Both intensified turbulence and lower-frequency motions associated with KH-type structures are critical factors in enhancing mixing. Consequently, the OML depth increases by up to a factor of two compared with cases without a canopy. Denser canopies and stronger geostrophic currents lead to more pronounced KH-type structures and internal waves, stronger turbulence and greater OML deepening. Additionally, vertical nutrient transport is enhanced as the OML deepens due to the presence of the canopy. Considering that the canopy density investigated in this study closely represents offshore macroalgal farms, these findings suggest a mechanism for passive nutrient entrainment conducive to sustainable farming. Overall, this study reveals the intricate interactions between the suspended canopy, turbulent mixing and stratification, underscoring their importance in reshaping OML characteristics.

Bo, Tong (ORCID:0000000260300561)↗

DOE OSTI2026

Shear-free, inhomogeneous turbulence in a stably stratified fluid

High-resolution large eddy simulations are conducted of locally forced, shear-free turbulence in the presence of an initially sharp density interface. The simulations are reminiscent of oscillating grid turbulence experiments used to isolate the effect of turbulent diffusion and entrainment from background shear. By simulating such a flow we avoid common challenges of the experiments such as secondary-flow contamination due to sidewall effects and the inevitable interaction of the stratifying agent and forcing region. To address the latter concern, we add a heating term (potential energy sink) to the governing equations in the forcing layer, thereby preventing a heat flux through the source region. This modification sets up a continuous stratification in the mixed layer that is often assumed to be negligible in experiments. Despite this difference, we are able to make meaningful comparisons in terms of the overall entrainment rate, which varies as a power law with a turbulent Richardson number. Two exponents, −2 and −1, are measured depending on the definition of the Richardson number and entrainment rate used. The definition leading to −1 is consistent with most experiments, and we argue it is the superior choice if one is able to measure the relevant quantities. We also verify the self-similar scaling of turbulence velocity and length scales in the homogeneous fluid and propose ‘inner’ and ‘outer’ scalings for the stratified cases based on a local Froude number. The detailed scaling results are useful for turbulence model validation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Aitken Mode Aerosols Buffer Decoupled Mid‐Latitude Boundary Layer Clouds Against Precipitation Depletion

Aerosol-cloud-precipitation interactions are a leading source of uncertainty in estimating climate sensitivity. Remote marine boundary layers where accumulation mode (~100–400 nm diameter) aerosol concentrations are relatively low are very susceptible to aerosol changes. These regions also experience heightened Aitken mode aerosol (~10–100 nm) concentrations associated with ocean biology. Aitken aerosols may significantly influence cloud properties and evolution by replenishing cloud condensation nuclei and droplet number lost through precipitation (i.e., Aitken buffering). We use a large-eddy simulation with an Aitken-mode enabled microphysics scheme to examine the role of Aitken buffering in a mid-latitude decoupled boundary layer cloud regime observed on 15 July 2017 during the Aerosol and Cloud Experiments in the Eastern North Atlantic flight campaign: cumulus rising into stratocumulus under elevated Aitken concentrations (~100–200 mg -1 ). In situ measurements are used to constrain and evaluate this case study. Our simulation accurately captures observed aerosol-cloud-precipitation interactions and reveals time-evolving processes driving regime development and evolution. Aitken activation into the accumulation mode in the cumulus layer provides a reservoir for turbulence and convection to carry accumulation aerosols into the drizzling stratocumulus layer above. Further Aitken activation occurs aloft in the stratocumulus layer. Together, these activation events buffer this cloud regime against precipitation removal, reducing cloud break-up and associated increases in heterogeneity. We examine cloud evolution sensitivity to initial aerosol conditions. With halved accumulation number, Aitken aerosols restore accumulation concentrations, maintain droplet number similar to original values, and prevent cloud break-up. Without Aitken aerosols, precipitation-driven cloud break-up occurs rapidly. In this regime, Aitken buffering sustains brighter, more homogeneous clouds for longer.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Unsteady Land-Sea Breeze Circulations in the Presence of a Synoptic Pressure Forcing

Unsteady land-sea breezes (LSBs) that result from time-varying surface temperature contrasts Δθ(t) are explored in the presence of a constant synoptic pressure forcing, M g , oriented from sea to land (α = 0°) or land to sea (α = 180°). Large eddy simulations reveal the development of four distinctive regimes, depending on the joint interaction between M g , α, and Δθ(t) in modulating the fine-scale dynamics. Time lags, computed as the shifts that maximize correlation coefficients of the velocity between the unsteady and the corresponding steady scenarios at Δθ = Δθ max , are found to be significant and to extend 2 hr longer for α = 0° compared to α = 180°. These diurnal dynamics result in nonequilibrium conditions that are significantly affected by the flow history, and that behave differently over the two patches for the different α’s. Turbulence is found to be out of equilibrium with the mean flow, and the mean itself is found to be out of equilibrium with the thermal forcing. The sea surface heat flux is consistently more sensitive than its land counterpart to the time-varying external forcing Δθ(t), and more so for synoptic forcing from land to sea (α = 180°). Hence, although the land reaches equilibrium faster, the sea patch is found to exert a stronger control on the turbulence-mean flow equilibrium response. Finally, the vertical velocity profile at the shore and shore-normal velocity transects at the first grid level are shown to encode the multiscale regimes of the LSBs evolution and can thus be used to identify these regimes using k-means clustering.

58 GEOSCIENCES↗