Results for “Reynolds Averaged Navier Stokes” — Page 2

NASA NTRS2012

Large-Eddy / Reynolds-Averaged Navier-Stokes Simulations of a Dual-Mode Scramjet Combustor

Numerical simulations of reacting and non-reacting flows within a scramjet combustor configuration experimentally mapped at the University of Virginia s Scramjet Combustion Facility (operating with Configuration A ) are described in this paper. Reynolds-Averaged Navier-Stokes (RANS) and hybrid Large Eddy Simulation / Reynolds-Averaged Navier-Stokes (LES / RANS) methods are utilized, with the intent of comparing essentially blind predictions with results from non-intrusive flow-field measurement methods including coherent anti-Stokes Raman spectroscopy (CARS), hydroxyl radical planar laser-induced fluorescence (OH-PLIF), stereoscopic particle image velocimetry (SPIV), wavelength modulation spectroscopy (WMS), and focusing Schlieren. NC State's REACTMB solver was used both for RANS and LES / RANS, along with a 9-species, 19- reaction H2-air kinetics mechanism by Jachimowski. Inviscid fluxes were evaluated using Edwards LDFSS flux-splitting scheme, and the Menter BSL turbulence model was utilized in both full-domain RANS simulations and as the unsteady RANS portion of the LES / RANS closure. Simulations were executed and compared with experiment at two equivalence ratios, PHI = 0.17 and PHI = 0.34. Results show that the PHI = 0.17 flame is hotter near the injector while the PHI = 0.34 flame is displaced further downstream in the combustor, though it is still anchored to the injector. Reactant mixing was predicted to be much better at the lower equivalence ratio. The LES / RANS model appears to predict lower overall heat release compared to RANS (at least for PHI = 0.17), and its capability to capture the direct effects of larger turbulent eddies leads to much better predictions of reactant mixing and combustion in the flame stabilization region downstream of the fuel injector. Numerical results from the LES/RANS model also show very good agreement with OH-PLIF and SPIV measurements. An un-damped long-wave oscillation of the pre-combustion shock train, which caused convergence problems in some RANS simulations, was also captured in LES / RANS simulations, which were able to accommodate its effects accurately.

Fulton, Jesse A.↗

NASA NTRS2021

Transition Analysis for the CRM-NLF Wind Tunnel Configuration

This paper reports the results of a comprehensive linear stability analysis of the boundary layer flow over the common research model with natural laminar flow (CRM-NLF) 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 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. However, a significant scatter ( N ∈ (4,9)) was observed in the values of the logarithmic amplification factor along the measured transition front. This previous analysis was based on an approximate basic state (based on a boundary layer code with conical flow approximation) and parallel stability computations without surface curvature effects. Here, we examine the effects of these various approximations with the goal of quantifying the resulting changes in the N-factor correlations. Specifically, both linear stability theory (LST) and the parabolized stability equations (PSE)are used in conjunction with an accurate definition of the laminar boundary layer flow as computed with a Navier-Stokes solver with a Reynolds-Averaged-Navier-Stokes (RANS) based turbulence model within the turbulent parts of the flow. Furthermore, the effects of instability wave propagation within a fully three-dimensional boundary layer are also evaluated by integrating the disturbance growth rates along suitably chosen, curvilinear (i.e., nonplanar) propagation trajectories. The results of this analysis are also used in an accompanying paper by Venkatachari et al. to develop improved, physics based transition predictions for the same CRM-NLF configuration.

Boundary layer transition↗

DOE OSTI2024

ECAR-7932 Rev 0 Large Eddy Simulation of MARVEL Reactor Core Subchannel to Evaluate Model Uncertainty of Reynolds-Averaged Navier-Stokes Equation Based Computational Fluid Dynamics Analysis

