Search NASA⌕ Search

SEARCH · Search NASA

Results for “Turbulence simulations”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Simulations of a hypersonic turbulent boundary layer over wavy surfaces

Here, we conduct large-eddy simulations of a Mach 5.84 cold wall turbulent boundary layer over one-dimensional wavy walls with varying amplitudes and wavelengths. Across all wall topologies, a series of alternating shock and expansion waves is shown to influence the entire boundary layer, and generate repeating wave patterns in the turbulent stresses, dispersive stresses, and turbulent kinetic energy budget. The series of alternating shocks and expansions imposes repeating adverse and favourable pressure gradients across the wavy wall, and at sufficient wall amplitude, triggers flow separation in the trough of the wave. Flow separation is demonstrated to influence the behaviour of wall pressure fluctuations over the wavy wall. In attached flows, the prominent frequencies are consistent with integral-scale boundary layer turbulence, whereas in separated flows, a two-decade frequency range is present, akin to two-dimensional shock–boundary layer interactions. Counter-rotating streamwise-oriented structures are observed on the windward side of the wave, which diminish over the wave crest. A conditional analysis demonstrates that these structures are present in the upstream boundary layer, and are amplified with increasing wall amplitude. An examination of the Görtler number and probability density function (PDF) of the fluctuating lateral wall shear stress demonstrates the strong correlation between a large Görtler number and growth of the PDF tail density, suggesting that the amplification of the counter-rotating streamwise-oriented structures are linked to centrifugal instabilities in regions of streamline concavity.

boundary layers↗
DOE OSTI2021

An LES study of secondary motion and wall shear stresses in a pipe bend

In this study, large-eddy simulations (LES) of a single-phase, turbulent flow in a 90 ° pipe bend are performed at three Reynolds numbers (5300, 27 000, and 45 000) to investigate the correlation between secondary flow motion and wall shear stresses, which is suspected to be a potential mechanism responsible for material erosion. The isothermal flows are validated against available experimental and numerical data first. The snapshot proper orthogonal decomposition (POD) is applied for the medium and high Reynolds number flows to identify the secondary flow motions and the oscillation of the Dean vortices that are found to cause swirl-switching. Distinguished frequencies of the POD time coefficients at Strouhal numbers of 0.25 and 0.28 are identified for Reynolds numbers at 27 000 and 45 000, respectively. Moreover, shear stress on the pipe wall and the associated power spectral density are obtained and shown to have the same oscillating frequency as the swirl-switching.

42 ENGINEERING↗
DOE OSTI2024

GX: a GPU-native gyrokinetic turbulence code for tokamak and stellarator design

GX is a code designed to solve the nonlinear gyrokinetic system for low-frequency turbulence in magnetized plasmas, particularly tokamaks and stellarators. In GX, our primary motivation and target is a fast gyrokinetic solver that can be used for fusion reactor design and optimization along with wide-ranging physics exploration. Here, this has led to several code and algorithm design decisions, specifically chosen to prioritize time to solution. First, we have used a discretization algorithm that is pseudospectral in the entire phase space, including a Laguerre–Hermite pseudospectral formulation of velocity space, which allows for smooth interpolation between coarse gyrofluid-like resolutions and finer conventional gyrokinetic resolutions and efficient evaluation of a model collision operator. Additionally, we have built GX to natively target graphics processors (GPUs), which are among the fastest computational platforms available today. Finally, we have taken advantage of the reactor-relevant limit of small $\rho _*$ by using the radially local flux-tube approach. In this paper we present details about the gyrokinetic system and the numerical algorithms used in GX to solve the system. We then present several numerical benchmarks against established gyrokinetic codes in both tokamak and stellarator magnetic geometries to verify that GX correctly simulates gyrokinetic turbulence in the small $\rho _*$. Moreover, we show that the convergence properties of the Laguerre–Hermite spectral velocity formulation are quite favourable for nonlinear problems of interest. Coupled with GPU acceleration, which we also investigate with scaling studies, this enables GX to be able to produce useful turbulence simulations in minutes on one (or a few) GPUs and higher fidelity results in a few hours using several GPUs. GX is open-source software that is ready for fusion reactor design studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Separability of microtearing mode and electron temperature gradient turbulence regimes

