Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transition to 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.

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 OSTI2025

Characterization of the fractal dimension of helium plumes in laminar to turbulent transition

Abstract Schlieren imaging in conjunction with a high-speed camera was used to observe the behavior of metered helium plumes as they transition from laminar to turbulent flow in an air environment. The plumes were visualized at twelve jet Reynolds numbers ranging from 200 to 2980. The fractal dimension of the flows was obtained by applying a box counting algorithm to the recorded schlieren images. The results were analyzed to determine the correlation between the Reynolds number of the flow and the fractal dimension of the observed turbulence. A trend of increasing fractal dimension with increasing Reynolds number was observed for several different types of schlieren cutoffs including horizontal cutoff, vertical cutoff, circular cutoff, focused shadowgraphy and de-focused shadowgraphy. The vertical cutoff and focused shadowgraphy imaging methods showed the most consistent results for the fractal dimension characterization during the laminar to turbulent transition. For transitional plumes, it was observed that fractal dimension increased with distance from the jet outlet.

42 ENGINEERING↗
DOE OSTI2024

Turbulence Transition in Magnetically Confined Hydrogen and Deuterium Plasmas

In this study, we discovered a turbulence transition in a large helical device. The turbulence level and turbulence-driven energy transport decrease to a specific transition density and increase above it. The ruling turbulences below and above the transition density were ion-temperature gradient (ITG) and resistive-interchange (RI) turbulences, consistent with the predictions of gyrokinetic theory and two-fluid MHD model, respectively. Isotope experiments on hydrogen (H) and deuterium (D) clarified the role of transitions. In the ITG regime, turbulence levels and energy transport were comparable in the H and D plasmas. In contrast, in the RI regime, they were clearly suppressed in the D plasma. The results provide crucial knowledge for understanding isotope effects and future optimization of stellarator and heliotron devices. Published by the American Physical Society 2024

Physics↗
DOE OSTI2023

Using Taylor Green Vortex to Assess Modeling Turbulence Transition

The Taylor Green Vortex problem is an interesting problem to study the transition from laminar flow to fully developed turbulence. The problem consists initially of large-scale vortices that interact and eventually form smaller vortices that create a cascade of energy from larger scales to smaller scales. This document describes a) when fully developed turbulence occurs, and b) the scale separation in the inertial range. In addition, the resolution needed to capture perturbations in transitional flow is also proposed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Coarse grained simulation of convectively driven turbulent mixing, transition, and turbulence decay

Accurate predictions with quantifiable uncertainties are essential to many practical turbulent flow applications exhibiting extreme geometrical complexity and broad ranges of length and time scales. Under-resolved computer simulations are typically unavoidable in such applications, and implicit large-eddy simulation (ILES) often becomes the effective strategy. In this work, we focus on assessing ILES initialized with well-characterized 25 6 3 homogeneous isotropic turbulence datasets generated with direct numerical simulation (DNS). ILES solutions based on the LANL xRAGE code are studied as function of resolution for 6 4 3 , 12 8 3 , 25 6 3 , and 51 2 3 grids. ILES performance of new directionally-unsplit high-order numerical hydrodynamics algorithms in xRAGE is examined in this context. Compared to the initial 25 6 3 DNS, we find longer inertial subranges and higher turbulence Reynolds number (Re) for 25 6 3 and 51 2 3 xRAGE — attributed to having linked DNS (Navier–Stokes based) solutions to nominally inviscid (higher Re) Euler based ILES solutions. For fixed 25 6 3 resolution, we find that significantly larger simulated turbulence Re can be achieved with the higher-order unsplit (vs. split) discretizations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Coupling of Laminar-Turbulent Transition with RANS Computational Fluid Dynamics

This project combines several new concepts to create a boundary layer transition prediction capability that is suitable for analyzing modern hypersonic flight vehicles. The first new concept is the use of ''optimization'' methods to detect the hydrodynamic instabilities that cause boundary layer transition; the use of this method removes the need for many limiting assumptions of other methods and enables quantification of the interactions between boundary layer instabilities and the flow field imperfections that generate them. The second new concept is the execution of transition analysis within a conventional hypersonics CFD code, using the same mesh and numerical schemes for the transition analysis and the laminar flow simulation. This feature enables rapid execution of transition analysis with less user oversight required and no interpolation steps needed.

42 ENGINEERING↗
DOE OSTI2022

Design of a high-resolution Rayleigh-Taylor experiment with the Crystal Backlighter Imager on the National Ignition Facility

