Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Turbulence simulation”

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

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

At least 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Interface capturing simulations of bubble population effects in PWR subchannels

As the computational power of high-performance computing (HPC) facilities grows, so too does the feasibility of using first principle based simulation to study turbulent two-phase flows within complex pressurized water reactor (PWR) geometries. Direct numerical simulation (DNS), integrated with an interface capturing method, allows for the collection of high-fidelity numerical data using advanced analysis techniques. The research presented here employs the massively parallel, finite-element based, unstructured mesh code, PHASTA, to simulate a set of two-phase bubbly flows through PWR subchannel geometries including auxiliary structures (spacer grids and mixing vanes). The main objective of the presented work is to analyze bubble dynamics and turbulence interactions at varying bubble concentrations to support the development of advanced two-phase flow closure models. Turbulent two-phase flows in PWR subchannels were simulated at hydraulic Reynolds numbers of 81,000 with bubble concentrations of 3%–15% by gas volume fraction (768–3928 resolved bubbles, respectively) and compared against a 1% void fraction case (262 bubbles) that had been previously simulated. The finite element mesh utilized for the study at higher bubble concentrations was composed of 1.55 billion elements, compared to the previous study which employed 1.11 billion elements, ensuring all turbulence scales and individual bubbles within the flow are fully resolved. For each case, the resolved initial bubble size was 0.65 mm in diameter (resolved with 25 grid points across the diameter). The simulations were analyzed to find flow features such as the mean velocity profile, bubble relative velocity and the effect of the bubbles on the turbulent conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRSDate not supplied

Wall-modeled LES of Turbulent Flow Over a Two-dimensional Gaussian Bump

We perform wall-modeled large eddy simulations (WMLES) of turbulent flow over a nominally two-dimensional Gaussian-shaped bump geometry to assess its performance in the accelerating and separation regions of the flow. The flow conditions are based on high-fidelity numerical simulation of Uzun & Malik (AIAA Journal 2022). The oncoming flow Mach number is 0.2, with the bump length-based Reynolds number (ReL) of 2 million. In our previous study, while WMLES with the constant coefficient Vreman subgrid scale model and equilibrium wall model performed satisfactorily at lower ReL = 1 million, it failed to predict flow separation at the higher Re, contrary to the experimental observations. We investigate the sensitivity of WMLES to different subgrid scale models, wall models, and grid resolution and topology on flow separation by comparing with available data.

Computational Fluid Dynamics↗
DOE OSTI2022

Evolution of the marker distribution in gyrokinetic $δf$ particle-in-cell simulations

The evolution of the particle weight in a δf particle-in-cell simulation depends on the marker distribution that can evolve in a turbulent field due to turbulent diffusion. When Monte Carlo methods are used to implement the test particle collision operator, or when the particle motion is not strictly Hamiltonian in a collisionless simulation, the marker distribution will evolve along the particle trajectory and, in general, cannot be known exactly. A two-dimensional numerical marker distribution is proposed as an approximation. It is shown to be advantageous over other common methods for evaluating the marker distribution in long-time turbulence simulations. A generalized two-weight δf-method is proposed to mitigate the marker evolution problem.

, Monte Carlo methods↗
DOE OSTI2023

Mixing across stable density interfaces in forced stratified turbulence

Understanding how turbulence enhances irreversible scalar mixing in density-stratified fluids is a central problem in geophysical fluid dynamics. While isotropic overturning regions are commonly the focus of mixing analyses, we here investigate whether significant mixing may arise in anisotropic statically stable regions of the flow. Focusing on a single forced direct numerical simulation of stratified turbulence, we analyse spatial correlations between the vertical density gradient $\partial \rho /\partial z$ and the dissipation rates of kinetic energy $\epsilon$ and scalar variance $\chi$ , the latter quantifying scalar mixing. The domain is characterized by relatively well-mixed density layers separated by sharp stable interfaces that are correlated with high vertical shear. While static instability is most prevalent within the mixed layers, much of the scalar mixing is localized to the intervening interfaces, a phenomenon not apparent if considering local static instability or $\epsilon$ alone. While the majority of the domain is characterized by the canonical flux coefficient $\varGamma \equiv \chi /\epsilon =0.2$ , often assumed in ocean mixing parametrizations, extreme values of $\chi$ within the statically stable interfaces, associated with elevated $\varGamma$ , strongly skew the bulk statistics. Our findings suggest that current parametrizations of turbulent mixing may be biased by undersampling, such that the most common, but not necessarily the most significant, mixing events are overweighted. Having focused here on a single simulation of stratified turbulence, it is hoped that our results motivate a broader investigation into the role played by stable density interfaces in mixing, across a wider range of parameters and forcing schemes representative of ocean turbulence.

