Search NASA⌕ Search

SEARCH · Search NASA

Results for “Turbulent flows”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

COMPUTATION FLUID DYNAMICS ANALYSIS FOR GENERIC SMALL MODULAR REACTOR CONTAINMENT SEPARATE EFFECTS TEST

It is desirable for fourth-generation Small Modular Reactors to be passively cooled in standard and accident operations. Passive Containment Cooling Systems can reject heat from the containment structure, without using pumps or blowers. The targeted design containment structure is a large, domed, stainless steel, cylindrical vessel. In a postulated Design Basis Accident, steam will flash inside containment. Steam condensation occurs on the inner containment wall and transfers heat through the steel containment into a large body of water known as the annular reservoir (AR) surrounding the vessel serving as the ultimate heat sink. Natural circulation drives the flow in the AR and heat will be released to the environment by evaporation of water. Unique containment geometry requires a separate effects test (SET) facility for the verification and validation of the computer code and evaluation model development and assessment for reactor licensing efforts. In this study, STAR-CCM+, a computational fluid dynamics (CFD) code was used to inform the decision-making process on the design of the SET. The CFD simulation modeled, a two-phase turbulent flow with fluid film development and heat transfer for different containment geometries. The Reactor Excursion and Leak Analysis Program will also be used in a code-to-code verification against the CFD results.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2022

Wind farm flow control: prospects and challenges

Abstract. Wind farm control has been a topic of research for more than two decades. It has been identified as a core component of grand challenges in wind energy science to support accelerated wind energy deployment and to transition to a clean and sustainable energy system for the 21st century. The prospect of collective control of wind turbines in an array, to increase energy extraction, reduce structural loads, improve the balance of systems, reduce operation and maintenance costs, etc. has inspired many researchers over the years to propose innovative ideas and solutions. However, practical demonstration and commercialization of some of the more advanced concepts has been limited by a wide range of challenges, which include the complex physics of turbulent flows in wind farms and the atmosphere, uncertainties related to predicting structural load and failure statistics, and the highly multi-disciplinary nature of the overall design optimization problem, among others. In the current work, we aim at providing a comprehensive overview of the state of the art and outstanding challenges, thus identifying the key research areas that could further enable commercial uptake and success of wind farm control solutions. To this end, we have structured the discussion on challenges and opportunities into four main areas: (1) insight in control flow physics, (2) algorithms and AI, (3) validation and industry implementation, and (4) integrating control with system design (co-design).

Meyers, Johan (ORCID:0000000228284397)↗
DOE OSTI2021

ExaWind: Predictive Wind Energy Simulations

This presentation describes the ExaWind project and the team's progress in creating a suite of performance-portable codes designed for predictive simulations of wind farms on next-generation exascale-class supercomputers. Such simulations will require the resolution of scales spanning many orders of magnitude, from blade boundary layers to wind farm flow structures. In the U.S., the first exascale systems will be GPU accelerated, and different GPU manufacturers have been chosen for the different systems. At the heart of the ExaWind software is a hybrid-solver approach based on the codes Nalu-Wind and AMR-Wind, which are computational fluid dynamics solvers for the incompressible Navier-Stokes equations. Nalu-Wind is an unstructured-grid code used to resolve wind turbine geometry and blade boundary layers, whereas AMR-Wind is a structured-grid background solver for atmospheric turbulent flow and turbine wake propagation. The models are coupled with overset meshes and global linear systems are approximated through a loose-coupling algorithm. Results will include validation-quality high-fidelity simulations and strong/weak scaling results from the Summit supercomputer.