The separability of microtearing mode (MTM)-dominated and electron temperature gradient (ETG)-driven turbulence regimes is studied with multiscale nonlinear gyrokinetic simulations. The simulations are based on National Spherical Torus Experiment-like, high-confinement mode pedestal parameters, where electromagnetic perturbations are large. Linear analysis indicates a wide scale-separation between the MTM and ETG modes in binormal wavenumber space (perpendicular to the magnetic field line), with no unstable modes at intermediate scales. Likewise, single-scale nonlinear analyses, retaining ion-only or electron-only spatio-temporal scales, produce seemingly well-converged transport states. Surprisingly, the multiscale simulation, which contains both the ion and electron scales, closely follows the transport from the electron-scale simulation. This trend is robust over a wide range of electron temperature gradient. Remarkably, compared to ion-scale simulations, MTM turbulence is significantly reduced at multiscale resolution even when ETG turbulence is low. In this case, traditional ion-scale resolution overestimates the electron energy flux, and it is not possible to accurately simulate the MTM turbulence with separable ion-scale simulations. Here, while the analysis confirms the validity of electron-scale simulations for predicting the electron transport, it also indicates that multiscale simulation may be required for reproducing the turbulence spectrum for systems with coupled MTM-ETG turbulence.

H-mode↗
DOE OSTI2021

The discrete Green's function paradigm for two-way coupled Euler–Lagrange simulation

We outline a methodology for the simulation of two-way coupled particle-laden flows. The drag force that couples fluid and particle momentum depends on the undisturbed fluid velocity at the particle location, and this latter quantity requires modelling. We demonstrate that the undisturbed fluid velocity, in the low particle Reynolds number limit, can be related exactly to the discrete Green's function of the discrete Stokes equations. In addition to hydrodynamics, the method can be extended to other physics present in particle-laden flows such as heat transfer and electromagnetism. The discrete Green's functions for the Navier–Stokes equations are obtained at low particle Reynolds number in a two-plane channel geometry. We perform verification at different Reynolds numbers for a particle settling under gravity parallel to a plane wall, for different wall-normal separations. Compared with other point-particle schemes, the Stokesian discrete Green's function approach is the most robust at low particle Reynolds number, accurate at all wall-normal separations. To account for degradation in accuracy away from the wall at finite Reynolds number, we extend the present methodology to an Oseen-like discrete Green's function. The extended discrete Green's function method is found to be accurate within 6% at all wall-normal separations for particle Reynolds numbers up to 24. Furthermore, the discrete Green's function approach is well suited to dilute systems with significant mass loading and this is highlighted by comparison against other Euler–Lagrange as well as particle-resolved simulations of gas–solid turbulent channel flow. Strong particle–turbulence coupling is observed in the form of turbulence modification and turbophoresis suppression, and these observations are placed in context of the different methods.

42 ENGINEERING↗
DOE OSTI2024

Turbulence modelling in neutron star merger simulations

Observations of neutron star mergers have the potential to unveil detailed physics of matter and gravity in regimes inaccessible by other experiments. Quantitative comparisons to theory and parameter estimation require nonlinear numerical simulations. However, the detailed physics of energy and momentum transfer between different scales, and the formation and interaction of small scale structures, which can be probed by detectors, are not captured by current simulations. This is where turbulence enters neutron star modelling. This review will outline the theory and current status of turbulence modelling for relativistic neutron star merger simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A Large Eddy Simulation Study of Flow Turbulence, Alumina Transport, and Bath Temperature Evolution in Conventional Aluminum-Smelting Cell Using OpenFOAM

