SEARCH · Search NASA

Results for “Mixing layer”

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 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2025

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

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

mixing layer height↗

DOE OSTI2022

Impact of Rainfall on Tropical Cyclone-Induced Sea Surface Cooling

Tropical cyclones (TCs) are often accompanied by strong winds and torrential rains. While the winds associated with TCs tend to enhance mixing in the upper-ocean, the freshwater input from rain can stratify the water column and limit mixing. However, the extent to which the stabilizing effect of rainfall can compete with wind-induced mixing, and to what degree it modulates TC-induced sea surface cooling, remains unknown. Here we show, using a suite of observations, that heavy rains under weak TCs can significantly reduce the magnitude of cold wakes induced by them. Additionally, when compared to storms with low rain rates, the ocean surface under TCs with high rain rates freshens significantly and cools less. High-resolution climate model simulations and idealized experiments with an ocean mixed layer model support these results and reveal that oceanic mixing processes are primarily responsible for the reduced cooling under TCs, with a lesser role for surface fluxes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Impact of surface hydrophilicity on the ordering and transport properties of bicontinuous microemulsions

Microemulsions (MEs) have many industrial applications, where recent developments have shown that MEs can be utilized for electrochemical applications, including potentially in redox flow batteries. However, understanding the structure and dynamics of these systems, including at a surface, is needed to direct and rationally control their electrochemical behavior. While bulk solution measurements have provided insight into their structure, their assembly at an interface also impacts the electron (to the electrode) and ion (across the surfactant) charge transfer processes in the system. To address this shortcoming, neutron reflectivity experiments and molecular simulations have been completed that document the near surface structure of a series of deuterated water (D 2 O)/polysorbate-20/toluene MEs on hydrophilic and amphiphilic surfaces. These results show that the microemulsions form complex layered structures near a hard electrode surface, where most layers are not purely one component. Decreasing the D 2 O concentration in the ME increases the number of and purity of the layers established on the solid surface. These lamellar-type layers transition from the surface to the bulk microemulsion as a series of mixed layers (i.e., containing oil, water, and surfactant) that are consistent with perforated lamellae. Additionally, these mixed lamellae appear to become more perforated with oil and water pathways on an amphiphilic surface. Furthermore, the purity and thickness of these layers will influence the accessibility of an electrode by redox active species, as well as ion transport required to satisfy the electroneutrality condition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Evaluating turbulence models at high energy densities

The RESHOCK campaign at Lawrence Livermore National Laboratory has been working on increasing our understanding of the evolution of turbulent, unstable plasma interfaces which are applicable in various HED applications including ICF implosions. A common approach to model such interfaces involves the use of Renolds-Averaged-Navier-Stokes (RANS) models whose parameters have been constrained by theory and experiments. These parameter sets are not unique and the applicability of one set of tuned model parameters to systems with differing shock strengths or material densities has not been well tested. Our new NIF experiments, along with related studies of turbulent shear flows, provide data and model validation for plasma interfaces in the HED regime. Our experiments measure mixing-layer width at an unstable interface in the plasma regime and are specifically designed to challenge mix models of the Reynolds-Averaged-Navier-Stokes (RANS) type where turbulence is assumed to be fully developed, i.e. the flow exhibits a broad spectrum of length scales without memory of the initial condition.4 We utilize precise control of the initial interface conditions (densities on both sides of the interface, an initial rippled perturbation pattern consisting of equally weighted wavelengths between 10 and 20 microns, and initial material compositions) along with a repeatable drive history

36 MATERIALS SCIENCE↗

DOE OSTI2022

Aerosols consistently suppress the convective boundary layer development

