Search NASA⌕ Search

SEARCH · Search NASA

Results for “turbulent flows”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Uniform momentum and temperature zones in unstably stratified turbulent flows

Wall-bounded turbulent flows exhibit a zonal arrangement, in which streamwise velocity organizes into uniform momentum zones (UMZs), separated by thin layers of elevated interfacial shear. While significant research efforts have focused on these structural features in neutrally stratified flows, the effects of unstable thermal stratification on UMZs and on analogous uniform temperature zones (UTZs) have not been considered previously. In this article, statistical properties of UMZs and UTZs are investigated using a suite of large eddy simulations of unstably stratified turbulent channel flow spanning weakly to highly convective conditions. When normalized by the friction velocity and stability-dependent mixing length, the mean velocity gradient based on UMZ interfacial velocity jumps and the vorticity thickness exhibits good collapse for all stabilities, establishing a link between UMZ properties and scaling predictions from Monin–Obukhov similarity theory. A similar relationship is found between UTZ properties and surface-layer scaling of the mean temperature gradient. In the mixed layer, mean UMZ depth is quasi-constant with wall-normal distance, while the deepest UTZs are found in the centre of the boundary layer. These instantaneous structures are found to be linked to the well-mixed velocity and temperature profiles in the convective mixed layer. Conditional averaging indicates that both UMZ and UTZ interfaces are associated with ejections of momentum and warm updrafts below the interface and sweeps of momentum and cool downdrafts above the interface. These results demonstrate a tangible connection between instantaneous structural features, mean properties and scaling laws in unstably stratified flows.

Mechanics↗
DOE OSTI2025

Multimode turbulent flow measurements using magnetic resonance imaging- and laser-based techniques and computational fluid dynamics simulations

We studied the flow field characteristics of a turbulent flow over a regularized cube array with a perpendicular injection flow through the floor between the second and third cubical elements, representing the complex flow interactions of a 3D jet and the wake flows behind cubical obstacles. Four different experimental measurements were performed: two magnetic resonance imaging-based measurements for three-dimensional three-component velocity (MRV) and concentration (MRC) and two laser-based techniques, particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF), for two-dimensional two-component velocity and concentration measurement, respectively. The mainstream Reynolds number is Re = 15 000⁠, based on the primary inlet velocity U m and channel height D H ⁠, whereas the injector Reynolds number is Re j = 3400⁠, based on the injector velocity U j and the injector's exit width D j ⁠. Numerical simulations were performed for the studied flow configuration of turbulent flow over a regularized cube array using Reynolds-averaged Navier–Stokes (RANS) and large-eddy simulation (LES) approaches. Results obtained from experimental measurements—including MRV, MRC, PIV, and PLIF—as well as RANS and LES simulations are discussed and compared along several horizontal and vertical planes of the studied configuration. In addition, 3D turbulent flow structures, such as leading-edge vortex, horseshoe vortex, and jet shear-layer vortex, and the isosurfaces of scalar concentration successfully revealed by MRV and MRC techniques were found to be in very good agreement with those 3D features extracted from RANS and LES simulations. In conclusion, the high-resolution experimental and numerical database obtained from this study could be useful for validation and verification of numerical codes.

Computational fluid dynamics↗
DOE OSTI2021

Applying Time-Parallelization to Turbulent Flows

Parallelization of the temporal domain is explored for the solution of turbulent flows. Multigrid reduction-in-time (MGRIT) is used to advance the large-scale fluid dynamics in time sequentially on the coarsest space-time grid but propagate the information in time parallel on all other levels. The goal of this process is to accurately and efficiently resolve the coarse-scale turbulence structure and use that to drive the fine-scales of the turbulent flow. The extra forcing from nonlinear multigrid facilitates the coupling and interaction between fine and coarse scales, through which the multiscale nonlinear physics is properly captured. Adaptive mesh refinement is employed to finely resolve only the regions with strong gradients, which provides further computational efficiency. The underlying computational fluid dynamics solver is a fourth-order finite-volume scheme with the standard 4-stage Runge-Kutta method. An advanced approach is devised and implemented to enable MGRIT to solve highly turbulent flows successfully. Furthermore, the method is applied to solve a Taylor-Green vortex problem and a doubleshear-layer turbulent mixing flow. Results are promising, validating that MGRIT with the filtering approach has the potential to efficiently solve general turbulent flows.

Computational Fluid Dynamics↗
DOE OSTI2026