In this study, a Large Eddy Simulation (LES) of the aluminum-smelting process is performed using OpenFOAM. To understand the coupled behavior of heat transfer, mass transfer, and flow of the smelting process, a multi-physics computational fluid dynamics (CFD) model based on the Eulerian–Eulerian multi-fluid approach is adopted. The model accounts for CO 2 bubble and magnetohydrodynamics (MHD)-driven flow, along with alumina dissolution, transport, and bath temperature evolution. The simulation predictions show small-scale turbulent vortical structures in the anode–cathode space caused by combined effect of MHD and CO 2 bubble-bath interactions and relatively large-scale asymmetric vortices in the inter-anode space caused by the CO 2 bubble-bath interactions. The vortex formation at the edges of the anodes evidently aids in transporting alumina from the central channel to the bottom of the anodes and prevents accumulation of gas bubbles in the periphery of the anode bottom. Symmetric bath cold spots are observed in the vicinity of the feeder. Cold spots are also observed in the anode–cathode distance space below the anode bottom due to the transport of undissolved solid to this region by the flow. The findings from the work are useful in developing and designing alumina-feeding strategy leading to reduced anode effects and smooth operation of the cell. The work also highlights the important flow structures in conventional aluminum-smelting cell.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Frame invariant neural network closures for Kraichnan turbulence

Numerical simulations of geophysical and atmospheric flows have to rely on parameterizations of subgrid scale processes due to their limited spatial resolution. Despite substantial progress in developing parameterization (or closure) models for subgrid scale (SGS) processes using physical insights and mathematical approximations, they remain imperfect and can lead to inaccurate predictions. In recent years, machine learning has been successful in extracting complex patterns from high-resolution spatio-temporal data, leading to improved parameterization models, and ultimately better coarse grid prediction. However, the inability to satisfy known physics and poor generalization hinders the application of these models for real-world problems. In this work, we put forth a frame invariant closure approach to improve the accuracy and generalizability of deep learning-based subgrid scale closure models by embedding physical symmetries directly into the structure of the neural network. Specifically, we utilized specialized layers within the convolutional neural network in such a way that desired constraints are theoretically guaranteed without the need for any regularization terms. We demonstrate our framework for a two-dimensional decaying turbulence test case mostly characterized by the forward enstrophy cascade. We show that our frame invariant SGS model (i) accurately predicts the subgrid scale source term, (ii) respects the physical symmetries such as translation, Galilean, and rotation invariance, and (iii) is numerically stable when implemented in coarse-grid simulation with generalization to different initial conditions and Reynolds number. This work opens up a possibility of connecting physics-based theories and data-driven modeling paradigms, and thus represents a promising step towards the development of physically consistent data-driven turbulence closure models.

42 ENGINEERING↗
DOE OSTI2025

Implementation and Evaluation of Physics-Driven Dynamic Entrainment-Mixing Parameterization in a Climate Model and Its Impact on Low-Cloud Simulation

The turbulent entrainment-mixing process in the Community Earth System Model version 1.2 (CESM1.2) is assumed to follow the extremely inhomogeneous entrainment-mixing. However, different entrainment-mixing scenarios can occur in real clouds. To address this deficiency, a unifying parameterization that represents different entrainment-mixing processes is implemented and evaluated in CESM1.2. The results indicate that the homogeneous mixing degree values simulated by the new parameterization in CESM1.2 are predominantly greater than 50%, suggesting a tendency toward homogeneous mixing. Compared to the extremely inhomogeneous mixing mechanism, the new parameterization increases the cloud droplet number concentration (Nc). More importantly, the new parameterization improves low-cloud fraction (CLDLOW) simulation in Northwest Pacific (NWP) and Southeast Pacific (SEP) regions, with relative improvements of 2.95% and 4.17%, respectively. Furthermore, the improvements reach up to 44.6% and 16.2% in the NWP and SEP regions, respectively, when considering the relationship between N c and CLDLOW. Further analysis reveals that the new parameterization enhances cloud optical depth, longwave radiative cooling effect, net condensation rate, cloud water mixing ratio, lower-troposphere stability, and CLDLOW by increasing N c . Additionally, these results underscore the importance of improving entrainment-mixing parameterization in climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Reynolds number dependence of Lagrangian dispersion in direct numerical simulations of anisotropic magnetohydrodynamic turbulence

