SEARCH · Search NASA

Results for “TWO EQUATION MODEL”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Three-Dimensional Navier-Stokes Method with Two-Equation Turbulence Models for Efficient Numerical Simulation of Hypersonic Flows

A new computational efficient 3-D compressible Reynolds-averaged implicit Navier-Stokes method with advanced two equation turbulence models for high speed flows is presented. All convective terms are modeled using an entropy satisfying higher-order Total Variation Diminishing (TVD) scheme based on implicit upwind flux-difference split approximations and arithmetic averaging procedure of primitive variables. This method combines the best features of data management and computational efficiency of space marching procedures with the generality and stability of time dependent Navier-Stokes procedures to solve flows with mixed supersonic and subsonic zones, including streamwise separated flows. Its robust stability derives from a combination of conservative implicit upwind flux-difference splitting with Roe's property U to provide accurate shock capturing capability that non-conservative schemes do not guarantee, alternating symmetric Gauss-Seidel 'method of planes' relaxation procedure coupled with a three-dimensional two-factor diagonal-dominant approximate factorization scheme, TVD flux limiters of higher-order flux differences satisfying realizability, and well-posed characteristic-based implicit boundary-point a'pproximations consistent with the local characteristics domain of dependence. The efficiency of the method is highly increased with Newton Raphson acceleration which allows convergence in essentially one forward sweep for supersonic flows. The method is verified by comparing with experiment and other Navier-Stokes methods. Here, results of adiabatic and cooled flat plate flows, compression corner flow, and 3-D hypersonic shock-wave/turbulent boundary layer interaction flows are presented. The robust 3-D method achieves a better computational efficiency of at least one order of magnitude over the CNS Navier-Stokes code. It provides cost-effective aerodynamic predictions in agreement with experiment, and the capability of predicting complex flow structures in complex geometries with good accuracy.

Bardina, J. E.↗

NASA NTRS1996

Prediction of Airfoil Characteristics With Higher Order Turbulence Models

This study focuses on the prediction of airfoil characteristics, including lift and drag over a range of Reynolds numbers. Two different turbulence models, which represent two different types of models, are tested. The first is a standard isotropic eddy-viscosity two-equation model, and the second is an explicit algebraic stress model (EASM). The turbulent flow field over a general-aviation airfoil (GA(W)-2) at three Reynolds numbers is studied. At each Reynolds number, predicted lift and drag values at different angles of attack are compared with experimental results, and predicted variations of stall locations with Reynolds number are compared with experimental data. Finally, the size of the separation zone predicted by each model is analyzed, and correlated with the behavior of the lift coefficient near stall. In summary, the EASM model is able to predict the lift and drag coefficients over a wider range of angles of attack than the two-equation model for the three Reynolds numbers studied. However, both models are unable to predict the correct lift and drag behavior near the stall angle, and for the lowest Reynolds number case, the two-equation model did not predict separation on the airfoil near stall.

Gatski, Thomas B.↗

NASA NTRS1999

Heat Transfer on a Film-Cooled Rotating Blade Using Different Turbulence Models

A three-dimensional Navier Stokes code has been used to compute the heat transfer coefficient on a film-cooled, rotating turbine blade. The blade chosen is the ACE rotor with five rows containing 93 film cooling holes covering the entire span. This is the only film-cooled rotating blade over which experimental data is available for comparison. Over 2.278 million grid points are used to compute the flow over the blade including the tip clearance region. using Wilcox's k-omega model, Coakley's q-omega model, and the zero-equation Baldwin-Lomax (B-L) model. A reasonably good comparison with the experimental data is obtained on the suction surface for all the turbulence models. At the leading edge, the B-L model yields a better comparison than tile two-equation models. On the pressure surface however the comparison between the experimental data and the prediction from the k-omega model is much better than from the other two models. Overall, the k-omega model provides the best comparison with the experimental data. However, the two-equation models require at least 40% more computational resources than the B-L model.

Garg, Vijay K.↗

NASA NTRS1989

Transonic viscous flow calculations for a turbine cascade with a two equation turbulence model

A numerical method for the study of steady, transonic, turbulent viscous flow through plane turbine cascades is presented. The governing equations are written in Favre-averaged form and closed with a first order model. The turbulent quantities are expressed according to a two-equation kappa-epsilon model where low Reynolds number and compressibility effects are included. The solution is obtained by using a pseudo-unsteady method with improved perturbation propagation properties. The equations are discretized in space by using a finite volume formulation. An explicit multistage dissipative Runge-Kutta algorithm is then used to advance the flow equations in the pseudo-time. First results of calculations compare fairly well with experimental data.

