Search NASA⌕ Search

SEARCH · Search NASA

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 271 records · Page 15

Numerical Study on the Aerodynamics of an Iced Airfoil With Scale-Resolving Simulations

Scale-resolving simulations of the NACA 23012 airfoil with horn ice accretion on the leading edge are conducted using the hybrid Reynolds-averaged Navier–Stokes/large-eddy simulation (hybrid RANS/LES) and wall-modeled large-eddy simulation (WMLES) approaches implemented in the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework. Aerodynamic results at the Reynolds number of 1.8 million show good comparison with the experimental measurements at different angles of attack from pre-stall to post-stall regimes. The pressure plateaus caused by the flow separation and the recovery of pressure inside the separation bubble around the iced leading edge are well predicted with the scale-resolving simulations when sufficient grid resolution is used around the accreted ice. The unsteadiness of the turbulent flows around the iced airfoil is also examined through the turbulent kinetic energy with the Reynolds normal stress anisotropy. Kelvin–Helmholtz instability (KHI) arises at the shear layer triggered by the upper ice horn and leads to rapid laminar-to-turbulent transition over a large range of angle of attack. With the increase of the angle of attack, the region with high turbulence intensity induced by the unstable shear layer spreads quickly over the entire upper surface of the airfoil. The coherent KHI modes from the upper and lower ice horns are extracted using the spectral proper orthogonal decomposition (SPOD) technique. The SPOD modes extracted from the upper shear layer have large-scale variations in the spanwise direction and low-rank behavior where the energy of the leading SPOD mode at each Strouhal number of the KHI largely represents the total energy when the mode number in the spanwise direction is small.

ARMD↗

Numerical Study on the Aerodynamics of an Iced Airfoil With Scale-Resolving Simulations

Scale-resolving simulations of the NACA 23012 airfoil with horn ice accretion on the leading edge are conducted using the hybrid Reynolds-averaged Navier–Stokes/large-eddy simulation (hybrid RANS/LES) and wall-modeled large-eddy simulation (WMLES) approaches implemented in the Launch, Ascent, and Vehicle Aerodynamics (LAVA) framework. Aerodynamic results at the Reynolds number of 1.8 million show good comparison with the experimental measurements at different angles of attack from pre-stall to post-stall regimes. The pressure plateaus caused by the flow separation and the recovery of pressure inside the separation bubble around the iced leading edge are well predicted with the scale-resolving simulations when sufficient grid resolution is used around the accreted ice. The unsteadiness of the turbulent flows around the iced airfoil is also examined through the turbulent kinetic energy with the Reynolds normal stress anisotropy. Kelvin–Helmholtz instability (KHI) arises at the shear layer triggered by the upper ice horn and leads to rapid laminar-to-turbulent transition over a large range of angle of attack. With the increase of the angle of attack, the region with high turbulence intensity induced by the unstable shear layer spreads quickly over the entire upper surface of the airfoil. The coherent KHI modes from the upper and lower ice horns are extracted using the spectral proper orthogonal decomposition (SPOD) technique. The SPOD modes extracted from the upper shear layer have large-scale variations in the spanwise direction and low-rank behavior where the energy of the leading SPOD mode at each Strouhal number of the KHI largely represents the total energy when the mode number in the spanwise direction is small.

ARMD↗

Effect of Distributed Patch of Smooth Roughness Elements on Transition in a High-Speed Boundary Layer

Surface roughness is known to have a substantial impact on the aerothermodynamic loading of hypersonic vehicles, particularly via its influence on the laminar-turbulent transition process within the boundary layer. Numerical simulations are performed to investigate the effects of a distributed region of densely packed, smooth-shaped roughness elements on the laminar boundary layer over a 7-degree half-angle, circular cone for flow conditions corresponding to a selected trajectory point from the ascent phase of the HIFiRE-1 flight experiment. For peak-to-valley roughness heights of 50 percent or less in comparison with the thickness of the unperturbed boundary layer, the computations converge to a stationary flow, suggesting that the flow is globally stable. Analysis of convective instabilities in the wake of the roughness patch indicates two dominant families of unstable disturbances, namely, a high frequency mode that corresponds to Mack mode waves modified by the wake and a lower frequency mode that corresponds to shear layer instabilities associated with the streaks in the roughness wake. Even though the peak growth rate of the later mode is more than 35 percent greater than the peak growth rates of the Mack modes, the latter modes achieve higher amplification ratios, and hence, are likely to dominate the onset of transition, which is estimated to occur slightly later than that in the unperturbed, i.e., smooth surface boundary layer. Additional computations are performed to investigate the effects of various roughness patch configurations on a Mach 3.5 flat plate boundary layer, to help guide an upcoming experiment in the Mach 3.5 Supersonic Low Disturbance Tunnel at NASA Langley Research Center. In this case, the cumulative reinforcement of basic state distortion over the length of the roughness patch is predicted to yield a significantly earlier transition than that over a smooth plate or a plate with a shorter length roughness patch.

