Search NASA⌕ Search

SEARCH · Search NASA

Results for “unsteady state”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

One-Parameter Analytical Derivation in Modified Budyko Framework for Unsteady-State Streamflow Elasticity in Humid Catchments

The changes in climate and catchment properties have altered the hydrological processes significantly at different spatiotemporal scales around the world. In particular for finer time scales, changes in water storage, which has been commonly neglected for long-term temporal scales, may play an important role on hydrological processes. Nevertheless, few studies addressed this question in modifying the Budyko framework, with most of them remaining a steady-state assumption on catchment characteristics. Here, in this work, we derive new analytical formulas of unsteady-state streamflow elasticity in a modified Budyko framework, incorporating both storage change and one specific parameter for catchment properties. We study 78 humid catchments in the USA with simulation data of daily hydrological processes from a probability-distribution-based hydrological model (the modified HyMOD). As indicated by results, the annual storage change ratio is linearly correlated with the annual aridity index, and this relationship can be used to estimate elasticity coefficients with our formulas. The estimated elasticity coefficients perform well in simulating the annual streamflow with the power-law model. For different catchments, variability of the unsteady-state elasticity is higher than that of the steady-state elasticity. Unsteady-state streamflow coefficients show significant linear correlation with catchment properties, such as the average slope, average elevation, and catchment area. This study provides a new analytical approach to investigate the interannual stability of catchments with varying climate and catchment properties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Equation-Free Coarse Control of Distributed Parameter Systems via Local Neural Operators

The control of high-dimensional distributed parameter systems (DPS) remains a challenge when explicit coarse-grained equations are unavailable. Classical equation-free (EF) approaches rely on fine-scale simulators treated as black-box timesteppers. However, repeated simulations for steady-state computation, linearization, and control design are often computationally prohibitive, or the microscopic timestepper may not even be available, leaving us with data as the only resource. We propose a data-driven alternative that uses local neural operators, trained on spatiotemporal microscopic/mesoscopic data, to obtain efficient short-time solution operators. These surrogates are employed within Krylov subspace methods to compute coarse steady and unsteady-states, while also providing Jacobian information in a matrix-free manner. Krylov-Arnoldi iterations then approximate the dominant eigenspectrum, yielding reduced models that capture the open-loop slow dynamics without explicit Jacobian assembly. Both discrete-time Linear Quadratic Regulator (dLQR) and pole-placement (PP) controllers are based on this reduced system and lifted back to the full nonlinear dynamics, thereby closing the feedback loop.

93B52, 93C20, 47N70, 65J15, 65M32, 68T07, 68T20, 6↗
DOE OSTI2025

An experimental study on gas-liquid phase fluid migration in hydrate-bearing sediments during hydrate dissociation

Natural gas hydrate production tests face problems such as severe sand blockage, poor gas-liquid phase separation, and significant land subsidence. This is because of the insufficient understanding of the complex phase transition and gas-liquid multi-phase fluid migration during hydrate dissociation. In hydrate-bearing sediment systems, hydrate phase transition couples with gas-liquid fluid migration. The phase transition causes changes in pore structure, which in turn modifies porous infiltration parameters and fluid flow capacity. Meanwhile, alterations in phase interfaces affect key parameters like surface tension and wettability. Gas-liquid fluid migration influences heat and mass transfer, thus affecting phase equilibrium and dissociation rates. To bridge the gap in describing gas-liquid fluid migration during hydrate dissociation in experiments, this research innovatively integrated an unsteady-state gas displacement by water and a quantitative hydrate dissociation process, independently developed a multi-phase seepage experimental system suitable for hydrate dissociation and determined the relationship between seepage parameters and hydrate saturation under different porosity. The results are as follows: a) Core samples with higher initial porosity show a greater recovery rate of fluid flow capacity. b) The retarding effect of multi-phase fluid has a more significant impact on the migration of the wetting phase fluid (water) than that of the non-wetting phase fluid (methane). c) During hydrate dissociation, the evolution of absolute permeability shows an “S-shaped” pattern, and the evolution of relative permeability shows a “wiring-harness” pattern. In conclusion, the findings can provide a theoretical basis for preventing geological disasters and for geotechnical engineering design during hydrate production.