Boretti, A. A.↗

NASA NTRS2000

Implementation of Advanced Two Equation Turbulence Models in the USM3D Unstructured Flow Solver

USM3D is a widely-used unstructured flow solver for simulating inviscid and viscous flows over complex geometries. The current version (version 5.0) of USM3D, however, does not have advanced turbulence models to accurately simulate complicated flow. We have implemented two modified versions of the original Jones and Launder k-epsilon "two-equation" turbulence model and the Girimaji algebraic Reynolds stress model in USM3D. Tests have been conducted for three flat plate boundary layer cases, a RAE2822 airfoil and an ONERA M6 wing. The results are compared with those from direct numerical simulation, empirical formulae, theoretical results, and the existing Spalart-Allmaras one-equation model.

Wang, Qun-Zhen↗

NASA NTRS1994

An analysis of supersonic flows with low-Reynolds number compressible two-equation turbulence models using LU finite volume implicit numerical techniques

A generalized flow solver using an implicit Lower-upper (LU) diagonal decomposition based numerical technique has been coupled with three low-Reynolds number kappa-epsilon models for analysis of problems with engineering applications. The feasibility of using the LU technique to obtain efficient solutions to supersonic problems using the kappa-epsilon model has been demonstrated. The flow solver is then used to explore limitations and convergence characteristics of several popular two equation turbulence models. Several changes to the LU solver have been made to improve the efficiency of turbulent flow predictions. In general, the low-Reynolds number kappa-epsilon models are easier to implement than the models with wall-functions, but require much finer near-wall grid to accurately resolve the physics. The three kappa-epsilon models use different approaches to characterize the near wall regions of the flow. Therefore, the limitations imposed by the near wall characteristics have been carefully resolved. The convergence characteristics of a particular model using a given numerical technique are also an important, but most often overlooked, aspect of turbulence model predictions. It is found that some convergence characteristics could be sacrificed for more accurate near-wall prediction. However, even this gain in accuracy is not sufficient to model the effects of an external pressure gradient imposed by a shock-wave/ boundary-layer interaction. Additional work on turbulence models, especially for compressibility, is required since the solutions obtained with base line turbulence are in only reasonable agreement with the experimental data for the viscous interaction problems.

Lee, J.↗

NASA NTRS1996

Transonic Turbulent Flow Predictions With Two-Equation Turbulence Models

Solutions of the Favre-averaged Navier-Stokes equations for two well-documented transonic turbulent flows are compared in detail with existing experimental data. While the boundary layer in the first case remains attached, a region of extensive flow separation has been observed in the second case. Two recently developed k-epsilon, two-equation, eddy-viscosity models are used to model the turbulence field. These models satisfy the realizability constraints of the Reynolds stresses. Comparisons with the measurements are made for the wall pressure distribution, the mean streamwise velocity profiles, and turbulent quantities. Reasonably good agreement is obtained with the experimental data.

Liou, William W.↗

NASA NTRS1994

Modification of the Two-equation Turbulence Model in NPARC to a Chien Low Reynolds Number K-epsilon Formulation

This report documents the changes that were made to the two-equation k-epsilon turbulence model in the NPARC (National-PARC) code. The previous model based on the low Reynolds number model of Speziale, was replaced with the low Reynolds number k-epsilon model of Chien. The most significant difference was in the turbulent Prandtl numbers appearing in the diffusion terms of the k and epsilon transport equations. A new inflow boundary condition and stability enhancements were also implemented into the turbulence model within NPARC. The report provides the rationale for making the change to the Chien model, code modifications required, and comparisons of the performances of the new model with the previous k-epsilon model and algebraic models used most often in PARC/NPARC. The comparisons show that the Chien k-epsilon model installed here improves the capability of NPARC to calculate turbulent flows.

Georgiadis, Nicholas J.↗

NASA NTRS1985

Diagnosis of ageostrophic circulations in a two-dimensional primitive equation model of frontogenesis

A two-dimensional primitive equation frontogenesis model is applied to the identification of dynamical mechanisms which contribute to transverse ageostrophic circulation in warm and cold advection situations. Attention is focused on conditions conducive to significant alongfront ageostrophic flow with divergence confined to the transverse plane. The circulation is partitioned as a function of forcing mechanisms accompanied by confluence, horizontal shear and the alongfront component of the ageostrophic wind. A series of diagnostic equations is defined for the transverse geostrophic circulation and instantaneous contribution from each partitioned component of the circulation are quantified. In the cold advection case, results are obtained which are consistent with positive feedback between the subsidence pattern and forcing by horizontal shear.

Keyser, D.↗

