Results for “turbulence methods” — Page 89

NASA NTRS1989

Numerical modeling of turbulent supersonic reacting coaxial jets

The paper considers the mixing and subsequent combustion within turbulent reacting shear layers. A computer program was developed to solve the axisymmetric Reynolds averaged Navier-Stokes equations. The numerical method integrates the Reynolds averaged Navier-Stokes equations using a finite volume approach while advancing the solution forward in time using a Runge-Kutta scheme. Three separate flowfields are investigated and it is found that no single turbulence model considered could accurately predict the degree of mixing for all three cases.

Eklund, Dean R.↗

NASA NTRS1988

Investigation of viscous/inviscid interaction in transonic flow over airfoils with suction

The viscous/inviscid interaction over transonic airfoils with and without suction is studied. The streamline angle at the edge of the boundary layer is used to couple the viscous and inviscid flows. The potential flow equations are solved for the inviscid flow field. In the shock region, the Euler equations are solved using the method of integral relations. For this, the potential flow solution is used as the initial and boundary conditions. An integral method is used to solve the laminar boundary-layer equations. Since both methods are integral methods, a continuous interaction is allowed between the outer inviscid flow region and the inner viscous flow region. To avoid the Goldstein singularity near the separation point the laminar boundary-layer equations are derived in an inverse form to obtain solution for the flows with small separations. The displacement thickness distribution is specified instead of the usual pressure distribution to solve the boundry-layer equations. The Euler equations are solved for the inviscid flow using the finite volume technique and the coupling is achieved by a surface transpiration model. A method is developed to apply a minimum amount of suction that is required to have an attached flow on the airfoil. The turbulent boundary layer equations are derived using the bi-logarithmic wall law for mass transfer. The results are found to be in good agreement with available experimental data and with the results of other computational methods.

Vemuru, C. S.↗

DOE OSTI2024

Langmuir turbulence in suspended kelp farms

This study investigates the influence of suspended kelp farms on ocean mixed layer hydrodynamics in the presence of currents and waves. We use the large eddy simulation method, where the wave effect is incorporated by solving the wave-averaged equations. Distinct Langmuir circulation patterns are generated within various suspended farm configurations, including horizontally uniform kelp blocks and spaced kelp rows. Intensified turbulence arises from the farm-generated Langmuir circulation, as opposed to the standard Langmuir turbulence observed without a farm. The creation of Langmuir circulation within the farm is attributed to two primary factors depending on farm configuration: (i) enhanced vertical shear due to kelp frond area density variability, and (ii) enhanced lateral shear due to canopy discontinuity at lateral edges of spaced rows. Both enhanced vertical and lateral shear of streamwise velocity, representing the lateral and vertical vorticity components, respectively, can be tilted into downstream vorticity to create Langmuir circulation. This vorticity tilting is driven by the Craik–Leibovich vortex force associated with the Stokes drift of surface gravity waves. In addition to the farm-generated Langmuir turbulence, canopy shear layer turbulence is created at the farm bottom edge due to drag discontinuity. The intensity of different types of turbulence depends on both kelp frond area density and the geometric configuration of the farm. The farm-generated turbulence has substantial consequences for nutrient supply and kelp growth. These findings also underscore the significance of the presence of obstacle structures in modifying ocean mixed layer characteristics.

Bo, Tong (ORCID:0000000260300561)↗

NASA NTRS1989

A review of high-speed, convective, heat-transfer computation methods

The objective of this report is to provide useful engineering formulations and to instill a modest degree of physical understanding of the phenomena governing convective aerodynamic heating at high flight speeds. Some physical insight is not only essential to the application of the information presented here, but also to the effective use of computer codes which may be available to the reader. A discussion is given of cold-wall, laminar boundary layer heating. A brief presentation of the complex boundary layer transition phenomenon follows. Next, cold-wall turbulent boundary layer heating is discussed. This topic is followed by a brief coverage of separated flow-region and shock-interaction heating. A review of heat protection methods follows, including the influence of mass addition on laminar and turbulent boundary layers. Also discussed are a discussion of finite-difference computer codes and a comparison of some results from these codes. An extensive list of references is also provided from sources such as the various AIAA journals and NASA reports which are available in the open literature.

Tauber, Michael E.↗

NASA NTRSDate not supplied

MEDLI2: MISP Inferred Aerothermal Environment and Flow Transition Assessment

