Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 OSTI2020

Deep learning for in situ data compression of large turbulent flow simulations

As the size of turbulent flow simulations continues to grow, in situ data compression is becoming increasingly important for visualization, analysis, and restart checkpointing. For these applications, single-pass compression techniques with low computational and communication overhead are crucial. In this paper we present a deep-learning approach to in situ compression using an autoencoder architecture that is customized for three-dimensional turbulent flows and is well suited for contemporary heterogeneous computing resources. The autoencoder is compared against a recently introduced randomized single-pass singular value decomposition (SVD) for three different canonical turbulent flows: decaying homogeneous isotropic turbulence, a Taylor-Green vortex, and turbulent channel flow. Our proposed fully convolutional autoencoder architecture compresses turbulent flow snapshots by a factor of 64 with a single pass, allows for arbitrarily sized input fields, is cheaper to compute than the randomized single-pass SVD for typical simulation sizes, performs well on unseen flow configurations, and has been made publicly available. The results reported here show that the autoencoder dramatically outperforms a randomized single-pass SVD with similar compression ratio and yields comparable performance to a higher-rank decomposition with an order of magnitude less compression in regard to preserving a number of important statistical quantities such as turbulent kinetic energy, enstrophy, and Reynolds stresses.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Effects of turbulent flow regime on perforated plate stabilized piloted lean premixed flames

An experimental study of the effects of turbulent flow regime on the flame structure is conducted by using perforated-plate-stabilized, hydrogen-piloted, lean premixed methane/air turbulent flames. The underlying turbulent flow field was investigated using two-dimensional three-component particle imaging velocimetry (2D3C-PIV) with two perforated plates of different blockage ratios. The flow data allowed a separation of the turbulent flow regime into axial velocity stream dominated and vortex dominated flows. A plate with 62% blockage ratio was used to establish the stream-dominated flow regime and another with 86% blockage ratio was used to establish the vortex-dominated flow regime. OH laser-induced fluorescence was used to study the effects of the turbulent flow regime on the mean progress variable, flame brush thickness, flame surface density, and global consumption speed. In comparison with the stream-dominated flow, the vortex-dominated flow makes a wider and shorter flame. Also, the vortex-dominated flow has a thicker horizontal flame brush and a thinner longitudinal flame brush. Especially, the spatial variation of the horizontal flame brush thickness for the vortex-dominated flow does not follow the turbulent diffusion theory. The vortex-dominated flow shows a relatively constant flame surface density and the stream-dominated flow shows a decreasing flame surface density along the streamwise direction. The flame surface density for the vortex-dominated flow is higher than the one for the stream-dominated flow. Lastly, the vortex-dominated turbulent flow shows a much higher consumption speed in comparison with that of the stream-dominated turbulent flow with identical velocity fluctuation levels.

42 ENGINEERING↗
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↗