Search NASA⌕ Search

SEARCH · Search NASA

Results for “turbulence model”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Deep modelling of plasma and neutral fluctuations from gas puff turbulence imaging

The role of turbulence in setting boundary plasma conditions is presently a key uncertainty in projecting to fusion energy reactors. To robustly diagnose edge turbulence, we develop and demonstrate a technique to translate brightness measurements of HeI line radiation into local plasma fluctuations via a novel integrated deep learning framework that combines neutral transport physics and collisional radiative theory for the $3^3 D - 2^3 P$ transition in atomic helium. The tenets for experimental validity are reviewed, illustrating that this turbulence analysis for ionized gases is transferable to both magnetized and unmagnetized environments with arbitrary geometries. Based upon fast camera data on the Alcator C-Mod tokamak, we present the first 2-dimensional time-dependent experimental measurements of the turbulent electron density, electron temperature, and neutral density revealing shadowing effects in a fusion plasma using a single spectral line.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

On analysis and stochastic modeling of the particle kinetic energy equation in particle-laden isotropic turbulent flows

In this study, we analyze three-dimensional particle-laden, isotropic turbulence to develop an understanding of inertial particle dynamics from a kinetic energy perspective. Data trends implying inhomogeneous sampling of the flow by particles are identified and used to support a proposed particle behavior: particles appear to accumulate in regions of low flow kinetic energy over time because they lose kinetic energy and slow down in such regions, ultimately causing them to spend more time there. To elucidate this behavior, we derive a particle kinetic energy equation from the particle momentum equation, which incorporates inertial effects through the Schiller–Naumann drag correlation. Upon extracting fundamental physics from this equation, hypotheses regarding the role of the Stokes number in the temporal change of particle kinetic energy and the previously proposed particle behavior are evaluated using simulation data considering three Stokes numbers. Finally, a Fokker–Planck equation is used to derive the steady-state probability density function of the particle kinetic energy. The model fits the simulation data well and provides a tool for further investigation into understanding preferential concentration, as well as a reduced order model for predicting particle kinetic energy in turbulent flows.

42 ENGINEERING↗
DOE OSTI2020

First implementation of gyrokinetic exact linearized Landau collision operator and comparison with models

Gyrokinetic simulations are fundamental to understanding and predicting turbulent transport in magnetically confined fusion plasmas. Previous simulations have used model collision operators with approximate field-particle terms of unknown accuracy and/or have neglected collisional finite Larmor radius (FLR) effects. We have implemented the linearized Fokker–Planck collision operator with exact field-particle terms and full FLR effects in a gyrokinetic code (GENE). The new operator, referred to as “exact” in this paper, allows the accuracy of model collision operators to be assessed. The conservative Landau form is implemented because its symmetry underlies the conservation laws and the H-theorem, and enables numerical methods to preserve this conservation, independent of resolution. The implementation utilizes the finite-volume method recently employed to discretize the Sugama collision model in GENE, allowing direct comparison between the two operators. Results show that the Sugama model appears accurate for the growth rates of trapped electron modes (TEMs) driven only by density gradients, but appreciably underestimates the growth rates as the collisionality and electron temperature gradient increase. The TEM turbulent fluxes near the nonlinear threshold using the exact operator are similar to the Sugama model for the n e = d ln T e /d ln n e = 0 case, but substantially larger than the Sugama model for the n e = 1 case. The FLR effects reduce the growth rates increasingly with wavenumber deepening a “valley” at the intermediate binormal wavenumber as the unstable mode extends from the TEM regime to the electron temperature gradient instability regime. Application to the Hinton–Rosenbluth problem shows that zonal flows decay faster as the radial wavenumber increases and the exact operator yields weaker decay rates.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Evaluation of RANS vs. LES simulation of fluid flow through 3 × 3 rod bundle with a simple spacer grid as a precursor to coupled fluid–structure interaction simulations