Aerosols with different vertical distribution and various optical properties induce diverse heating rates and thereby affecting convective boundary layer (CBL) development. Our results showed consistent CBL-suppression of aerosols during daytime with numerical experiments, in which aerosols were specified at different heights with synthesized single scattering albedo from 64 studies and asymmetry factor from 20 studies globally. Absorbing aerosols concentrated below but close to the CBL top had the strongest suppression effect on CBL development relative to that concentrated near surface or above CBL. Aerosol cooling effect by attenuating incident solar radiation and surface heat flux exceeded its warming effect by reheating the atmosphere layer with absorbed shortwave radiation, and eventually declined net heating rate, which inhibited CBL development, lowered mixed-layer potential temperature and stabilized atmospheric stratification. Stove effect of absorbing aerosols (CBL enhancement) under a zero background aerosol extinction coefficient is negligible for dominant dome effect (CBL suppression) which consistently suppresses CBL development regardless of aerosol vertical height and background aerosol extinction coefficient. Our study also highlighted the importance of specifying background aerosol extinction coefficient in numerical experiments for accurate assessment of aerosol radiative forcing and CBL-aerosol interactions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Hydrodynamic Modeling of Stratification and Mixing in a Shallow, Tropical Floodplain Lake

Abstract Floodplain lakes are widespread and ecologically important throughout tropical river systems, however data are rare that describe how temporal variations in hydrological, meteorological and optical conditions moderate stratification and mixing in these shallow lakes. Using time series measurements of meteorology and water‐column temperatures from 17 several day campaigns spanning two hydrological years in a representative Amazon floodplain lake, we calculated surface energy fluxes and thermal stratification, and applied and evaluated a 3‐dimensional hydrodynamic model. The model successfully simulated diel cycles in thermal structure characterized by buoyancy frequency, depth of the actively mixing layer, and other terms associated with the surface energy budget. Diurnal heating with strong stratification and nocturnal mixing were common; despite considerable heat loss at night, the strong stratification during the day meant that mixing only infrequently extended to the bottom at night. Simulations indicated that the diurnal thermocline up and downwelled creating lake‐wide differences in near‐surface temperatures and mixing depths. Infrequent full mixing creates conditions conducive to anoxia in these shallow lakes given their warm temperatures.

Environmental Sciences & Ecology↗

DOE OSTI2022

Atomic-Scale Insights into the Lateral and Vertical Epitaxial Growth in Two-Dimensional Pd 2 Se 3 –MoS 2 Heterostructures

Two-dimensional (2D) materials form heterostructures in both the lateral and vertical directions when two different materials are interfaced, but with totally different bonding mechanisms of covalent in-plane to van der Waal’s layered interactions. Understanding how the competition between lateral and vertical forces influences the epitaxial growth is important for future materials development of complex mixed layered heterostructures. In this study, we use atomic-resolution annular dark-field scanning transmission electron microscopy to study the detailed atomic arrangements at mixed 2D heterostructure interfaces composed of two semiconductors with distinctly different crystal symmetry and elemental composition, Pd 2 Se 3 :MoS 2 , in order to understand the role of different chemical bonds on the resultant epitaxy. Pd 2 Se 3 is grown off the step edge in bilayer MoS 2 , and the vertical and lateral epitaxial relationships of the Pd 2 Se 3 –MoS 2 heterostructures are investigated. We find that the similarity of geometry at the interface with one metal (Pd or Mo) atoms bonded with two chalcogens (S or Se) are the crucial factors to make the atomically stitched lateral junction of 2D heterostructures. In addition, the vertical van der Waal interactions that are normally dominant in layered materials can be overcome by in-plane forces if the interfacial atomic stitching is high in quality and low in defect density. This knowledge should help guide the approaches for improving the epitaxy in mixed 2D heterostructures and seamless stitching of in-plane 2D heterostructures with various complex monolayer structures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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

Boundary Layer Structures Over the Northwest Atlantic Derived From Airborne High Spectral Resolution Lidar and Dropsonde Measurements During the ACTIVATE Campaign

