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.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Development of MHD for NekRS

A spectral-element-based formulation of incompressible MHD is presented in the context of the open-source fluid-thermal code, Nek5000/RS. The formulation supports magnetic fields in a solid domain that surrounds the fluid domain. Several steady-state and time-transient model problems are presented as part of the code verification process. Nek5000/RS is designed for large-scale turbulence simulations, which will be the next step with this new MHD capability.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Magnetic Field Strength from Turbulence Theory. I. Using Differential Measure Approach

The mean plane-of-sky magnetic field strength is traditionally obtained from the combination of polarization and spectroscopic data using the Davis–Chandrasekhar–Fermi (DCF) technique. However, we identify the major problem of the DCF technique to be its disregard of the anisotropic character of MHD turbulence. On the basis of the modern MHD turbulence theory we introduce a new way of obtaining magnetic field strength from observations. Unlike the DCF technique, the new technique uses not the dispersion of the polarization angle and line-of-sight velocities, but increments of these quantities given by the structure functions. To address the variety of astrophysical conditions for which our technique can be applied, we consider turbulence in both media with magnetic pressure higher than the gas pressure, corresponding, e.g., to molecular clouds, and media with gas pressure higher than the magnetic pressure, corresponding to the warm neutral medium. We provide general expressions for arbitrary admixtures of Alfvén, slow, and fast modes in these media and consider in detail particular cases relevant to diffuse media and molecular clouds. We successfully test our results using synthetic observations obtained from MHD turbulence simulations. We demonstrate that our differential measure approach, unlike the DCF technique, can be used to measure the distribution of magnetic field strengths, can provide magnetic field measurements with limited data, and is much more stable in the presence of induced large-scale variations of nonturbulent nature. Furthermore, our study uncovers the deficiencies of earlier DCF research.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

JET-ILW Nonlinear Pedestal ETG Data

Contains data for paper on nonlinear pedestal electron-temperature-gradient turbulence simulations: *.py # python plotting scripts for generating paper plots. *.txt, *.npy # post-processed data required for plotting scripts. *.png, *.eps, *.pdf # paper plots. *.in # gs2 and stella input files for simulations in this paper. stella used for nonlinear simulations (https://github.com/mabarnes/stella)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Plasma edge and scrape-off layer turbulence in gyrokinetic simulations of negative triangularity plasmas

Gyrokinetic simulations in the long-wavelength or drift-kinetic limit are carried out of DIII-D inner-wall-limited (IWL) plasmas to investigate the effect of triangularity on edge and scrape-off layer (SOL) turbulence. The effect of neutral interactions and triangularity on plasma blobs is explored due to the impact blobs can have in setting the SOL width or introducing impurities through interactions with plasma-facing components. Seeded blob simulations with neutrals in shaped SOL scenarios demonstrate that increasing elongation, triangularity, or Shafranov shift decreases radial blob velocities, but neutral interactions have a minor effect. Fully turbulent simulations of DIII-D IWL plasmas include both open- and closed-field-line regions. The negative triangularity (NT) simulation has lower average core T e , lower normalized T e fluctuations, and lower fluxes, but a greater number of coherent structures (blobs) identified with increased size and velocity, on average. Density and electron temperature profiles are within a factor of 2 of experimental values. In conclusion, the increased trapped electron particle fraction in NT simulations is consistent with previous studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Improved Representation of Horizontal Variability and Turbulence in Mesoscale Simulations of an Extended Cold-Air Pool Event

Abstract Cold-air pools (CAPs), or stable atmospheric boundary layers that form within topographic basins, are associated with poor air quality, hazardous weather, and low wind energy output. Accurate prediction of CAP dynamics presents a challenge for mesoscale forecast models in part because CAPs occur in regions of complex terrain, where traditional turbulence parameterizations may not be appropriate. This study examines the effects of the planetary boundary layer (PBL) scheme and horizontal diffusion treatment on CAP prediction in the Weather Research and Forecasting (WRF) Model. Model runs with a one-dimensional (1D) PBL scheme and Smagorinsky-like horizontal diffusion are compared with runs that use a new three-dimensional (3D) PBL scheme to calculate turbulent fluxes. Simulations are completed in a nested configuration with 3-km/750-m horizontal grid spacing over a 10-day case study in the Columbia River basin, and results are compared with observations from the Second Wind Forecast Improvement Project. Using event-averaged error metrics, potential temperature and wind speed errors are shown to decrease both with increased horizontal grid resolution and with improved treatment of horizontal diffusion over steep terrain. The 3D PBL scheme further reduces errors relative to a standard 1D PBL approach. Error reduction is accentuated during CAP erosion, when turbulent mixing plays a more dominant role in the dynamics. Last, the 3D PBL scheme is shown to reduce near-surface overestimates of turbulence kinetic energy during the CAP event. The sensitivity of turbulence predictions to the master length-scale formulation in the 3D PBL parameterization is also explored. Significance Statement In this article, we demonstrate how a new framework for modeling atmospheric turbulence improves cold pool predictions, using a case study from January 2017 in the Columbia River basin (U.S. Pacific Northwest). Cold pools are regions of cold, stagnant air that form within valleys or basins, and improved forecasts could help to mitigate the risks they pose to air quality, transportation, and wind energy production. For the chosen case study, our tests show a reduction in temperature and wind speed errors by up to a factor of 2–3 relative to standard model options. These results strongly motivate continued development of the framework as well as its application to other complex weather events.