The Rayleigh-Taylor (RT) instability affects a vast range of High Energy Density (HED) length scales, spanning from supernova explosions (10 13 m) to inertial confinement fusion (10 -6 m). In inertial confinement fusion, the RT instability is known to induce mixing or turbulent transition, which in turn cools the hot spot and hinders ignition. The fine-scale features of the RT instability, which are difficult to image in HED physics, may help determine if the system is mixing or is transitioning to turbulence. Earlier diagnostics lacked the spatial and temporal resolution necessary to diagnose the dynamics that occur along the RT structure. A recently developed diagnostic, the Crystal Backlighter Imager (CBI), [1,2] can now produce an x-ray radiograph capable of resolving the fine-scale features expected in these RT unstable systems. This paper describes an experimental design that adapts a well-characterized National Ignition Facility (NIF) platform to accommodate the CBI diagnostic. Simulations and synthetic radiographs highlight the resolution capabilities of the CBI in comparison to previous diagnostics. Furthermore, the improved resolution of the system can provide new observations to study the RT instability's involvement in mixing and the transition to turbulence in the HED regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Impact of Numerical Hydrodynamics in Turbulent Mixing Transition Simulations

Underresolved simulations are unavoidable in high Reynolds (Re) and Mach (Ma) number turbulent flow applications at scale. Implicit large-Eddy simulation (ILES) often becomes the effective strategy to capture the dominating effects of convectively driven flow instabilities. We evaluate the impact of three distinct numerical strategies in simulations of transition and turbulence decay with ILES: the Harten–Lax–van Leer (HLL) Riemann solver applying Strang splitting and a Lagrange-plus-Remap formalism to solve the directional sweep—denoted split; the Harten–Lax–Van Leer-Contact (HLLC) Riemann solver using a directionally unsplit strategy and parabolic reconstruction—denoted unsplit; and the HLLC Riemann solver using unsplit and a low-Ma correction (LMC)—denoted unsplit*. Three case studies are considered: (1) a shock tube problem prototyping shock-driven turbulent mixing, (2) the Taylor–Green Vortex (TGV) prototyping transition to turbulence, and, (3) an homogeneous isotropic turbulence (HIT) case, focusing on the impact of discretization on transition and decay from fixed well-characterized initial conditions. Significantly more accurate predictions are provided by the unsplit schemes, in particular, when augmented with the LMC. For given resolution, only the unsplit schemes predict the turbulent mixing transition after reshock observed in the shock tube experiments. Relevant comparisons of ILES based on Euler and Navier–Stokes equations addressing potential occurrence of low-Re regimes in the applications are presented. Unsplit schemes are instrumental in allowing to capture the spatial development of the TGV flow and its validation at prescribed Re with significantly less resolution. HIT analysis confirms higher simulated turbulence Re and increased small-scale content associated with the unsplit discretizations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Laminar-to-turbulent flame transition and cycle-to-cycle variations in large eddy simulation of spark-ignition engines

This paper investigates the effect of laminar-to-turbulent flame transition modeling on the prediction of cycle-to-cycle variations (CCVs) in large eddy simulation (LES) of spark-ignition (SI) engines. A laminar-to-turbulent flame transition model that describes the non-equilibrium sub-filter flame speed evolution during an early stage of flame kernel growth is developed. In the present model, the flame transition is characterized by the flame kernel size at which the flame transition ends, defined here as the flame transition scale. The proposed model captures the effects that variations in a turbulent flow field have on the evolution of early-stage burning rates, through variations in the flame transition scale. The proposed flame transition model is combined with the front propagation formulation (FPF) method and a spark-ignition model to predict CCVs in a gasoline direct injection SI engine. It is found that multi-cycle LES with the proposed flame transition model reproduces experimentally-observed CCVs satisfactorily. When the transition model is not considered or when variations in the transition process are neglected, CCVs are significantly under-predicted for the case considered here. These results indicate the importance of modeling the laminar-to-turbulent flame transition and the effect of turbulence on the transition process, when predicting CCVs, under certain engine conditions. The LES results are also used to analyze sources for variations in the flame transition. It is found, for the present engine case, that the most important source is the cycle-to-cycle variation in the turbulence dissipation rate, which is used to measure the strength of turbulence in the proposed model, near a spark plug. The large-scale velocity field and the variations of the laminar flame speed due to the mixture composition and thermal stratification are also found to be important factors to contribute to the variations in the flame transition.

42 ENGINEERING↗
DOE OSTI2023

Transition to turbulence in viscoelastic channel flow of dilute polymer solutions