49 EE - Wind and Water Power Program - Wind (EE-4W↗
DOE OSTI2022

Long-Term Static and Dynamic Corrosion Stability of Nonwetting Surfaces

Superhydrophobic surfaces (SHSs) and lubricant-infused surfaces (LISs) are two classes of nonwetting surfaces that have drawn attention due to their advanced functional properties including corrosion inhibition. Yet there is a conspicuous lack of corrosion study of SHSs and LISs with respect to their fabrication and material parameters, especially at high temperatures and under dynamic flow conditions over long durations, which is sought to be addressed in this article. Considering copper SHSs and LISs, a full factorial combinatorial study of two facile texturing processes, electrodeposition and etching, two different functionalization agents, stearic acid and mercaptan, and two types of infused lubricants, Krytox 104 and DOWSIL 510, is presented, encompassing over 650 measurements on 90 tested surfaces. All fabricated surfaces demonstrated water repellency with a contact angle above 150° and a sliding angle below 7°. For the first time, the study examines high-temperature corrosion stability and long-term corrosion durability of the nonwetting surfaces in both static fluid and dynamic turbulent flow conditions over a period of 30 days. Here, LISs and SHSs are shown to provide excellent corrosion inhibition over all tested corrosion conditions, with negligible presence of corrosion species on the surfaces and no deterioration of the texturing. The surfaces are also shown to rejuvenate easily to the initial wettability and corrosion resistance values. This study provides valuable insights into the selection of materials and processing parameters for the fabrication of nonwetting surfaces for the application of interest.

Copper↗
DOE OSTI2020

Laminar-to-turbulent flame transition and cycle-to-cycle variations in large eddy simulation of spark-ignition engines

This paper investigates the effect of laminar-to-turbulent flame transition modeling on the prediction of cycle-to-cycle variations (CCVs) in large eddy simulation (LES) of spark-ignition (SI) engines. A laminar-to-turbulent flame transition model that describes the non-equilibrium sub-filter flame speed evolution during an early stage of flame kernel growth is developed. In the present model, the flame transition is characterized by the flame kernel size at which the flame transition ends, defined here as the flame transition scale. The proposed model captures the effects that variations in a turbulent flow field have on the evolution of early-stage burning rates, through variations in the flame transition scale. The proposed flame transition model is combined with the front propagation formulation (FPF) method and a spark-ignition model to predict CCVs in a gasoline direct injection SI engine. It is found that multi-cycle LES with the proposed flame transition model reproduces experimentally-observed CCVs satisfactorily. When the transition model is not considered or when variations in the transition process are neglected, CCVs are significantly under-predicted for the case considered here. These results indicate the importance of modeling the laminar-to-turbulent flame transition and the effect of turbulence on the transition process, when predicting CCVs, under certain engine conditions. The LES results are also used to analyze sources for variations in the flame transition. It is found, for the present engine case, that the most important source is the cycle-to-cycle variation in the turbulence dissipation rate, which is used to measure the strength of turbulence in the proposed model, near a spark plug. The large-scale velocity field and the variations of the laminar flame speed due to the mixture composition and thermal stratification are also found to be important factors to contribute to the variations in the flame transition.

42 ENGINEERING↗
DOE OSTI2022

A Hybrid Bulk Algorithm to Predict Turbulent Fluxes over Dry and Wet Bare Soils

Measurements made in the Columbia River basin (Oregon) in an area of irregular terrain during the second Wind Forecast Improvement Project (WFIP2) field campaign are used to develop an optimized hybrid bulk algorithm to predict the surface turbulent fluxes from readily measured or modeled quantities over dry and wet bare or lightly vegetated soil surfaces. The hybrid (synthetic) algorithm combines (i) an aerodynamic method for turbulent flow, which is based on the transfer coefficients (drag coefficient and Stanton number), roughness lengths, and Monin–Obukhov similarity; and (ii) a modified Priestley–Taylor (P-T) algorithm with physically based ecophysiological constraints, which is essentially based on the surface energy budget (SEB) equation. Soil heat flux in the latter case was estimated from measurements of soil temperature and soil moisture. In the framework of the hybrid algorithm, bulk estimates of the momentum flux and the sensible heat flux are derived from a traditional aerodynamic approach, whereas the latent heat flux (or moisture flux) is evaluated from a modified P-T model. Direct measurements of the surface fluxes (turbulent and radiative) and other ancillary atmospheric/soil parameters made during WFIP2 for different soil conditions (dry and wet) are used to optimize and tune the hybrid bulk algorithm. The bulk flux estimates are validated against the measured eddy-covariance fluxes. We also discuss the SEB closure over dry and wet surfaces at various time scales based on the modeled and measured fluxes. Although this bulk flux algorithm is optimized for the data collected during the WFIP2, a hybrid approach can be used for similar flux-tower sites and field campaigns.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

The effect of dynamic near-wake modulation on utility-scale wind turbine wake development

High-resolution field-scale experiments using flow visualization with natural snowfall and high-fidelity large eddy simulations are combined to investigate the effect of dynamic turbine operation and atmospheric conditions on wind turbine wake mixing and recovery in the wake of a 2.5 MW wind turbine. Instantaneous near-wake expansion and deflection in response to changes in blade pitch and wind direction, termed dynamic wake modulation, is quantified using both techniques, demonstrating excellent agreement. The simulations are used to extend these results by calculating the energy flux into the wake 7 rotor diameters downstream, showing that dynamic turbine-atmospheric interactions enhance mixing in the far-wake. This finding is exhibited under both uniform and turbulent inflow conditions. Under turbulent flow, a synergistic relationship is also observed between dynamic wake modulation and wake meandering, as wake recovery can be further accelerated when the two phenomena occur together. The results of this study have implications for the development of more realistic far-wake models that include the significant impact of dynamic wake modulation on wake mixing and development. Additionally, the findings from the current study can be used to develop advanced control algorithms to speed up wake breakdown and recovery, further improving wind farm efficiency.

17 WIND ENERGY↗
DOE OSTI2020

Energy Transfer from Large to Small Scales in Turbulence by Multiscale Nonlinear Strain and Vorticity Interactions

An intrinsic feature of turbulent flows is an enhanced rate of mixing and kinetic energy dissipation due to the rapid generation of small-scale motions from large-scale excitation. The transfer of kinetic energy from large to small scales is commonly attributed to the stretching of vorticity by the strain rate, but strain self-amplification also plays a role. Previous treatments of this connection are phenomenological or inexact, or cannot distinguish the contribution of vorticity stretching from that of strain self-amplification. In this Letter, an exact relationship is derived which quantitatively establishes how intuitive multiscale mechanisms such as vorticity stretching and strain self-amplification together actuate the interscale transfer of energy in turbulence. Further, numerical evidence verifies this result and uses it to demonstrate that the contribution of strain self-amplification to energy transfer is higher than that of vorticity stretching, but not overwhelmingly so.

42 ENGINEERING↗
DOE OSTI2025

Gas-Gap Calorimeter Sizing and Rating Tools for Heat Pipe Experimentation

Effective gas-gap calorimeter sizing and rating tools are required to design calorimeters that enable well-defined operating conditions, high cooling powers, and accurate calorimetry. The present report describes two tools utilizing Microsoft Excel and MathWorks MATLAB that enable the sizing and rating of gas-gap calorimeters with a He-Ar binary gas mixture in the gap. The Excel tools are two easy-to-use spreadsheets for calculating heat pipe operating temperature or He-Ar mole fractions given the required inputs. The MATLAB tool includes similar models with more robust solution algorithms along with sub-options for full-vacuum, static gas-gap, and forced circulation in the gas-gap. The accuracy of the developed tools were demonstrated by comparing with existing work in literature. The MATLAB tool was used for a parametric study to determine the gas-gap thickness, coolant gap thickness, coolant flow rate, and coolant inlet temperatures for the testing of a 3/4 in outer diameter heat pipe up to powers of 10 kW and temperatures of 1,000°C. The effects of heat pipe and calorimeter inner shell surface emissivities were investigated, considering the inability to control or accurately measure surface emissivities in most applications. It was found that controlling the He-Ar mole fractions provides a wide operating range for gas-gap thicknesses around ~ 0.042--0.090 in (~ 1.07--2.29 mm) for a surface emissivity range of 0.4--0.8, since conduction heat transfer is a significant fraction of the heat transfer across the gas-gap. In addition, it was found that turbulent flow in the coolant gap is needed at high input powers to prevent shell temperatures from approaching the boiling temperature of the water coolant. The parametric study resulted in choosing stainless steel tubes with an outer diameter and thickness of 1 x 0.065 in as the inner shell, and a 1-1/2 x 0.156 in as the outer shell of the calorimeter. Overall, this report presents the necessary information for the sizing and rating of gas-gap calorimeters with He-Ar mixtures for the testing of high-temperature heat pipes (~ 500--1,000°C).

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Filtering, averaging, and scale dependency in homogeneous variable density turbulence

We investigate relationships between statistics obtained from filtering and from ensemble or Reynolds-averaging turbulence flow fields as a function of length scale. Generalized central moments in the filtering approach are expressed as inner products of generalized fluctuating quantities, q ' ( ξ , x ) = q ( ξ ) - q ¯ ( x ) , representing fluctuations of a field q ( ξ ) , at any point ξ, with respect to its filtered value at x. For positive-definite filter kernels, these expressions provide a scale-resolving framework, with statistics and realizability conditions at any length scale. In the small-scale limit, scale-resolving statistics become zero. In the large-scale limit, scale-resolving statistics and realizability conditions are the same as in the Reynolds-averaged description. Using direct numerical simulations (DNS) of homogeneous variable density turbulence, we diagnose Reynolds stresses, T i j , resolved kinetic energy, kr, turbulent mass-flux velocity, a i , and density-specific volume covariance, b, defined in the scale-resolving framework. These variables, and terms in their governing equations, vary smoothly between zero and their Reynolds-averaged definitions at the small and large scale limits, respectively. At intermediate scales, the governing equations exhibit interactions between terms that are not active in the Reynolds-averaged limit. For example, in the Reynolds-averaged limit, b follows a decaying process driven by a destruction term; at intermediate length scales, it is a balance between production, redistribution, destruction, and transport, where b grows as the density spectrum develops and then decays when mixing becomes strong enough. This work supports the notion of a generalized, length-scale adaptive model that converges to DNS at high resolutions and to Reynolds-averaged statistics at coarse resolutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Xcompact3D: An open-source framework for solving turbulence problems on a Cartesian mesh

Xcompact3D is a Fortran 90–95 open-source framework designed for fast and accurate simulations of turbulent flows, targeting CPU-based supercomputers. It is an evolution of the flow solver Incompact3D which was initially designed in France in the mid-90’s for serial processors to solve the incompressible Navier–Stokes equations. Incompact3D was then ported to parallel High Performance Computing (HPC) systems in the early 2010’s. Very recently the capabilities of Incompact3D have been extended so that it can now tackle more flow regimes (from incompressible flows to compressible flows at low Mach numbers), resulting in the design of a new user-friendly framework called Xcompact3D. The present manuscript presents an overview of Xcompact3D with a particular focus on its functionalities, its ready-to-run simulations and a few case studies to demonstrate its impact.

17 WIND ENERGY↗
DOE OSTI2021

Adversarial sampling of unknown and high-dimensional conditional distributions

Many engineering problems require the prediction of realization-to-realization variability or a refined description of modeled quantities. In that case, it is necessary to sample elements from unknown high-dimensional spaces with possibly millions of degrees of freedom. While there exist methods able to sample elements from probability density functions (PDF) with known shapes, several approximations need to be made when the distribution is unknown. In this paper the sampling method, as well as the inference of the underlying distribution, are both handled with a data-driven method known as generative adversarial networks (GAN), which trains two competing neural networks to produce a network that can effectively generate samples from the training set distribution. In practice, it is often necessary to draw samples from conditional distributions. When the conditional variables are continuous, only one (if any) data point corresponding to a particular value of a conditioning variable may be available, which is not sufficient to estimate the conditional distribution. This work handles this problem using an a priori estimation of the conditional moments of a PDF. Herein, two approaches, stochastic estimation, and an external neural network are compared for computing these moments; however, any preferred method can be used. The algorithm is demonstrated in the case of the deconvolution of a filtered turbulent flow field. It is shown that all the versions of the proposed algorithm effectively sample the target conditional distribution with minimal impact on the quality of the samples compared to state-of-the-art methods. Additionally, the procedure can be used as a metric for the diversity of samples generated by a conditional GAN (cGAN) conditioned with continuous variables.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A numerical investigation of coal particle modeling for the inlet section of a pressurized oxy-combustion burner

Concerns over climate change have led to numerous efforts to develop low-carbon technologies, and pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emissions in the power industry. Designing an effective burner for such technologies is vital for the development of these and other coal combustion technologies. Because of the high pressure, the volume fraction of coal particles at the fuel inlet of the pressurized oxy-combustor burner is close to or even higher than the maximum limit for commercial CFD codes (e.g., ANSYS FLUENT). At these high particle volume fractions, the interactions among particles, fluid flow and wall need to be re-evaluated. The present computational work is a first step in a systematic analysis of the influence of various parameters, like method of particle release, release location, and particle size, in a pilot-scale POC combustor, developed at Washington University in St. Louis (WUSTL). In the POC process, pulverized coal is burned under elevated pressure in an O2-CO2 environment. Specifically, a 15-bar, 100 kWth, POC combustor is modeled employing Ansys FLUENT commercial sodtware using the Reynolds-Averaged Navier-Stokes (RANS) approach. It is revealed that for this pilot-scale, pressurized burner, velocity profiles in the near-wall region exhibit some anomalies. In order to investigate the influence of particle loading in the near-wall region, particle release location was investigated. The numerical simulations incorporate the coupling between the turbulent flow and the particles. A sensitivity investigation of particle release location found that by tuning the release location, the velocity profile can be consistent with the pure gas flow velocity profile. More importantly, the particle releasing location also affects flame stability. Particle size was also found to have a significant impact on particle trajectory, flame stability and temperature. Finally, Large Eddy Simulations (LES) were performed and compared with the results from two-dimensional RANS.

Li, Lei↗
DOE OSTI2024

Validation of new and existing methods for time-domain simulations of turbulence and loads

We seek to obtain a second-by-second match between the simulated and measured structural loads of a utility-scale wind turbine. To obtain the one-to-one load simulations, we start with the furthest upstream component of the modeling chain: the turbulent inflow. We consider new and existing methods to generate constrained-turbulence flow fields. The new method is based on large-eddy simulations (LES) and machine learning (ML). The existing methods include Kaimal-based TurbSim and the superstatistical wind field model. The inflow measurements used to constrain these simulations are obtained with a nacelle-mounted scanning lidar. We compare the flow fields for the different inflow simulation approaches and validate their associated load predictions against measurements collected in the Rotor Aero-dynamics, Aeroelastics, and Wake (RAAW) field campaign. We find that the rotor-position control developed for this study is key in enabling the time match between measurements and simulations. When this control approach is used, the load simulation performance tracks with the inflow simulation fidelity, with LES+ML yielding errors ≤ 4% for the damage-equivalent loads of flapwise bending moment, and tower fore-aft bending moments.

17 WIND ENERGY↗
DOE OSTI2021

Convection in Scaled Turbine Internal Cooling Passages With Additive Manufacturing Roughness

Abstract Additive manufacturing processes, such as direct metal laser sintering (DMLS), enable the creation of novel turbine cooling internal passages and systems. However, the DMLS method produces a significant and unique surface roughness. Previous work in scaled passages analyzed pressure losses and friction factors associated with the rough surfaces, as well as investigated the velocity profiles and turbulent flow characteristics within the passage. In this study, the heat transfer characteristics of scaled additively manufactured surfaces were measured using infrared (IR) thermography. Roughness panels were CNC machined from plates of aluminum 6061 to create near isothermal roughness elements when heated. Fluid resistance differences between the aluminum roughness panels and roughness panels constructed from ABS plastic using the same roughness patterns from McClain et al. (2020, “Flow in a Simulated Turbine Blade Cooling Channel With Spatially Varying Roughness Caused by Additive Manufacturing Orientation,” ASME Turbo Expo 2020, Turbomachinery Technical Conference and Exposition, Virtual Conference, Sept. 21–25, GT2020-16069) were investigated. Finally, the overall thermal performance enhancements and friction losses were assessed through the calculation of surface averaged “global thermal performance” ratios. The global thermal performance characterizations indicate results in-line with those found for traditional commercial roughness and slightly below traditional internal passage convection enhancement methods such as swirl chambers, dimples, and ribs. The passages investigated in this study do not include compressibility effects or the long-wavelength artifacts and channel geometric deviations observed by Wildgoose et al. (2020, “Impact of Additive Manufacturing on Internal Cooling Channels with Varying Diameters and Build Directions,” ASME Turbo Expo 2020, Turbomachinery Technical Conference and Exposition, Virtual Conference, Sept. 21–25, GT2020-15049). However, the results of this study indicate that, based on the roughness augmentation alone, artificial convective cooling enhancers such as turbulators or dimples may still be required for additively manufactured turbine component cooling.

Engineering↗
DOE OSTI2021

Experimental investigation of elemental and isotopic evaporation processes by laser heating in an aerodynamic levitation furnace

In this report we carried out evaporation experiments on a B-type calcium–aluminium-rich inclusion (CAI) melt in a gas-mixing aerodynamic levitation laser furnace, at 1873 K and an oxygen partial pressure of 10 -9.1 atm, for durations ranging from 60 to 600 s. Evaporation of SiO 2 and MgO follow the same trend as those observed in vacuum furnace experiments at the same temperature and starting composition, showing that their evaporation relative to one another from the melt is independent of pressure, oxygen fugacity, and hydrodynamical regime specific to the furnace. Isotopic ratios of Mg and Si in evaporation residues are used to derive fractionation factors of α 26/24 Mg vap-liq = 0.9906± 0.0004 and α 30/28 Si vap-lip = 0.9943±0.0003, which are both significantly closer to unity than those found for evaporation in a vacuum, which translates to less isotope fractionation. The residues are also less isotopically fractionated than expected for cases in which transport of the gas species away from the melt is diffusion-controlled at 1-atm. By analysing the flow regimes in our furnace, we find that advection by the levitating gas is the primary mode of mass transport away from the melt surface, as opposed to diffusion-limited transport in a vacuum or 1-atm tube furnace. A modified Hertz–Knudsen–Langmuir formulation accounts for this process, and shows that isotopic fractionation of both Si and Mg reflect a saturation factor (ratio of the pressure of the evaporating species to vapour saturation pressure) equal to 0.75. This is in perfect accord with recent measurements of Cu isotopic fractionation using a similar furnace. The fact that three elements (Mg, Si, Cu) with varying equilibrium vapour pressures, activity coefficients in the liquid, and diffusion coefficients in the gas have the same scaling behaviour to saturation pressure is a strong indication that the mechanism controlling evaporation is driven by the hydrodynamical regime imposed in the furnace. Therefore, this class of experiments can be used to constrain processes in which advection dominates over diffusion, such as (but not limited to) planetary ejecta, tektites, giant impacts, nebular condensation in a turbulent flow, or nuclear fallout material. Finally, the possibility to reach high temperatures (in excess of 3500 K) in this furnace allows it to be used to evaluate the activity coefficients of melt components in extreme conditions relevant to molten planetary interiors (i.e., magma oceans), with a specific focus on refractory elements.

58 GEOSCIENCES↗
DOE OSTI2024

Properties of molecular clumps and cores in colliding magnetized flows

ABSTRACT We simulate the formation of molecular clouds in colliding flows of warm neutral medium with the adaptive mesh refinement code flash in eight simulations with varying initial magnetic field strength, between 0.01–5 μG. We include a chemical network to treat heating and cooling and to follow the formation of molecular gas. The initial magnetic field strength influences the fragmentation of the forming cloud because it prohibits motions perpendicular to the field direction and hence impacts the formation of large-scale filamentary structures. Molecular clump and core formation occurs anyhow. We identify 3D clumps and 3D cores, which are defined as connected, CO-rich regions. Additionally, 3D cores are heavily shielded. While we do not claim those 3D objects to be directly comparable to observations, this enables us to analyse their full virial state. With increasing field strength, we find more fragments with a smaller average mass; yet the dynamics of the forming clumps and cores only weakly depends on the initial magnetic field strength. The molecular clumps are mostly unbound, probably transient objects, which are weakly confined by ram pressure or thermal pressure, indicating that they are swept up by the turbulent flow. They experience significant fluctuations in the mass flux through their surface, such that the Eulerian reference frame shows a dominant time-dependent term due to their indistinct nature. We define the cores to encompass highly shielded molecular gas. Most cores are in gravitational-kinetic equipartition and are well described by the common virial parameter $\alpha _\mathrm{vir}$, while some undergo minor dispersion by kinetic surface effects.

Weis, M. (ORCID:0000000256838860)↗
DOE OSTI2024

Machine-learned closure of URANS for stably stratified turbulence: connecting physical timescales & data hyperparameters of deep time-series models

Stably stratified turbulence (SST), a model that is representative of the turbulence found in the oceans and atmosphere, is strongly affected by fine balances between forces and becomes more anisotropic in time for decaying scenarios. Moreover, there is a limited understanding of the physical phenomena described by some of the terms in the Unsteady Reynolds-Averaged Navier–Stokes (URANS) equations—used to numerically simulate approximate solutions for such turbulent flows. Rather than attempting to model each term in URANS separately, it is attractive to explore the capability of machine learning (ML) to model groups of terms, i.e. to directly model the force balances. We develop deep time-series ML for closure modeling of the URANS equations applied to SST. We consider decaying SST which are homogeneous and stably stratified by a uniform density gradient, enabling dimensionality reduction. We consider two time-series ML models: long short-term memory and neural ordinary differential equation. Both models perform accurately and are numerically stable in a posteriori (online) tests. Furthermore, we explore the data requirements of the time-series ML models by extracting physically relevant timescales of the complex system. We find that the ratio of the timescales of the minimum information required by the ML models to accurately capture the dynamics of the SST corresponds to the Reynolds number of the flow. The current framework provides the backbone to explore the capability of such models to capture the dynamics of high-dimensional complex dynamical system like SST flows.

97 MATHEMATICS AND COMPUTING↗