17 WIND ENERGY↗

DOE OSTI2023

Grid resolution requirement of chemical explosive mode analysis for large eddy simulations of premixed turbulent combustion

Full Article Figures & data References Citations Metrics Reprints & Permissions Read this article Abstract The grid resolution requirement for trustworthy Chemical Explosive Mode Analysis (CEMA) in Large Eddy Simulation (LES) of premixed turbulent combustion is proposed. Explicit filtering, to emulate the effect of the LES filter, is applied to one-dimensional laminar flame and three-dimensional planar turbulent flames across a wide range of Karlovitz numbers (5 - 239). The identification of the flame front by CEMA is found relatively insensitive to the cell size (Δ), while the combustion mode identification shows more significant sensitivity. Specifically, increasing Δ falsely enhances the auto-ignition and local extinction modes and suppresses the diffusion-assisted mode. Limited dependence of the CEMA performance on the turbulent combustion regime (Karlovitz number) is observed. A simple grid size criterion for reliable CEMA mode identification in LES is proposed as Δ ≲ δ L /2; The criterion can be relaxed to Δ ≲ δ L in the laminar flame limit. Furthermore, theoretical analysis is conducted on an idealised chemistry-diffusion system. The effects of the filtering process and turbulence on the local combustion mode are demonstrated, which is consistent with the numerical observations. Further, by incorporating turbulent combustion models in CEMA, potential improvement in identifying local combustion modes can be expected.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Direct numerical simulation of hypersonic turbulent boundary layers: effect of spatial evolution and Reynolds number

Direct numerical simulations (DNS) are performed to investigate the spatial evolution of flat-plate zero-pressure-gradient turbulent boundary layers over long streamwise domains ( >300δ i , with δ i the inflow boundary-layer thickness) at three different Mach numbers, 2.5, 4.9 and 10.9, with the surface temperatures ranging from quasiadiabatic to highly cooled conditions. The settlement of turbulence statistics into a fully developed equilibrium state of the turbulent boundary layer has been carefully monitored, either based on the satisfaction of the von Kármán integral equation or by comparing runs with different inflow turbulence generation techniques. The generated DNS database is used to characterize the streamwise evolution of multiple important variables in the high-Mach-number, cold-wall regime, including the skin friction, the Reynolds analogy factor, the shape factor, the Reynolds stresses, and the fluctuating wall quantities. The data confirm the validity of many classic and newer compressibility transformations at moderately high Reynolds numbers (up to friction Reynolds number Re τ ≈1200 ) and show that, with proper scaling, the sizes of the near-wall streaks and superstructures are insensitive to the Mach number and wall cooling conditions. The strong wall cooling in the hypersonic cold-wall case is found to cause a significant increase in the size of the near-wall turbulence eddies (relative to the boundary-layer thickness), which leads to a reduced-scale separation between the large and small turbulence scales, and in turn to a lack of an outer peak in the spanwise spectra of the streamwise velocity in the logarithmic region.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Validation of turbulence intensity as simulated by the Weather Research and Forecasting model off the US northeast coast