Couchman, Miles M. P. (ORCID:0000000246676829)↗
DOE OSTI2023

Statistics of pressure fluctuations in turbulent kinetic plasmas

ABSTRACT In this study, we explore the statistics of pressure fluctuations in kinetic collisionless turbulence. A 2.5D kinetic particle-in-cell simulation of decaying turbulence is used to investigate pressure balance via the evolution of thermal and magnetic pressure in a plasma with β of order unity. We also discuss the behaviour of thermal, magnetic, and total pressure structure functions and their corresponding wavenumber spectra. The total pressure spectrum exhibits a slope of −7/3 extending for about a decade in the ion-inertial range. In contrast, shallower −5/3 spectra are characteristic of the magnetic pressure and thermal pressure. The steeper total pressure spectrum is a consequence of cancellation caused by density-magnetic field magnitude anti-correlation. Further, we evaluate higher order total pressure structure functions in an effort to discuss intermittency and compare the power exponents with higher order structure functions of velocity and magnetic fluctuations. Finally, applications to astrophysical systems are also discussed.

Astronomy & Astrophysics↗
NASA NTRS1988

Aerodynamic applications of an efficient incompressible Navier-Stokes solver

A robust and accurate upwind scheme based on flux-difference splitting for three-dimensional, incompressible flows in general coordinates is described. Viscous effects are included through a thin-layer approximation to the Navier-Stokes equations. Turbulent flows are simulated with extensions to the Baldwin-Lomax (1978) turbulence model. The effects of accuracy and gridding on the computed flow field results are assessed. The results obtained from the code for vortical flows over a delta wing, a double delta wing, and a tangent-ogive forebody compare well with experimental data.

Hartwich, Peter M.↗
NASA NTRS1985

Experimental aerodynamic heating to simulated space shuttle tiles in laminar and turbulent boundary layers with variable flow angles at a nominal Mach number of 7

The heat transfer to simulated shuttle thermal protection system tiles was investigated experimentally by using a highly instrumented metallic thin wall tile arranged with other metal tiles in a staggered tile array. Cold wall heating rate data for laminar and turbulent flow were obtained in the Langley 8 foot high Temperature Tunnel at a nominal Mach number of 7, a nominal total temperature of 3300R, a free stream unit Reynolds number from 3.4 x 10 sup 5 to 2.2 10 sup 6 per foot, and a free stream dynamic pressure from 2.1 to 9.0 psia. Experimental data are presented to illustrate the effects of flow angularity and gap width on both local peak heating and overall heating loads. For the conditions of the present study, the results show that localized and total heating are sensitive to changes in flow angle only for the test conditions of turbulent boundary layer flow with high kinetic energy and that a flow angle from 30 deg to 50 deg will minimize the local heating.

Avery, D. E.↗
DOE OSTI2023

Electrostatic turbulence in EAST plasmas with internal transport barrier