In the previous work (ECAR-7210), the peak cladding temperature of the MARVEL microreactor has been evaluated by steady-state Reynolds-Averaged Navier-Stokes (RANS) based computational fluid dynamics (CFD) simulations. Although numerical uncertainties of RANS-based CFD simulations has been assessed in ECAR-7210, the model uncertainty of RANS turbulence models must be investigated to resolve the issues related to inaccurate prediction of turbulent heat flux and flow pulsation in a tight lattice rod bundle using the steady-state RANS simulations. Consequently, this ECAR conducted a high-fidelity CFD analysis utilizing Large Eddy Simulation (LES) to generate reference data and investigated the model uncertainty of RANS-based CFD simulations.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

NASA NTRS2006

Reynolds-Averaged Navier-Stokes Analysis of Zero Efflux Flow Control Over a Hump Model

The unsteady flow over a hump model with zero efflux oscillatory flow control is modeled computationally using the unsteady Reynolds-averaged Navier-Stokes equations. Three different turbulence models produce similar results, and do a reasonably good job predicting the general character of the unsteady surface pressure coefficients during the forced cycle. However, the turbulent shear stresses are underpredicted in magnitude inside the separation bubble, and the computed results predict too large a (mean) separation bubble compared with experiment. These missed predictions are consistent with earlier steady-state results using no-flow-control and steady suction, from a 2004 CFD validation workshop for synthetic jets.

Rumsey, Christopher L.↗

NASA NTRS2006

Reynolds-Averaged Navier-Stokes Analysis of Zero Efflux Flow Control over a Hump Model

The unsteady flow over a hump model with zero efflux oscillatory flow control is modeled computationally using the unsteady Reynolds-averaged Navier-Stokes equations. Three different turbulence models produce similar results, and do a reasonably good job predicting the general character of the unsteady surface pressure coefficients during the forced cycle. However, the turbulent shear stresses are underpredicted in magnitude inside the separation bubble, and the computed results predict too large a (mean) separation bubble compared with experiment. These missed predictions are consistent with earlier steady-state results using no-flow-control and steady suction, from a 2004 CFD validation workshop for synthetic jets.

Rumsey, Christopher L.↗

NASA NTRS2006

Reynolds-Averaged Navier-Stokes Analysis of Zero Efflux Flow Control over a Hump Model

The unsteady flow over a hump model with zero efflux oscillatory flow control is modeled computationally using the unsteady Reynolds-averaged Navier-Stokes equations. Three different turbulence models produce similar results, and do a reasonably good job predicting the general character of the unsteady surface pressure coefficients during the forced cycle. However, the turbulent shear stresses are underpredicted in magnitude inside the separation bubble, and the computed results predict too large a (mean) separation bubble compared with experiment. These missed predictions are consistent with earlier steady-state results using no-flow-control and steady suction, from a 2004 CFD validation workshop for synthetic jets.

Rumsey, Christopher L.↗

NASA NTRSDate not supplied

A Reynolds-Averaged Navier-Stokes Perspective for the High Lift Common Research Model Using the LAVA Framework

An in-depth assessment of the High Lift-Common Research Model (HL-CRM) aerodynamic performance is made using a Reynolds-Averaged Navier-Stokes (RANS) methodology. Results presented here were produced in support of the 4th AIAA High Lift Prediction Workshop(HLPW-4), with additional simulation results and analysis used to support findings from the workshop. Mesh resolution, RANS closure models, and tunnel effects were among the sources of simulation sensitivity addressed, all in the pursuit of evaluating best-practice simulation methods for complex high-lift configurations. In general, free-air simulations using RANS yield accurate predictions of integrated loads in the pre-stall regime for the HL-CRM, while strong sensitivity to simulation method is observed near CLmax. Additionally, a study assessing changes in aerodynamic performance due to flap deflection angle indicated strong sensitivity to turbulence model closure and mesh resolution, particularly at the highest deflection setting. Several available corrections for the Spalart-Allmaras (SA) turbulence model were tested and generally introduced more error relative to the experiment in comparison to the baseline SA closure model. The Standard Menter Baseline Two-Equation (k-ωBSL) turbulence model was also explored, which arguably yielded results highly consistent with experiment, but at a substantially higher computational cost. Therefore, as a compromise between computational cost and accuracy, the baseline SA closure model was selected for additional HL-CRM studies that assess alternative solution methods. These strategies include “warm start” approaches, which initialize the simulation from a previously converged solution at a nearby angle of attack, unsteady calculations, and simulations that model the experimental wind tunnel. Results from these studies would reinforce the shortcomings of RANS in predicting aircraft performance at high-lift conditions, namely at high-αand highly loaded flap conditions, while also answering questions about appropriate modeling techniques to minimize the errors associated with RANS where possible