The Mars Entry, Descent, and Landing Instrumentation 2 (MEDLI2) sensor suite on the Mars2020 mission contained multiple sensors on the aeroshell to measure the aerothermal environment during entry into the Martian atmosphere. These sensors performed superbly and successfully returned forebody and aftbody heating measurements. Analysis of MEDLI2 data indicated flow transitioning from a laminar to turbulent state on the heatshield. No evidence of flow transition was observed on the backshell. Two methods were used to estimate flow transition times on the heatshield: (1) temperature gradient of near-surface thermocouple data and (2) heat flux gradient from an inverse reconstruction approach using thermocouple data and material response modeling. Both methods produced similar transition times with an estimated accuracy of ±1 s. To assess various transition criteria, transition parameters were evaluated at each sensor location using flow field solutions from computational fluid dynamics (CFD) simulations. The idea was to use conservative values inferred from MEDLI2 data as transition criteria for other Mars missions. To test this hypothesis, MEDLI data from the Mars Science Laboratory (MSL) mission was used to compare predicted vs. actual flow transition times. The comparisons suggest smooth wall transition criteria are not well-suited in modeling the rapid progression of a turbulent transition front. Transition criteria containing a roughness element parameter agreed better with the flight data. In summary, critical transition values derived from MEDLI2 data may be used as a starting point in constructing a flow transition model for future Mars missions.

Chun Y Tang↗

DOE OSTI2022

Challenge Problem 1: Preliminary Model Development and Assessment of Flexible Heat Transfer Modeling Approaches

This report presents the modeling progress of a university consortium studying flexible modeling for turbulent heat transfer. In the context of advanced nuclear reactor design, where challenges include non-unity Prandtl fluids, free and mixed convection, and transitional flow, computational fluid dynamics methods are underutilized due to low confidence in modeling approaches and high computational cost. This report evaluates the performance of engineering CFD models in diabatic, buoyant turbulent flow. It finds that all turbulence models including several variants of the k-ε and k-ω models struggle to predict accurate turbulent momentum and heat transfer in such flows. The Nusselt numbers have been compared between the models and the DNS data, where calculations have been performed for each case and trends have shown a good agreement between DNS estimated Nusselt numbers and available correlations as well as experimental data. A novel DNS correlation for high Pr cases as they are transitioning from mixed convection to natural convection has been proposed. While Nusselt number errors relative to DNS range from 20% to 50%, the models capture similar trends to DNS with respect to Nusselt suppression and amplification under varying levels of buoyancy effect. This report also highlights model form error as a significant contributor to CFD predictions and proposes a framework for quantifying model error and improving confidence in CFD calculations. To improve the predictive capability of engineering CFD models, data-driven approaches for turbulence models are investigated. Theoretical frameworks based on the invariant tensor / vector basis neural networks for prediction of Reynolds stresses and turbulent heat fluxes are employed. The models are developed using direct numerical simulations data for forced convection flows of different fluids in vertical planar channel domain. The framework is implemented in spectral element solvers Nek5000 / nekRS and has shown a potential for future development and consideration of mixed convection flows.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1987

An aerothermodynamic study of slender conical vehicles

A numerical study was performed to assess the applicability of some current techniques which can be used for aerothermal predictions over slender spherically blunted cones. Predictions using a viscous-shock-layer method and several engineering approaches were compared with experimental results from flight and ground-based tests, with each other, and with other detailed results. Good agreement was obtained in comparisons with laminar and turbulent heating data from the Reentry F flight vehicle and with the wind-tunnel data. In particular, the viscous-shock-layer method was shown to yield excellent comparisons and should be useful in providing detailed flowfield and surface values for slender blunted cones. Additional predictions were obtained with these methods for two 5-deg half-angle cones with different nose radii to illustrate the effects of nose bluntness and angle of attack on drag and heat transfer. These results demonstrate the benefit of nose blunting with respect to heating and drag for laminar and transitional flow at zero-lift conditions and the benefit of heating reduction at angle of attack. Detailed comparison of the engineering code predictions with the viscous-shock-layer results for these additional cases generally showed good agreement except for the laminar prodictions at angle of attack on the forward cone surface. Nonequilibrium calculations at 0-deg angle of attack showed that substantial benefits of low surface catalysis existed only in the nose-dominated region.

Thompson, R. A.↗

DOE OSTI2020

Evaluation of a nonlinear variational multiscale method for fluid transport problems