Meelan Choudhari↗

Hypersonic Boundary Layer Transition Measurements Using NO2 approaches NO Photo-dissociation Tagging Velocimetry

Measurements of instantaneous and mean streamwise velocity profiles in a hypersonic laminar boundary layer as well as a boundary layer undergoing laminar-to-turbulent transition were obtained over a 10-degree half-angle wedge model. A molecular tagging velocimetry technique consisting of a NO2 approaches?NO photo-dissociation reaction and two subsequent excitations of NO was used. The measurement of the transitional boundary layer velocity profiles was made downstream of a 1-mm tall, 4-mm diameter cylindrical trip along several lines lying within a streamwise measurement plane normal to the model surface and offset 6-mm from the model centerline. For laminar and transitional boundary layer measurements, the magnitudes of streamwise velocity fluctuations are compared. In the transitional boundary layer the fluctuations were, in general, 2-4 times larger than those in the laminar boundary layer. Of particular interest were fluctuations corresponding to a height of approximately 50% of the laminar boundary layer thickness having a magnitude of nearly 30% of the mean measured velocity. For comparison, the measured fluctuations in the laminar boundary layer were approximately 5% of the mean measured velocity at the same location. For the highest 10% signal-to-noise ratio data, average single-shot uncertainties using a 1 ?Es and 50 ?Es interframe delay were ~115 m/s and ~3 m/s, respectively. By averaging single-shot measurements of the transitional boundary layer, uncertainties in mean velocity as low as 39 m/s were obtained in the wind tunnel. The wall-normal and streamwise spatial resolutions were 0.14-mm (2 pixel) and 0.82-mm (~11 pixels), respectively. These measurements were performed in the 31-inch Mach 10 Air Wind Tunnel at the NASA Langley Research Center.

Bathel, Brett F.↗

Flow regime and Reynolds number variation effects on the mixing behavior of parallel flows

The hydraulic single-phase mixing of three parallel rectangular channels is experimentally investigated at various Reynolds numbers (Re) and flow regime combinations. Particle Image Velocimetry results for seven mixing cases are presented and discussed with varying Re combinations ranging from 1,824 to 20,844. While all cases result in the same Re ratio of ~0.69 between the inner and outer flows, two cases represent multi-regime mixing with the inner-outer regime pair of laminar-transitional and transitional-turbulent, while the other 5 cases are all characteristic of turbulent mixing with varying levels of turbulence. The outer channels initially share characteristics with a backward facing step. The center channel is found to initially behave like a slot jet, but then sees a significant increase in velocity decay. This inner flow velocity decay increased dramatically in the laminar-transitional mixing case, whose centerline velocity decay was ~6 times larger than the decay in the turbulent mixing cases. Second order statistics revealed a consistent mixing layer thickness of ~0.1 hydraulic diameters for all the cases but showed more intense shearing in the multi-regime mixing cases. The combined point and thereby the mixing layer length is determined using centerline velocity decay profiles, which show a much more aggressive mixing in multi-regime flows. Multi-regime mixing demonstrated superior characteristics relative to turbulent mixing due to a more dramatic velocity decay in the inner flow and a shorter mixing length. The contributions of this work include communicating the benefits of multi-regime mixing and providing detailed characterization efforts that can serve future efforts for validating computational models. Here this research also lays the groundwork for future studies aimed at achieving high levels of mixing without a severe penalty in pressure drop.

42 ENGINEERING↗