We report turbulence intensity (TI) is often used to quantify the strength of turbulence in wind energy applications and serves as the basis of standards in wind turbine design. Thus, accurately characterizing the spatiotemporal variability in TI should lead to improved predictions of power production. Nevertheless, turbulence measurements over the ocean are far less prevalent than over land due to challenges in instrumental deployment, maintenance, and operation. Atmospheric models such as mesoscale (weather prediction) and large-eddy simulation (LES) models are commonly used in the wind energy industry to assess the spatial variability of a given site. However, the TI derivation from atmospheric models has not been well examined. An algorithm is proposed in this study to realize online calculation of TI in the Weather Research and Forecasting (WRF) model. Simulated TI is divided into two components depending on scale, including sub-grid (parameterized based on turbulence kinetic energy (TKE)) and grid resolved. The sensitivity of sea surface temperature (SST) on simulated TI is also tested. An assessment is performed by using observations collected during a field campaign conducted from February to June 2020 near the Woods Hole Oceanographic Institution Martha's Vineyard Coastal Observatory. Results show that while simulated TKE is generally smaller than the lidar-observed value, wind speed bias is usually small. Overall, this leads to a slight underestimation in sub-grid-scale estimated TI. Improved SST representation subsequently reduces model biases in atmospheric stability as well as wind speed and sub-grid TI near the hub height. Large TI events in conjunction with mesoscale weather systems observed during the studied period pose a challenge to accurately estimating TI from models. Due to notable uncertainty in accurately simulating those events, this suggests summing up sub-grid and resolved TI may not be an ideal solution. Efforts in further improving skills in simulating mesoscale flow and cloud systems are necessary as the next steps.

17 WIND ENERGY↗

DOE OSTI2025

Gyrokinetic Particle Simulations of Core Turbulence in Fusion Plasmas

In this project, we have further developed GTC gyrokinetic particle simulations to study plasma confinement properties in the stellarators and tokamaks with 3D magnetic fields, where both neoclassical and turbulent transport are self-consistently treated on an equal footing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

On the convergence of statistics in simulations of stationary incompressible turbulent flows

When reporting statistics from simulations of statistically stationary chaotic phenomenon, it is important to verify that the simulations are time-converged. This condition is connected with the statistical error or number of digits with which statistics can be reliably reported. In this work we consider numerical experiments of low Reynolds number incompressible homogeneous and isotropic turbulence as a model problem to investigate statistical convergence over finite simulation times. Specifically, we investigate the time integration requirements that allow meaningful reporting of the statistical error associated with finiteness of the temporal domain. We address two key questions: (1) How long should a simulation be performed in terms of large eddy time, and (2) How should the simulation time be divided among temporal windows over which a quantity of interest is estimated so that its statistical error could be reliably reported? We find that reliable reporting of statistical errors requires simulations on the order of 10 4 large eddy times, which is orders of magnitude longer than typically performed. Additionally, data post-processing should employ windows of at least ten times the large eddy time scale, with the most robust computation of statistical error of the mean requiring window sizes of an additional factor of ten. For practical simulations, we demonstrate that it is possible to estimate the statistical error within a factor of two under a less stringent condition in which a minimum of four windows with size at least ten large eddy times are used. In conclusion, our observations for homogeneous isotropic turbulence are also shown to hold in turbulent channel flow.

42 ENGINEERING↗

DOE OSTI2024

Scale Separation Effects on Simulations of Plasma Turbulence

Abstract Understanding plasma turbulence requires a synthesis of experiments, observations, theory, and simulations. In the case of kinetic plasmas such as the solar wind, the lack of collisions renders the fluid closures such as viscosity meaningless and one needs to resort to higher-order fluid models or kinetic models. Typically, the computational expense in such models is managed by simulating artificial values of certain parameters such as the ratio of the Alfvén speed to the speed of light ( v A / c ) or the relative mass ratio of ions and electrons ( m i / m e ). Although, typically care is taken to use values as close as possible to realistic values within the computational constraints, these artificial values could potentially introduce unphysical effects. These unphysical effects could be significant at sub-ion scales, where kinetic effects are the most important. In this paper, we use the 10-moment fluid model in the Gkeyll framework to perform controlled numerical experiments, systematically varying the ion–electron mass ratio from a small value down to the realistic proton–electron mass ratio. We show that the unphysical mass ratio has a significant effect on the kinetic range dynamics as well as the heating of both plasma species. The dissipative process for both ions and electrons becomes more compressive in nature, although the ions remain nearly incompressible in all cases. The electrons move from being dominated by incompressive viscous-like heating/dissipation to very compressive heating/dissipation dominated by compressions/rarefactions. While the heating change is significant for the electrons, a mass ratio of m i / m e ∼ 250 captures the asymptotic behavior of electron heating.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Experiments and Simulations on the Turbulent, Rarefaction Wave Driven Rayleigh–Taylor Instability