58 GEOSCIENCES↗
DOE OSTI2024

Electrochemical leaching of spent LIBs: Kinetics, novel reactor, and modeling

The use of electrons as main reagent for the recovery and recycling of critical metals from spent lithium-ion batteries (LIBs) is a process electrification strategy that can be used to close the life-cycle loop of LIBs through more sustainable methods. Electrochemical leaching, a process that uses a reductant that is constantly regenerated electrochemically for the leaching of lithium-ion battery black mass (LIBBM), has shown high extraction efficiencies and sustainable scores. However, slow kinetics, reactor design challenges and lack of deeper understanding of the underlying processes are barriers to the optimization, scale-up, and market adoption of this technology. In this paper, a kinetic study and mathematical model for dissolving LIBBM is presented to better understand the underlying mechanisms aiming to reduce the processing time and make predictions for future design and scale-up. The effect of acid and electrochemically mediated reductant concentrations, LIBBM loading, and cathode/reactor designs were explored. As a result, the leaching time was reduced from 7h to under 1h at a pulp density of 73 g/L, without external heating. A novel reactor with parallel baffle electrodes (PBE) was developed, which significantly reduced the leaching time by improving convection in a stirred slurry electrochemical reactor. Dimensionless numbers were deduced from an unsteady state model, which can be used in dimensional analysis for future process design and scale-up.

25 ENERGY STORAGE↗
DOE OSTI2023

Developing Renewable Energy Applications for Water Treatment Technologies and Transmission

The need for smaller treatment facilities for more isolated or rural municipal entities and industrial water treatment facilities continues to grow. From environmental justice and energy needs in remote communities that have never been serviced by larger scale utilities to smaller onsite treatment requirements to meet discharge regulations, energy efficient treatment systems that can leverage renewable energy sources can be an important asset to owners. Climate change and population dynamics have changed the needs for distributed water treatment systems and the manner in which they function. New treatment technologies that can operate under dynamic conditions and leverage variable energy sources (e.g., renewable energy, RE) can provide substantial efficiencies and cost savings. The National Renewable Energy Laboratory (NREL) is embarking on research projects that will perform multi-scale studies to investigate treatment technologies on varying source water quality that can operate under unsteady state flow conditions (timing) in response to energy availability, supply, and demand. A pilot-scale system will be built to assess modularity and scaling of treatment system operating conditions and performance. Modeling will also be used to determine feasibility and integration into grid emulation platforms to understand the use and value of operation with RE. Optimization will be completed for the treatment and distribution systems to document strengths and weaknesses of using RE power sources. Results of this effort will advance dynamic water treatment understanding under variable conditions and integration of RE sources to power these systems. This will lead to further decarbonization of municipal and private treatment facilities while supporting a paradigm shift in water treatment process operation and system design to achieve net-zero goals and move toward energy-sensitive resilient water treatment systems.

ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

Mapping Spatiotemporal Solvent Velocity from Measured Concentration Gradients in a Polarized Electrolyte

The electric-field induced motion of neutral species impedes the efficacy of electrochemical devices. By combining operando X-ray transmission measurements with continuum mechanics, we have developed a methodology for determining the velocity of neutral solvent molecules under an applied field. The X-ray transmission experiments were used to determine ion concentration profiles as a function of space and time in a polymer electrolyte. The unsteady state solvent mass balance equation was solved numerically with experimental concentration profiles to map spatiotemporal solvent velocities. We compare our experimentally derived results with predictions made with concentrated solution theory. We use the cation transference number as the only adjustable parameter to match experimental measurements of both concentration and solvent velocity. Our approach may be used to determine solvent velocity with any operando technique used to measure time-dependent ion concentration profiles.

Abdo, Emily E↗
DOE OSTI2023

Time-Gated Raman Spectroscopy (Intern Final Report)