TTT↗

NASA NTRS1982

Reynolds averaged Navier-Stokes computations of transonic flows The-state-of-the-art

The state of the art in computer simulations for transonic flowfields requiring solutions for the Navier-Stokes equations is assessed. It is noted that current simulations of transonic flowfields require comparisons with experimental results because the simulations are not free from discretization errors. Cases of turbulence are treated with weighted variables in a time-averaged scheme to yield Reynolds averaged Navier-Stokes equations. The turbulence is modeled in a first-order approach with a Reynolds stress tensor or a second-order approach where the tensor is obtained from the Navier-Stokes equations. The choice of coordinate systems and dependent variables for simulating the flow around airfoils is discussed, with particular attention to curvilinear coordinates. The determinations of boundary conditions is examined, along with numerical methods related to physical phenomena

Mehta, U.↗

NASA NTRS2007

Parametric Study of Flow Control Over a Hump Model Using an Unsteady Reynolds- Averaged Navier-Stokes Code

This is an expanded version of a limited-length paper that appeared at the 5th International Symposium on Turbulence and Shear Flow Phenomena by the same authors. A computational study was performed for steady and oscillatory flow control over a hump model with flow separation to assess how well the steady and unsteady Reynolds-averaged Navier-Stokes equations predict trends due to Reynolds number, control magnitude, and control frequency. As demonstrated in earlier studies, the hump model case is useful because it clearly demonstrates a failing in all known turbulence models: they under-predict the turbulent shear stress in the separated region and consequently reattachment occurs too far downstream. In spite of this known failing, three different turbulence models were employed to determine if trends can be captured even though absolute levels are not. Overall the three turbulence models showed very similar trends as experiment for steady suction, but only agreed qualitatively with some of the trends for oscillatory control.

Rumsey, Christopher L.↗

NASA NTRS2012

Performance of Reynolds Averaged Navier-Stokes Models in Predicting Separated Flows: Study of the Hump Flow Model Problem

Separation can be seen in most aerodynamic flows, but accurate prediction of separated flows is still a challenging problem for computational fluid dynamics (CFD) tools. The behavior of several Reynolds Averaged Navier-Stokes (RANS) models in predicting the separated ow over a wall-mounted hump is studied. The strengths and weaknesses of the most popular RANS models (Spalart-Allmaras, k-epsilon, k-omega, k-omega-SST) are evaluated using the open source software OpenFOAM. The hump ow modeled in this work has been documented in the 2004 CFD Validation Workshop on Synthetic Jets and Turbulent Separation Control. Only the baseline case is treated; the slot flow control cases are not considered in this paper. Particular attention is given to predicting the size of the recirculation bubble, the position of the reattachment point, and the velocity profiles downstream of the hump.

Cappelli, Daniele↗

DOE OSTI2022

Large-eddy simulation and Reynolds-averaged Navier-Stokes modeling of three Rayleigh-Taylor mixing configurations with gravity reversal

In this work, high-fidelity large-eddy simulation (LES) is performed of Rayleigh-Taylor (RT) mixing in three different configurations involving gravity reversal. In each configuration, LES results are compared with one-dimensional Reynolds-averaged Navier-Stokes (RANS) results, and a deficiency in a commonly used transport equation for the mass-flux velocity, $a_j$, is identified. In the first configuration, a classical two-component RT mixing layer is allowed to develop before it is subjected to rapid acceleration reversal. In the second configuration, a three-component RT mixing layer with an intermediate density layer is allowed to develop before being subjected to rapid acceleration reversal. Finally, in the third configuration, a light layer is interposed between two heavy layers; in this configuration, only one interface is RT-unstable at a time as it undergoes rapid acceleration reversal. In all cases, a commonly used buoyancy production closure in the $a_j$ transport equation is shown to lead to significant over-prediction of mixing layer growth after gravity reversal. An alternative formulation for this closure is then presented which is shown to more accurately capture the stabilization effect of gravity reversal.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Large Eddy Simulation of Low-Reynolds-Number Turbulent Flow of Low-Prandtl-Number Fluid in a Tight Lattice Bundle for Assessment of Reynolds-Averaged Navier-Stokes Turbulence Model

