SEARCH · Engineering Papers

Results for “turbulence”

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 217 records · Page 12

All sources NASA NTRS DOE OSTI

DOE OSTI2022

On the Injection Scale of the Turbulence in the Partially Ionized Very Local Interstellar Medium

The cascade of magnetohydrodynamic (MHD) turbulence is subject to ion–neutral collisional damping and neutral viscous damping in the partially ionized interstellar medium. By examining the damping effects in the warm and partially ionized local interstellar medium, we find that the interstellar turbulence is damped by neutral viscosity at ~261 au and cannot account for the turbulent magnetic fluctuations detected by Voyager 1 and 2. The MHD turbulence measured by Voyager in the very local interstellar medium (VLISM) should be locally injected in the regime where ions are decoupled from neutrals for its cascade to survive the damping effects. With the imposed ion–neutral decoupling condition and the strong turbulence condition for the observed Kolmogorov magnetic energy spectrum, we find that the turbulence in the VLISM is sub-Alfvénic, and its largest possible injection scale is ~194 au.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Anisotropic Turbulent Flows Observed in Above-the-loop-top Regions during Solar Flares

Abstract Solar flare above-the-loop-top (ALT) regions are vital for understanding solar eruptions and fundamental processes in plasma physics. Recent advances in three-dimensional (3D) magnetohydrodynamic (MHD) simulations have revealed unprecedented details on turbulent flows and MHD instabilities in flare ALT regions. Here, for the first time, we examine the observable anisotropic properties of turbulent flows in ALT by applying a flow-tracking algorithm on narrow-band extreme-ultraviolet images that are observed from the face-on viewing perspective. First, the results quantitatively confirm the previous observation that vertical motions dominate and that the anisotropic flows are widely distributed in the entire ALT region with the contribution from both upflows and downflows. Second, the anisotropy shows height-dependent features, with the most substantial anisotropy appearing at a certain middle height in ALT, which agrees well with the MHD modeling results where turbulent flows are caused by Rayleigh–Taylor-type instabilities in the ALT region. Finally, our finding suggests that supra-arcade downflows (SADs), the most prominently visible dynamical structures in ALT regions, are only one aspect of turbulent flows. Among these turbulent flows, we also report the antisunward-moving underdense flows that might develop due to MHD instabilities, as suggested by previous 3D flare models. Our results indicate that the entire flare fan displays group behavior of turbulent flows where the observational bright spikes and relatively dark SADs exhibit similar anisotropic characteristics.

Xie, Xiaoyan (ORCID:0009000705827807)↗

DOE OSTI2022

Tearing Instability in Alfven and Kinetic-Alfven Turbulence

Recently, it has been realized that magnetic plasma turbulence and magnetic field reconnection are inherently related phenomena. Turbulent fluctuations generate regions of sheared magnetic field that become unstable to the tearing instability and reconnection, thus modifying turbulence at the corresponding scales. In this contribution, we give a brief review of some recent results on tearing-mediated magnetic turbulence.We illustrate the main ideas of this rapidly developing field of study by concentrating on two important examples—magnetohydrodynamic Alfven turbulence and small-scale kinetic-Alfven turbulence. Due to various potential applications of these phenomena in space physics and astrophysics, we specifically try not to overload the text by heavy analytical derivations but rather present a qualitative discussion accessible to a non-expert in the theories of turbulence and reconnection.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Near wake evolution of a tidal stream turbine due to asymmetric sheared turbulent inflow with different integral length scales

Tidal stream turbines deployed at highly energetic open water sites are subjected to sheared inflow in the rotor plane. The inflow shear is expected to cause asymmetric loading on the rotor blades and affect the downstream wake. In the current study, two different turbulent inflow conditions, static-high shear and dynamic shear, were generated via an active-grid turbulence generator. A 1:20 scaled three-bladed horizontal axis tidal turbine model was tested in those conditions. The results were compared to a quasi-laminar case with no imposed turbulence or shear. The results show that the high shear reduces the average performance, with a drop of up to 16% in the optimal power coefficient. Besides, the shear profiles increase torque fluctuations and induce significant differences in wake hydrodynamics between the high-speed (upper) and low-speed (lower) regions. The large integral length scales further enhance the load fluctuations perceived by the rotor but have a negligible effect on the mean wake field quantities and the wake recovery. The lower half region featured a faster breakdown of tip vortex structure and a rapid drop of swirl number, a phenomenon conjectured to be a consequence of the strong turbulence intensities and Reynolds stresses in the lower half region. Furthermore, the sheared turbulent inflow also results in a very intensive energy redistribution process towards large-scale, low-frequency motions, which is important to the downstream turbines.