MRV Challenge 3: velocity and passive scalar comparison in a 3D turbulent flow

The third iteration of a magnetic resonance velocimetry (MRV) challenge comparison activity consisting of 3D flow measurements for a turbulent flow in a water channel moving past centrally positioned cubic obstacles is reported. In this challenge iteration, MRV measurements are coupled with either temperature or concentration measurements to extend the diagnostic utility. Five research teams from around the world conducted the measurements on a single shared apparatus. The water channel included partial elements along the channel sidewalls that varied in height and precluded easy optical accessibility. In addition to a turbulent mainstream flow, a secondary flow entered the channel between the second and third cubic elements from a square injector hole at the channel bottom wall. The injector jet interacts with the mainstream flow and mixes turbulently as it advects downstream. For the selected flow regime and water solutions used as working fluids, the mixing of the higher temperature or concentration through the secondary flow with the mainstream flow satisfies the same dynamics so that both temperature and concentration fields behave as passive scalars and can be directly compared. The measurements are explained in detail by each participating team, and the results are interpolated onto a common grid and compared using line profiles, contour plots, and isosurfaces.

Benson, Mike [ORNL] (ORCID:000000023210116X)↗
DOE OSTI2022

An approximate inertial manifold (AIM) based closure for turbulent flows

A closure model for turbulent flows is developed based on a dynamical system theory. An appropriately discretized formulation of the governing equations is considered for this process. The key ingredient is an approximation of the system’s attractor, where all the trajectories in phase space are confined. This approximate inertial manifold based approach provides a path to track trajectories of the system in a lower-dimensional subspace. Unlike conventional coarse-graining approaches, the turbulent field is decomposed into resolved and unresolved dynamics using the properties of the governing equations. The novelty of the approach relies on the reconstruction of the unresolved field constrained by the governing equations. A posteriori tests for homogeneous isotropic turbulence and the Kuramoto–Sivashinsky equation show promising results for considerable dimension reduction with strong convergence properties. The proposed model outperforms the dynamic Smagorinsky model, and the computational overhead is competitive with similar approaches.

97 MATHEMATICS AND COMPUTING↗
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

Direct Numerical Simulations of Turbulent Flow over Misaligned Traveling Waves

Direct numerical simulations of turbulent flow over prescribed traveling waves in a half-channel flow setup subject to a combined streamwise and spanwise pressure gradient are performed at a friction Reynolds number of 180. The simulations undertaken in this study consider both flow-aligned (including opposing and following) and misaligned waves as we attempt to quantify the effect of wave misalignment on turbulence statistics, the pressure drag force, and the wave attenuation rate. For the simulations, we consider three characteristic wave age values corresponding to slow-, intermediate-, and fast-moving waves. Wave misalignment is taken into account by applying a spanwise pressure gradient vertical to the traveling waves, which results in a three-dimensional turbulent flow field above the moving waves. Key flow quantities such as the mean velocity, velocity variances, and momentum fluxes are found to vary with the wave parameters, confirming the findings of previous studies.

direct numerical simulation↗
DOE OSTI2021

Can Bump Arrays Separate Particles from Turbulent Flows?

In this paper, we evaluate the hypothesis that bump arrays can be used to separate particles from turbulent flows entering the array. Microfluidic bump arrays are known for separating particles by size from laminar inlet flows. However, turbulent inlet flows have not been explored but become important as microfluidic bump arrays are scaled up to mesofluidic bump arrays. We find experimentally that particle separation is indeed effective at higher Reynolds numbers. These experimental findings portend industrial scale particle separation due to the higher flow rates they facilitate.

Pease, Leonard F.↗
DOE OSTI2022

An optimization method for chaotic turbulent flow

Evidence indicates that quantities-of-interest in some turbulent flows can be controlled despite the overall chaotic dynamics. It is typically thought that this is via relatively deterministic, larger-scale components of the turbulence. However, finding such controls, if they exist, is challenging because chaos causes sensitivity gradients to explode and the search space to become intractably non-convex. This challenge is analyzed, and a penalty method is introduced to cope with it. In the new approach, the time domain is broken into segments approximately matching the chaos time scales, so that the solution within each segment is both physical and relatively deterministic. The initial condition of each segment is included in an adjoint-based gradient optimization, which temporarily introduces artificial Δq discontinuities in the overall solution. The optimization then proceeds in stages with increasing penalization of Δq. The method is developed and illustrated for a logistic map, the Lorenz Equation, and an advection augmented Kuramoto–Sivashinsky Equation. These examples show how the Δq temporarily increases the search scale prior to the strong Δq → 0 penalization that recovers a physical solution. It is then applied to turbulent Kolmogorov flow, for which it also far outperforms a standard adjoint-based gradient search. Finally, the utility of such an optimized chaotic solution is discussed.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2023