Time-gated Raman spectroscopy is a powerful tool for rapid non-destructive analysis of chemical compounds, particularly when combined with other investigative methods. Time-gating the signal obtained from a sample during Raman spectroscopy allows for the removal of background emission and black body radiation, improving clarity and accuracy of scans. The main objectives of this project were to prove the lab’s time-gated Raman spectroscopy system works as proof of concept and to make select improvements to the in-house LabVIEW software used to run the time-gated system. This time-gated system will be added to the CTK’s TAP reactors in the future, but this is out of the scope of this internship project. This research is important, as understanding the mechanisms behind catalysts leads to designing better systems/materials and the perfection of catalysts, and energy-efficient chemical manufacturing is key when designing a greener future. The research for this project was done with samples at both ambient conditions in sample holders and at operando conditions (i.e., high temperature, unsteady-state) inside a small-volume chemical reactor. Numerous spectro-kinetic scans of the catalysts were obtained, but the advanced dynamics of these systems are still being interpreted. There is data that demonstrates that the time-gated system is correctly rejecting emission from samples, and one can observe real-time changes in the signal from the sample, due to coking or changes induced by redox chemistry. Chemical manufacturers and the planet are the main beneficiaries of catalysis research, as better catalysts will reduce carbon emissions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Unsteady- and Steady-State Relative Permeability Study with X-ray and Acoustic Monitoring for CO 2 Storage in Deep Saline Aquifers

In this study, we link the multiphase flow measurements with controls of sedimentary structures (e.g., heterogeneity and anisotropy) on relative permeability to variations in ultrasonic velocities for two deep saline aquifer formations (Entrada and Bluff Sandstones) in the San Juan Basin of the Southwestern USA. The rock specimens were extracted from outcrop sites near Durango, CO, USA. They have distinct differences in grain size, cementation composition, and individual chemical amounts, despite both formations being eolian sandstones. We performed a series of unsteady- and steady-state CO 2 -brine relative permeability experiments under capillary-controlled displacement rates. Unsteady-state experiments were conducted at 71 °C and 9.65 MPa; steady-state experiments were conducted at 85 °C and 22.8 MPa and 89 °C and 24.1 MPa for the Bluff and Entrada Sandstones, respectively. During the unsteady experiments, X-ray computed tomography was used to visualize multiphase flow in porous media and quantify saturations during brine drainage under various flow rates. Scan images and saturation profiles indicate that the CO 2 distribution in the pore volume was strongly impacted by the presence of high-angle cross-laminations, heterogeneous rock structure, and direction of bedding orientation. Those factors contribute to dramatic and quick initial breakthroughs and affect the overall saturation dynamics. Steady-state relative permeability tests were conducted at net flow rates of 1 mL/min for both brine drainage and imbibition. During the experimental steps, the CO 2 fractional flow was increased and decreased for both drainage and imbibition scenarios to mimic the front when CO 2 contacts brine and behind the front when brine enters space previously occupied by CO 2 . It was found that compressional velocity decreased, while shear waves slightly increased as brine saturation decreased. The hysteresis effects for the relative permeability and acoustic velocities were distinct. Furthermore, the CO 2 /brine front stability is quantified by applying a mobility ratio approach to spot saturations at which the boundary line between fluids becomes uneven. The results presented in this work can potentially boost the quality and precision of forecasts for the CO 2 storage projects in which the vertical and horizontal core-scale heterogeneity and anisotropy impact the plume migration within host reservoirs.

carbon dioxide (CO2)↗
DOE OSTI2025

Fuel reid vapor pressure level and ethanol content on stochastic preignition, effects at steady and unsteady engine operation