Transition Delay via Vortex Generators in a Hypersonic Boundary Layer at Flight Conditions

The potential of realizable, stationary streaks undergoing non-modal growth to stabilize a hypersonic boundary-layer flow and, subsequently, delay the laminar-turbulent transition onset, is studied via numerical computations. The geometry and flow conditions are selected to match a relevant trajectory location from the ascent phase of the HIFiRE-1 flight experiment, namely, a 7-degree half-angle cone with 2.5 mm nose radius, freestream Mach number of 5.30, freestream unit Reynolds number equal to 13.42 x 10(exp 6)/m, and wall-to-adiabatic temperature ratio of approximately 0.35 over most of the test article. This paper investigates flow modifications induced by wall-mounted vortex generators (VGs), followed by an analysis of the modal instability of the perturbed, streaky boundary-layer flow. Results are presented both for a single array of VGs that was designed on the basis of optimal growth theory and for a VG configuration involving two separate arrays with opposite orientations that ware designed to provide staged control of flow instabilities while simultaneously reducing the amplification of streak instabilities resulting from the control devices. Earlier research had shown that the onset of transition during the HIFiRE-1 flight experiment, which did not include any control devices, correlated with an amplification factor of N = 14.7 for the planar Mack modes. If one assumes that the transition N -factor is not affected by the introduction of the VGs, then the control configurations based on a single array of VGs and two separate arrays would result in a transition delay of 17% and 40%, respectively. These findings suggest a passive flow control s to induce streaks that would delay transition in hypersonic boundary dominated by Mack-mode instabilities.

Incompressible boundary layer↗

Nature of transonic sub-Alfvénic turbulence and density fluctuations in the near-Sun solar wind

Context. Recent Parker Solar Probe (PSP) measurements have revealed that solar wind (SW) turbulence transits from a subsonic to a transonic regime near the Sun, while remaining sub-Alfvénic. These observations call for a revision of the existing SW models, where turbulence is considered to be both subsonic and sub-Alfvénic. Aims. In this work, we introduce a new magnetohydrodynamic (MHD) model of transonic sub-Alfvénic turbulence (TsAT). Methods. We used 3D MHD simulations initialized with parameters measured by PSP to investigate the properties of the new near-Sun SW transonic turbulent regime. We then derived a reduced set of MHD equations in the transonic sub-Alfvénic limit to interpret our numerical results. Results. Our TsAT model shows that turbulence is effectively nearly incompressible (NI) and has a 2D + slab (quasi-2D) geometry not only in the subsonic limit, but also in the transonic regime, as long as it remains sub-Alfvénic, a condition essentially enforced everywhere in the heliosphere by the strong local magnetic field. These predictions are consistent with 3D MHD simulations, showing that transonic turbulence is dominated by low-frequency quasi-2D incompressible structures, while compressible fluctuations are a minor component corresponding to low-frequency slow modes and high-frequency fast modes. Conclusions. Our new TsAT model extends existing NI theories of turbulence, and is potentially relevant for the theoretical and numerical modeling of space and astrophysical plasmas, including the near-Sun SW, the solar corona, and the interstellar medium.

79 ASTRONOMY AND ASTROPHYSICS↗

A three-dimensional simulation of transition and early turbulence in a time-developing mixing layer

The physics of the transition and early turbulence regimes in the time developing mixing layer was investigated. The sensitivity of the mixing layer to the disturbance field of the initial condition is considered. The growth of the momentum thickness, the mean velocity profile, the turbulence kinetic energy, the Reynolds stresses, the anisotropy tensor, and particle track pictures of computations are all examined in an effort to better understand the physics of these regimes. The amplitude, spectrum shape, and random phases of the initial disturbance field were varied. A scheme of generating discrete orthogonal function expansions on some nonuniform grids was developed. All cases address the early or near field of the mixing layer. The most significant result shows that the secondary instability of the mixing layer is produced by spanwise variations in the straining field of the primary vortex structures.

Cain, A. B.↗

DeepONet-Assisted Optimization of Surface Topography for Transition Delay in A Mach 4.5 Boundary Layer