Based on first-principles nonlinear gyrokinetic simulations, the electrostatic turbulence properties in the internal transport barrier (ITB) region of an Experimental Advanced Superconducting Tokamak discharge (#93890) are investigated. Specifically, ITBs with steep density and temperature gradients are located in the weakly negative magnetic shear region at the plasma center. In the linear stage, the growth rate and frequency of the ion temperature gradient (ITG) mode increase significantly due to resonant excitation by trapped electrons. That is, the resonance between trapped electrons and the ITG becomes strong due to the precession drift reversal of trapped electrons by the negative magnetic shear and Shafranov shift. Meanwhile, the trapped electron mode is stable in the ITB region due to only a very small fraction of electrons precessing in the direction of the electron diamagnetic drift. Nonlinear simulations show that, after considering the non-adiabatic effect of trapped electrons, the heat conductivity of ions and the turbulence intensity increase by at least a factor of 7 compared with the results only considering the adiabatic effect of electrons. The zonal charge density of trapped electrons can partially cancel that of ions, which weakens the intensity of the zonal flow, and consequently reduces the zonal flow regulation and enhances the turbulent transport.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Shaping effects on scrape-off layer plasma turbulence: A rigorous validation of three-dimensional simulations against TCV measurements

The effect of plasma shaping on scrape-off layer (SOL) plasma turbulence is investigated through a rigorous validation exercise. Two- and three-dimensional simulations of the SOL plasma dynamics in three TCV limited discharges are carried out with the GBS code. These discharges realize an almost circular magnetic equilibrium, an elongated equilibrium, and an elongated equilibrium with negative triangularity. For the three plasma discharges, three simulations are performed, considering (i) a three-dimensional model with an explicit dependence on elongation, triangularity, and the inverse aspect ratio; (ii) a circular three-dimensional model in the infinite aspect ratio limit; and (iii) a two-dimensional model, which assumes a circular magnetic equilibrium in the infinite aspect ratio limit, cold ions, and interchange driven turbulence in the sheath limited regime. Ten validation observables common to simulations and experimental measurements from a reciprocating probe located at the TCV outer midplane are identified, and the agreement between experimental and numerical results relative to each observable is evaluated. The composite metric introduced by Ricci et al. [Phys. Plasmas 18, 032109 (2011)] is then used to assess the overall agreement between simulations and experimental measurements. It is found that the shaping model implemented in GBS improves the description of SOL plasma turbulence, taking into account the impact of elongation and triangularity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Dynamic Subgrid Turbulence Modeling for Shallow Cumulus Convection Simulations beyond LES Resolutions

A scale-dependent dynamic Smagorinsky model is implemented in the Met Office/NERC Cloud (MONC) model using two averaging flavors, along Lagrangian pathlines and local moving averages. The dynamic approaches were compared against the conventional Smagorinsky–Lilly scheme in simulating the diurnal cycle of shallow cumulus convection. The simulations spanned from the LES to the near-gray-zone and gray-zone resolutions and revealed the adaptability of the dynamic model across the scales and different stability regimes. The dynamic model can produce a scale- and stability-dependent profile of the subfilter turbulence length scale across the chosen resolution range. At gray-zone resolutions the adaptive length scales can better represent the early precloud boundary layer leading to temperature and moisture profiles closer to the LES compared to the standard Smagorinsky. As a result, the initialization and general representation of the cloud field in the dynamic model is in good agreement with the LES. In contrast, the standard Smagorinsky produces a less well-mixed boundary layer, which fails to ventilate moisture from the boundary layer, resulting in the delayed spinup of the cloud layer. Moreover, strong downgradient diffusion controls the turbulent transport of scalars in the cloud layer. However, the dynamic approaches rely on the resolved field to account for nonlocal transports, leading to overenergetic structures when the boundary layer is fully developed and the Lagrangian model is used. Introducing the local averaging version of the model or adopting a new Lagrangian time scale provides stronger dissipation without significantly affecting model behavior.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Global gyrokinetic simulations of electrostatic microturbulent transport using kinetic electrons in LHD stellarator

Global gyrokinetic simulations of ion temperature gradient (ITG) and trapped electron mode (TEM) in the LHD stellarator are carried out using the gyrokinetic toroidal code (GTC) with kinetic electrons. ITG simulations show that kinetic electron effects increase the growth rate by more than 50% and more than double the turbulent transport levels compared with simulations using adiabatic electrons. Zonal flow dominates the saturation mechanism in the ITG turbulence. Nonlinear simulations of the TEM turbulence show that the main saturation mechanism is not the zonal flow but the inverse cascade of high to low toroidal harmonics. Further nonlinear simulations with various pressure profiles indicate that the ITG turbulence is more effective in driving heat conductivity whereas the TEM turbulence is more effective for particle diffusivity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

On the Two-Dimensionalization of Homogeneous Rotating Turbulence

Large-eddy simulation of the incompressible Navier-Stokes equations has been used to examine the long-time development of initially isotropic turbulence subjected to solid-body rotation. The simulations were carried out using a pseudo-spectral method with 128 x 128 x 512 collocation points in a computational domain that is four times larger along the rotation axis than in the other directions; subgrid-scale motions were parameterized using a spectral eddy viscosity model modified for system rotation. Simulation results show that the correlation length along the rotation am's of velocities orthogonal to the rotation vector exhibits rapid growth while the integral length-scale of velocities aligned with the rotation axis is relatively unaffected by rotation. Examination of the energy spectrum of two-dimensional, two-component motions indicates the presence of an inverse cascade of energy. System rotation also causes an alignment of vorticity along the rotation axis with relatively stronger cyclonic vorticity than anticyclonic. The onset of anisotropic effects are well characterized by Rossby numbers defined in terms of both macroscopic and microscopic quantities.

Squires, K. D.↗
DOE OSTI2021

Role of Parallel Solenoidal Electric Field on Energy Conversion in 2.5D Decaying Turbulence with a Guide Magnetic Field

We perform 2.5D particle-in-cell simulations of decaying turbulence in the presence of a guide (out-of-plane) background magnetic field. The fluctuating magnetic field initially consists of Fourier modes at low wavenumbers (long wavelengths). With time, the electromagnetic energy is converted to plasma kinetic energy (bulk flow+thermal energy) at the rate per unit volume of J· E for current density J and electric field E. Such decaying turbulence is well known to evolve toward a state with strongly intermittent plasma current. Here we decompose the electric field into components that are irrotational, E ir , and solenoidal (divergence-free), E so . E ir is associated with charge separation, and J · E ir is a rate of energy transfer between ions and electrons with little net change in plasma kinetic energy. Therefore, the net rate of conversion of electromagnetic energy to plasma kinetic energy is strongly dominated by J · E so , and for a strong guide magnetic field, this mainly involves the component E so,∥ parallel to the total magnetic field B. We examine various indicators of the spatial distribution of the energy transfer rate J ∥ · E so,∥ , which relates to magnetic reconnection, the best of which are (1) the ratio of the out-of-plane electric field to the in-plane magnetic field, (2) the out-of-plane component of the nonideal electric field, and (3) the magnitude of the estimate of current helicity.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1987

Turbulence statistics in fully developed channel flow at low Reynolds number

A direct numerical simulation of a turbulent channel flow is performed. The unsteady Navier-Stokes equations are solved numerically at a Reynolds number of 3300, based on the mean centerline velocity and channel half-width, with about 4 million grid points. All essential turbulence scales are resolved on the computational grid and no subgrid model is used. A large number of turbulence statistics are computed and compared with the existing experimental data at comparable Reynolds numbers. Agreements as well as discrepancies are discussed in detail. Particular attention is given to the behavior of turbulence correlations near the wall. A number of statistical correlations which are complementary to the existing experimental data are reported for the first time.

Kim, John↗
NASA NTRS1985

Heat transfer from a simulated shuttle external tank in the naturally turbulent outdoor environment

To estimate the local heat-transfer coefficients around an external-tank (ET) surface, a 9.8-m-high by 8.4-m-diameter ET simulator is constructed, instrumented, and tested in the naturally turbulent outdoor environment. Two different configurations under which the testing is conducted are described. The time-average, local values of Nusselt number at time-average Reynolds numbers of 2.2 x 10 to the 6th, 2.6 x 10 to the 6th, and 4.5 x 10 to the 6th are presented. It is shown that, at the same Reynolds numbers, the heat transfer coefficients due to the atmospheric airflow are higher than those due to airflow in a low-turbulent-intensity wind tunnel.

Lin, F. N.↗
NASA NTRS1999

Direct-Numerical and Large-Eddy Simulations of a Non-Equilibrium Turbulent Kolmogorov Flow

A non-equilibrium form of turbulent Kolmogorov flow is set up by making an instantaneous change in the amplitude of the spatially-periodic forcing. It is found that the response of the flow to this instantaneous change becomes more dramatic as the wavenumber of the forcing is increased, and, at the same time, that the faithfulness with which the large-eddy-simulation results agree with the direct-numerical results decreases.

Woodruff, S. L.↗
NASA NTRS1982

Influences of horizontal hypokinesia on performance and circadian physiological rhythms in female humans

Eight females from 35-45 yr of age were subjected to seven days of ambulatory control, seven days of bed rest, and a five day recovery period, with 30 min of centrifugation on day seven of bedrest to determine the effects of weightlessness on the circadian rhythms of females in that age group. Heart rate and rectal temperature (RT) were monitored and each subject was tested in a flight simulator twice a day in conditions of varying levels of turbulence. The flight simulations were run during the morning and acrophase of the circadian RT and performance errors wery monitored for 6 min. No significant differences were detected in the group performance data pre-, during, and post-bedrest, although better performance in the simulator was observed after the centrifuge exposure. An RT phase shift was statistically significant between pre- and during bedrest stages.

Winget, C. M.↗
NASA NTRS2002

Filtered Mass Density Function for Design Simulation of High Speed Airbreathing Propulsion Systems

The objective of this research is to improve and implement the filtered mass density function (FDF) methodology for large eddy simulation (LES) of high speed reacting turbulent flows. NASA is interested in the design of various components involved in air breathing propulsion systems such as the scramjet. There is a demand for development of robust tools that can aid in the design procedure. The physics of high speed reactive flows is rich with many complexities. LES is regarded as one of the most promising means of simulating turbulent reacting flows.

Givi, P.↗