The MARVEL (Microreactor Applications Research Validation and Evaluation) microreactor utilizes natural circulation as core cooling mechanism and liquid metal as a primary coolant. Moreover, the reactor core has a pitch-to-diameter ratio of 1.056, which is considered a tight lattice configuration. Numerous studies have widely reported that Reynolds-Averaged Navier-Stokes (RANS) turbulence models inaccurately predict heat transfer in liquid metals and fail to capture flow pulsations that can occur within tight lattices, leading to further inaccuracies in simulation results. Therefore, evaluating the accuracy of RANS turbulence models in the thermal-hydraulic analysis of the MARVEL microreactor core is crucial for assessing reactor safety. In this study, a Large Eddy Simulation (LES) of the MARVEL microreactor core subchannel was conducted and compared with a RANS simulation to evaluate the accuracies and conservatism of the RANS model. The flow pulsation in a tight lattice predicted by LES enhanced the heat transfer, whereas the RANS model underpredicted it. Consequently, the RANS model predicted the peak cladding temperature higher than the LES model, but the discrepancy between the two approaches was not significant due to the good thermal characteristics of the liquid metal. It can be concluded that the steady-state RANS model is effective for the thermal analysis of liquid-metal-cooled MARVEL microreactor core and can provide conservative predictions from a safety analysis perspective.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

DOE OSTI2021

Scalar mixing in a Kelvin-Helmholtz shear layer and implications for Reynolds-averaged Navier-Stokes modeling of mixing layers

Large-eddy simulation of a temporally evolving Kelvin-Helmholtz (KH) mixing layer is performed with the tenth-order compact difference code miranda to examine the steady-state behavior of a passive scalar in a shear-driven mixing layer. It is shown that the integral behavior of scalar variance in a KH mixing layer behaves similarly to the integral behavior of scalar variance in a Rayleigh-Taylor (RT) mixing layer, and mixedness of the simulated KH shear layer tends towards a value of about 0.8. It is further shown that if the k-L-a-V Reynolds-averaged Navier-Stokes (RANS) model [B. E. Morgan et al., Phys. Rev. E 98, 033111 (2018)], calibrated to reproduce steady-state mixing in an RT layer, is applied to simulate a KH mixing layer, the RANS model will significantly overpredict the magnitude of scalar variance in the KH layer. A straightforward addition to the k-L-a-V model is then suggested, and self-similarity analysis is applied to determine constraints on model coefficients. Furthermore, it is shown that with the addition of a buoyancy production term in the model equation for scalar variance, it becomes possible to eliminate the model deficiency and match steady-state mixedness in simulations of both RT and KH mixing layers with a single model calibration.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Simulation and Reynolds-averaged Navier-Stokes modeling of a three-component Rayleigh-Taylor mixing problem with thermonuclear burn

Rayleigh-Taylor mixing in the presence of a third component with intermediate density is investigated through three-dimensional large-eddy simulation (LES) with a high-order compact finite-difference code. Two configurations are considered: (1) a symmetric configuration in which the Atwood number between the heavy and intermediate components matches the Atwood number between the intermediate and light components and (2) an asymmetric configuration in which the Atwood number between the heavy and intermediate components is an order of magnitude greater than the Atwood number between the intermediate and light components. Mass fraction covariances are extracted, and proposed Reynolds-averaged Navier-Stokes (RANS) closures for density-specific-volume and density-mass-fraction covariances are evaluated in an a priori fashion. Additionally, a multicomponent extension of the k - Φ - L - a - V RANS model [Morgan, Phys. Rev. E 104, 015107 (2021)] is presented which includes model equations for the upper-triangular elements of the mass fraction covariance matrix. This model, referred to as the k - Φ - L - a - C model, is compared against results from LES and against other RANS models. Profiles of average mass fraction, mass-fraction covariance, and density-specific-volume covariance obtained with the k - Φ - L - a - C model are found to agree well with LES data. Finally, the impact of three-component turbulent mixing on average reaction rate is investigated in both premixed and nonpremixed cases by heating the mixing layer and allowing it to undergo thermonuclear (TN) burn. A closure model for average reaction rate is proposed for use with the k - Φ - L - a - C model, and when this model is applied, improved agreement is obtained between LES and RANS in total TN neutron production.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2023