The research literature on Computational Fluid Dynamics (CFD) of coolant flow through rod bundles with spacer-grids and mixing vanes is replete, ranging from high fidelity Large Eddy Simulation (LES)/Direct Numerical Simulation (DNS) simulations to Reynolds-Averaged Navier–Stokes (RANS) modeled studies. The mixing of flow between subchannels and the pressure drop through the bundle are fundamental quantities useful for comparing and evaluating CFD methods. Less commonly observed and compared are the forces exerted onto the structure by the fluid. The present study seeks to evaluate the use of RANS simulations for predicting the structural response to fluid flow. Wall resolved RANS simulations are benchmarked against LES simulations of fluid flow at a Reynolds number of 15,000 through a 3 × 3 fuel rod bundle with a simple spacer grid. Velocity line-plots are compared showing good agreement between RANS and LES results, ascertaining that the former is capable of capturing the essential time-averaged velocity profile. Additionally, the distribution of forces on the spacer grid and fuel rods are collected as a function of time and space. The RANS methods are evaluated using the frequency and magnitude of the fluctuating forces on various portions of the structure as compared to LES. In conclusion, the power spectral density evaluation of the models reveal underprediction of force amplitude on the rod walls by RANS and also discrepancy in the prediction of high frequency spectra, especially in the immediate vicinity of spacer-grid structure, which may be attributed to the lack of random turbulence fluctuation or insufficient modeling of small-scale eddies in RANS simulation.

FIV↗
DOE OSTI2025

Learning Physically Interpretable Atmospheric Models From Data With WSINDy

The multiscale and turbulent nature of Earth's atmosphere has historically rendered accurate weather modeling a hard problem. Recently, there has been an explosion of interest surrounding data-driven approaches to weather modeling, which in many cases show improved forecasting accuracy and computational efficiency when compared to traditional methods. However, many of the current data-driven approaches employ highly parameterized neural networks, often resulting in uninterpretable models and limited gains in scientific understanding. In this work, we address the interpretability problem by explicitly discovering partial differential equations governing atmospheric phenomena, identifying symbolic mathematical models with direct physical interpretations. The purpose of this paper is to demonstrate that, in particular, the weak-form sparse identification of nonlinear dynamics (WSINDy) algorithm can learn effective atmospheric models from both simulated and assimilated data. Our approach adapts the standard WSINDy algorithm to work with high-dimensional fluid data of arbitrary spatial dimension.

58 GEOSCIENCES↗
DOE OSTI2021

A model investigation of the impact of lower hybrid wave scattering angle on current drive profile in EAST and Alcator C-Mod

Lower hybrid current drive (LHCD) is beneficial for developing a steady-state operation scenario in a tokamak. This paper conducts a modeling investigation to identify an optimum rotation angle of the initial lower hybrid perpendicular (to the background magnetic field $\overrightarrow{B}$) wavevector for best matching the experimental RF current profile. It is hypothesized that central RF power deposition widely observed in the present-day LHCD experiments arises from wave scattering by turbulence. In a standard model without considering such interactions, the predicted power deposition profile is generally broad with off-axis peaking, not in agreement with experimental observations. A heuristic approach is adopted by introducing a spectral broadening mechanism by modifying the initial orientation of the perpendicular wavevector. The ray-tracing/Fokker–Planck solver GENRAY/CQL3D is utilized within the python-based π-scope framework. A focus is given to identify the perpendicular wavenumber orientation angle with respect to the magnetic surface normal vector at the initial ray location. Our modeling study shows that rotating the perpendicular wavevector in such a way as to increase the initial poloidal component is effective in reproducing the centrally peaked current profile observed in normal shear plasmas on both EAST and C-Mod. These waves can readily be absorbed to the central plasma, which reduces the sensitivity of the power deposition profile to a slight change of the plasma condition. The same approach is also found to help broaden the off-axis power deposition profile in a reverse-shear EAST plasma, leading to a better agreement with the experiment. The results presented here suggest that spectral modification arising from edge density fluctuations in a tokamak may need to be considered in understanding wave propagation and absorption. A further experimental and theoretical/modeling study is vital as a reverse approach is adopted in this study. Finally, our work suggests that mitigation or control measures are critical for parasitic effects occurring on the first pass in a reactor regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 OSTI2020

What can the GRACE seasonal cycle tell us about lake-aquifer interactions?