Large-scale magnetic fields thread through the electrically conducting matter of the interplanetary and interstellar medium, stellar interiors and other astrophysical plasmas, producing anisotropic flows with regions of high-Reynolds-number turbulence. It is common to encounter turbulent flows structured by a magnetic field with a strength approximately equal to the root-mean-square magnetic fluctuations. In this work, direct numerical simulations of anisotropic magnetohydrodynamic (MHD) turbulence influenced by such a magnetic field are conducted for a series of cases that have identical resolution, and increasing grid sizes up to $2048^3$ . The result is a series of closely comparable simulations at Reynolds numbers ranging from 1400 up to 21 000. We investigate the influence of the Reynolds number from the Lagrangian viewpoint by tracking fluid particles and calculating single-particle and two-particle statistics. The influence of Alfvénic fluctuations and the fundamental anisotropy on the MHD turbulence in these statistics is discussed. Single-particle diffusion curves exhibit mildly superdiffusive behaviours that differ in the direction aligned with the magnetic field and the direction perpendicular to it. Competing alignment processes affect the dispersion of particle pairs, in particular at the beginning of the inertial subrange of time scales. Scalings for relative dispersion, which become clearer in the inertial subrange for a larger Reynolds number, can be observed that are steeper than indicated by the Richardson prediction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A review of tidal energy—Resource, feedbacks, and environmental interactions

The ocean contains a variety of renewable energy resources, little of which has been exploited. Here, we review both tidal range and tidal stream energy, with a focus on the resource, feedbacks, and environmental interactions. The review covers a wide range of timescales of relevance to tidal energy, from fortnightly (spring-neap) and semi-diurnal variability, down to array, and device-scale turbulence. When simulating the regional tidal energy resource, and to assess environmental impacts, it is necessary to account for feedbacks between the tidal array and the resource itself. We critically review various methods for simulating energy extraction, from insights gained through theoretical studies of “tidal fences” in idealized channels, to realistic three-dimensional model studies with complex geometry and arrays of turbines represented by momentum sinks and additional turbulence due to the presence of rotors and support structures. We discuss how variability can be reduced by developing multiple (aggregated) sites with a consideration of the enhanced phase diversity offered by exploiting less energetic tidal currents. This leads to future research questions that have not yet been explored in depth at first-generation tidal sites in relatively sheltered channels (e.g., the interaction of waves with currents). Such enhanced understanding of real sea conditions, including the effects of wind and waves, leads to our other identified primary future research direction—reduced uncertainties in turbulence predictions, including the development of realistic models that simulate the interaction between ambient turbulence and the turbulence resulting from multiple wakes, and changes to system-wide hydrodynamics, water quality, and sedimentation.

Environmental impacts↗
DOE OSTI2021

A physics-based model for frost buildup under turbulent flow using direct numerical simulations

We present a new model for frost buildup under turbulent (and laminar) flow using direct numerical simulations. The physical model consists of two layers, the air and the frost. The air layer is fully resolved and consists of solving for the velocity, temperature, and vapor mass fraction fields. The frost layer thickness is resolved using conservation of mass and energy. Both phases are dynamically coupled using the immersed boundary method. Three-dimensional simulations are conducted in an open-channel configuration. A number of challenges need to be overcome to make these simulations feasible. First, to enforce far-field conditions of zero gradient and prescribed mean temperature and humidity, a source term is added to the energy and transport equations in the flow solver. Second, the mean frost thickness is subtracted after each time step to ensure a constant mean flow thickness and level of turbulence in the numerical domain. Third, a slow-time acceleration approach, which accelerates the frost buildup by a predetermined factor, is employed to bridge the gap between the fast turbulent and slow frost buildup time scales. Finally, a frost densification scheme is used to overcome the difficulties of vertically varying frost properties. The model is validated by comparing the frost thickness and frost thickness buildup rate over a period of one hour from a cooled flat plate experiment. As a result, both quantities compare favorably with experiments.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Investigation of turbulent inflow specification in Euler–Lagrange simulations of mid-field spray