16 TIDAL AND WAVE POWER↗

DOE OSTI2023

New scaling law for turbulent boundary layers with high surface mass transfer

A fundamental understanding of the phenomena occurring in the turbulent boundary layer in the presence of surface blowing is limited, and considerable disagreements persist even in describing primary quantities, such as the boundary layer profile. The theories based on the linear boundary layer equations show that the thickness of the sublayer increases in the presence of surface blowing; therefore, the viscous sublayer and law of the wall modify. In this study, direct numerical simulations (DNS) of turbulent boundary layers with uniform surface mass transfers are carried out in order to scale the velocity profile. Emphasis is placed on moderate to high mass transfer rates, which are relevant to the most common hybrid rockets configuration. DNS data are used to establish a functional law of the wall and a law of wake by means of the relation between the wall shear stress and surface mass transfer. Analysis of the mean kinetic energy budget shows that the magnitude of turbulent kinetic energy increases by surface mass transfer, and the production rate extends significantly in the inner layer as the injection rate increases. DNS data of various surface blowing are used to complete the closure of turbulence kinetic energy equation and develop an eddy viscosity model. The predicted turbulent kinetic energy and eddy viscosity agree with DNS data for moderate to high blowing rates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Modeling of convective cells, turbulence, and transport induced by a radio-frequency antenna in the tokamak boundary plasma

The edge turbulence model Hermes (Dudson et al 2017 Plasma Phys. Control. Fusion 59 05401) is set up for plasma boundary simulations with an radiofrequency (RF) antenna, using parameters characteristic of a tokamak edge. Cartesian slab geometry is used with thin plate limiters representing the ion cyclotron range of frequency (ICRF) antenna side-wall limiters. Ad-hoc DC electric biasing of the limiters, motivated by calculations with VSim (Nieter et al 2004 J. Comput. Phys. 196 448), represents an induced RF sheath rectified potential in the plasma turbulence model. Flux-driven turbulence simulations demonstrate a realistic distribution of plasma profiles and fluctuations. There is a clear effect of the antenna sheath voltage leading to formation of convective cells; bias-induced convective transport flattens the scrape-off layer density profile and fluctuations penetrate into the shadow region of the limiters as the bias voltage increases. Turbulent transport for impurity ions is inferred by following ion trajectories in the simulated plasma turbulence fields, showing Bohm-like effective diffusion rates. All in all, the model elucidates the key physical phenomena governing the effects of ICRF-induced antenna biasing on the tokamak boundary plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Spatiotemporal Energy Cascade in Three-Dimensional Magnetohydrodynamic Turbulence

We present a new scale decomposition method to investigate turbulence in wavenumber-frequency space. Using 3D magnetohydrodynamic turbulence simulations, we show that magnetic fluctuations with time scales longer than the nonlinear time exhibit an inverse cascade toward even smaller frequencies. Low frequency magnetic fluctuations support turbulence, acting as an energy reservoir that is converted into plasma kinetic energy, the latter cascading toward large wavenumbers and frequencies, where it is dissipated. Our results shed new light on the spatiotemporal properties of turbulence, potentially explaining the origin and role of low frequency turbulent fluctuations in the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Dynamic Phase Alignment in Inertial Alfvén Turbulence

In weakly collisional plasma environments with sufficiently low electron beta, Alfvénic turbulence transforms into inertial Alfvénic turbulence at scales below the electron skin depth, k ⊥ d e ≳ 1. We argue that, in inertial Alfvénic turbulence, both energy and generalized kinetic helicity exhibit direct cascades. We demonstrate that the two cascades are compatible due to the existence of a strong scale dependence of the phase alignment angle between velocity and magnetic field fluctuations, with the phase alignment angle scaling as cos α k ∝ $k$$^{–1}_{⊥}$. The kinetic and magnetic energy spectra scale as ∝ $k$$^{–5/3}_{⊥}$ and ∝ $k$$^{–11/3}_{⊥}$, respectively. As a result of the dual direct cascade, the generalized helicity spectrum scales as ∝ $k$$^{–5/3}_{⊥}$, implying progressive balancing of the turbulence as the cascade proceeds to smaller scales in the k ⊥ d e >> 1 range. Turbulent eddies exhibit a phase-space anisotropy k ∥ ∝ $k$$^{5/3}_{⊥}$, consistent with critically balanced inertial Alfvén fluctuations. Furthermore, our results may be applicable to a variety of geophysical, space, and astrophysical environments, including the Earth’s magnetosheath and ionosphere, solar corona, and nonrelativistic pair plasmas, as well as to strongly rotating nonionized fluids.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Physics-informed Machine Learning for Modeling Turbulence in Supernovae