The present work investigates relations between fuel Reid vapor pressure (RVP) and biofuel (ethanol) content on stochastic preignition (SPI) at both sustained steady-state engine operation and following load transients. This work stems from in-field observations that automotive original equipment manufacturers have observed consistent seasonal increases in United States customer drivability complaints and warranty claims during September and October where SPI is suspected to be responsible. The seasonal timing of these events coincides with the United States seasonal fuel property changeover initiating on September 15 each year, where fuel RVP increases. To explore potential linkage between fuel RVP and SPI the present study employs engine SPI experiments coupled with laboratory spray measurements of fuels with RVPs of 8, 12, and 16 psi in both E10 (10% ethanol) and E25 (25% ethanol) fuels. Engine results are partitioned into fuel RVP and ethanol content effects on SPI in steady-state, sustained high-load engine operation and unsteady-state low- to high-load transitions, where off-engine spray vessel patternation and tip penetration results help to elucidate the observed fuel effects on SPI. A boosted direct-injected, spark-ignition engine was fueled with three market relevant E10 and E25 fuels with RVPs of 8, 12, and 16 to characterize the interplay between winter fuels and abnormal combustion behavior. The steady-state work shows that for high-load, steady-state engine operation, SPI is directly linked to fuel retention, which was found to be dependent on fuel distillation. The unsteady-state engine operation work shows that following low-to high-load transitions, SPI can occur from a memory of fuel property effects at low-load operation. Specifically, the fuel RVP effect on fuel spray collapse at low loads was found to correlate with SPI with a more than 95% confidence interval following low- to high-engine-load transitions. Results suggest that fuel-wall impingement at low-load operation could carry over into high-load transitions and generate SPI events following low- to high-load transitions.

09 BIOMASS FUELS↗
DOE OSTI2023

Capturing Time-Varying Wake Dynamics Using Hybridized Actuator Disks in Steady-State Simulations for Improved Optimization Efficiency

When optimizing turbine blade properties, using a time varying, unsteady simulation allows for high-fidelity representations of the wake dynamics. Unfortunately these types of optimizations are prohibitively expensive in terms of computation time and memory requirements. This presentation aims to show a method of transferring the wake dynamics of an unsteady simulation to a much more computationally tractable steady state simulation using hybrid actuator disks informed by unsteady, actuator line model dynamics.

actuator disk model↗
DOE OSTI2024

A theoretical/computational framework to measure SiO2 and MgO viscosity at high pressure

The convection of the mantle of Earth and super-Earths is important for many terrestrial phenomena, from plate tectonics to outgassing. Rheological properties, such as viscosity, regulate the transport of thermal energy and mass. However, the viscosity of mantle-relevant materials, such as MgO, at relevant pressures (>120 GPa) are not well constrained. The objective of this work was to develop a computational platform to simulate novel experiments aiming to measure the viscosity of MgO at high pressures. Experiments performed by our collaborators at Johns Hopkins University and Lawrence Livermore National Laboratory use the OMEGA EP laser facility to shock a corrugated MgO interface to 170 GPa, with resulting velocity evolution governed by the viscous Richtmyer-Meshkov instability. We used an in-house hydrocode to simulate this process and thus provide a bound on viscosity by comparing our simulations results to these experiments, as well as to examine the physical processes at play. We simulated the experiments with different values of MgO viscosity (from inviscid to 10,000 Pa∙s) while taking into account the unsteady laser pulse, the material's equation of state, and the rate-dependent constitutive relation for MgO, and the materials' equations of state. Our results suggest that MgO at these conditions has a viscosity within an order of magnitude of 5000 Pa∙s.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Performance characterization and modeling of an oscillating surge wave energy converter

Abstract Testing wave energy converters in the ocean could be expensive and complex, which necessitates the use of numerical modeling. However, accurately modeling the response of wave energy converters with high-fidelity simulations can be computationally intensive in the design stage where different configurations must be considered. Reduced-order models based on simplified equations of motion can be very useful in the design, optimization, or control of wave energy converters. Given the complex dynamics of wave energy converters, accurate representation, and evaluation of relative contributions by different forces are required. This effort is concerned with a performance characterization of the hydrodynamic response of an oscillating surge wave energy converter that is based on a reduced-order model. A state-space model is used to represent the radiation damping term. Morison’s representation of unsteady forces is used to account for the nonlinear damping. Wave tank tests are performed to validate simulations. A free response simulation is used to determine the coefficients of the state-space model. Torque-forced simulations are used to identify the coefficients of the nonlinear damping term for different amplitudes and wave frequencies. The impact of varying these coefficients on the response is investigated. An assessment of the capability of the model in predicting the hydrodynamic response under irregular forcing is performed. The results show that the maximum error is 3% when compared with high-fidelity simulations. It is determined that the nonlinear damping is proportional to the torque amplitude and its effects are more pronounced as the amplitude of the flap oscillations increases.