Diverse transport problems, especially those based on fluid flow models, are intrinsically multiscale and nonlinear, characteristics that often lead to intricate dynamics such as the development of instabilities and turbulence. Computational simulations that resolve all scales in these problems are often unfeasible, prompting to coarse-grained simulation strategies in which small-scale features are modeled instead of resolved. Variational Multiscale (VMS) methods, and particularly residual-based Large-Eddy Simulation (LES) approaches, have proven effective and robust for the coarse-grained simulation of complex transport problems. VMS methods avoid the assumption of separable nonlinearity and the reliance on empirical small-scale models by using a variational decomposition of scales together with a residual-based approximation of the small-scales. Evaluation of a nonlinear VMS approach, denoted as VMS n , is presented for the coarse-grained simulation of transient-advective-diffusive-reactive (TADR) transport problems arising from fluid flow models. In contrast to classical VMS approaches that neglect the effect of the small scales on the transport operator, VMS n treats the inter-dependence between large- and small-scales upfront. The treatment of inter-scale coupling involves the solution of a local algebraic nonlinear system describing the evolution of the small-scales. The VMS n approach is complemented with two algebraic approximations of the small-scales: one based on the main diagonal of the transport matrices and another that preserves transport fluxes and is suitable for generic TADR systems. The suitability of the VMS n approach for handling general TADR problems and regimes is evaluated with benchmark incompressible, compressible, and magnetohydrodynamic laminar flow problems, the incompressible Taylor-Green vortex flow, the turbulent free jet, and the two-temperature arc in crossflow. Furthermore, simulation results show that VMS n leads to minor improvements in accuracy with respect to the classical VMS for the laminar flow problems, but to significantly greater accuracy for the turbulent flows and the unsteady plasma flow problems, while using the same cohesive numerical formulation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1984

Aeroacoustics of turbulent shear flows

Recent analytical, numerical-simulation and experimental studies of sound generation by high-Reynolds-number turbulent shear flows are reviewed, with a focus on the application of linear rapid-distortion theory to the calculation of the unsteady flow field producing the sound. This approach is considered the most important alternative to acoustic-analogy methods. Topics surveyed include the linear theory of solid-surface interactions, the jet-noise problem, extensions to more complex turbulent flows, and supersonic flows. Graphs comparing theoretical and experimental results are shown.

Goldstein, M. E.↗

NASA NTRS1992

Comments on the present state and future directions of PDF methods

The one point probability density function (PDF) method is examined in light of its use in actual engineering problems. The PDF method, although relatively complicated, appears to be the only format available to handle the nonlinear stochastic difficulties caused by typical reaction kinetics. Turbulence modeling, if it is to play a central role in combustion modeling, has to be integrated with the chemistry in a way which produces accurate numerical solutions to combustion problems. It is questionable whether the development of turbulent models in isolation from the peculiar statistics of reactant concentrations is a fruitful line of development as far as propulsion is concerned. There are three issues for which additional viewgraphs are prepared: the one point pdf method; the amplitude mapping closure; and a hybrid strategy for replacing a full two point pdf treatment of reacting flows by a single point pdf and correlation functions. An appeal is made for the establishment of an adequate data base for compressible flow with reactions for Mach numbers of unity or higher.

Obrien, E. E.↗

NASA NTRS2023

Fan Noise Predictions of the NASA Source Diagnostic Test Using Unsteady Simulations with LAVA Part I: Near-Field Aerodynamics and Turbulence

A sliding mesh technique within the Launch, Ascent, and Vehicle Aerodynamics (LAVA) computational framework is validated using the experimental dataset collected as part of the NASA Source Diagnostic Test (SDT) campaign. Two modeling approaches are explored: the unsteady Reynolds-Averaged Navier Stokes (URANS) with Spalart-Allmaras (SA) turbulence model closure, and a hybrid Reynolds-Averaged Navier Stokes/Large Eddy Simulation (RANS/LES) paradigm employing a Zonal Detached Eddy Simulation (ZDES) closure with enhanced shielding protection. Fan stage performance metrics, aerodynamic quantities and turbulent flow structures are analyzed in this work. Initial studies focusing on grid and time-step sensitivity are presented. Sensitivity to different variants of the SA turbulence model is analyzed, supporting the use of the baseline SA model in the production runs. Two conditions are analyzed in detail using URANS and hybrid RANS/LES (HRLES). Mean flow quantities are well-captured by both methods in the low-speed (approach) regime. While URANS misses all the upstream-propagating noise in the inlet due to the rotor-locked tones being evanescent in nature at subsonic fan tip speeds, HRLES captures this broadband component in its pressure field. At the high-speed (sideline) condition, URANS shows better agreement with the SDT data than HRLES in the interstage flow-field. In this regime, URANS captures the tonal content propagating through the inlet, since the tones are now cut-on. Both methods are suitable to capture fan stage performance metrics and mean flow quantities, but only HRLES is able to resolve the fine turbulent structures responsible for broadband noise. The results support the use of the sliding mesh technique implemented in this work for future turbomachinery applications within the LAVA solver framework.