Experiments were performed to observe the growth of the turbulent, Rayleigh–Taylor unstable mixing layer generated between air and SF 6 , with an Atwood number of A=(ρ 2 –ρ 1 )/(ρ 2 +ρ 1 ) = 0.64, where ρ 1 and ρ 2 are the densities of air and SF 6 , respectively. A nonconstant acceleration with an average value of 2300g 0 , where g 0 is the acceleration due to gravity, was generated by interaction of the interface between the two gases with a rarefaction wave. Three-dimensional, multimode perturbations were generated on the diffuse interface, with a diffusion layer thickness of δ = 3.6 mm, using a membraneless vertical oscillation technique, and 20 experiments were performed to establish a statistical ensemble. The average perturbation from this ensemble was extracted and used as input for a numerical simulation using the Lawrence Livermore National Laboratory (LLNL) Miranda code. Good qualitative agreement between the experiment and simulation was observed, while quantitative agreement was best at early to intermediate times. Several methods were used to extract the turbulent growth constant α from experiments and simulations while accounting for time varying acceleration. Experimental, average bubble and spike asymptotic self-similar growth rate values range from α = 0.022 to α = 0.032 depending on the method used, and accounting for variable acceleration. Values found from the simulations range from α = 0.024 to α = 0.041. Here, the values of α measured in the experiments are lower than what are typically measured in the literature but are more in line with those found in recent simulations.

42 ENGINEERING↗

DOE OSTI2022

Nonresonant particle acceleration in strong turbulence: Comparison to kinetic and MHD simulations

Collisionless, magnetized turbulence offers a promising framework for the generation of nonthermal high-energy particles in various astrophysical sites. Yet, the detailed mechanism that governs particle acceleration has remained subject to debate. By means of 2D and 3D particle-in-cell, as well as 3D (incompressible) magnetohydrodynamic (MHD) simulations, we test here a recent model of nonresonant particle acceleration in strongly magnetized turbulence, which ascribes the energization of particles to their continuous interaction with the random velocity flow of the turbulence, in the spirit of the original Fermi model. To do so, we compare, for a large number of particles that were tracked in the simulations, the predicted and the observed histories of particles momenta. The predicted history is that derived from the model, after extracting from the simulations, at each point along the particle trajectory, the three force terms that control acceleration: the acceleration of the field line velocity projected along the field line direction, its shear projected along the same direction, and its transverse compressive part. Overall, we find a clear correlation between the model predictions and the numerical experiments, indicating that this nonresonant model can successfully account for the bulk of particle energization through Fermi-type processes in strongly magnetized turbulence. Additionally we also observe that the parallel shear contribution tends to dominate the physics of energization in the particle-in-cell simulations, while in the magnetohydrodynamic incompressible simulation, both the parallel shear and the transverse compressive term provide about equal contributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Direct numerical simulation of a turbulent methane/air flame impinged by a sub-breakdown electric field

Here, the effects of an external electric field on a turbulent methane/air diffusion flame are analyzed in this work using direct numerical simulations. The analyzed configuration consists of a temporally evolving mixing layer of air and a mixture of methane and nitrogen at 1 atm that is impinged by an electric field in the direction normal to the mean mixing plane. The combustion and chemi-ionization reactions involved in the flow are modeled using finite-rate chemistry and a reduced reaction scheme consisting of 26 species and 134 reactions. The mass diffusion and ion-wind effects are modeled using a detailed description of the diffusion coefficients and electric mobilities based on kinetic theory. The presented calculations show that the turbulence generated within the mixing layer is mostly unaffected by the applied electric field for the configuration under exam. In fact, the electric body force is developed away from the mixing region where the flow is uniform. Conversely, the turbulence is able to introduce very high intermittency in the electrically charged species concentration and, consequently, in the electric body force. Such intermittency will constitute a challenge in the future formulation of combustion models that take into account ion-wind effects.

09 BIOMASS FUELS↗

DOE OSTI2023

RANS Simulation of Variable Density Turbulent Round Jets with Coflow using xRAGE Hydrodynamic Code and BHR Turbulence Models