Measurements and Computations of Natural Transition on the NASA Juncture-Flow Model with a Symmetric Wing

Experiments were performed in the 14- by 22-Foot Subsonic Tunnel to assess natural transition on the symmetric-airfoil wings of the NASA Juncture-Flow Model. Infrared thermography was used to visualize the heating on the upper surface of both wings of the full-span model, and on the fuselage, for angles of incidence ranging from -10° to 10° at a fixed Reynolds number of 2.4E6 based on the chord length at the wing planform break. The fuselage boundary layer transitioned well upstream of the wing-root leading edge for all conditions. Transition fronts were identified by a steep rise in the surface temperature, and the transition coordinates were transformed from an image-based to a body-fixed system. Additionally, the state of the boundary layer was estimated at pressure ports distributed on the wings through observation of the pressure coefficient as a function of the angle of incidence. For increasing angles of incidence, the transition front was observed to advance upstream, in a mostly spanwise-uniform fashion, from near midchord at α = 0°; however, for increasingly negative angles of incidence, the transition front first receded and then advanced in a nonuniform jagged manner that is typically observed with stationary crossflow. The transition wedges first appeared inboard of the wing break and then spread outboard to near the tip by α = -6°. The upstream shift in transition at positive angles of incidence and the outboard progression of crossflow-dominated transition at increasingly negative angles of incidence are consistent with trends identified in a computational assessment of the boundary-layer transition based on both linear stability analysis and Reynolds-averaged-Navier-Stokes-based transition models. The stability results obtained from the Langley Stability and Transition Analysis Code were used to recalibrate a dual N-factor criterion, which allowed for the prediction of transition fronts that showed excellent agreement with the experiment. The Reynolds-averaged-Navier-Stokes-based models, from the NASA OVERFLOW 2.3 solver, that accounted for the crossflow instability showed mixed results in comparison with the experiment, with the helicity-based Langtry-Menter model performing the best. The experimental data, particularly the cases involving strong influence from both Tollmien-Schlichting and crossflow instabilities, will be valuable for the continued validation and improvement of transition models.

Boundary layer transition↗

NASA NTRS2023

Measurements and Computations of Natural Transition on the NASA Juncture-Flow Model with a Symmetric Wing

Experiments were performed in the 14- by 22-Foot Subsonic Tunnel to assess natural transition on the symmetric-airfoil wings of the NASA Juncture-Flow Model. Infrared thermography was used to visualize the heating on the upper surface of both wings of the full-span model, and on the fuselage, for angles of incidence ranging from -10° to 10° at a fixed Reynolds number of 2.4E6 based on the chord length at the wing planform break. The fuselage boundary layer transitioned well upstream of the wing-root leading edge for all conditions. Transition fronts were identified by a steep rise in the surface temperature, and the transition coordinates were transformed from an image-based to a body-fixed system. Additionally, the state of the boundary layer was estimated at pressure ports distributed on the wings through observation of the pressure coefficient as a function of the angle of incidence. For increasing angles of incidence, the transition front was observed to advance upstream, in a mostly spanwise-uniform fashion, from near midchord at α = 0°; however, for increasingly negative angles of incidence, the transition front first receded and then advanced in a nonuniform jagged manner that is typically observed with stationary crossflow. The transition wedges first appeared inboard of the wing break and then spread outboard to near the tip by α = -6°. The upstream shift in transition at positive angles of incidence and the outboard progression of crossflow-dominated transition at increasingly negative angles of incidence are consistent with trends identified in a computational assessment of the boundary-layer transition based on both linear stability analysis and Reynolds-averaged-Navier-Stokes-based transition models. The stability results obtained from the Langley Stability and Transition Analysis Code were used to recalibrate a dual N-factor criterion, which allowed for the prediction of transition fronts that showed excellent agreement with the experiment. The Reynolds-averaged-Navier-Stokes-based models, from the NASA OVERFLOW 2.3 solver, that accounted for the crossflow instability showed mixed results in comparison with the experiment, with the helicity-based Langtry-Menter model performing the best. The experimental data, particularly the cases involving strong influence from both Tollmien-Schlichting and crossflow instabilities, will be valuable for the continued validation and improvement of transition models.