Abstract Turbulence plays an important role in astrophysical phenomena, including core-collapse supernovae (CCSNe), but current simulations must rely on subgrid models, since direct numerical simulation is too expensive. Unfortunately, existing subgrid models are not sufficiently accurate. Recently, machine learning (ML) has shown an impressive predictive capability for calculating turbulence closure. We have developed a physics-informed convolutional neural network to preserve the realizability condition of the Reynolds stress that is necessary for accurate turbulent pressure prediction. The applicability of the ML subgrid model is tested here for magnetohydrodynamic turbulence in both the stationary and dynamic regimes. Our future goal is to utilize this ML methodology (available on GitHub) in the CCSN framework to investigate the effects of accurately modeled turbulence on the explosion of these stars.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Investigation of Main Bearing Fatigue Estimate Sensitivity to Synthetic Turbulence Models Using a Novel Drivetrain Model Implemented in OpenFAST

ABSTRACT A coupled medium‐fidelity drivetrain model is developed and implemented in OpenFAST for a 10‐MW land‐based reference turbine. The implementation is verified against a fully coupled multibody wind turbine model, including a detailed drivetrain. The new model can simultaneously and accurately estimate main bearing loads and represent elastic bending of the drivetrain. It has low computational cost and is useful for early design phases, sensitivity analyses and complex systems like wind farms (where computational expense must be expended elsewhere). Here, the model is implemented for a monopile offshore wind turbine and used to investigate the sensitivity of main bearing basic rating life to different synthetic turbulence models. Large‐eddy simulations (LES) targeting stable, neutral, and unstable atmospheric conditions at below‐, near‐ and above‐rated wind speeds are used as a reference. The turbulence models recommended by the International Electrotechnical Commission, the Mann spectral tensor model, and the Kaimal spectral model with exponential coherence are fitted to the LES data. Additionally, a constrained turbulence generator, PyConTurb (short for Python Constrained Turbulence ), based on LES data, is applied in the aero‐hydro‐servo‐elastic simulations. Taking PyConTurb as the baseline, the Kaimal model significantly underestimates fatigue of the downwind main bearing, with between 10% and 40% less damage. The Mann model also underestimates the downwind main bearing fatigue by up to 30%. The upwind main bearing damage is driven by mean loads, and differences between models are less significant, although the trends are similar. Reasons for these discrepancies are investigated and attributed to differences in spatial and temporal variations among the turbulence models.

17 WIND ENERGY↗

DOE OSTI2023

Turbulence regimes in the nocturnal roughness sublayer: Interaction with deep convection and tree mortality in the Amazon

We investigated the influence of seasonality and proximity to the forest canopy on nocturnal turbulence regimes in the roughness sublayer of a Central Amazon forest. Since convective systems of different scales are common in this region, we also analyzed the effect of extreme wind gusts (propagated from convective downdrafts) on the organization of the turbulence regimes, and their potential to cause the mortality of canopy trees. Our data include high-frequency winds, temperature and ozone concentration at different heights during the dry and wet seasons of 2014. In addition, we used critical wind-speed data derived from a tree-winching experiment and a modeling study conducted in the same study site. Two different turbulence regimes were identified at three heights above the canopy: a weakly stable (WS) and a very stable regime (VS). The threshold wind speeds that mark the transition between turbulence regimes were larger during the dry season and increased as a function of the height above the canopy. The turbulent fluxes of sensible heat and momentum during the WS accounted for 88% of the entire nighttime flux. Downdrafts occurred only in the WS and favored a fully coupled state of wind flow along the canopy profile. The destructive potential of winds was four times higher than on nights without downdrafts.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

On the use of a multigrid-reduction-in-time algorithm for multiscale convergence of turbulence simulations

Simulations of turbulent flow present challenges in terms of accuracy and affordability on modern highly-parallel computer architectures. A multigrid-reduction-in-time algorithm is used to provide a framework for separately evolving different scales of turbulence and for parallelizing the temporal domain, thereby increasing the concurrency. It is hypothesized that the space–time locality of the small scales of turbulence can be used to circumvent difficulties in applying temporal multigrid to flows dominated by inertial physics. For algorithms that fall well short of spectral accuracy (fourth-order is used in this work) attention must be paid to the accuracy of features on scales transferred between multigrid levels. Numerical experiments were performed using implicit large-eddy simulation. Results from applying the approach to an infinite-Reynolds number Taylor–Green flow and a double-shear flow at a Reynolds number of 11650 provide strong evidence that the approach has merit. The multigrid-reduction-in-time framework can be used to parallelize the temporal domain of a high-Reynolds-number turbulent flow and permit independent convergence of different scales. Establishing this foundation allows for future research in reducing the wall-clock time to solve turbulent flows while retaining the same accuracy as sequential solvers. In conclusion, current performance results from parallelizing the temporal domain are not competitive with those from sequential-in-time methods.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