Strongly Heated Turbulent Flow in a Channel with Pin Fins

Large-eddy simulations (LES) were performed to study the turbulent flow in a channel of height H with a staggered array of pin fins with diameter D = H/2 as a function of heating loads that are relevant to the cooling of turbine blades and vanes. The following three heating loads were investigated—wall-to-coolant temperatures of Tw/Tc = 1.01, 2.0, and 4.0 - where the Reynolds number at the channel inlet was 10,000 and the back pressure at the channel outlet was 1 bar. For the LES, two different subgrid-scale models—the dynamic kinetic energy model (DKEM) and the wall-adapting local eddy-viscosity model (WALE)—were examined and compared. This study was validated by comparing with data from direct numerical simulation and experimental measurements. The results obtained show high heating loads to create wall jets next to all heated surfaces that significantly alter the structure of the turbulent flow. Results generated on effects of heat loads on the mean and fluctuating components of velocity and temperature, turbulent kinetic energy, the anisotropy of the Reynolds stresses, and velocity-temperature correlations can be used to improve existing RANS models.

42 ENGINEERING↗
DOE OSTI2021

DNS/LES Study of Representative Wall-Bounded Turbulent Flows using SIERRA/Fuego

This report summarizes a series of SIERRA/Fuego validation efforts of turbulent flow models on canonical wall-bounded configurations. In particular, direct numerical simulations (DNS) and large eddy simulations (LES) turbulence models are tested on a periodic channel, a periodic pipe, and an open jet for which results are compared to the velocity profiles obtained theoretically or experimentally. Velocity inlet conditions for channel and pipe flows are developed for application to practical simulations. To show this capability, LES is performed over complex terrain in the form of two natural hills and the results are compared with other flow solvers. The practical purpose of the report is to document the creation of inflow boundary conditions of fully developed turbulent flows for other LES calculations where the role of inflow turbulence is critical.

42 ENGINEERING↗
DOE OSTI2021

Application of Scalar Filtered Density Function to Turbulent Flows Under Supercritical Condition

Abstract The scalar filtered density function (FDF) methodology is extended and employed for large eddy simulation (LES) of turbulent flows under supercritical condition. To describe real fluid behavior, the extended methodology incorporates the generalized heat and mass diffusion models along with real fluid thermodynamic relations which are derived using the cubic Peng–Robinson equation of state. These models are implemented within the stochastic differential equations comprising the scalar FDF transport. Simulations are conducted of a temporally developing mixing layer under supercritical condition and the results are assessed by comparing with data generated by direct numerical simulation (DNS) of the same layer. The consistency of the proposed FDF methodology is assessed. The LES-FDF predictions are shown to agree favorably with the DNS data and exhibit several key features pertaining to supercritical turbulent flows.

Energy & Fuels↗
DOE OSTI2023

Direct Numerical Simulation of Transitional and Turbulent Flows Over Multi-Scale Surface Roughness—Part I: Methodology and Challenges

Abstract High-fidelity simulation of transitional and turbulent flows over multi-scale surface roughness presents several challenges. For instance, the complex and irregular geometrical nature of surface roughness makes it impractical to employ conforming structured grids, commonly adopted in large-scale numerical simulations due to their high computational efficiency. One possible solution to overcome this problem is offered by immersed boundary methods, which allow wall boundary conditions to be enforced on grids that do not conform to the geometry of the solid boundary. To this end, a three-dimensional, second-order accurate boundary data immersion method (BDIM) is adopted. A novel mapping algorithm that can be applied to general three-dimensional surfaces is presented, together with a newly developed data-capturing methodology to extract and analyze on-surface flow quantities of interest. A rigorous procedure to compute gradient quantities such as the wall shear stress and the heat flux on complex non-conforming geometries is also introduced. The new framework is validated by performing a direct numerical simulation (DNS) of fully developed turbulent channel flow over sinusoidal egg-carton roughness in a minimal-span domain. For this canonical case, the averaged streamwise velocity profiles are compared against results from the literature obtained with a body-fitted grid. General guidelines on the BDIM resolution requirements for multi-scale roughness simulation are given. Momentum and energy balance methods are used to validate the calculation of the overall skin friction and heat transfer at the wall. The BDIM is then employed to investigate the effect of irregular homogeneous surface roughness on the performance of an LS-89 high-pressure turbine blade at engine-relevant conditions using DNS. This is the first application of the BDIM to realize multi-scale roughness for transitional flow in transonic conditions in the context of high-pressure turbines. The methodology adopted to generate the desired roughness distribution and to apply it to the reference blade geometry is introduced. The results are compared to the case of an equivalent smooth blade.