boundary layer transition↗

NASA NTRS2023

Measurements and Computations of Natural Transition on the NASA Juncture-Flow Model with a Symmetric Wing

Experiments were performed in the 14- by 22-Foot Subsonic Tunnel to assess natural transition on the symmetric-airfoil wings of the NASA Juncture-Flow Model. Infrared thermography was used to visualize the heating on the upper surface of both wings of the full-span model, and on the fuselage, for angles of incidence ranging from -10° to 10° at a fixed Reynolds number of 2.4E6 based on the chord length at the wing planform break. The fuselage boundary layer transitioned well upstream of the wing-root leading edge for all conditions. Transition fronts were identified by a steep rise in the surface temperature, and the transition coordinates were transformed from an image-based to a body-fixed system. Additionally, the state of the boundary layer was estimated at pressure ports distributed on the wings through observation of the pressure coefficient as a function of the angle of incidence. For increasing angles of incidence, the transition front was observed to advance upstream, in a mostly spanwise-uniform fashion, from near midchord at α = 0°; however, for increasingly negative angles of incidence, the transition front first receded and then advanced in a nonuniform jagged manner that is typically observed with stationary crossflow. The transition wedges first appeared inboard of the wing break and then spread outboard to near the tip by α = -6°. The upstream shift in transition at positive angles of incidence and the outboard progression of crossflow-dominated transition at increasingly negative angles of incidence are consistent with trends identified in a computational assessment of the boundary-layer transition based on both linear stability analysis and Reynolds-averaged-Navier-Stokes-based transition models. The stability results obtained from the Langley Stability and Transition Analysis Code were used to recalibrate a dual N-factor criterion, which allowed for the prediction of transition fronts that showed excellent agreement with the experiment. The Reynolds-averaged-Navier-Stokes-based models, from the NASA OVERFLOW 2.3 solver, that accounted for the crossflow instability showed mixed results in comparison with the experiment, with the helicity-based Langtry-Menter model performing the best. The experimental data, particularly the cases involving strong influence from both Tollmien-Schlichting and crossflow instabilities, will be valuable for the continued validation and improvement of transition models.

boundary layer transition↗

NASA NTRS2014

On Parametric Sensitivity of Reynolds-Averaged Navier-Stokes SST Turbulence Model: 2D Hypersonic Shock-Wave Boundary Layer Interactions

Examined is sensitivity of separation extent, wall pressure and heating to variation of primary input flow parameters, such as Mach and Reynolds numbers and shock strength, for 2D and Axisymmetric Hypersonic Shock Wave Turbulent Boundary Layer interactions obtained by Navier-Stokes methods using the SST turbulence model. Baseline parametric sensitivity response is provided in part by comparison with vetted experiments, and in part through updated correlations based on free interaction theory concepts. A recent database compilation of hypersonic 2D shock-wave/turbulent boundary layer experiments extensively used in a prior related uncertainty analysis provides the foundation for this updated correlation approach, as well as for more conventional validation. The primary CFD method for this work is DPLR, one of NASA's real-gas aerothermodynamic production RANS codes. Comparisons are also made with CFL3D, one of NASA's mature perfect-gas RANS codes. Deficiencies in predicted separation response of RANS/SST solutions to parametric variations of test conditions are summarized, along with recommendations as to future turbulence approach.

Hypersonic↗