NASA NTRS1992

A compressible near-wall turbulence model for boundary layer calculations

A compressible near-wall two-equation model is derived by relaxing the assumption of dynamical field similarity between compressible and incompressible flows. This requires justifications for extending the incompressible models to compressible flows and the formulation of the turbulent kinetic energy equation in a form similar to its incompressible counterpart. As a result, the compressible dissipation function has to be split into a solenoidal part, which is not sensitive to changes of compressibility indicators, and a dilational part, which is directly affected by these changes. This approach isolates terms with explicit dependence on compressibility so that they can be modeled accordingly. An equation that governs the transport of the solenoidal dissipation rate with additional terms that are explicitly dependent on the compressibility effects is derived similarly. A model with an explicit dependence on the turbulent Mach number is proposed for the dilational dissipation rate. Thus formulated, all near-wall incompressible flow models could be expressed in terms of the solenoidal dissipation rate and straight-forwardly extended to compressible flows. Therefore, the incompressible equations are recovered correctly in the limit of constant density. The two-equation model and the assumption of constant turbulent Prandtl number are used to calculate compressible boundary layers on a flat plate with different wall thermal boundary conditions and free-stream Mach numbers. The calculated results, including the near-wall distributions of turbulence statistics and their limiting behavior, are in good agreement with measurements. In particular, the near-wall asymptotic properties are found to be consistent with incompressible behavior; thus suggesting that turbulent flows in the viscous sublayer are not much affected by compressibility effects.

So, R. M. C.↗

NASA NTRS1992

Hypersonic flows as related to the national aerospace plane

The study in the last 6 months has observed a clear evidence that the current two-equation models tend to under-predict flow separation and over-predict heat transfer rate near flow re-attachment regions. In hypersonic flow calculations, these model deficiencies appear to be even more pronounced. This is particularly true in the incapability of the model to predict the extent of the flow separation. Two major deficiencies of the current two-equation models in predicting complex hypersonic flows have been reported, i.e., under-prediction of flow separation and over-prediction of peak heat transfer rate. Two modifications to the k - epsilon model were reported and tested over a range of flows. Based on our limited study, the modified models have been found to give better agreements in both surface pressure and heat transfer predictions for several complex shock-wave boundary-layer interaction flows. However, in order to confirm our observation, more calculations will be performed in the future study covering a wider range of flows and conditions than reported here.

Kussoy, Marvin↗

NASA NTRS1993

Numerical simulation of crossing/turbulent boundary layer interaction at Mach 8.3 comparison of zero and two-equation turbulence models

A 3D hypersonic crossing shock wave/turbulent boundary layer interaction is examined numerically. The test geometry consists of a pair of opposing sharp fins of angle alpha = 15 deg mounted on a flat plate. The freestream Mach number is 8.28. Two theoretical models are evaluated. The full 3D Reynolds-averaged Navier-Stokes equations are solved using the Baldwin-Lomax algebraic turbulent eddy viscosity model and the Rodi turbulence model. Computed results for both cases show good agreement with experiment for flat plate surface pressure and for pitot pressure and yaw angle profiles in the flowfield. General agreement is obtained for surface flow direction. Fair to poor agreement is obtained for surface heat transfer, indicating a need for more accurate turbulence models. The overall flowfield structure is similar to that observed in previous crossing shock interaction studies.

Narayanswami, N.↗

NASA NTRS1986

Incompressible viscous flow simulations of the NFAC wind tunnel

The capabilities of an existing 3-D incompressible Navier-Stokes flow solver, INS3D, are extended and improved to solve turbulent flows through the incorporation of zero- and two-equation turbulence models. The two-equation model equations are solved in their high Reynolds number form and utilize wall functions in the treatment of solid wall boundary conditions. The implicit approximate factorization scheme is modified to improve the stability of the two-equation solver. Applications to the 3-D viscous flow inside the 80 by 120 feet open return wind tunnel of the National Full Scale Aerodynamics Complex (NFAC) are discussed and described.

Champney, Joelle Milene↗

NASA NTRS2007

Small Scale Response and Modeling of Periodically Forced Turbulence

The response of the small scales of isotropic turbulence to periodic large scale forcing is studied using two-point closures. The frequency response of the turbulent kinetic energy and dissipation rate, and the phase shifts between production, energy and dissipation are determined as functions of Reynolds number. It is observed that the amplitude and phase of the dissipation exhibit nontrivial frequency and Reynolds number dependence that reveals a filtering effect of the energy cascade. Perturbation analysis is applied to understand this behavior which is shown to depend on distant interactions between widely separated scales of motion. Finally, the extent to which finite dimensional models (standard two-equation models and various generalizations) can reproduce the observed behavior is discussed.