The process of atomization of a liquid jet by a parallel high-speed gas stream results in a spray, whose downstream development is of considerable interest to several applications. The round jet spray can be spatially divided into (i) a near-field (near-nozzle) region of liquid atomization and (ii) a downstream mid-field region of fully-dispersed droplets. In order to accurately model mid-field droplet dispersion, this work aims at developing a rigorous and robust injection model for Euler–Lagrange spray simulations. Results from experiments are used to obtain the relevant droplet number density, size distribution, and mean and standard deviation velocity distributions of the injection model, systematically in a step-by-step process. Two-phase large eddy simulations are performed by stochastically generating the Lagrangian droplets at the inlet of the mid-field region. Number flux, diameter distribution, mean velocity, and other time-averaged statistics at several downstream locations are shown to agree well with the corresponding experimental data.

42 ENGINEERING↗
DOE OSTI2022

Adaptive Mesh Refinement Large Eddy Simulation of the Supercritical Carbon Dioxide Round Turbulent Jet

Supercritical carbon dioxide (sCO2) is of interest to a range of engineering problems, including carbon capture, utilization, and storage (CCUS) as well as advanced cycles for power generation. Non-ideal variations in physical properties of sCO2 impact the physics of these systems. In this study, we simulate turbulent sCO2 jets to gain a better understanding of these physics.We use a second order finite volume method with adaptive mesh refinement as implemented in the first-principles simulation code PeleC to perform a Large Eddy Simulation (LES) of three turbulent jets of sCO2. Additionally, we use the Soave-Redlich-Kwong equation of state to close the system and examine the impact of a cubic equation of state on the turbulent flow physics. We look at velocity and Reynolds stress profiles at different downstream locations for three cases in which the temperature of the jet andthat of the ambient fluid differ in order to capture the effects of widely varying thermal properties in the pseudocritical region. These results are then contrasted with established theory for ideal gas jets.

adaptive mesh refinement↗
DOE OSTI2022

Impact of Cloud-Base Turbulence on CCN Activation: Single-Size CCN

Abstract This paper examines the impact of cloud-base turbulence on activation of cloud condensation nuclei (CCN). Following our previous studies, we contrast activation within a nonturbulent adiabatic parcel and an adiabatic parcel filled with turbulence. The latter is simulated by applying a forced implicit large-eddy simulation within a triply periodic computational domain of 64 3 m 3 . We consider two monodisperse CCN. Small CCN have a dry radius of 0.01 μ m and a corresponding activation (critical) radius and critical supersaturation of 0.6 μ m and 1.3%, respectively. Large CCN have a dry radius of 0.2 μ m and feature activation radius of 5.4 μ m and critical supersaturation 0.15%. CCN are assumed in 200-cm −3 concentration in all cases. Mean cloud-base updraft velocities of 0.33, 1, and 3 m s −1 are considered. In the nonturbulent parcel, all CCN are activated and lead to a monodisperse droplet size distribution above the cloud base, with practically the same droplet size in all simulations. In contrast, turbulence can lead to activation of only a fraction of all CCN with a nonzero spectral width above the cloud base, of the order of 1 μ m, especially in the case of small CCN and weak mean cloud-base ascent. We compare our results to studies of the turbulent single-size CCN activation in the Pi chamber. Sensitivity simulations that apply a smaller turbulence intensity, smaller computational domain, and modified initial conditions document the impact of specific modeling assumptions. The simulations call for a more realistic high-resolution modeling of turbulent cloud-base activation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Automated and efficient local adaptive regression for principal component-based reduced-order modeling of turbulent reacting flows