The Planetary Boundary Layer height (PBLH) is essential for studying PBL and ocean-atmosphere interactions. Marine PBL is usually defined to include a mixed layer (ML) and a capping inversion layer. The ML height (MLH) estimated from the measurements of aerosol backscatter by a lidar was usually compared with PBLH determined from radiosondes/dropsondes in the past, as the PBLH is usually similar to MLH in nature. However, PBLH can be much greater than MLH for decoupled PBL. Here, in this study, we evaluate the retrieved MLH from an airborne lidar (HSRL-2) by utilizing 506 co-located dropsondes during the ACTIVATE field campaign over the Northwest Atlantic from 2020 to 2022. First, we define and determine the MLH and PBLH from the temperature and humidity profiles of each dropsonde, and find that the MLH values from HSRL-2 and dropsondes agree well with each other, with a coefficient of determination of 0.66 and median difference of 18 m. In contrast, the HSRL-2 MLH data do not correspond to dropsonde-derived PBLH, with a median difference of -47 m. Therefore, we modify the current operational and automated HSRL-2 wavelet-based algorithm for PBLH retrieval, decreasing the median difference significantly to -8 m. Further data analysis indicates that these conclusions remain the same for cases with higher or lower cloud fractions, and for decoupled PBLs. These results demonstrate the potential of using HSRL-2 aerosol backscatter data to estimate both marine MLH and PBLH and suggest that lidar-derived MLH should be compared with radiosonde/dropsonde-determined MLH (not PBLH) in general.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

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

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

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Optimizing feed modulation for coupled methane and NO x conversion over Pd-Pt/Mn 0.5 Fe 2.5 O 4 /Al 2 O 3 monolith catalyst

Here the impacts of feed modulation (frequency, amplitude) and catalyst design (composition and architecture) parameters are reported for the conversion of methane and NO x over a dual-layer Pt+Pd/Al 2 O 3 + Mn 0.5 Fe 2.5 O 4 /Al 2 O 3 monolith. CH 4 and NO x conversion data show that the dual-layer catalyst outperforms single-layer samples having the same catalyst loadings, with and without spinel. Close proximity of the PGM and MFO functions in the mixed-layer catalyst lowers the CH 4 conversion at high temperature while separating the PGM and spinel layers with an intermediate Al 2 O 3 layer does not. Methane conversion enhancement is linked to its nonmonotonic dependence on O 2 . The performance gains are tied to a transient activity spike that occurs during the lean-to-rich feed transition when water is present in the feed. The transient spike is attributed to the removal of CO and H 2 products via reactions with stored O 2 in the spinel, eliminating inhibition of methane steam reforming.

03 NATURAL GAS↗

DOE OSTI2025

Evaluation of probability density function descriptions for three-component Rayleigh–Taylor mixing

Results from simulations of a three-component Rayleigh–Taylor (RT) mixing problem are presented. These simulations are conducted in heavy–light–heavy and heavy–intermediate–light configurations, and each of these configurations are further considered in high- and low-Reynolds-number regimes. This results in RT-unstable flow with one or both interfaces initially unstable, permitting the influence of problem configuration on the statistical description of three-component RT-driven mixing to be considered. Mass fraction covariances are observed to undergo a sign change through the mixing layer in all four configurations considered. This appears to be unique to the multi-component case and represents another way in which multi-component RT mixing differs from the two-component case. Qualitative and quantitative comparisons of joint and marginal probability density function (PDF) descriptions of species concentration are made. Three-, five-, and six-parameter model PDFs are compared against simulation data to assess how accurately they describe the mixing, and it is found that three-component mixing requires at least a five-parameter model PDF to accurately describe the mixing. Notably, the marginal distributions of three-component mixing do not appear to conform to a beta distribution, representing a departure from the classical two-component RT case. In conclusion, statistical neutrality also appears to influence the optimal choice of model PDF, which is found to be a function of problem configuration.

Large-eddy simulation↗

DOE OSTI2025

Advanced Model Development for Large Eddy Simulation of Oxy-Combustion and Supercritical Carbon Dioxide Power Cycles