Bos, Wouter↗

NASA NTRS1993

Turbulence modeling for complex hypersonic flows

The paper presents results of calculations for a range of 2D turbulent hypersonic flows using two-equation models. The baseline models and the model corrections required for good hypersonic-flow predictions will be illustrated. Three experimental data sets were chosen for comparison. They are: (1) the hypersonic flare flows of Kussoy and Horstman, (2) a 2D hypersonic compression corner flow of Coleman and Stollery, and (3) the ogive-cylinder impinging shock-expansion flows of Kussoy and Horstman. Comparisons with the experimental data have shown that baseline models under-predict the extent of flow separation but over-predict the heat transfer rate near flow reattachment. Modifications to the models are described which remove the above-mentioned deficiencies. Although we have restricted the discussion only to the selected baseline models in this paper, the modifications proposed are universal and can in principle be transferred to any existing two-equation model formulation.

Huang, P. G.↗

NASA NTRS2018

Development and Documentation of kL-Based Linear, Nonlinear, and Full Reynolds Stress Turbulence Models

The development and implementation of kL-based Reynolds Average Navier-Stokes (RANS) turbulence models are reported herein. The kL is based on Abdol-Hamid's closure and Menter's modi cation to Rotta's two-equation model. Rotta shows that a reliable transport equation can be formed from the turbulent length scale L, and the turbulent kinetic energy k. Rotta's kL equation is well suited for term-by-term modeling and displays useful features compared to other scale formulation. One of the important di erences is the inclusion of higher order velocity derivatives in the source terms of the scale equation. This can enhance the ability of RANS solvers to simulate unsteady ows in URANS mode. The present report documents the formulation of three model levels of turbulence models as implemented in the CFD code FUN3D. Methodology and calibration examples are shown in detail. The levels are the linear k-kL and the two-equation algebraic Reynolds stress model (ARSM) as well as the full Reynolds Stress Model (RSM). Attached, separated and corner ow cases are documented and compared with experimental, theoretical and other turbulence model data. The results show generally very good comparisons with canonical and experimental data. The results from this formulation are similar or better than results using the SST two- equation turbulence model. ARSM shows great promise with similar level of computational resources as general two equation turbulence models.

Abdol-Hamid, Khaled S.↗

NASA NTRS2006

PAB3D: Its History in the Use of Turbulence Models in the Simulation of Jet and Nozzle Flows

This is a review paper for PAB3D s history in the implementation of turbulence models for simulating jet and nozzle flows. We describe different turbulence models used in the simulation of subsonic and supersonic jet and nozzle flows. The time-averaged simulations use modified linear or nonlinear two-equation models to account for supersonic flow as well as high temperature mixing. Two multiscale-type turbulence models are used for unsteady flow simulations. These models require modifications to the Reynolds Averaged Navier-Stokes (RANS) equations. The first scheme is a hybrid RANS/LES model utilizing the two-equation (k-epsilon) model with a RANS/LES transition function, dependent on grid spacing and the computed turbulence length scale. The second scheme is a modified version of the partially averaged Navier-Stokes (PANS) formulation. All of these models are implemented in the three-dimensional Navier-Stokes code PAB3D. This paper discusses computational methods, code implementation, computed results for a wide range of nozzle configurations at various operating conditions, and comparisons with available experimental data. Very good agreement is shown between the numerical solutions and available experimental data over a wide range of operating conditions.

Abdol-Hamid, Khaled S.↗

NASA NTRS1991

Near-wall modeling of the dissipation rate equation

Near-wall modeling of the dissipation rate equation is investigated and its asymptotic behavior is studied in detail using a k-epsilon model. It is found that all existing modeled dissipation rate equations predict an incorrect behavior for the dissipation rate near a wall. An improvement is proposed and the resulting near-wall dissipation rate distribution is found to be similar to that given by numerical simulation data. To further validate the improved k-epsilon model, it is used to calculate flat-plate turbulent boundary-layer flows at high- as well as low-turbulence Reynolds numbers, and the results are compared with measurements, numerical simulation data, and the calculations of three different two-equation models. These comparisons show that all the models tested give essentially the same flow properties away from the wall; significant differences only occur in a region very close to the wall. In this region, the calculations of the improved k-epsilon model are in better agreement with measurements and numerical simulation data. In particular, the modeled distribution of the dissipation rate is significantly improved and a maximum is predicted at the wall instead of away from the wall. Furthermore, the improved k-epsilon model is found to be the most asymptotically consistent among the four different two-equation models examined.

So, R. M. C.↗