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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Pressure–strain interaction. III. Particle-in-cell simulations of magnetic reconnection

How energy is converted into thermal energy in weakly collisional and collisionless plasma processes, such as magnetic reconnection and plasma turbulence, has recently been the subject of intense scrutiny. The pressure–strain interaction has emerged as an important piece, as it describes the rate of conversion between bulk flow and thermal energy density. In two companion studies, we presented an alternate decomposition of the pressure–strain interaction to isolate the effects of converging/diverging flow and flow shear instead of compressible and incompressible flow, and we derived the pressure–strain interaction in magnetic field-aligned coordinates. Here, we use these results to study pressure–strain interaction during two-dimensional anti-parallel magnetic reconnection. We perform particle-in-cell simulations and plot the decompositions in both Cartesian and magnetic field-aligned coordinates. We identify the mechanisms contributing to positive and negative pressure–strain interaction during reconnection. Furthermore, this study provides a roadmap for interpreting numerical and observational data of the pressure–strain interaction, which should be important for studies of reconnection, turbulence, and collisionless shocks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Langmuir turbulence in suspended kelp farms

This study investigates the influence of suspended kelp farms on ocean mixed layer hydrodynamics in the presence of currents and waves. We use the large eddy simulation method, where the wave effect is incorporated by solving the wave-averaged equations. Distinct Langmuir circulation patterns are generated within various suspended farm configurations, including horizontally uniform kelp blocks and spaced kelp rows. Intensified turbulence arises from the farm-generated Langmuir circulation, as opposed to the standard Langmuir turbulence observed without a farm. The creation of Langmuir circulation within the farm is attributed to two primary factors depending on farm configuration: (i) enhanced vertical shear due to kelp frond area density variability, and (ii) enhanced lateral shear due to canopy discontinuity at lateral edges of spaced rows. Both enhanced vertical and lateral shear of streamwise velocity, representing the lateral and vertical vorticity components, respectively, can be tilted into downstream vorticity to create Langmuir circulation. This vorticity tilting is driven by the Craik–Leibovich vortex force associated with the Stokes drift of surface gravity waves. In addition to the farm-generated Langmuir turbulence, canopy shear layer turbulence is created at the farm bottom edge due to drag discontinuity. The intensity of different types of turbulence depends on both kelp frond area density and the geometric configuration of the farm. The farm-generated turbulence has substantial consequences for nutrient supply and kelp growth. These findings also underscore the significance of the presence of obstacle structures in modifying ocean mixed layer characteristics.

Bo, Tong (ORCID:0000000260300561)↗
DOE OSTI2025

Incorporating long-range dependence and fractal features in turbulence spectra

We introduce an advanced turbulence spectrum model developed from mathematical foundations from a covariance function class and empirically validated using extensive field data. This model captures the complex dynamics of long-range dependence, and fractal characteristics prevalent in riverine and atmospheric boundary layer (ABL) flows that are ignored by classical spectrum models, such as IEC (International Electrotechnical Commission) von Kármán and Kaimal model. The model delineates scaling behaviors across distinct frequency bands and offers substantial flexibility through five well-defined parameters each characterizing a distinct physical aspect of the velocity time series. A detailed procedure for obtaining each parameter from time series data is outlined. The comprehensive validations with field data from tidal currents and ABL flows substantiate the model’s fidelity in accurately replicating observed phenomena. This validation establishes the reliability of the proposed model and, when incorporated into stochastic full-field simulators such as TurbSim, demonstrates its potential to advance the predictive modeling and analysis of turbulent flows in environmental science and engineering contexts.