This work for the fiscal year 2023 (FY23) is a continuation of previous efforts to evaluate the BHR turbulence models for their ability to accurately simulate variable density turbulent round jets with coflow. As before, RANS simulations are carried out using the xRAGE hydrodynamic code. The following are some of the previous findings. Israel showed that i) the symmetry boundary conditions for the BHR models in axisymmetric simulations were in error, ii) three grids of different resolutions did not lead to converging solutions, and iii) BHR 3.1 simulation exhibited instabilities and did not reach a steady state. Saenz and Rauenzahn derived and implemented into xRAGE the correct BHR boundary conditions at the symmetry axis. Cline conducted sensitivity studies with various parameters including the BHR models (versions 2, 3.1 and 4), gravity, material pressure, specific heat, initial turbulent kinetic energy and initial turbulent length scale, and found that the largest factor impacting on simulation results was the BHR model version, followed by the initial turbulent length scale. In addition, freeze boundary conditions at the exit and the side wall of the computational domain were used to remove anomalous flow behavior. Cline adjusted the inlet jet velocity, initial turbulent kinetic energy and initial turbulent length scale to obtain the best reasonable match with the experimental data of Charonko and Prestridge. It was found that the BHR 2 and 3.1 models performed in a similar manner, but the BHR 2 model produced much lower levels of density-specific-volume covariance and turbulent kinetic energy. The main focus for the FY23 is to study the effects of computational parameters related to the boundary conditions, mesh, domain size and timestep size. The reasoning behind this is that, unless simulation results can be shown to be reasonably independent from the aforementioned computational parameters, it would be difficult to attribute any discrepancies between simulation and experimental results to turbulence models. This important aspect has largely been overlooked in the previous years, and therefore will be studied comprehensively here. Additionally, effects of varying the initial turbulent length scale will be examined because it was previously identified as a major factor affecting the flow fields.

42 ENGINEERING↗

DOE OSTI2022

Fully Kinetic Shearing-box Simulations of Magnetorotational Turbulence in 2D and 3D. I. Pair Plasmas

The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e., kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic particle-in-cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from 3D MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e., MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron–ion simulations, on which we will focus in future work.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Gyrokinetic simulations of core turbulence and thermal transport in the high-β P discharge on EAST

The properties of core turbulence and thermal transport are investigated for EAST high-β P (β P ~ 3.1) plasmas with dominant electron heating (Te/Ti > 1) via gyrokinetic simulation with the NLT code. Linear simulations identify that the electrostatic η e -driven trapped electron mode (η e -TEM) dominates in the core region (ρ< 0.7) and the ion temperature gradient (ITG) mode dominates in the out region (ρ≥ 0.7), consistent with the linear threshold analysis of micro-instabilities. Sensitivity analysis shows that the normalized electron density gradient (R/L ne ) and ITG (R/L Ti ) are two effective parameters to stabilize TEM instability. Nonlinear simulations are also carried out and compared with the experimental results, which show that the electron thermal internal transport barrier (ITB) in EAST high-β P plasma is determined by the TEM-induced turbulence. A higher zonal flow shearing rate is observed in the ITB region (0.2 < ρ < 0.34), which can regulate energy transport induced by TEM turbulence and facilitate the formation of e-ITB. A plausible positive feedback mechanism to mitigate turbulence transport and improve energy confinement via enhanced ion heating is proposed for future experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Validation of edge turbulence codes against the TCV-X21 diverted L-mode reference case

Self-consistent full-size turbulent-transport simulations of the divertor and scrape-off-layer (SOL) of existing tokamaks have recently become feasible. This enables the direct comparison of turbulence simulations against experimental measurements. In this work, we perform a series of diverted ohmic L-mode discharges on the tokamak à configuration variable (TCV) tokamak, building a first-of-a-kind dataset for the validation of edge turbulence models. This dataset, referred to as TCV-X21, contains measurements from five diagnostic systems from the outboard midplane (OMP) to the divertor targets—giving a total of 45 one- and two-dimensional comparison observables in two toroidal magnetic field directions. The experimental dataset is used to validate three flux-driven 3D fluid-turbulence models—GBS, GRILLIX and TOKAM3X. With each model, we perform simulations of the TCV-X21 scenario, individually tuning the particle and power source rates to achieve a reasonable match of the upstream separatrix value of density and electron temperature. We find that the simulations match the experimental profiles for most observables at the OMP—both in terms of profile shape and absolute magnitude—while a comparatively poorer agreement is found towards the divertor targets. The match between simulation and experiment is seen to be sensitive to the value of the resistivity, the heat conductivities, the power injection rate and the choice of sheath boundary conditions. Additionally, despite targeting a sheath-limited regime, the discrepancy between simulations and experiment also suggests that the neutral dynamics should be included. The results of this validation show that turbulence models are able to perform simulations of existing devices and achieve reasonable agreement with experimental measurements. Where disagreement is found, the validation helps to identify how the models can be improved. By publicly releasing the experimental dataset and validation analysis, this work should help to guide and accelerate the development of predictive turbulence simulations of the edge and SOL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