16 TIDAL AND WAVE POWER↗
DOE OSTI2026

Statorless mixed flow turbine for transonic pulsating inflow

Effective harvesting of power from high speed highly transient inflows such as the outflow of rotating detonation combustors (RDCs) is key to achieving their promised cycle efficiency step jump. To increase power density and efficiency simultaneously, a concept that can directly ingest transonic outflow, addressing the choking is needed, i.e. without additional transition elements. A new statorless design was assessed using a comprehensive approach to quantify all contributions to loss generation in transient flows, locally and globally. The turbine rotor was designed under steady flow conditions with a genetic algorithm. The comparison of power and loss generation between steady and unsteady flow results shows that a design methodology under steady-state conditions is suitable to characterize the performance of different designs. Three design families were further assessed under highly transient transonic conditions including traveling shock waves with a relative total pressure amplitude of 149.6% around the mean value and inflow angle variations from -26.5 deg to +51.8 deg. The oblique shock impinging on the pressure side (PS) of the turbine augments the shaft power extraction. When the oblique shock reflects between the pressure side and suction side (SS), its strength diminishes. This paper provides design guidelines on efficient turbine work extraction from the shocks emanating from detonation combustors.

rotating detonation engines↗
DOE OSTI2025

Chemical timescale effects on detonation convergence