We use deep learning, an ensemble variationaltechnique (EnVar), and direct numerical simulations(DNS) to design an optimal topography for a two-dimensional roughness element that delays the on-set of laminar-turbulent transition in a Mach 4.5 flat-plate boundary layer. Deep operator networks (Deep-ONets), which have the known ability to learn com-plex nonlinear operators within dynamical systems,are used for machine learning. For the baseline config-uration of a smooth flat plate, the second-mode wavesat the DNS inflow cause a quick nonlinear breakdownof the high-speed boundary layer within the computa-tional domain. Results reported in the present studyvalidate the ability of DeepONets to model the tran-sition delay via a given topography of the roughnesselement. The computing cost to optimize the rough-ness element for minimal skin-friction drag is substan-tially lowered by the DeepONets-based reduced-ordermodel. In comparison to the baseline method of EnVaroptimization based on DNS alone, the DeepONets-based EnVar optimizer is able to delay transition pastthe outflow boundary of the computational domainwhile utilizing almost 5–6 times fewer DNS.

Machine Learning↗

Transition Analysis for the CRM-NLF Wind Tunnel Configuration

This paper presents the results of an ongoing study into the linear stability characteristics of the boundary layer flow over the common research model with natural laminar flow (CRMNLF) aircraft configuration. The flow conditions match selected test conditions from a recent wind tunnel experiment in the National Transonic Facility at the NASA Langley Research Center. Previous work involving parallel stability computations of a boundary layer flow based on the conical flow approximation has shown that the measured onset of laminar-turbulent transition during the experiments can be correlated with the linear amplification of Tollmien- Schlichting (TS) and stationary crossflow (CF) instabilities in the swept wing boundary layer. Here, we examine the effects of the simplifying approximations in both basic state computation and the stability analysis, with the goal of quantifying the resulting changes in the N-factor correlations. Specifically, the basic states are computed by using full Navier-Stokes equations and the stability analysis is performed by using a nonorthogonal coordinate system that allows a clear distinction between planar TS and CF instabilities. Furthermore, the effects of curvature and nonparallel mean flow have been included in the stability computations based on the parabolized stability equations (PSE). The fully turbulent Reynolds-Averaged-Navier-Stokes (RANS) mean flow solutions show good agreement with the measured wall pressure distribution. Viscous-inviscid interactive effects are observed to be important because the shock fronts along the suction surface are influenced by the imposed transition front. The stability results confirm the previous findings related to TS amplification within the inboard region of the wing and the dominance of stationary CF modes in the outboard region. However, given the close proximity of the measured transition front and the dual shock system within the outer part of the wing, the onset of transition may well be shock limited within the outboard region. In general, the transition criterion based on the dual N-factor method with N TS = N CF = 6 is reasonably successful at correlating with the measured transition fronts at Re MAC = 15 million and AoA = 1.5, 2 degrees; however, the low values of the correlating N-factors at Re MAC = 17.5 million support the hypothesis that the measured transition at the higher Reynolds number may have been strongly influenced by the merging of turbulent wedges that originate from surface imperfections near the leading edge.

Boundary layer transition↗

DeepONet-Assisted Optimization of Surface Topography for Transition Delay in a Mach 4.5 Boundary Layer

We use deep learning, an ensemble variational technique (EnVar), and direct numerical simulations(DNS) to design an optimal topography for a two-dimensional roughness element that delays the on-set of laminar-turbulent transition in a Mach 4.5 flat-plate boundary layer. Deep operator networks (DeepONets), which have the known ability to learn complex nonlinear operators within dynamical systems, are used for machine learning. For the baseline configuration of a smooth flat plate, the second-mode waves at the DNS inflow cause a quick nonlinear breakdown of the high-speed boundary layer within the computational domain. Results reported in the present study validate the ability of DeepONets to model the transition delay via a given topography of the roughness element. The computing cost to optimize the rough-ness element for minimal skin-friction drag is substantially lowered by the DeepONets-based reduced-order model. In comparison to the baseline method of EnVar optimization based on DNS alone, the DeepONets-based EnVar optimizer is able to delay transition past the outflow boundary of the computational domain while utilizing almost 5–6 times fewer DNS.

Machine Learning↗

Numerical Investigation of Two-Phase Flows With Charged Droplets in Electrostatic Field