A joint experimental and numerical study is performed to observe the characteristics of a supercritical carbon dioxide turbulent mixing layer in the presence of strong nonlinearities in the thermodynamic and transport properties. A bespoke experimental setup is designed and employed for this purpose and provides insight into macroscopic mixing behavior. The mixing is experimentally observed using two techniques: shadowgraphy and spontaneous Raman scattering. Qualitative and quantitative intensity fields obtained via these techniques yield instantaneous and mean density data. Spanwise temperature data is also collected using analogue resistance temperature detectors. These measurements are used to quantify the level of mixed material within the field. The experimental data are supplemented by a companion high-fidelity numerical study. The numerical results are obtained through fully resolved, three-dimensional direct numerical simulation. The numerical dataset permits observation of the near-field mixing characteristics, which are difficult to measure experimentally due to the rapid dynamics and sharp thermophysical gradients in this area. Qualitative field visualizations are presented, followed by quantitative mixed material results and observations regarding thermodynamic property trends at select locations within the field. One-dimensional spectra of the turbulent kinetic energy and solenoidal dissipation are provided to observe the spectral characteristics of the flow. Reynolds stress anisotropy is analyzed graphically through anisotropy invariance maps (Lumley triangles). The mixing quantification, spectral data and anisotropy analysis of a flow at these thermodynamic conditions represent the main outcomes of the work.

20 FOSSIL-FUELED POWER PLANTS↗

DOE OSTI2024

Profiling the molecular destruction rates of temperature and humidity as well as the turbulent kinetic energy dissipation in the convective boundary layer

A simultaneous deployment of Doppler, temperature, and water-vapor lidars is able to provide profiles of molecular destruction rates and turbulent kinetic energy (TKE) dissipation in the convective boundary layer (CBL). Horizontal wind profiles and profiles of vertical wind, temperature, and moisture fluctuations are combined, and transversal temporal autocovariance functions (ACFs) are determined for deriving the dissipation and molecular destruction rates. These are fundamental loss terms in the TKE as well as the potential temperature and mixing ratio variance equations. These ACFs are fitted to their theoretical shapes and coefficients in the inertial subrange. Error bars are estimated by a propagation of noise errors. Sophisticated analyses of the ACFs are performed in order to choose the correct range of lags of the fits for fitting their theoretical shapes in the inertial subrange as well as for minimizing systematic errors due to temporal and spatial averaging and micro- and mesoscale circulations. We demonstrate that we achieve very consistent results of the derived profiles of turbulent variables regardless of whether 1 or 10 s time resolutions are used. We also show that the temporal and spatial length scales of the fluctuations in vertical wind, moisture, and potential temperature are similar with a spatial integral scale of ≈160 m at least in the mixed layer (ML). The profiles of the molecular destruction rates show a maximum in the interfacial layer (IL) and reach values of ϵ m ≃ 7 x 10 -4 g 2 kg −2 s −1 for mixing ratio and ϵ θ ≃ 1.6 x 10 -3 K 2 s −1 for potential temperature. In contrast, the maximum of the TKE dissipation is reached in the ML and amounts to ≃ 10 -2 m 2 s −3 . We also demonstrate that the vertical wind ACF coefficient k w ∝ $\overline{w'^{2}}$ and the TKE dissipation ϵ ∝ ($\overline{w'^{2}}$) 3/2 . For the molecular destruction rates, we show that ϵ m ∝ $\overline{m'^{2}}$ ($\overline{w'^{2}}$) 1/2 and ϵ θ ∝ $\overline{θ'^{2}}$ ($\overline{w'^{2}}$) 1/2 . These equations can be used for parameterizations of ϵ, ϵ m , and ϵ θ . All noise error bars are derived by error propagation and are small enough to compare the results with previous observations and large-eddy simulations. The results agree well with previous observations but show more detailed structures in the IL. Consequently, the synergy resulting from this new combination of active remote sensors enables the profiling of turbulent variables such as integral scales, variances, TKE dissipation, and the molecular destruction rates as well as deriving relationships between them. The results can be used for the parameterization of turbulent variables, TKE budget analyses, and the verification of large-eddy simulations.