The transition to turbulence in a plane Poiseuille flow of dilute polymer solutions is studied by direct numerical simulations of a finitely extensible nonlinear elastic fluid with the Peterlin closure. The range of Reynolds number ($Re$)$2000 \le Re \le 5000$is studied but with the same level of elasticity in viscoelastic flows. The evolution of a finite-amplitude perturbation and its effects on the transition dynamics are investigated. A viscoelastic flow begins transition at an earlier time than its Newtonian counterparts, but the transition time appears to be insensitive to polymer concentration in the dilute and semi-dilute regimes studied. Increasing polymer concentration, however, decreases the maximum attainable energy growth during the transition process. The critical or minimum perturbation amplitude required to trigger transition is computed. Interestingly, both Newtonian and viscoelastic flows follow almost the same power-law scaling of$Re^\gamma$with the critical exponent$\gamma \approx -1.25$, which is in close agreement with previous studies. However, a shift downward is observed for viscoelastic flow, suggesting that smaller perturbation amplitudes are required for the transition. A mechanism of the early transition is investigated by the evolution of wall-normal and spanwise velocity fluctuations and flow structure. The early growth of these fluctuations and the formation of quasi-streamwise vortices around low-speed streaks are promoted by polymers, hence causing an early transition. These vortical structures are found to support the critical exponent$\gamma \approx -1.25$. Once the transition process is completed, polymers play a role in dampening the wall-normal and spanwise velocity fluctuations and vortices to attain a drag-reduced state in viscoelastic turbulent flows.

Mechanics↗
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 OSTI2024

Pressure-Gradient-Based RANS Model for Predicting Separation in Transitional and Turbulent Flows

Predicting flow separation poses a significant challenge for RANS models, particularly in transitional flows over airfoils. We propose a novel improvement to RANS models to predict incipient separation in both transitional and fully turbulent flows. Our approach modifies the eddy viscosity model in regions indicated by a pressure-gradient criterion that accounts for intermittency - determining whether the boundary layer is laminar or turbulent. This model demonstrates robust generalization across various airfoil shapes and Reynolds numbers. Applied to the NREL Phase VI wind turbine rotor, our model shows improved aerodynamic performance predictions compared to the baseline RANS model.

k-omega SST↗
DOE OSTI2023

Direct Numerical Simulation of Transitional and Turbulent Flows Over Multi-Scale Surface Roughness—Part II: The Effect of Roughness on the Performance of a High-Pressure Turbine Blade

Abstract Turbine blades generally present surface roughness introduced in the manufacturing process or caused by in-service degradation, which can have a significant impact on aero-thermal performance. A better understanding of the fundamental physical mechanisms arising from the interaction between the roughness and the turbine flow at engine-relevant conditions can provide insights for the design of blades with improved efficiency and longer operational life. To this end, a high-fidelity numerical framework combining a well-validated solver for direct numerical simulation and a second-order accurate immersed boundary method is employed to predict roughness-induced aero-thermal effects on an LS89 high-pressure turbine (HPT) blade at engine-relevant conditions. Different amplitudes and distributions of surface roughness are investigated and a reference smooth-blade simulation under the same flow conditions is conducted for comparison. Roughness of increasing amplitude progressively shifts the blade suction side boundary layer transition upstream, producing larger values of the turbulent kinetic energy and higher total wake losses. The on-surface data-capturing capabilities of the numerical framework provide direct measurements of the heat flux and the skin friction coefficient, hence offering quantitative information between the surface topology and engineering-relevant performance parameters. This work may provide a benchmark for future numerical studies of turbomachinery flows with roughness.

Engineering↗
DOE OSTI2024

A method for examining ensemble averaging forms during the transition to turbulence in HED systems for application to RANS models