A numerical method to solve two-phase turbulent flows with charged droplets in an electrostatic field is presented. The ensemble-averaged Navier-Stokes equations and the electrostatic potential equation are solved using a finite volume method. The transitional turbulence field is described using multiple-time-scale turbulence equations. The equations of motion of droplets are solved using a Lagrangian particle tracking scheme, and the inter-phase momentum exchange is described by the Particle-In-Cell scheme. The electrostatic force caused by an applied electrical potential is calculated using the electrostatic field obtained by solving a Laplacian equation and the force exerted by charged droplets is calculated using the Coulombic force equation. The method is applied to solve electro-hydrodynamic sprays. The calculated droplet velocity distributions for droplet dispersions occurring in a stagnant surrounding are in good agreement with the measured data. For droplet dispersions occurring in a two-phase flow, the droplet trajectories are influenced by aerodynamic forces, the Coulombic force, and the applied electrostatic potential field.

Kim, Sang-Wook↗

Receptivity of flat-plate boundary layer in a non-uniform free stream (vorticity normal to the plate)

Recent progress in both the linear and nonlinear aspects of stability theory has highlighted the importance of the receptivity problem. One of the most unclear aspects of receptivity study is the receptivity of boundary-layer flow normal to vortical disturbances. Some experimental and theoretical results permit the proposition that quasi-steady outer-flow vortical disturbances may trigger by-pass transition. In present work such interaction is investigated for vorticity normal to a leading edge. The interest in these types of vortical disturbances arise from theoretical work, where it was shown that small sinusoidal variations of upstream velocity along the spanwise direction can produce significant variations in the boundary-layer profile. In the experimental part of this work, such non-uniform flow was created and the laminar-turbulent transition in this flow was investigated. The experiment was carried out in a low-turbulence direct-flow wind tunnel T-361 at the Central Aerohydrodynamic Institute (TsAGI). The non-uniform flow was produced by laminar or turbulent wakes behind a wire placed normal to the plate upstream of the leading edge. The theoretical part of the work is devoted to studying the unstable disturbance evolution in a boundary layer with strongly non-uniform velocity profiles similar to that produced by outer-flow vorticity. Specifically, the Tollmien-Schlichting wave development in the boundary layer flow with spanwise variations of velocity is investigated.

Kogan, M. N.↗

CFD Vision 2030 Study: A Path to Revolutionary Computational Aerosciences

This report documents the results of a study to address the long range, strategic planning required by NASA's Revolutionary Computational Aerosciences (RCA) program in the area of computational fluid dynamics (CFD), including future software and hardware requirements for High Performance Computing (HPC). Specifically, the "Vision 2030" CFD study is to provide a knowledge-based forecast of the future computational capabilities required for turbulent, transitional, and reacting flow simulations across a broad Mach number regime, and to lay the foundation for the development of a future framework and/or environment where physics-based, accurate predictions of complex turbulent flows, including flow separation, can be accomplished routinely and efficiently in cooperation with other physics-based simulations to enable multi-physics analysis and design. Specific technical requirements from the aerospace industrial and scientific communities were obtained to determine critical capability gaps, anticipated technical challenges, and impediments to achieving the target CFD capability in 2030. A preliminary development plan and roadmap were created to help focus investments in technology development to help achieve the CFD vision in 2030.

Jeffrey P Slotnick↗

Measured Boundary Layer Transition and Rotor Hover Performance at Model Scale

An experiment involving a Mach-scaled, 11:08 f t: diameter rotor was performed in hover during the summer of 2016 at NASA Langley Research Center. The experiment investigated the hover performance as a function of the laminar to turbulent transition state of the boundary layer, including both natural and fixed transition cases. The boundary layer transition locations were measured on both the upper and lower aerodynamic surfaces simultaneously. The measurements were enabled by recent advances in infrared sensor sensitivity and stability. The infrared thermography measurement technique was enhanced by a paintable blade surface heater, as well as a new high-sensitivity long wave infrared camera. The measured transition locations showed extensive amounts, x=c>0:90, of laminar flow on the lower surface at moderate to high thrust (CT=s > 0:068) for the full blade radius. The upper surface showed large amounts, x=c > 0:50, of laminar flow at the blade tip for low thrust (CT=s < 0:045). The objective of this paper is to provide an experimental data set for comparisons to newly developed and implemented rotor boundary layer transition models in CFD and rotor design tools. The data is expected to be used as part of the AIAA Rotorcraft SimulationWorking Group