Principal Component Analysis can be used to reduce the cost of Computational Fluid Dynamics simulations of turbulent reacting flows by reducing the dimensionality of the transported variables through projection of the thermochemical state onto a lower-dimensional manifold. However, because of the nonlinearity of the principal component source terms, nonlinear regression techniques must be utilized for the source terms in terms of the principal components. Unfortunately, widely available and utilized nonlinear regression techniques can have prohibitive computational requirements and/or accuracy that is highly dependent on user experience in ad hoc tuning of model architecture and hyperparameters. Here, in this work, a new nonlinear regression approach is proposed that is both computationally efficient and automated so does not require any user input. The approach is evaluated through a priori prediction of principal component source terms using data from a Direct Numerical Simulation of a turbulent nonpremixed n-heptane/air jet flame. In particular, the proposed framework consists of local regressions whose complexity is adapted according to the local nonlinearity of the data: local linear regression when accurate enough and local Artificial Neural Networks when nonlinear regression is required. The number of local clusters for local regression is determined automatically using the Davies-Bouldin index. In addition, Bayesian optimization is utilized for model training (i.e., to select the best architectures and hyperparameters of the nonlinear regressions in an unsupervised fashion), eliminating ad hoc hand-tuning and/or expensive grid searches. Overall, compared to a single, global neural network, the new local adaptive regression approach is shown to have comparable accuracy but 69% less training time due to the utilization of local linear regression and faster training of local neural networks.

42 ENGINEERING↗
DOE OSTI2021

Assessment of Turbulent Combustion Models for Simulating Prechamber Ignition in a Natural Gas Engine

Lean combustion is a promising strategy to increase thermal efficiency in an internal combustion engine, by exploiting a favorable specific heat ratio of the fresh mixture while simultaneously suppressing the heat losses to the cylinder wall. However, unstable ignition and slow flame propagation at fuel-lean conditions lead to large cycle-to-cycle variability and limit the high-efficiency engine operating range. Prechamber ignition is considered an effective concept to extend the lean operating limit, by providing spatially distributed ignition with multiple turbulent flame-jets and enabling a faster combustion rate compared to the conventional spark ignition approach. From a numerical modeling standpoint to date science base and available simulation tools are inadequate to properly understand and predict the combustion processes in prechamber ignited engines. In this paper, conceptually different Reynolds-averaged Navier–Stokes (RANS) combustion models widely adopted in the engine modeling community are used to simulate the ignition and combustion processes in a medium-duty natural gas engine with a prechamber spark-ignition system. Here, a flamelet-based turbulent combustion model, i.e., G-equation, and a multizone well-stirred reactor model are employed for this modeling study. Simulation results are compared with experimental data in terms of in-cylinder pressure and heat release rate. Finally, the analysis of the performance of the two models is carried out to highlight the strengths and limitations of the two evaluated approaches.

42 ENGINEERING↗
DOE OSTI2020

Assessment of turbulent combustion models for simulating pre-chamber ignition in a natural gas engine

Lean combustion in an internal combustion engine is a promising strategy to increase thermal efficiency by leveraging a more favorable specific heat ratio of the fresh mixture and simultaneously suppressing the heat losses to the cylinder wall. However, unstable ignition events and slow flame propagation at fuel-lean condition lead to high cycle-to-cycle variability and hence limit the high-efficiency engine operating range. Pre-chamber ignition is considered an effective concept to extend the lean operating limit, by providing spatially distributed ignition with multiple turbulent flame-jets and enabling faster combustion rate compared to the conventional spark ignition approach. From a numerical modeling perspective, to date, still the science base and available simulation tools are inadequate for understanding and predicting the combustion processes in pre-chamber ignited engines. In this paper, conceptually different RANS combustion models widely adopted in the engine modeling community were used to simulate the ignition and combustion processes in a medium-duty natural gas engine with a pre-chamber spark-ignition system. A flamelet-based turbulent combustion model, i.e., G-equation, and a multi-zone well-stirred reactor model were employed for the multi-dimensional study. Simulation results were compared with experimental data in terms of in-cylinder pressure and heat release rate. Finally, the analysis of the performance of the two models is carried out to highlight the strengths and limitations of the two formulations respectively.

combustion model evaluation↗