Numerical simulations of detonation-containing flows have emerged as crucial tools for designing next- generation power and propulsion devices. As these tools mature, it is important for the combustion community to properly understand and isolate grid resolution effects when simulating detonations. To this end, the objective of this work is to provide a comprehensive analysis of the numerical convergence of unsteady detonation simulations, with focus on isolating the impacts of chemical timescale modifications on convergence characteristics in the context of operator splitting. With the aid of an AMReX-based adaptive mesh refinement flow solver-which enables resolutions up to ($\mathcal{O}$ (1000) cells-per-induction length-the convergence analysis is conducted using two kinetics configurations: (1) the simplified three-step Arrhenius-based model mechanism of Short and Quirk (1997), where chemical timescales in the detonation are modified by adjusting activation energies in the initiation and branching reactions, and (2) a detailed hydrogen- air mechanism, where the chemical timescales are adjusted by varying the ambient pressure. The convergence of unsteady self-sustained detonations in one-dimensional channels is then analyzed with reference to steady-state theoretical baseline solutions using these mechanisms. The goal of the analysis is to provide a detailed comparison of the effects of grid resolution on both macroscopic (peak pressures and wave speeds) and microscopic (wave structure) quantities of interest, drawing connections between the deviations from steady-state baselines and minimum chemical timescales. In particular, chemical timescale reductions were found to have minimal impact on the convergence of macroscopic properties. Furthermore, analyses of microscopic convergence trends, particularly in the reaction front location, revealed a key insight: maintaining the induction time while eliminating prohibitive chemical timescales through mechanism simplifications and combustion modeling can significantly enhance detonation convergence properties. Ultimately, this work uncovers resolution-dependent unsteady detonation convergence regimes and highlights the important role played by not only the chemical timescales, but also the ratio between the chemical timescale and induction time on the numerical convergence of the detonation wave structure.

Adaptive mesh refinement↗
DOE OSTI2026

Multiphysics Analysis of Li Cooled Divertor Substrate During Loss of Coolant Accident (LOCA)

In the ongoing study of potential designs for liquid metal (LM) plasma-facing components (PFCs), so-called “slow” and “fast” Li flow divertor concepts are under investigation. In the previous studies on design and analysis of the slow Li flow divertor and comparison with the fast Li flow divertor, the magnetohydrodynamics (MHD)/heat transfer effects of the Li flowing inside the substrate as a second coolant were comprehensively investigated under the normal steady-state operation conditions. Here, in the present study, the multiphysics analysis is extended to the unsteady abnormal divertor scenario where the Li layer on top of the substrate does not provide full coverage or even totally disappears for a certain period of time, so that the substrate becomes directly exposed to the incident high plasma heat flux. Such an unwanted event may happen regardless of the concept of the divertor and is worth detailed investigations, typically referred to as a loss of coolant accident (LOCA). To address this situation, a simplified scoping analysis is conducted first in 2-D, and then an integrated 3-D modeling is performed using a time-dependent multiphysics model in COMSOL Multiphysics that integrates LM MHD, heat transfer, and solid mechanics. The main goal is to evaluate conditions under which the major material limits, such as the maximum allowable temperature, stress, and displacement of the substrate, can still be met. It was shown that the maximum time over which the substrate of RAFM steel can retain structural integrity during the LOCA is around 0.2 ∼ 0.3 s. Any divertor concept that utilizes RAFM steel as a substrate material and liquid Li as a second coolant should take such a permitted time into consideration.

divertor↗
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 OSTI2024

A generalized wind turbine cross section as a reduced-order model to gain insights in blade aeroelastic challenges

In this work, we present an approach to study the aeroelastic stability of a wind turbine by focusing on the dynamics of a blade cross section. We present a methodology to obtain a reduced-order model of the blade dynamics in the form of generalized cross-sectional quantities that approximates the aerodynamic and structural properties of the full blade. The motivation for the work is to gain a physical understanding of the influence of aerodynamic models such as dynamic wake and dynamic stall on the frequency and damping of the structure using a reduced-order model with low computational cost. The model may be coupled to two-dimensional computational fluid dynamics softwares or engineering unsteady airfoil aerodynamics models accounting for dynamic wake and dynamic stall. In the latter case, we can obtain monolithic state-space forms of the aeroelastic system of equations, which simplifies the determination of the modal parameters and therefore the study of stability. The work investigates wind turbines in operation or at standstill, where vortex-induced vibrations and stall-induced vibrations, respectively, might be an issue. The implementation is made available as part of the open-source Python package WELIB and as part of the open-source unsteady aerodynamic driver of OpenFAST.

17 WIND ENERGY↗
DOE OSTI2023

Deep learning closure models for large-eddy simulation of flows around bluff bodies

Near-wall flow simulation remains a central challenge in aerodynamics modelling: Reynolds-averaged Navier–Stokes predictions of separated flows are often inaccurate, and large-eddy simulation (LES) can require prohibitively small near-wall mesh sizes. A deep learning (DL) closure model for LES is developed by introducing untrained neural networks into the governing equations and training in situ for incompressible flows around rectangular prisms at moderate Reynolds numbers. The DL-LES models are trained using adjoint partial differential equation (PDE) optimization methods to match, as closely as possible, direct numerical simulation (DNS) data. They are then evaluated out-of-sample – for aspect ratios, Reynolds numbers and bluff-body geometries not included in the training data – and compared with standard LES models. The DL-LES models outperform these models and are able to achieve accurate LES predictions on a relatively coarse mesh (downsampled from the DNS mesh by factors of four or eight in each Cartesian direction). We study the accuracy of the DL-LES model for predicting the drag coefficient, near-wall and far-field mean flow, and resolved Reynolds stress. A crucial challenge is that the LES quantities of interest are the steady-state flow statistics; for example, a time-averaged velocity component $\langle {u}_i\rangle (x) = \lim _{t \rightarrow \infty } ({1}/{t}) \int _0^t u_i(s,x)\, {\rm d}s$ . Calculating the steady-state flow statistics therefore requires simulating the DL-LES equations over a large number of flow times through the domain. It is a non-trivial question whether an unsteady PDE model with a functional form defined by a deep neural network can remain stable and accurate on $t \in [0, \infty )$ , especially when trained over comparatively short time intervals. Our results demonstrate that the DL-LES models are accurate and stable over long time horizons, which enables the estimation of the steady-state mean velocity, fluctuations and drag coefficient of turbulent flows around bluff bodies relevant to aerodynamics applications.

Mechanics↗