Search NASA⌕ Search

SEARCH · Search NASA

Results for “Direct Numerical Simulation”

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 1,207 records · Page 67

The Intermediate Wake of a Thick Flat Plate with a Circular Trailing Edge

The intermediate wake region of a thick flat plate with a circular trailing edge (TE) is investigated with a direct numerical simulation (DNS). The upper and lower separating boundary layers are both turbulent and are statistically identical; the resulting wake is symmetric in the mean. Earlier research dealt with the near/very-near wake of the same plate (x/D < 13.0, x is the streamwise distance from the center of the circular TE and D is the plate-thickness/TE-diameter). In the present investigation the emphasis is on the evolution of shed-vortex structure and turbulence intensity distributions with increasing x; the focus is on the region 20.0 < x/D < 40.0. Profile similarity in wake velocity statistics is explored.

Rai, Man Mohan↗

Stabilization of a Swept-Wing Boundary Layer by Discrete Roughness Elements at High Reynolds Numbers

Direct numerical simulations (DNS) are performed to study potential stabilizing ef- fect of spanwise periodic discrete roughness elements (DREs) on cross ow instabilities in a spatially developing three-dimensional boundary layer over an in nite-swept natural- laminar- ow wing at a freestream Mach number of 0:75 and a chord Reynolds number of approximately 25 million. In the DNS, both the spanwise periodic DREs and distributed roughness in the leading-edge region are implemented to simulate a typical experimen- tal scenario in which multiple steady cross ow modes including the most unstable mode (i.e., the \target" mode) emerge because of the presence of naturally distributed surface roughness in the leading edge region and spanwise periodic control cylinders of subcritical wavelength are used to force small-wavelength disturbances (i.e., the control mode) for damping the target mode. The DNS results show that the e ectiveness of DRE control is sensitive to roughness diameter, height, and chordwise placement. For the DRE parame- ters considered in this study, the stabilizing e ect on the target mode is small within the computational domain that ended at about 35% of the chord.

Nicholson, Gary↗

DNS of Laminar to Turbulent Transition on NACA 0012 Airfoil with Sand Grain Roughness

The Lattice-Boltzmann-based solver PowerFLOW is used to perform direct numerical simulations of the transitional flow over an airfoil at Reynolds number equal to 0.657 million. The leading edge of the airfoil is covered with sand particles, represented by polyhedra, to mimic the grit used in experiments. The sensitivity of the laminar to turbulent transition to the size of these particles, grid resolution, spanwise length is evaluated and rectangular trips are also tested.

Ribeiro, A. F. P.↗

Effects of Wall Cooling on Supersonic Modes in High-Enthalpy Hypersonic Boundary Layers over a Cone

High-enthalpy hypersonic boundary-layer transition plays an important role in many entry/descent vehicles. Aerothermodynamic performance of these vehicles strongly depends on the transition location on the surface. However, detailed transition flow physics in these chemically reacting boundary layers are poorly understood and transition estimates during the design phase rely heavily on empirically-derived transition criteria such as 𝑹𝒆(𝜽)/𝑴(𝒆). One of the most intriguing characteristics in hypersonic boundary layers is the presence of unstable supersonic modes, first identified in the 1990s. Due to a recent surge in hypersonic applications, there has been renewed interest in studying the flow physics of supersonic modes using either theory or direct numerical simulations. This paper investigates the rise of supersonic modes in a high-enthalpy hypersonic flow over a 𝟓° half-angle cone at various wall temperatures using both quasiparallel linear stability theory and linear PSE. It was found that supersonic modes exist in all wall temperature conditions including the adiabatic wall case. Cooler wall temperature causes the second Mack mode to become an unstable supersonic mode naturally downstream of the upper-branch neutral location when the wall is sufficiently cooled. In terms of the integrated growth, the second mode is still the dominant mode. Nonetheless, supersonic modes can cause additional series of relatively weaker growth than that of the second mode beyond the peak amplitude location. According to the present linear PSE results, contrary to what was speculated in the literature, supersonic-mode pressure disturbance structure radiated into the freestream is nonacoustic in nature and the formation of unstable supersonic modes is mainly associated with the synchronization of phase speed between the instability and acoustic waves in a nonp ayer, not due to nonlinear modal interaction as suggested in the literature.

Chang, Chau-Lyan↗

Laminar to Turbulence Transition in Boundary Layers due to Tripping Devices

Direct numerical simulations (DNS)and linear stability analysis were performed to identify the mechanisms by which two-and three-dimensional trips promote early transition in subsonic boundary layers over a flat plate. Investigations were conducted for two boundary layer transition trip configurations: one is a two-dimensional rectangular-shaped rod and the other is a three-dimensional zigzag tape. Reynolds numbers based on the trip heights are 683 and 300 for the two-and three-dimensional trips, respectively. Simulations with the two-dimensional trip showed that the flow separates at the top of the trip and forms a long recirculation zone downstream. However, the flow remains steady in the simulation without any external disturbances. Linear stability analysis showed that the inflectional profile in the recirculation zone is strongly unstable to frequencies in the range of 100-300 Hz. Simulations with the external disturbances forced at the frequency of 200 Hz showed that disturbances grow strongly inside the recirculation zone as observed in the Klebanoff1 experiment.

Transition↗

Laminar to Turbulence Transition in Boundary Layers due to Tripping Devices

Direct numerical simulations (DNS)and linear stability analyses were performed to identify the mechanisms by whichtwo-and three-dimensional trips promote early transition in subsonic boundary layers over a flat plate. Investigations were conducted for threeboundary-layer transition trip configurations:a two-dimensional smooth hump,a two-dimensional rectangular-shapedrod and a three-dimensional zig-zagtape. The smooth hump was used by Park et al.1in the linear and nonlinear parabolized stability (PSE) computations. Reynolds numbersbased on the tripheights are 710, 683,and 300 for the two-and three-dimensional trips, respectively. Simulations withthe smooth hump showed that the flow becomes stable upstream and strongly unstable downstream of the hump. Simulations with the two-dimensional rectangular trip showed that the flow separates at the top of the trip and forms a long recirculation zonedownstream. However, the flow remains steady in the simulation without any external disturbances. Linear stability analysis showed that the inflectional profile in the recirculation zone is strongly unstable to frequencies in the range of 100-300 Hz. Simulations with the external disturbances forced at the frequency of 200 Hz showed that disturbances grow strongly inside the recirculation zone as observed in the Klebanoff2 experiment. Three-dimensional simulations with secondary disturbances of fundamental, subharmonic, and oblique types produced turbulent flows shortlydownstream. The three-dimensional zig-zag trip with a small roughness Reynolds number lead to early breakdown to a turbulent flow.

Stability↗

A Priori Subgrid Analysis of Temporal Mixing Layers with Evaporating Droplets

Subgrid analysis of a transitional temporal mixing layer with evaporating droplets has been performed using three sets of results from a Direct Numerical Simulation (DNS) database, with Reynolds numbers (based on initial vorticity thickness) as large as 600 and with droplet mass loadings as large as 0.5.

subgrid↗

Turbulent mixing in supercritical jets: effect of compressibility factor and inflow condition

Fuel injection and turbulent mixing at supercritical pressures determines ignition and combustion in numerous engineering applications. Flow evolution under such conditions is characterized by strong non-linear coupling between dynamics, transport coefficients, and thermodynamics. Experimental studies observe that the jets injected at supercritical pressures exhibit significantly different dynamics from the jets at subcritical conditions, owing to the lack of distinct liquid and gas phases in supercritical state. Thus, the averaged flow quantities such as the potential core length, jet spatial growth rate and velocity decay profiles differ in the two conditions, resulting in different mixed-fluid distributions. In this study, turbulent jet direct numerical simulations (DNS) are performed to examine the variations in flow statistics between injection of Nitrogen (N₂) in Nitrogen (N₂) at both subcritical (perfect-gas) and supercritical conditions. In all cases, isothermal round jets at Reynolds number (Re_{D}), based on jet diameter (D) and jet orifice velocity (U₀), of 5000 are considered. For mixing analyses, a passive scalar transported with the flow is examined.

Sharan, Nek↗

Budgets of Reynolds Stresses and Turbulent Heat Flux for Hypersonic Turbulent Boundary Layers Subject to Pressure Gradients

In this paper, the budget terms in the transport equations of the Reynolds stresses, the turbulent kinetic energy (TKE), the internal energy, and the turbulent heat flux are computed using a direct numerical simulation (DNS) database of hypersonic turbulent boundary layers (TBLs) subject to favorable pressure gradients (FPG). Study of the Reynolds stress and TKE budgets shows that FPG has a damping effect on turbulence production, transport, and destruction. The semilocal scaling fails to collapse the Reynolds stress and TKE budget terms between the zero pressure gradient (ZPG) and the strong FPG cases, suggesting that the mechanical nonequilibrium due to FPG leads to changes in Reynolds-stress and TKE transport that cannot be captured by Morkovin’s hypothesis. Minimal change is seen in the internal energy budget in the presence of FPG, indicating a minimal impact due to FPG on the mean thermal field. While the near-wall transport of wall-normal turbulent heat flux is largely governed by the balance between the turbulent viscous-thermal transport term and the turbulent viscous-thermal dissipation term, the transport of the streamwise turbulent heat flux is influenced by contributions from almost all of the budget terms. For the strong FPG case in particular, production due to both the mean and the fluctuating strain rates begins to play a significant role in the near-wall transport of the streamwise turbulent heat flux, and their inclusion or modeling may be required to accurately predict turbulent heat fluxes for a strong FPG.

Reynolds stress↗

Streaming Instability with Multiple Dust Species – I. Favourable Conditions for the Linear Growth

A recent study suggests that the streaming instability, one of the leading mechanisms for driving the formation of planetesimals, may not be as efficient as previously thought. Under some disc conditions, the growth time-scale of the instability can be longer than the disc lifetime when multiple dust species are considered. To further explore this finding, we use both linear analysis and direct numerical simulations with gas fluid and dust particles tomutually validate and study the unstablemodes of the instability in more detail.We extend the previously studied parameter space by one order of magnitude in both the range of the dust-size distribution [Ts, min, Ts, max] and the total solid-to-gas mass ratio ε and introduce a third dimension with the slope q of the size distribution. We find that the fast-growth regime and the slow-growth regime are distinctly separated in the ε–Ts, max space, while this boundary is not appreciably sensitive to q or Ts, min.With a wide range of dust sizes present in the disc (e.g. Ts, min  10−3), the growth rate in the slow-growth regime decreases as more dust species are considered.With a narrow range of dust sizes (e.g. Ts, max/Ts, min = 5), on the other hand, the growth rate in most of the ε–Ts, max space is converged with increasing dust species, but the fast and the slow growth regimes remain clearly separated. Moreover, it is not necessary that the largest dust species dominate the growth of the unstable modes, and the smaller dust species can affect the growth rate in a complicated way. In any case, we find that the fast-growth regime is bounded by ε  1 or Ts, max  1, which may represent the favourable conditions for planetesimal formation.

Zhaohuan Zhu↗

Interaction of a Tunnel-like Acoustic Disturbance Field with a Blunt Cone Boundary Layer at Mach 8

The existing measurements of laminar-to-turbulent transition over circular cones in conventional (i.e., “noisy”) hypersonic wind tunnels have established that the transition location moves downstream when the nose radius is increased from zero. However, this initially downstream movement slows down and ultimately reverses beyond a critical value of the nose radius, and may be related to external forcing in the form of freestream disturbances and/or surface roughness. To understand the effects of freestream acoustic disturbances on transition reversal over a blunt body, hypersonic boundary-layer receptivity to broadband freestream acoustic disturbances from the nozzle wall of a digital conventional wind tunnel is investigated by both direct numerical simulations (DNS) and modal and nonmodal stability analysis. A Mach 8 flow over a 7 deg half-angle cone with a nose radius of 𝑅 𝑛 = 5.2mm and freestream Reynolds number of 12.2 × 10 6 m -1 is considered. The results show that the broadband tunnel noise in the free stream of a convectional hypersonic wind tunnel (i.e., outside of the nozzle-wall turbulent boundary layer) can be well represented by an acoustic model with an ansatz of slow acoustic waves. With successful calibration of the model parameters against the precursor tunnel DNS, such an acoustic ansatz can successfully reproduce both the frequency-wavenumber spectra and the temporal evolution of the broadband tunnel noise radiated from the nozzle wall. Additionally, the DNS of the Mach 8 blunt cone with tunnel-like acoustic input above the bow shock showed that the spectra of wall-pressure and heat-transfer fluctuations recovers the signature of the axisymmetric waves predicted by the nonmodal analysis. Furthermore, the azimuthal wavenumber and frequency spectrum of the temperature fluctuations as a function of the wall-normal distance show higher amplitudes for three-dimensional waves above the boundary-layer edge. The numerical schlieren contours show the inclined structures commonly observed in blunt cone experiments, demonstrating that they correspond to three-dimensional structures due to freestream disturbances in the presence of an entropy layer.

Boundary layer transition↗

Numerical Study of the Turbulent Prandtl Number in Supersonic Plane-Channel Flow – the Effect of Thermal Boundary Conditions

Direct numerical simulation is used to investigate the turbulent Prandtl number Prt above cold (isothermal) and hot (adiabatic) walls in a family of supersonic channel flows. A range of Reynolds numbers and mean temperature/density variations is considered. The value of Prt away from the wall approaches 0.85 above both the isothermal and adiabatic walls, with the width of the Prt = 0.85 region increasing with wall-variable-based Reynolds number Reτw (the latter strongly affected by the thermal boundary condition). The variable, near-wall Prt profiles from both the present and previous, passive-scalar simulations collapse as a function of the semilocal y∗ wall scaling proposed by Huang et al., with the y∗-dependence weakly parameterized by Reτw.

Turbulence↗

Computational Study of Low-Speed Boundary-Layer Laminar-Turbulent Transition over a 30º Swept Plate in a Channel

Reliable boundary-layer laminar-turbulent transition predictions are critical for the design of aircraft that use natural or controlled laminar flow to achieve low drag and better fuel efficiency. During cruise flight through low turbulence, the three-dimensional boundary layer over swept surfaces can transition prematurely due to stationary crossflow instabilities induced by surface roughness. Parabolized stability equations analyses have been successfully used to explain the different phases of the transition process under well controlled perturbations such as steps or discrete roughness elements. However, methodologies are still lacking to model the boundary layer receptivity to distributed random roughness, and the following non-linear instability stages leading to transition. For engineering applications, several semi-empirical methods based on linear amplification factors or local helicity have found growing adoption, but their reliability is limited by their level of simplification and the sparsity of data they are based on. To address this gap, a new study of the effect of distributed random surface roughness on crossflow transition was started, using a 30◦ swept plate as a test vehicle. The plate, which is mounted in the 2 Foot × 3 Foot Low-Speed Boundary-Layer Channel at the NASA Langley Research Center, has been tested with micro-scale roughness of varied geometry near the leading edge. In order to analyze the flow in the test section and prepare Direct Numerical Simulations with resolved roughness, Reynolds-Averaged Navier-Stokes (RANS) simulations of the flow over the plate installed in the wind tunnel test section were performed using the Langtry-Menter transition model with stationary crossflow extension. This paper presents the results of these simulations, with comparison to the experiment.

ARMD↗

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↗

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↗

Hypersonic Boundary-Layer Transition over a Blunt Circular Cone in a Mach 8 Digital Wind Tunnel

To understand the effects of freestream acoustic disturbances on transition reversal over a blunt body, three-dimensional direct numerical simulations (DNS) were conducted to investigate the interaction of a tunnel-like acoustic disturbance field with a Mach 8 laminar boundary layer over a 7 deg half-angle blunt cone with a nose radius of 𝑅 𝑛 = 5.2 mm. The flow conditions and the axisymmetric cone geometry matched those measured in the Sandia Hypersonic Wind Tunnel at Mach 8 (Sandia HWT-8), and the incident freestream acoustic disturbances in the DNS were generated by acoustic radiation from the nozzle-wall turbulent boundary layer in the same tunnel. Consistent with the predictions of the nonmodal instability analysis, the study shows the emergence of a second spectral peak for moderately oblique waves between 150 kHz and 200 kHz in the aft part of the cone. The spectral energy of temperature and streamwise velocity fluctuations is initially concentrated within the entropy layer. These fluctiations are eventually “swallowed” by the boundary layer as the entropy-layer thickness 𝛿 𝑆 becomes similar to the boundary-layer thickness 𝛿h𝑡 at downstream locations. The numerical schlieren contours show the same inclined structures commonly observed in blunt cone experiments, and these structures are found to be three dimensional and advect with a constant speed approximately equal to the mean flow speed near the edge of the boundary layer.

Hypersonic Flow↗

Hypersonic Boundary-Layer Transition over a Blunt Circular Cone in a Mach 8 Digital Wind Tunnel

To understand the effects of freestream acoustic disturbances on transition reversal over a blunt body, three-dimensional direct numerical simulations (DNS) were conducted to investigate the interaction of a tunnel-like acoustic disturbance field with a Mach 8 laminar boundary layer over a 7 deg half-angle blunt cone with a nose radius of 𝑅 𝑛 = 5.2 mm. The flow conditions and the axisymmetric cone geometry matched those measured in the Sandia Hypersonic Wind Tunnel at Mach 8 (Sandia HWT-8), and the incident freestream acoustic disturbances in the DNS were generated by acoustic radiation from the nozzle-wall turbulent boundary layer in the same tunnel. Consistent with the predictions of the nonmodal instability analysis, the study shows the emergence of a second spectral peak for moderately oblique waves between 150 kHz and 200 kHz in the aft part of the cone. The spectral energy of temperature and streamwise velocity fluctuations is initially concentrated within the entropy layer. These fluctiations are eventually “swallowed” by the boundary layer as the entropy-layer thickness 𝛿 𝑆 becomes similar to the boundary-layer thickness 𝛿h𝑡 at downstream locations. The numerical schlieren contours show the same inclined structures commonly observed in blunt cone experiments, and these structures are found to be three dimensional and advect with a constant speed approximately equal to the mean flow speed near the edge of the boundary layer.

Hypersonic Flow↗

New Capabilities and Improvements to the High-Order Glenn Flux Reconstruction Code

The Glenn Flux Reconstruction (GFR) code is a computational fluid dynamics (CFD) code under development at NASA Glenn Research Center. GFR is based on the high-order flux reconstruction (FR) method and provides a large-eddy simulation (LES) capability that is both accurate and efficient for complex aeropropulsion flows. Three significant new capabilities have been added to the code that improve its performance and functionality. First, a variety of explicit Runge-Kutta methods, including some with adaptive time stepping, were added to GFR with two methods offering a 33% improvement in time-to-solution. Second, GFR can now utilize fully unstructured, mixed-element meshes to more easily facilitate the grid generation process for complex geometries. Finally, a rotating reference frame capability has been added to GFR for solving rotating turbomachinery problems. A selection of results demonstrating these new capabilities are presented in this work. The Taylor-Green vortex problem is used to verify the new unstructured capability by showing similar accuracy and resolution for all element types. LES of the Turbulent Heat Flux Phase III (THX3) experiment with comparison to another high-order LES code and a popular Reynolds-averaged Navier-Stokes (RANS) code demonstrates the accuracy of the code for complex aeropropulsion flows. Finally, LES of a spacecraft cabin ventilation fan shows the ability of GFR to efficiently establish a fan performance map and identify operating points for further analysis at high orders of accuracy.

High-Order Methods↗