AATT↗

NASA NTRS2023

Fan Noise Predictions of the NASA Source Diagnostic Test Using Unsteady Simulations with LAVA Part I: Near-Field Aerodynamics and Turbulence

A sliding mesh technique within the Launch, Ascent, and Vehicle Aerodynamics (LAVA) computational framework is validated using the experimental dataset collected as part of the NASA Source Diagnostic Test (SDT) campaign. Two modeling approaches are explored: the unsteady Reynolds-Averaged Navier Stokes (URANS) with Spalart-Allmaras (SA) turbulence model closure, and a hybrid Reynolds-Averaged Navier Stokes/Large Eddy Simulation (RANS/LES) paradigm employing a Zonal Detached Eddy Simulation (ZDES) closure with enhanced shielding protection. Fan stage performance metrics, aerodynamic quantities and turbulent flow structures are analyzed in this work. Initial studies focusing on grid and time-step sensitivity are presented. Sensitivity to different variants of the SA turbulence model is analyzed, supporting the use of the baseline SA model in the production runs. Two conditions are analyzed in detail using URANS and hybrid RANS/LES (HRLES). Mean flow quantities are well-captured by both methods in the low-speed (approach) regime. While URANS misses all the upstream-propagating noise in the inlet due to the rotor-locked tones being evanescent in nature at subsonic fan tip speeds, HRLES captures this broadband component in its pressure field. At the high-speed (sideline) condition, URANS shows better agreement with the SDT data than HRLES in the interstage flow-field. In this regime, URANS captures the tonal content propagating through the inlet, since the tones are now cut-on. Both methods are suitable to capture fan stage performance metrics and mean flow quantities, but only HRLES is able to resolve the fine turbulent structures responsible for broadband noise. The results support the use of the sliding mesh technique implemented in this work for future turbomachinery applications within the LAVA solver framework.

AATT↗

NASA NTRSDate not supplied

Comparison of Corcos-based and experimentally-derived coherence factors for BFFs estimation

In this paper, high-spatial-resolution unsteady Pressure Sensitive Paint (uPSP) data are utilized to compare two methods for panel Buffet Forcing Functions (BFF) estimation for the Space Launch System (SLS). Such methods are based on discrete pressure measurements within a panel but employ coherence factors to account for partially-correlated fluctuating pressures across the whole panel. In one method, coherence factors are derived based on the Corcos model, whereas the second method utilizes experimentally-derived coherence factors. To simulate discrete measurements using uPSP data, suitable subsets of the data are extracted. When full uPSP resolution is retained, uPSP data provide a benchmark to assess discrete-measurements-based methods. The analysis focuses on the peak SLS buffet environment located downstream of the Forward Attachment Hardware (FAH) between the core stage and solid rocket boosters. Trends of Corcos-based and experimentally-derived coherence factors are in reasonable agreement with the benchmark. However, at certain frequencies, experimentally-derived coherence factors are sensitive to the separation distance between pressure measurements utilized to compute coherence lengths. Such sensitivity originates from deviation of the experimental-based coherence function from an exponential decay assumption. On the other hand, the present implementation of the Corcos model fails to capture certain non-turbulent boundary layer related environments, such as a subharmonic of FAH vortex-shedding. For all methods presented in this paper, at near transonic conditions, increased pressure coherence and spatial nonuniformity lead to BFF overestimation and sensitivity to the pressure measurement location within the panel.

transonic buffet↗

NASA NTRS1990

Flow separation in shock wave boundary layer interactions at hypersonic speeds