Overmeyer, Austin D.↗

Secondary Instability of Second Modes in Hypersonic Boundary Layers

Second mode disturbances dominate the primary instability stage of transition in a number of hypersonic flow configurations. The highest amplification rates of second mode disturbances are usually associated with 2D (or axisymmetric) perturbations and, therefore, a likely scenario for the onset of the three-dimensionality required for laminar-turbulent transition corresponds to the parametric amplification of 3D secondary instabilities in the presence of 2D, finite amplitude second mode disturbances. The secondary instability of second mode disturbances is studied for selected canonical flow configurations. The basic state for the secondary instability analysis is obtained by tracking the linear and nonlinear evolution of 2D, second mode disturbances using nonlinear parabolized stability equations. Unlike in previous studies, the selection of primary disturbances used for the secondary instability analysis was based on their potential relevance to transition in a low disturbance environment and the effects of nonlinearity on the evolution of primary disturbances was accounted for. Strongly nonlinear effects related to the self-interaction of second mode disturbances lead to an upstream shift in the upper branch neutral location. Secondary instability computations confirm the previously known dominance of subharmonic modes at relatively small primary amplitudes. However, for the Purdue Mach 6 compression cone configuration, it was shown that a strong fundamental secondary instability can exist for a range of initial amplitudes of the most amplified second mode disturbance, indicating that the exclusive focus on subharmonic modes in the previous applications of secondary instability theory to second mode primary instability may not have been fully justified.

Li, Fei↗

Heat Transfer Through Turbulent Friction Layers

The "general Prandtl number" Pr(exp 1) - A(sub q)/A Pr, aside from the Reynolds number determines the ratio of turbulent to molecular heat transfer, and the temperature distribution in turbulent friction layers. A(sub q) = exchange coefficient for heat; A = exchange coefficient for momentum transfer. A formula is derived from the equation defining the general Prandtl number which describes the temperature as a function of the velocity. For fully developed thermal boundary layers all questions relating to heat transfer to and from incompressible fluids can be treated in a simple manner if the ratio of the turbulent shear stress to the total stress T(sub t)/T in the layers near the wall is known, and if the A(sub q)/A can be regarded as independent of the distance from the wall. The velocity distribution across a flat smooth channel and deep into the laminar sublayer was measured for isothermal flow to establish the shear stress ratio T(sub t)/T and to extend the universal wall friction law. The values of T(sub t)/T which resulted from these measurements can be approximately represented by a linear function of the velocity in the laminar-turbulent transition zone. The effect of the temperature relationship of the material values on the flow near the wall is briefly analyzed. It was found that the velocity at the laminar boundary (in contrast to the thickness of the laminar layer) is approximately independent of the temperature distribution. The temperature gradient at the wall and the distribution of temperature and heat flow in the turbulent friction layers were calculated on the basis of the data under two equations. The derived formulas and the figures reveal the effects of the Prandtl number, the Reynolds number, the exchange quantities and the temperature relationship of the material values.

Reichardt, H.↗

Numerical and physical aspects of aerodynamic flows III; Proceedings of the Third Symposium, California State University, Long Beach, CA, January 21-24, 1985

The present symposium considers hydrodynamic stability and turbulent transition, transition calculations in three-dimensional flows, a quasi-simultaneous finite difference approach for strongly interacting flows, the significance of the thin layer Navier-Stokes approximation, unsteady airfoil boundary layers, predictions and experiments on airfoils at low Reynolds numbers, and a comparison of interactive boundary layer and thin layer Navier-Stokes procedures. Also discussed are a viscous-inviscid interaction method for computing unsteady transonic separation, massive separation and dynamic stall on a cusped trailing edge airfoil, the computation of turbulent separated flows over wings, an iterative scheme for three-dimensional transonic flows, and the computation of three-dimensional flows with shock wave/boundary layer interaction. Many of the papers in this symposium were abstracted previously (cf., A85-42951).

Cebeci, T.↗