Lake-aquifer interactions have been the subject of investigation and debate for decades. Traditional investigation methods include measurement of water flux across the groundwater-surface water interface, application of heat and environmental tracer methods, conducting numerical simulations of the water flow, and mass balance-based approaches. We first review, evaluate and describe the different approaches that have been applied for examining lake-aquifer interactions and propose an additional complementary approach. While the leakage from lakes and its slow movement through porous media (laminar flow) is well characterized by existing methods, modeling the rapid and turbulent flow through preferred pathways (faults/karst) remains a challenging task. To better understand the nature, and the full scale, of the lake aquifer interactions, and using Lake Nasser and the underlying Dakhla subbasin of the Nubian Sandstone Aquifer System (NSAS) in the Western Desert of Egypt (area: 0.66 x 10 6 km 2 ) as a test site, we examine, model, and correlate temporally and spatially, the variations in Gravity Recovery and Climate Experiment terrestrial water storage (GRACE TWS ) solutions, with precipitation, lake levels, area, and water volume. We review current GRACE applications in hydrology and present our novel approach. Findings include: firstly, large seasonal fluctuations (peak: Nov./Dec.; trough: July/Aug.) in Lake Nasser's surface water levels (average annual fluctuations: from 3.7 m to 7.5 m), area (average area: 3622 km 2 but up to 4530 km 2 ), and volume (annual average: 13.4 km 3 , and up to 34 km 3 ) are observed throughout years 2006 to 2015. These fluctuations are accompanied by an increase in GRACE TWS (average: 50 ± 13 mm/yr, up to 77 ± 18 mm/yr) over Lake Nasser and by a progression of a front of increasing GRACE TWS values (> 50 ± 13 mm) away from the lake reaching distances of up to 700 km some 3 to 5 months following peak lake level periods. The areas witnessing the seasonal increase in GRACE TWS display a progressive increase in phase and decrease in amplitude with distance from the lake. Secondly, the negligible precipitation over the Dakhla subbasin cannot account for the observed seasonal GRACE TWS patterns and neither can the leakage signal from Lake Nasser. Thirdly, overall similarities in the distribution of GRACE TWS seasonal spatial patterns are observed. During high lake level periods (e.g., Dec. 2007; Feb. 2015) additional, and more distant, areas from Lake Nasser saw increases in GRACE TWS and vice versa during low lake level periods (e.g., July 2006; Aug. 2010). These observations are consistent with Lake Nasser being the main source of modern recharge for the Dakhla subbasin and suggest a new conceptual model for the subbasin: a slow groundwater flow through a porous matrix and a rapid groundwater flow along a network of faults, fractures, and karst topography across the matrix. We suggest that the proposed conceptual model for the interaction between Lake Nasser and the Dakhla subbasin could be applicable to aquifers of similar geologic, climatic, and hydrologic settings worldwide and that approaches similar to those advocated here could be used to investigate the validity of this suggestion.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Ultrasonic oscillatory two-phase flow in microchannels

Experimental and numerical investigations are performed to provide an assessment of the transport behavior of an ultrasonic oscillatory two-phase flow in a microchannel. The work is inspired by the flow observed in an innovative ultrasonic fabric drying device using a piezoelectric bimorph transducer with microchannels, where a water-air two-phase flow is transported by harmonically oscillating microchannels. The flow exhibits highly unsteady behavior as the water and air interact with each other during the vibration cycles, making it significantly different from the well-studied steady flow in microchannels. Here, the computational fluid dynamics (CFD) modeling is realized by combing the turbulence Reynolds-averaged Navier-Stokes (RANS) k – ω model with the phase-field method to resolve the dynamics of the two-phase flow. The numerical results are qualitatively validated by the experiment. Through parametric studies, we specifically examined the effects of vibration conditions (i.e., frequency and amplitude), microchannel taper angle, and wall surface contact angle (i.e., wettability) on the flow rate through the microchannel. The results will advance the potential applications where oscillatory or general unsteady microchannel two-phase flows may be present.

42 ENGINEERING↗
DOE OSTI2023

The myth of URANS