Engineering↗
DOE OSTI2022

Clustering of inertial particles in turbulent flow through a porous unit cell

We report direct numerical simulation is used to investigate effects of turbulent flow in the confined geometry of a face-centred cubic porous unit cell on the transport, clustering and deposition of fine particles at different Stokes numbers (St = 0.01, 0.1, 0.5, 1, 2) and at a pore Reynolds number of 500. Particles are advanced using one-way coupling and the collision of particles with pore walls is modelled as perfectly elastic with specular reflection. Tools for studying inertial particle dynamics and clustering developed for homogeneous flows are adapted to take into account the embedded, curved geometry of the pore walls. The pattern and dynamics of clustering are investigated using the volume change of Voronoi tesselation in time to analyse the divergence and convergence of the particles. Similar to the case of homogeneous, isotropic turbulence, the cluster formation is present at large volumes, while cluster destruction is prominent at small volumes and these effects are amplified with the Stokes number. However, unlike homogeneous, isotropic turbulence, the formation of a large number of very small volumes was observed at all Stokes numbers and attributed to the collision of particles with the pore wall. Multiscale wavelet analysis of the particle number density indicates that the peak of the energy density spectrum, representative of enhanced particle clustering, shifts towards larger scales with an increase in the Stokes number. Scale-dependent skewness and flatness quantify the intermittent void and cluster distribution, with cluster formation observed at small scales for all Stokes numbers, and void regions at large scales for large Stokes numbers.

42 ENGINEERING↗
DOE OSTI2023

Reduced Sediment Settling in Turbulent Flows Due To Basset History and Virtual Mass Effects

The behavior of suspended particles in turbulent flows is a recalcitrant problem spanning wide-ranging fields including geomorphology, hydrology, and dispersion of particulate matter in the atmosphere. One key mechanism underlying particle suspension is the difference between particle settling velocity (w s ) in turbulence and its still water counterpart (w so ). This difference is explored here for a range of particle-to-fluid densities (1 – 10) and particle diameter to Kolmogorov micro-eddy sizes (0.1 – 10). Conventional models of particle fluxes that equate w s to w so result in eddy diffusivities and turbulent Schmidt numbers contradictory to laboratory experiments. Incorporating virtual mass and Basset history forces resolves these inconsistencies, providing clarity as to why w s /w so is sub-unity for the aforementioned conditions. The proposed formulation can be imminently used to model particle settling in turbulence, especially when sediment distribution outcomes over extended time scales far surpassing turbulence time scales are sought.

Geology↗
DOE OSTI2024

Coherent structures in stably stratified wall-bounded turbulent flows

To date, a growing body of literature has documented the existence and impacts of coherent structures known as large- and very-large-scale motions within wall-bounded turbulent flows under neutral and unstable thermal stratification. These coherent structures can account for a considerable fraction of the overall turbulent transport and have been found to modulate small-scale turbulent fluctuations near the wall. In the context of stably stratified flows, however, the examination of such coherent structures has garnered relatively little attention. Stable stratification limits vertical transport and turbulent mixing within flows, which makes it unclear the extent to which previous findings on coherent structures under unstable and neutral stratification are applicable to stably stratified flows. In this study, we investigate the existence and characteristics of coherent structures under stable stratification with a wide range of statistical and spectral analyses. Outer peaks in premultiplied spectrograms under weak stability indicate the presence of large-scale motions, but these peaks become weaker and eventually vanish with increasing stability. Quadrant analysis of turbulent transport efficiencies (the ratio of net fluxes to their respective downgradient components) demonstrates dependencies on both stability and height above ground, which is evidence of morphological differences in the coherent structures under increasing stability. Amplitude modulation by large-scale streamwise velocity was found to decrease with increasing gradient Richardson number, whereas modulation by large-scale vertical velocity was approximately zero across all stability ranges. For sufficiently stable stratification, large eddies are suppressed enough to limit any inner–outer scale interactions.

Greene, Brian R. (ORCID:0000000343766818)↗