GPU-enabled extreme-scale turbulence simulations: Fourier pseudo-spectral algorithms at the exascale using OpenMP offloading

Fourier pseudo-spectral methods for nonlinear partial differential equations are of wide interest in many areas of advanced computational science, including direct numerical simulation of three-dimensional (3-D) turbulence governed by the Navier-Stokes equations in fluid dynamics. This paper presents a new capability for simulating turbulence at a new record resolution up to 35 trillion grid points, on the world's first exascale computer, Frontier, comprising AMD MI250x GPUs with HPE's Slingshot interconnect and operated by the US Department of Energy's Oak Ridge Leadership Computing Facility (OLCF). Key programming strategies designed to take maximum advantage of the machine architecture involve performing almost all computations on the GPU which has the same memory capacity as the CPU, performing all-to-all communication among sets of parallel processes directly on the GPU, and targeting GPUs efficiently using OpenMP offloading for intensive number-crunching including 1-D Fast Fourier Transforms (FFT) performed using AMD ROCm library calls. With 99% of computing power on Frontier being on the GPU, leaving the CPU idle leads to a net performance gain via avoiding the overhead of data movement between host and device except when needed for some I/O purposes. Memory footprint including the size of communication buffers for MPI_ALLTOALL is managed carefully to maximize the largest problem size possible for a given node count. Detailed performance data including separate contributions from different categories of operations to the elapsed wall time per step are reported for five grid resolutions, from 2048 3 on a single node to 32768 3 on 4096 or 8192 nodes out of 9408 on the system. Both 1D and 2D domain decompositions which divide a 3D periodic domain into slabs and pencils respectively are implemented. The present code suite (labeled by the acronym GESTS, GPUs for Extreme Scale Turbulence Simulations) achieves a figure of merit (in grid points per second) exceeding goals set in the Center for Accelerated Application Readiness (CAAR) program for Frontier. The performance attained is highly favorable in both weak scaling and strong scaling, with notable departures only for 2048 3 where communication is entirely intra-node, and for 32768 3 , where a challenge due to small message sizes does arise. Communication performance is addressed further using a lightweight test code that performs all-to-all communication in a manner matching the full turbulence simulation code. Performance at large problem sizes is affected by both small message size due to high node counts as well as dragonfly network topology features on the machine, but is consistent with official expectations of sustained performance on Frontier. Overall, although not perfect, the scalability achieved at the extreme problem size of 32768 3 (and up to 8192 nodes — which corresponds to hardware rated at just under 1 exaflop/sec of theoretical peak computational performance) is arguably better than the scalability observed using prior state-of-the-art algorithms on Frontier's predecessor machine (Summit) at OLCF. New science results for the study of intermittency in turbulence enabled by this code and its extensions are to be reported separately in the near future.

3D fast Fourier transform↗

DOE OSTI2026

Hybrid combustion modeling approach for turbulent jet ignition in natural-gas pre-chamber spark-ignition engines at high EGR

Here, this study presented a hybrid modeling approach for simulating turbulent jet ignition and combustion processes in a natural-gas pre-chamber spark-ignition engine operating under exhaust gas recirculation (EGR) diluted conditions. In-depth analyses of experimental data and simulation results from previous work [Chinnathambi et al., ICEF2021-67836; Kim et al., Fuel 409: 137815, 2026] revealed two key findings: (i) the magnitude of pressure difference between the pre-chamber and main chamber ($∆P_{PC-MC}$) was positively correlated with the combustion duration from the moment of $∆P_{PC-MC}=0$ to the point of 5% mass fraction burned, with larger $∆P_{PC-MC}$ associated with longer duration; and (ii) the turbulent combustion regime in the main chamber transitioned from the broken reaction zone to the corrugated flamelet regime, with the Karlovitz number exceeding 100 immediately after turbulent hot jets were ejected from nozzles, coinciding with observed local extinction events. To accurately simulate the entire combustion process, a hybrid approach was developed under Reynolds-Averaged Navier Stokes framework, combining the G-equation model for pre-chamber combustion with the multi-zone well-stirred reactor approach and a turbulence-chemistry interaction (TCI) submodel for main chamber combustion. The TCI submodel accounted for the attenuation of reaction rates due to turbulent strain and modeled local extinction by suppressing reaction rates under certain flow and flame conditions. When applied to three EGR rate conditions toward the dilution limit, the hybrid modeling approach accurately reproduced experimental data in terms of cylinder pressure, apparent heat release rate, and the observed positive correlation, including the delayed onset of main chamber combustion—a feature not captured by existing combustion models.