Since the 1990s, RANS practitioners have observed spontaneous unsteadiness in RANS simulations. Some have suggested deliberately using this as a method of resolving large turbulent structures. However, to date, no one has produced a theoretical justification for this unsteady RANS (URANS) approach. Here, we extend the dynamical system fixed point analysis to create a theoretical model for URANS dynamics. The results are compared to URANS simulations for homogeneous isotropic decaying turbulence. The model shows that URANS can predict incorrect decay rates and that the solution tends towards steady RANS over time. Similar analysis for forced turbulence shows a fixed modelled energy of about 30% of total energy, regardless of the model parameters. The same analysis can be used to show how hybrid type models can begin to address these issues.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Large Eddy Simulations of Turbulent and Buoyant Flows in Urban and Complex Terrain Areas Using the Aeolus Model

Fast and accurate predictions of the flow and transport of materials in urban and complex terrain areas are challenging because of the heterogeneity of buildings and land features of different shapes and sizes connected by canyons and channels, which results in complex patterns of turbulence that can enhance material concentrations in certain regions. To address this challenge, we have developed an efficient three-dimensional computational fluid dynamics (CFD) code called Aeolus that is based on first principles for predicting transport and dispersion of materials in complex terrain and urban areas. The model can be run in a very efficient Reynolds average Navier–Stokes (RANS) mode or a detailed large eddy simulation (LES) mode. The RANS version of Aeolus was previously validated against field data for tracer gas and radiological dispersal releases. As a part of this work, we have validated the Aeolus model in LES mode against two different sets of data: (1) turbulence quantities measured in complex terrain at Askervein Hill; and (2) wind and tracer data from the Joint Urban 2003 field campaign for urban topography. As a third set-up, we have applied Aeolus to simulate cloud rise dynamics for buoyant plumes from high-temperature explosions. For all three cases, Aeolus LES predictions compare well to observations and other models. These results indicate that Aeolus LES can be used to accurately simulate turbulent flow and transport for a wide range of applications and scales.

42 ENGINEERING↗
DOE OSTI2024

Clumpy structures within the turbulent primordial cloud

ABSTRACT The primordial clouds in the mini-haloes hatch the first generation stars of the Universe, which play a crucial role in cosmic evolution. In this paper, we investigate how turbulence impacts the structure of primordial star-forming clouds. Previous cosmological simulations of the first star formation predicted a typical mass of around $\mathrm{ 100 \, M_\odot }$. This conflicts with recent observations of extremely metal-poor stars, suggesting a lower mass scale of about $\mathrm{25 \, M_\odot }$. The discrepancy may arise from unresolved turbulence in the star-forming cloud, driven by primordial gas accretion during mini-halo formation in the previous simulations. To quantitatively examine the turbulence effect on the primordial cloud formation, we employ the adaptive mesh refinement code Enzo to model the gas cloud with primordial composition, including artificially driven turbulence on the cloud scale and relevant gas physics. This artificially driven turbulence utilizes a stochastic forcing model to mimic the unresolved turbulence inside mini-haloes. Our results show that the turbulence with high Mach number and compressional mode effectively fragments the cloud into several clumps, each with dense cores of $\mathrm{22.7 - 174.9 \, M_\odot }$ that undergo Jeans instability to form stars. Fragmentation caused by intense and compressive turbulence prevents a runaway collapse of the cloud. The self-bound clumps with smaller masses in the turbulent primordial clouds suggest a possible pathway to decrease the theoretical mass scale of the first stars, further reconciling the mass discrepancy between simulations and observations.

Astronomy & Astrophysics↗
DOE OSTI2023

Quantitative comparison of impurity transport in turbulence reduced and enhanced scenarios at Wendelstein 7-X

We assess the turbulent particle transport being responsible for the limitation of the confinement and, thus, the overall performance of the neoclassically optimized stellarator Wendelstein 7-X. The radial particle transport is experimentally inferred from the evaluation of impurity injection into turbulence reduced and enhanced plasma scenarios revealing a completely different confinement behavior. The impact of the density gradient on the turbulent ion transport is theoretically estimated using large-scale non-linear gyro-kinetic simulations enabling, for the first time in Wendelstein 7-X, a quantitative comparison to the experimentally assessed impurity transport properties. We demonstrate that impurity transport in most of the Wendelstein 7-X discharges, up to now impossible to cover only with neoclassical estimations, is dominated by turbulence and can be modelled via gyro-kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Collisionality driven turbulent particle transport changes in DIII-D H-mode plasmas