dissipation rate↗

DOE OSTI2024

Heat transfer and transport property contrast effects on the compressible Rayleigh-Taylor instability

In extreme environments, the Rayleigh-Taylor instability (RTI) may occur under large variations in density and temperature and with fluid transport properties strongly dependent on temperature. Direct numerical simulations of the 3D fully compressible RTI are conducted, examining the idealized configuration of a hotter, less dense fluid pushing against a colder, denser fluid. Various temperature ratios and transport property configurations are explored to examine how heat conduction, large variations in transport properties, and sudden changes in transport properties can affect the evolution of the mixing layer. Nonuniform fluid expansion and contraction induced by heat transfer can significantly affect local density differences and overall instability growth, causing profile asymmetries about the initial interface for flow and mixing statistics. The departures from classical self-similar development of the instability along with misalignment between regions of mixing and regions of most intense turbulent activity caused by both heat transfer and transport property contrasts are examined. After sudden changes in fluid transport properties, which may occur as a result of rapid heating (e.g., in inertial confinement fusion), the flow quickly responds and begins to relax towards quasi-self-similar late-time evolution. For many dynamical quantities such as vorticity and dissipation, this late-time evolution resembles that of the configuration that already started with the final transport property magnitudes, suggesting that these quantities depend only on the transport properties and not on past flow history, provided that the density field distributions for the flows remain similar. On the other hand, the mixing evolution after the transport property change is unique, implying that both property magnitudes and previous history are impactful on the mixing. In conclusion, these simulations demonstrate how various temperature-related effects are extremely important to consider in compressible RTI flows with large temperature variations.

42 ENGINEERING↗

DOE OSTI2022

Ice-shelf ocean boundary layer dynamics from large-eddy simulations

Abstract. Small-scale turbulent flow below ice shelves is regionally isolated and difficult to measure and simulate. Yet these small-scale processes, which regulate heat and salt transfer between the ocean and ice shelves, can affect sea-level rise by altering the ability of Antarctic ice shelves to “buttress” ice flux to the ocean. In this study, we improve our understanding of turbulence below ice shelves by means of large-eddy simulations at sub-meter resolution, capturing boundary layer mixing at scales intermediate between laboratory experiments or direct numerical simulations and regional or global ocean circulation models. Our simulations feature the development of an ice-shelf ocean boundary layer through dynamic ice melting in a regime with low thermal driving, low ice-shelf basal slope, and strong shear driven by the geostrophic flow. We present a preliminary assessment of existing ice-shelf basal melt parameterizations adopted in single component or coupled ice-sheet and ocean models on the basis of a small parameter study. While the parameterized linear relationship between ice-shelf melt rate and far-field ocean temperature appears to be robust, we point out a little-considered relationship between ice-shelf basal slope and melting worthy of further study.

54 ENVIRONMENTAL SCIENCES↗

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 OSTI2022

Deep Learning for Subgrid‐Scale Turbulence Modeling in Large‐Eddy Simulations of the Convective Atmospheric Boundary Layer

Abstract In large‐eddy simulations, subgrid‐scale (SGS) processes are parameterized as a function of filtered grid‐scale variables. First‐order, algebraic SGS models are based on the eddy‐viscosity assumption, which does not always hold for turbulence. Here we apply supervised deep neural networks (DNNs) to learn SGS stresses from a set of neighboring coarse‐grained velocity from direct numerical simulations of the convective boundary layer at friction Reynolds numbers Re τ up to 1243 without invoking the eddy‐viscosity assumption. The DNN model was found to produce higher correlation between SGS stresses compared to the Smagorinsky model and the Smagorinsky‐Bardina mixed model in the surface and mixed layers and can be applied to different grid resolutions and various stability conditions ranging from near neutral to very unstable. The DNN model can capture key statistics of turbulence in a posteriori (online) tests when applied to large‐eddy simulations of the atmospheric boundary layer.

54 ENVIRONMENTAL SCIENCES↗