This paper discusses a strategy to initialize a two-dimensional (2D) Reynolds-averaged Navier–Stokes model [LANL's Besnard–Harlow–Rauenzahn (BHR) model] in order to describe an unsteady transitional Richtmyer–Meshkov (RM)-induced flow observed in on-going high-energy-density ensemble experiments performed on the OMEGA-EP facility. The experiments consist of a nominal single-mode perturbation (initial amplitude a 0 ≈ 10 and wavelength $λ$ = 100μm) with target-to-target variations in the surface roughness subjected to the RM instability with delayed Rayleigh–Taylor in a heavy-to-light configuration. Our strategy leverages high-resolution three-dimensional (3D) implicit large eddy simulations (ILES) simulations to initialize BHR-relevant parameters and subsequently validate the 2D BHR results against the 3D ILES simulations. A suite of five 3D ILES simulations corresponding to five experimental target profiles is undertaken to generate an ensemble dataset. Using ensemble averages from the 3D simulations to initialize the turbulent kinetic energy in the BHR model ( K 0 ) demonstrates the ability of the model to predict the time evolution of the interface as well as the density-specific-volume covariance, b . To quantify the sensitivity of the BHR results to the choice of K 0 and the initial turbulent length scale, S 0 , we execute a parameter sweep spanning four orders of magnitude for both S 0 and K 0 , generating a parameter space consisting of 26 simulations. The Pearson's correlation coefficient is used as a measure of discrepancy between the 2D BHR and 3D ILES simulations and reveals that the ranges 8≲S 0 ≲20 μm and 10 9 ≲K 0 ≲10 10 cm 2 /s 2 produce predictions that agree best with the 3D ILES results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Phenomena Identification and Ranking Table (PIRT) for Heat Pipes

This Phenomena Identification and Ranking Table (PIRT) report provides an evaluation of key phenomena affecting the performance and operational regimes of heat pipes, particularly in the context of heat pipe microreactors (HPMRs). Heat pipes are advanced passive thermal management devices that utilize phase change and capillary action to achieve efficient heat transfer. However, due to the complexity of the phenomena coupled in the heat pipe, including phase change, turbulent transition, and compressibility effects, among others, there is high uncertainty in identifying and ranking the important phenomena affecting the operation of heat pipes and the current knowledge for their modeling and simulation and experimental measurements and instrumentation. This PIRT exercise, conducted as a collaborative effort involving the Department of Energy (DOE) Microreactor Program (MRP), the Nuclear Regulatory Commission (NRC), and university partners systematically identifies, reviews, and prioritizes critical phenomena affecting the operation of heat pipes based on their importance and knowledge levels. The report analyzes phenomena with high importance and low knowledge, such as wick de-wetting, critical heat flux, contact angles, and pressure dynamics, discussing challenges and future research directions for improving their modeling and simulation and experimental measurements. Additionally, the report addresses phenomena with low knowledge that could impact heat pipe operation during non-normal or transient operation, including frozen startup, laminar to turbulent transition, geysering, wick priming, underfilling conditions, surface roughness of the wick, NCGs trapped in the wick, and the timescales of startup and shutdown. This comprehensive evaluation serves as a valuable resource for guiding future research and development efforts, supporting the successful integration of heat pipes into critical applications such as nuclear reactors, and contributing to the advancement of heat pipe technologies in safety-critical industries.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2024

A Study of the Transition to Turbulence in a Bed of 67 Spherical Pebbles

Packed beds are commonly found in many engineering systems and have been widely studied for decades. A relatively new packed bed system is the Pebble Bed Reactor, a type of generation-IV nuclear reactor. Unlike many of the packed beds encountered in chemical and process engineering applications, Pebble Bed Reactors are larger and operate at significantly higher Reynolds numbers. As a result of these differences, there is a very limited amount of information on the detailed flow physics that exist in these complex geometries. This work seeks to contribute to a growing database of flow data for Pebble Bed Reactor systems by performing Direct Numerical Simulations of the flow in an experimental bed of 67 pebbles for a range of conditions. Simulations are performed at a Prandtl number of 0.66 and Reynolds numbers from 300–600. These Reynolds numbers are chosen to gain additional knowledge on the spatial development of turbulence in these systems. Analysis of the Turbulent Kinetic Energy, turbulence anisotropy, and Turbulent Heat Flux is performed. Results demonstrate significant development of the TKE across the tested range of Reynolds numbers. Examination of both the TKE and THF reveal that development first occurs near the center of the bed and propagates radially as the flow moves further into the bed. Notable regions of negative production of turbulent kinetic energy are observed in regions where flow accelerates around pebble contact points. Furthermore, these regions are found to coincide with regions of 1-component turbulence.Kindly check and confirm, all authors email id is correctly identified.These are correct

Direct numberical simulation↗
DOE OSTI2024

A Data-Driven Study of Rayleigh-Taylor Transition to Turbulence [Slides]

A temporal dynamics perspective of forgetting ICs; Autonomous system of neural ODEs of 0D dynamics; Steep growth of backward integration error quantifies forgetting; Does such a system exist in reality? (a low dim. dyn. sys. shadowing a projection of a high dim. dyn. sys.); Next: minimal QOIs to identify system; Extension to PDEs of 1D dynamics of DNS.

97 MATHEMATICS AND COMPUTING↗