An assessment is presented for the experimental data on separated flow in shock wave turbulent boundary layer interactions at hypersonic and supersonic speeds. The data base consists mainly of two dimensional and axisymmetric interactions in compression corners or cylinder-flares, and externally generated oblique shock interactions with boundary layers over flat plates or cylindrical surfaces. The conditions leading to flow separation and the subsequent changes in the flow empirical correlations for incipient separation are reviewed. The effects of the Mach number, Reynolds number, surface cooling and the methods of detecting separation are discussed. The pertinent experimental data for the separated flow characteristics in separated turbulent boundary layer shock interaction are also presented and discussed.

Hamed, A.↗

NASA NTRS1992

Hypersonic flow separation in shock wave boundary layer interactions

An assessment is presented for the experimental data on separated flow in shock wave turbulent boundary layer interactions at hypersonic and supersonic speeds. The data base consists mainly of two dimensional and axisymmetric interactions in compression corners or cylinder-flares, and externally generated oblique shock interactions with boundary layers over flat plates or cylindrical surfaces. The conditions leading to flow separation and the subsequent changes in the flow empirical correlations for incipient separation are reviewed. The effects of the Mach number, Reynolds number, surface cooling and the methods of detecting separation are discussed. The pertinent experimental data for the separated flow characteristics in separated turbulent boundary layer shock interaction are also presented and discussed.

Hamed, A.↗

NASA NTRS2006

A Correlation-Based Transition Model using Local Variables: Test Cases and Industrial Applications - Part 2

A new correlation-based transition model has been developed, which is built strictly on local variables. As a result, the transition model is compatible with modern computational fluid dynamics (CFD) methods using unstructured grids and massive parallel execution. The model is based on two transport equations, one for the intermittency and one for the transition onset criteria in terms of momentum thickness Reynolds number. The proposed transport equations do not attempt to model the physics of the transition process (unlike, e.g., turbulence models), but form a framework for the implementation of correlation-based models into general-purpose CFD methods.

Langtry, R. B.↗

NASA NTRSDate not supplied

Large-Eddy Simulations of a Single-Injector Cooling Flow Using the High-Order Flux Reconstruction Method

A single-injector cooling flow into a heated crossflow was used as a validation case for large eddy simulations (LES) from two high-order CFD codes with comparisons to a state-of-the-art Reynolds averaged Navier-Stokes (RANS) turbulence model. The focus of this paper is on the validation of LES for one of the high-order CFD codes, namely the GFR (Glenn Flux Reconstruction) code that is being developed at NASA Glenn Research Center. Fourth-order LES were performed for the single-injector cooling flow configuration at blowing ratios of both 1.0 and 2.0, and a fifth-order LES was also performed with a blowing ratio of 1.0. The GFR simulations agree very well with the experiment mean quantities and reasonably well with the experiment turbulence quantities. The LES solutions from GFR and the other high-order CFD code, FDL3DI, agreed very closely for nearly every flow quantity examined, despite the fact that these two codes use completely different numerical methods, different grid geometries and domains, and different inflow boundary conditions. Finally, both LES show a significant accuracy improvement over RANS turbulence models for flows with large temperature gradients, where the standard gradient diffusion approximation is insufficient for temperature predictions.

Large Eddy Simulations↗

NASA NTRSDate not supplied

Large Eddy Simulations of a Single-Injector Cooling Flow Using the High-Order Flux Reconstruction Method

A single-injector cooling flow into a heated crossflow was used as a validation case for large eddy simulations (LES) from two high-order CFD codes with comparisons to a state-of-the-art Reynolds averaged Navier-Stokes (RANS) turbulence model. The focus of this paper is on the validation of LES for one of the high-order CFD codes, namely the GFR (Glenn Flux Reconstruction) code that is being developed at NASA Glenn Research Center. Fourth-order LES were performed for the single-injector cooling flow configuration at blowing ratios of both 1.0 and 2.0, and a fifth-order LES was also performed with a blowing ratio of 1.0. The GFR simulations agree very well with the experiment mean quantities and reasonably well with the experiment turbulence quantities. The LES solutions from GFR and the other high-order CFD code, FDL3DI, agreed very closely for nearly every flow quantity examined, despite the fact that these two codes use completely different numerical methods, different grid geometries and domains, and different inflow boundary conditions. Finally, both LES show a significant accuracy improvement over RANS turbulence models for flows with large temperature gradients, where the standard gradient diffusion approximation is insufficient for temperature predictions.

Large Eddy Simulation↗