computational fluid dynamics simulation↗

DOE OSTI2020

Assessment of Low-Altitude Atmospheric Turbulence Models for Aircraft Aeroelasticity

We investigate the dynamic aeroelastic response of large but slow aircraft in low-altitude atmospheric turbulence. To this end, three turbulence models of increasing fidelity, namely, the one-dimensional von Kármán model, the two-dimensional Kaimal model and full three-dimensional wind fields extracted from large-eddy simulations (LES) are used to simulate ambient turbulence near the ground. Load calculations and flight trajectory predictions are conducted for a representative high-aspect-ratio wing aircraft, using a fully coupled nonlinear flight dynamics/aeroelastic model, when it operates in background atmospheric turbulence generated by the aforementioned models. Comparison of load envelopes and spectral content, on vehicles of varying flexibility, shows strong dependency between the selected turbulence model and aircraft aeroelastic response (e.g. 58% difference in the predicted magnitude of the wing root bending moment between LES and von Kármán models). This is mainly due to the presence of large flow structures at low altitudes that have comparable dimensions to the vehicle, and which despite the relatively small wind speeds within the Earth boundary layer, result in overall high load events for slow-moving vehicles. Results show that one-dimensional models that do not capture those effects provide fairly non-conservative load estimates and are unsuitable for very flexible airframe design.

atmospheric turbulence↗

DOE OSTI2022

An adaptive knowledge-based data-driven approach for turbulence modeling using ensemble learning technique under complex flow configuration: 3D PWR sub-channel with DNS data

This work describes a new approach to increase the accuracy of Reynolds-averaged Navier–Stokes (RANS) in modeling turbulence flow leveraging the machine learning technique. Traditionally, different turbulence models for Reynolds stress are developed for different flow patterns based on human knowledge. Each turbulence model has a certain application domain and prediction uncertainty. In recent years, with the rapid improvements of machine learning techniques, researchers start to develop an approach to compensate for the prediction discrepancy of traditional turbulence models with statistical models and data. However, the approach has deficiencies in several aspects. For example, the amount of human knowledge introduced to the statistical model couldn’t be controlled, which makes the statistical model learn from a very naïve stage and limits its application. In this work, a new approach is developed to address those deficiencies. Here, the new approach uses the “ensemble learning” technique to control the amount of human knowledge introduced into the statistical model. Therefore, the new approach could be adaptive to the multiple application domains. In conclusion, according to the results of case study, the new approach shows higher accuracy than both traditional turbulence models and the previous machine learning approach.

42 ENGINEERING↗

DOE OSTI2019

Fitting local, low-dimensional parameterizations of optical turbulence modeled from optimal transport velocity vectors

This work exploits a connection between optimal transport theory and the physics of image propagation to yield a locally low-dimensional model of turbulence-corrupted imagery. Optimal transport produces an invertible, pixel-wise linear trajectories to approximate the globally nonlinear turbulence between a clean and turbulence corrupted image pair. We use the low-dimensional model to fit subsets of the optimal transport vector fields and stitch the local models into a surrogate for the global map to be used for image cleaning. Experiments are performed on laboratory generated data of beam propagation using different values of the Fried parameter (a scale measuring turbulence coherence) as well as a toy data set. The results suggest this is a fruitful direction, and first step, towards using multiple realizations of turbulence corrupted images to learn a blind surrogate for the optimal transport vector field for image cleaning.

47 OTHER INSTRUMENTATION↗

DOE OSTI2020

Structure formation in turbulence as an instability of effective quantum plasma

Structure formation in turbulence can be understood as an instability of “plasma” formed by fluctuations serving as effective particles. These “particles” are quantumlike in the sense that their wavelengths are non-negligible compared to the sizes of background coherent structures. The corresponding “kinetic equation” describes the Wigner matrix of the turbulent field, and the coherent structures serve as collective fields. This formalism is usually applied to manifestly quantumlike or scalar waves. Here, we show how to systematically extend it to more complex systems using compressible Navier–Stokes turbulence as an example. In this case, the fluctuation Hamiltonian is a five-dimensional matrix operator and diverse modulational modes are present. As an illustration, we calculate these modes for a sinusoidal shear flow and find two modulational instabilities. One of them is specific to supersonic flows, and the other one is a Kelvin–Helmholtz-type instability that is a generalization of the known zonostrophic instability. Our calculations are readily extendable to other types of turbulence, for example, magnetohydrodynamic turbulence in plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