The results of the experimental dimensionless scan in this paper confirm that there is an increase in density peaking towards lower collisionality and that this can be partly linked to a shift in the turbulence regime from ITG towards TEM. However at the lowest collisionality, the changes in turbulence and transport are much more pronounced than expected from direct collisionality effect on the turbulence. In this paper, the collisionality, ν * is varied by a factor 5, while keeping ρ * , q , β , M , fixed. Additionally, a 3 Hz gas puff modulation is applied to modulate the electron density profile and extract the perturbed transport coefficients using two diagnostics. The transport analysis shows that the increase in density peaking at low ν * is linked to an increase in the inward particle pinch and not an increase in core fueling. These observations are not only in agreement with prior modeling scans of how turbulence changes as a function of collisionality and its impact upon the particle fluxes, but also with the multi-machine database (Fable E. et al 2010 Plasma Phys. Control. Fusion 52 015007) (Angioni C. et al 2003 Phys. Rev. Lett. 90 205003). The changes in turbulence across the collisionality scan were captured at large scale by the BES and at smaller scale by the DBS. A comparison with gradient-driven GENE simulations showed similar trends at both scales. Moreover, the changes observed in overall transport are in agreement with gradient-driven TGLF particle flux simulations. This indicates that TGLF/GENE when given the gradients as input, are able to reproduce the experimentally observed turbulence changes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Turbulent entrainment in finite-length wind farms

In this article, we present an entrainment-based model for predicting the flow and power output of finite-length wind farms. The model is an extension of the three-layer approach of Luzzatto-Fegiz & Caulfield ( Phys. Rev. Fluids , vol. 3, 2018, 093802) for wind farms of infinite length, and assumes dependence of key flow quantities, such as the wind farm bulk velocity, on the streamwise distance from the farm entrance. To assist our analysis and validate the proposed model, we undertake a series of large-eddy simulations with different turbine spacing arrangements and layouts. Comparisons are also made with the top-down model with entrance effects of Meneveau ( J. Turbul. , vol. 13, 2012, N7) and data from the literature. The finite-length entrainment model is shown to be capable of capturing the power drop between contiguous rows of turbines as well as describing the advection and turbulent transport of kinetic energy in both the entrance and fully developed regions. The fully developed regime is approximated only deep in the wind farm, after approximately 15 rows of turbines. Our data suggest that for the cases considered in this study, the empirical coefficients that can be used to describe turbulent entrainment and transfers above the wind farm exhibit little dependence on the farm layout and may be considered constant for modelling purposes. However, the flow field within the wind farm layer can be strongly modulated by the turbine density (spacing) as well as the array layout, and to that extent it can be argued that they are both primary factors determining the wind farm power output.

17 WIND ENERGY↗
DOE OSTI2022

Implications of inertial subrange scaling for stably stratified mixing

In this work, we investigate the effects of the turbulent dynamic range on active scalar mixing in stably stratified turbulence by adapting the theoretical passive scalar modelling and demonstrating their usefulness through consideration of the results of direct numerical simulations of statistically stationary homogeneous stratified and sheared turbulence. By analysis of inertial and inertial–convective subrange scalings, we show that the relationship between the active scalar and turbulence time scales is predicted by the ratio of the Kolmogorov and Obukhov–Corrsin constants, provided mean flow parameters permit the two subrange scalings to be appropriate approximations. We use the resulting relationship between time scales to parameterise an appropriate turbulent mixing coefficient $\varGamma \equiv \chi /\epsilon$, defined here as the ratio of available potential energy ($E_p$) and turbulent kinetic energy ($E_k$) dissipation rates. With the analysis presented here, we show that $\varGamma$ can be estimated by $E_p,E_k$ and a universal constant provided an appropriate Reynolds number is sufficiently high. This large Reynolds number regime appears here to occur at $ {{Re_b}} \equiv \epsilon / \nu N^{2} \gtrapprox 300$ where $\nu$ is the kinematic viscosity and $N$ is the characteristic buoyancy frequency. We propose a model framework for irreversible diapycnal mixing with robust theoretical parametrisation and asymptotic behaviour in this high-$ {{Re_b}}$ limit.

42 ENGINEERING↗