Cheng, Shyuan [Univ. of Illinois at Urbana-Champai↗
DOE OSTI2022

Scientific challenges to characterizing the wind resource in the marine atmospheric boundary layer

Abstract. With the increasing level of offshore wind energy investment, it is correspondingly important to be able to accurately characterize the wind resource in terms of energy potential as well as operating conditions affecting wind plant performance, maintenance, and lifespan. Accurate resource assessment at a particular site supports investment decisions. Following construction, accurate wind forecasts are needed to support efficient power markets and integration of wind power with the electrical grid. To optimize the design of wind turbines, it is necessary to accurately describe the environmental characteristics, such as precipitation and waves, that erode turbine surfaces and generate structural loads as a complicated response to the combined impact of shear, atmospheric turbulence, and wave stresses. Despite recent considerable progress both in improvements to numerical weather prediction models and in coupling these models to turbulent flows within wind plants, major challenges remain, especially in the offshore environment. Accurately simulating the interactions among winds, waves, wakes, and their structural interactions with offshore wind turbines requires accounting for spatial (and associated temporal) scales from O(1 m) to O(100 km). Computing capabilities for the foreseeable future will not be able to resolve all of these scales simultaneously, necessitating continuing improvement in subgrid-scale parameterizations within highly nonlinear models. In addition, observations to constrain and validate these models, especially in the rotor-swept area of turbines over the ocean, remains largely absent. Thus, gaining sufficient understanding of the physics of atmospheric flow within and around wind plants remains one of the grand challenges of wind energy, particularly in the offshore environment. This paper provides a review of prominent scientific challenges to characterizing the offshore wind resource using as examples phenomena that occur in the rapidly developing wind energy areas off the United States. Such phenomena include horizontal temperature gradients that lead to strong vertical stratification; consequent features such as low-level jets and internal boundary layers; highly nonstationary conditions, which occur with both extratropical storms (e.g., nor'easters) and tropical storms; air–sea interaction, including deformation of conventional wind profiles by the wave boundary layer; and precipitation with its contributions to leading-edge erosion of wind turbine blades. The paper also describes the current state of modeling and observations in the marine atmospheric boundary layer and provides specific recommendations for filling key current knowledge gaps.

17 WIND ENERGY↗
DOE OSTI2024

A two-dimensional numerical study of ion-acoustic turbulence

We investigate the linear and nonlinear evolution of the current-driven ion-acoustic instability in a collisionless plasma via two-dimensional (2-D) Vlasov–Poisson numerical simulations. We initialise the system in a stable state and gradually drive it towards instability with an imposed, weak external electric field, thus avoiding physically unrealisable super-critical initial conditions. A comprehensive analysis of the nonlinear evolution of ion-acoustic turbulence (IAT) is presented, including the detailed characteristics of the evolution of the particles’ distribution functions, (2-D) wave spectrum and the resulting anomalous resistivity. Our findings reveal the dominance of 2-D quasi-linear effects around saturation, with nonlinear effects, such as particle trapping and nonlinear frequency shifts, becoming pronounced during the later stages of the system's nonlinear evolution. Remarkably, the Kadomtsev–Petviashvili (KP) spectrum is observed immediately after the saturation of the instability. Another crucial and noteworthy result is that no steady saturated nonlinear state is ever reached: strong ion heating suppresses the instability, which implies that the anomalous resistivity associated with IAT is transient and short-lived, challenging earlier theoretical results. Towards the conclusion of the simulation, electron-acoustic waves are triggered by the formation of a double layer and strong modifications to the particle distribution induced by IAT.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Relativistic wind farm effect: Possibly turbulent flow of a charged, massless relativistic fluid in graphene

At low Reynolds numbers, the wind flow in the wake of a single wind turbine is generally not turbulent. However, turbines in wind farms affect each other's wakes so that a turbulent flow can arise. In the present work, an analog of this effect for the massless charge carrier flow around obstacles in graphene is outlined. We use a relativistic hydrodynamic simulation to analyze the flow in a sample containing impurities. Depending on the density of impurities in the sample, we indeed find evidence for potentially turbulent flow and discuss experimental consequences.

Watson, Mark (ORCID:0000000258959949)↗
DOE OSTI2024

Understanding changes in cloud simulations from E3SM version 1 to version 2

Abstract. This study documents clouds simulated by the Energy Exascale Earth System Model (E3SM) version 2 (E3SMv2) and attempts to understand what causes the model behavior change in clouds relative to E3SMv1. This is done by analyzing the last 30-year (1985–2014) data from the 165-year historical simulations using E3SMv1 and v2 and four sensitivity tests to isolate the impact of changes in model parameter choices in its turbulence, shallow convection, and cloud macrophysics parameterization (Cloud Layers Unified By Binormals, CLUBB); microphysical parameterization (MG2); and deep-convection scheme (ZM), as well as model physics changes in convective triggering. It is shown that E3SMv2 significantly improves the simulation of subtropical coastal stratocumulus clouds and clouds with optical depth larger than 3.6 over the stratocumulus-to-cumulus transition regimes, where the shortwave cloud radiative effect (SWCRE) is also improved, and the Southern Ocean (SO) while seeing an overall slight degradation in low clouds over other tropical and subtropical oceans. The better performance in E3SMv1 over those regions is partially due to error compensation between its simulated optically thin and intermediate low clouds for which E3SMv2 actually improves simulation of optically intermediate low clouds. Sensitivity tests indicate that the changes in low clouds are primarily due to the tuning done in CLUBB. The impact of the ZM tuning is mainly on optically intermediate and thick high clouds, contributing to an improved SWCRE and longwave cloud radiative effect (LWCRE). The impact of the MG2 tuning and the new convective trigger is primarily on the high latitudes and the SO. They have a relatively smaller impact on clouds than CLUBB tuning and ZM tuning do. This study offers additional insights into clouds simulated in E3SMv2 by utilizing multiple data sets and the Cloud Feedback Model Intercomparison Project (CFMIP) Observation Simulator Package (COSP) diagnostic tool as well as sensitivity tests. The improved understanding will benefit future E3SM developments.

58 GEOSCIENCES↗
DOE OSTI2024

Direct microstability optimization of stellarator devices

Turbulent transport is regarded as one of the key issues in magnetic confinement nuclear fusion, both for tokamaks and stellarators. Here, in this work, we show that a significant decrease in a microstability-based proxy, as opposed to a geometric one, for the turbulent heat flux, namely the quasilinear heat flux, can be obtained in an efficient manner by coupling stellarator optimization with linear gyrokinetic simulations. This is accomplished by computing the quasilinear heat flux at each step of the optimization process, as well as the deviation from quasisymmetry, and minimizing their sum, leading to a balance between neoclassical and the turbulent transport proxy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Zonal flow excitation in electron-scale tokamak turbulence

Abstract The derivation of an intermediate-scale gyrokinetic-electron theory in nonuniform tokamak plasmas (Chen et al 2021 Nucl. Fusion 61 066017 ) has shown that a Navier–Stokes type nonlinearity couples electron-temperature-gradient (ETG) modes and zonal flow (ZF) modes with wavelengths much shorter than the ion gyroradius but much longer than the electron gyroradius. This intermediate-scale ETG-ZF coupling is typically stronger than the Hasegawa–Mima type nonlinearity characteristic of the fluid approximation and is predicted to lead to relevant ZF generation and ETG mode regulation. Electron-scale, continuum, gyrokinetic simulation results are presented here which include both single-mode ETG and full-spectrum ETG turbulence. The ZF generation due to single ETG modes is investigated and the single-mode intermediate-scale results are found to be in agreement with theory. The full-spectrum results are then presented and explained qualitatively in terms of the single-mode results. It is found that the ETG-driven ZFs regulate intermediate-scale electron heat flux transport to levels in the predicted range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Expanded verification and validation studies of hypersonic aerodynamics with multiple physics-fidelity models

Hypersonic aerothermodynamics is an important domain of modern multiphysics simulation. The Multi-Fidelity Toolkit is a simulation tool being developed at Sandia National Laboratories to predict aerodynamic properties for compressible flows from a range of physics fidelities and computational speeds. These models include the Reynolds-averaged Navier–Stokes (RANS) equations, the Euler equations with momentum-energy integral technique (MEIT), and modified Newtonian aerodynamics with flat-plate boundary layer (MNA+FPBL) equations, and they can be invoked independently or coupled with hierarchical Kriging to interpolate between high-fidelity simulations using lower-fidelity data. However, as with any new simulation capability, verification and validation are necessary to gather credibility evidence. This work describes formal code- and solution-verification activities, as well as model validation with uncertainty considerations. Code verification activities on the MNA+FPBL model build on previous work by focusing on the viscous portion of the model. Viscous quantities of interest are compared against those from an analytical solution for flat-plate, inclined-plate, and cone geometries. The code verification methodology for the MEIT model is also presented. Test setup and results of code verification tests on the laminar and turbulent models within MEIT are shown. Solution-verification activities include grid-refinement studies on simulations that model the HIFiRE-1 wind tunnel experiments. These experiments are used for validation of all model fidelities. A thorough validation comparison with prediction error and uncertainty is also presented. Three additional HIFiRE-1 experimental runs are simulated in this study, and the solution verification and validation work examines the effects of the associated parameter changes on model performance. Finally, a study is presented that compares the computational costs and fidelities from each of the different models.

42 ENGINEERING↗
DOE OSTI2025

Applying a Compact Porous Media Model to Numerically Derive Resistance Coefficients for Lattice Structures

Additive Manufacturing allows for exploring various geometries to achieve specific engineering criteria. Lattices are one geometry with unique properties, including being periodically repeating structures which allow flow through them to be represented as a porous media according to Darcy-Forchheimer equations. These equation’s coefficients are generally experimentally derived, but this work demonstrates the ability to numerically derive them with CFD. Simulations were performed using three-dimensional stead state Reynolds-averaged Navier-Stokes with a k-ω Shear Stress Transport turbulence model using Ansys Fluent. Three lattice geometries were investigated and drag coefficients were derived. The method was validated against externally published data for similar geometries demonstrating strong agreement, and grid convergence for all simulations was calculated with a Grid Convergence Index method. Wall roughness is demonstrated to have a non-negligible impact on results and roughness values are considered for the primary focus Octahedral geometry where both smooth wall and rough wall coefficients were derived. The porosity coefficients for the Octahedral geometry at 1.0 [m/s] were found to be 2.89×10 6 and 2.90×10 6 [1/(Pa*m*s)] for the permeability coefficients, 6.37×10 1 and 5.44×10 1 [m 2 /kg] for the inertial resistance coefficients, and with a max pressure drop of 5116.7 [Pa] and 4429.5 [Pa] for the smooth walls and rough walls, respectively. The derived numerical method enables rapid exploration and optimization of new lattice designs for diverse engineering applications.

42 ENGINEERING↗
DOE OSTI2022

Universal Anomalous Diffusion of Quantized Vortices in Ultraquantum Turbulence

In classical viscous fluids, turbulent eddies are known to be responsible for the rapid spreading of embedded particles. However, in inviscid quantum fluids where the turbulence is induced by a chaotic tangle of quantized vortices, dispersion of the particles can be achieved via a nonclassical mechanism, i.e., their binding to the evolving vortices. However, knowledge on how the vortices diffuse and spread in quantum-fluid turbulence is very limited, especially for the so-called ultraquantum turbulence (UQT) generated by a random tangle of vortices. Here we report a systematic numerical study of the apparent diffusion of vortices in UQT in superfluid helium-4 using the full Biot-Savart simulation. We reveal that the vortices in the superfluid exhibit a universal anomalous diffusion (superdiffusion) at small times, which transits to normal diffusion at large times. This behavior is found to be the result of a generic scaling property of the vortex velocity. Our simulation at finite temperatures also nicely reproduces recent experimental observations. Lastly, the knowledge obtained from this study may form the base for understanding turbulent transport and universal vortex dynamics in various quantum fluids.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

An a priori evaluation of a principal component and artificial neural network based combustion model in diesel engine conditions

A principal component analysis (PCA) and artificial neural network (ANN) based chemistry tabulation approach is presented. ANNs are used to map the thermochemical state onto a low-dimensional manifold consisting of five control variables that have been identified using PCA. Three canonical configurations are considered to train the PCA-ANN model: a series of homogeneous reactors, a nonpremixed flamelet, and a two-dimensional lifted flame. The performance of the model in predicting the thermochemical manifold of a spatially-developing turbulent jet flame in diesel engine thermochemical conditions is a priori evaluated using direct numerical simulation (DNS) data. The PCA-ANN approach is compared with a conventional tabulation approach (tabulation using ad hoc defined control variables and linear interpolation). The PCA-ANN model provides higher accuracy and requires several orders of magnitude less memory. Here, these observations indicate that the PCA-ANN model is superior for chemistry tabulation, especially for modelling complex chemistries that present multiple combustion modes as observed in diesel combustion. The performance of the PCA-ANN model is then compared to the optimal estimator, i.e. the conditional mean from the DNS. The results indicate that the PCA-ANN model gives high prediction accuracy, comparable to the optimal estimator, especially for major species and the thermophysical properties. Higher errors are observed for the minor species and reaction rate predictions when compared to the optimal estimator. It is shown that the prediction of minor species and reaction rates can be improved by using training data that exhibits a variation of parameters as observed in the turbulent flame. The output of the ANN is analysed to assess mass conservation. It is observed that the ANN incurs a mean absolute error of 0.05% in mass conservation. Furthermore, it is demonstrated that this error can be reduced by modifying the cost function of the ANN to penalise for deviation from mass conservation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

A discontinuous Galerkin spectral element method for compressible reacting flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large-eddy simulations because of their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reacting Navier-Stokes equations. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of the DG approach. The framework, implemented in the spectral element code Nek5000, is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. An entropy-residual based artificial viscosity is added to smooth shocked regions of flow, and a positivity-preserving limiter is implemented to suppress non-physical oscillations. These enhancements support the numerical stability of the hydrodynamic sub-step, which is decoupled from the chemistry integration through a second-order operator splitting method. Here, a series of smooth and discontinuous validation cases are presented in increasing physical and computational complexity for both inviscid and viscous flows. In particular, simulations of canonical one-dimensional and two-dimensional detonations are performed, and the high-order numerical results are validated against available literature data. Additional validation studies are carried out for classical three-dimensional numerical simulations of incompressible and compressible turbulent flows.

Compressible reacting flows↗
DOE OSTI2024

DRDMannTurb: A Python package for scalable, data-driven synthetic turbulence

Synthetic turbulence models (STMs) are used in wind engineering to generate realistic flow fields and are employed as inputs to industrial wind simulations. Examples include prescribing inlet conditions in large eddy simulations that model loads on wind turbines and tall buildings. We are interested in STMs capable of generating fluctuations based on prescribed second-moment statistics since such models can simulate environmental conditions that closely resemble on-site observations. To this end, the widely used Mann model (see Mann, 1994, 1998) is the inspiration for DRDMannTurb. The Mann model is described by three physical parameters: a magnitude parameter influencing the global variance of the wind field and corresponding to the Kolmogorov constant multiplied by the rate of viscous dissipation of the turbulent kinetic energy to the two-thirds, αϵ 2/3 , a turbulence length scale parameter L, and a nondimensional parameter Γ related to the lifetime of the eddies. A number of studies, as well as international standards (e.g., those by the International Electrotechnical Commission (IEC)), include recommended values for these three parameters with the goal of standardizing wind simulations according to observed energy spectra. Yet, having only three parameters, the Mann model faces limitations in accurately representing the diversity of observable spectra. This Python package enables users to extend the Mann model and more accurately fit field measurements through flexible neural network models of the eddy lifetime function. Following Keith et al. (2021), we refer to this class of models as Deep Rapid Distortion (DRD) models. DRDMannTurb also includes a general module implementing an efficient method for synthetic turbulence generation based on a domain decomposition technique. This technique is also described in Keith et al. (2021).

17 WIND ENERGY↗
DOE OSTI2021

Advances in physics understanding of high poloidal beta regime toward steady-state operation of CFETR

Experimental and modeling investigations of the high βp scenarios on the DIII-D and EAST tokamaks show advantages in high energy confinement, avoidance of n = 1 MHD, and core-edge integration with reduced heat flux, making this scenario an attractive option for CFETR steady-state operation. Experiments show that plasmas with high confinement and high density can be achieved with neutral beam injection on DIII-D (β p ~ 2.2, β N ~ 3.5, f BS ~ 50%, f Gw ~ 1.0, H 98y2 ~ 1.5) and pure RF power on EAST (β P ~ 2.0, β N ~ 1.6, f BS ~ 50%, f Gw ~ 0.8, H 98y2 > 1.3). By tailoring the current density profile, a q-profile with local (off-axis) negative shear is achieved that yields improved confinement and MHD stability. Transport analysis and simulation suggest that the combination of high density gradient and high Shafranov shift allows turbulence stabilization and higher confinement. Using on-axis ECH injection, tungsten accumulation is avoided on EAST, and this is reproduced in modeling. Reduced heat flux (by >40%) and maintenance of high core confinement is achieved with active feedback control of the radiated divertor, an important result for long pulse operation in tokamaks. In conclusion, the improved physics understanding and validated modeling tools are used to design a 1GW steady-state scenario for CFETR.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Transport and chemistry of isoprene and its oxidation products in deep convective clouds

Deep convective clouds can transport trace gases from the planetary boundary layer into the upper troposphere where subsequent chemistry may impact aerosol particle formation and growth. In this modelling study, we investigate processes that affect isoprene and its oxidation products injected into the upper troposphere by an isolated deep convective cloud in the Amazon. We run a photochemical box model with coupled cloud microphysics along hundreds of individual air parcel trajectories sampled from a cloud-resolving model simulation of a convective event. The box model simulates gas-phase chemical reactions, gas scavenging by liquid and ice hydrometeors, and turbulent dilution inside a deep convective cloud. The results illustrate the potential importance of gas uptake to anvil ice in regulating the intensity of the isoprene oxidation and associated low volatility organic vapour concentrations in the outflow. Isoprene transport and fate also depends on the abundance of lightning-generated nitrogen oxide radicals (NO x = NO + NO 2 ). If gas uptake on ice is efficient and lightning activity is low, around 30% of the boundary layer isoprene will survive to the cloud outflow after approximately one hour of transport, while all the low volatile oxidation products will be scavenged by the cloud hydrometeors. If lightning NO x is abundant and gas uptake by ice is inefficient, then all isoprene will be oxidised during transport or in the immediate outflow region, while several low volatility isoprene oxidation products will have elevated concentrations in the cloud outflow. Reducing uncertainties associated with the uptake of vapours on ice hydrometeors, especially HO 2 and oxygenated organics, is essential to improve predictions of isoprene and its oxidation products in deep convective outflows and their potential contribution to new particle formation and growth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Theory of zonal flow growth and propagation in toroidal geometry

The toroidal geometry of tokamaks and stellarators is known to play a crucial role in the linear physics of zonal flows (ZFs), leading to e.g. the Rosenbluth–Hinton residual and geodesic acoustic modes. However, descriptions of the nonlinear ZF growth from a turbulent background typically resort to simplified models of the geometry. We present a generalised theory of the secondary instability to model the ZF growth from turbulent fluctuations in toroidal geometry, demonstrating that the radial magnetic drift substantially affects the nonlinear ZF dynamics. In particular, the toroidicity gives rise to a new branch of propagating ZFs, the toroidal secondary mode, which is nonlinearly supported by the turbulence. We present a theory of this mode and compare the theory against gyrokinetic simulations of the secondary mode. The connection with other secondary modes—the ion-temperature-gradient and Rogers–Dorland–Kotschenreuther secondary modes—is also examined.

Nies, Richard [Princeton Univ., NJ (United States)↗