Search NASA⌕ Search

SEARCH · Search NASA

Results for “Wall functions”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Computational Analysis of Water Film Thickness During Rain Events for Assessing Hydroplaning Risk Part 2: Rough Road Surfaces

Hydroplaning occurs when a water film forms on a road during a rain event and vehicles are traveling at a speed that does not provide sufficient time for the tires to push the water film out of the tire path. Under these conditions, the tire loses contact with the road and the driver may lose control over the vehicle. The present study is a part of a multiyear research effort aimed at gaining a better understanding of the factors that contribute to hydroplaning risk in order to minimize the occurrence of hydroplaning accidents. Three-dimensional computational fluid dynamics (CFD) analysis was used to investigate the water film thickness on multilane roads, one of the key parameters in evaluating hydroplaning risk. In the first part of the study, water films forming on wide roadways, i.e. multilane highways with 2, 4, and 6 lanes per side, were analyzed, with varying cross slope, longitudinal slope, and rainfall rate. Roadways with and without curbs and drainage were included in the analysis. The analysis in the first part was limited to nearly smooth asphalt or concrete surfaces because the maximum roughness height that can be specified in a CFD model using wall functions to determine shear stress at the road surface boundary is one half the thickness of computational cell layer adjacent to the boundary and the need to resolve the vertical velocity distribution in a thin water film requires cell layers that are a fraction of a millimeter thick. In this part of the study, the analysis was extended to consider the influence of pavement roughness over a range extending up 3.6 mm on the development of water film thickness. Various methods of modeling roughness were investigated and are presented: the roughness height model, a meshed-out geometry of the macrotexture of the pavement, and a porous region model. The computational results were analyzed and compared with experimental measurements performed by Gallaway et al.

42 ENGINEERING↗
DOE OSTI2025

Outflow Boundary Conditions for Turbine-Integrated Rotating Detonation Combustors

This study examines outflow boundary conditions (BCs) in computational fluid dynamics (CFD) simulations of a transition duct with and without guide vanes that converts supersonic flow exiting a rotating detonation combustor (RDC) to subsonic flow to drive a turbine. Since the flow exiting the transition duct has swirling shock waves with significant spatial and temporal variations in pressure, temperature, and Mach number, imposing proper BCs poses a challenge. To ensure all swirling shock waves exit the transition duct without creating non-physical reflected waves at its outlet, this study examined three outflow BCs: (1) the average pressure imposed at the duct’s outlet, (2) a nonreflecting BC (NRBC) with a specified average pressure imposed at the duct’s outlet, (3) the average pressure imposed at the outlet of an extension duct made up of a buffer layer and a sponge layer. This study is based on the three-dimensional, unsteady density-weighted-ensemble-averaged continuity, Navier–Stokes, and energy equations for a thermally perfect gas closed by the realizable k–ε model and “enhanced” wall functions. The results obtained show that imposing an average pressure at the transition duct’s outlet produces spurious waves that degrade the physical meaningfulness of the solution. When the NRBC was applied, swirling shock waves exited the duct’s outlet without creating spurious waves. However, its usage requires the gas to be thermally, as well as calorically, perfect, which this study shows could be a concern. By imposing the average pressure at the outlet of an extension duct, the gas does not need to be calorically perfect. The results obtained show the effects of the sponge layer’s length and coarsening ratio on damping nonuniformities in non-physical reflected waves to ensure the flow exiting the transition duct’s outlet can do so as if there are no boundaries present and has the desired average pressure—even though the BC is applied at the extension duct’s outlet.

gas turbines↗
NASA NTRS1982

A viscous-inviscid interaction algorithm for three-dimensional turbulent subsonic aerodynamic juncture region flow

An order-of magnitude analysis of the subsonic, three-dimensional, steady time-averaged Navier-Stokes equations for semibounded aerodynamic juncture geometries yields the parabolic Navier-Stokes simplification. The numerical solution of the resultant pressure Poisson equation is cast into complementary and particular parts, yielding an iterative interaction algorithm with an exterior three-dimensional potential flow solution. A parabolic transverse momentum equation set is constructed, wherein robust enforcement first-order continuity effects is accomplished using a penalty differential constraint concept within a finite element solution algorithm. A Reynolds stress constitutive equation, with low turbulence Reynolds number wall functions, is employed for closure, using parabolic forms of the two-equation turbulent kinetic energy-dissipation equation system. The algorithm is employed to predict the three-dimensional turbulent flowfield in the juncture region formed by two intersecting parabolic arcs at a free-stream Mach-number of 0.08, and Re/C = 600,000/m.

Baker, A. J.↗
NASA NTRS1982

The CMC:3DPNS computer program for prediction of three-dimensional, subsonic, turbulent aerodynamic juncture region flow. Volume 1: Theoretical

An order-of-magnitude analysis of the subsonic three dimensional steady time averaged Navier-Stokes equations, for semibounded aerodynamic juncture geometries, yields the parabolic Navier-Stokes simplification. The numerical solution of the resultant pressure Poisson equation is cast into complementary and particular parts, yielding an iterative interaction algorithm with an exterior three dimensional potential flow solution. A parabolic transverse momentum equation set is constructed, wherein robust enforcement of first order continuity effects is accomplished using a penalty differential constraint concept within a finite element solution algorithm. A Reynolds stress constitutive equation, with low turbulence Reynolds number wall functions, is employed for closure, using parabolic forms of the two-equation turbulent kinetic energy-dissipation equation system. Numerical results document accuracy, convergence, and utility of the developed finite element algorithm, and the CMC:3DPNS computer code applied to an idealized wing-body juncture region. Additional results document accuracy aspects of the algorithm turbulence closure model.

Baker, A. J.↗
NASA NTRS1983

Effects of expansion ratio on the calculation of parallel-walled backward-facing step flows - Comparison of four models of turbulence

This paper presents a numerical study of the effects of expansion ratio on two-dimensional separating and reattaching flows in plane backward-facing step geometries with parallel walls. Closure of the Reynolds equations was achieved by four different turbulence models: k-epsilon, 'modified' k-epsilon, algebraic stress, and 'modified' algebraic stress models. The k-epsilon model relates the Reynolds stresses to the mean rate of strain through the definition of an isotropic turbulent viscosity. The more advanced algebraic stress model calculates the stresses from implicit algebraic relationships containing the stresses themselves, the mean rate of strain, and the turbulent kinetic energy and its dissipation rate. 'Modified' versions of the models employ a new dissipation rate equation whose production term was made more sensitive to streamwise curvature effects. A new nonequilibrium wall function treatment proposed by Chieng and Launder (1980) was also incorporated into each model.

Sindir, M. M.↗
NASA NTRS1984

A computational study of complex three-dimensional compressible turbulent flow fields

Solutions of the time dependent, Reynolds-averaged, Navier-Stokes equations are presented and are compared with a family of experimental results for the three-dimensional interaction of a shock wave with a turbulent boundary layer. The solutions correctly predict the major features of the flow field independent of the shock strength and the extent of separation when using a two-equation turbulence model with wall functions. The experimentally observed boundary between cylindrical and conical flow regimes is also predicted. However, for the flow fields with large separated zones, the details of the measured pressure distributions are not accurately computed, which indicates a need for improved turbulence modeling and/or grid resolution.

Horstman, C. C.↗
NASA NTRS1984

SSME thrust chamber simulation using Navier-Stokes equations

The capability of the PHOENICS fluid dynamics code in predicting two-dimensional, compressible, and reacting flow in the combustion chamber and nozzle of the space shuttle main engine (SSME) was evaluated. A non-orthogonal body fitted coordinate system was used to represent the nozzle geometry. The Navier-Stokes equations were solved for the entire nozzle with a turbulence model. The wall boundary conditions were calculated based on the wall functions which account for pressure gradients. Results of the demonstration test case reveal all expected features of the transonic nozzle flows. Of particular interest are the locations of normal and barrel shocks, and regions of highest temperature gradients. Calculated performance (global) parameters such as thrust chamber flow rate, thrust, and specific impulse are also in good agreement with available data.

Przekwas, A. J.↗
NASA NTRS1986

A computational study of multiple jet and wall interaction

In the fuel preburner of the Space Shuttle Main Engine, face plate, injector, and baffle erosion have been observed. The observed patterns of erosion suggest that flame attachment to the walls is a contributing factor. To better understand the physical phenomena involved, a portion of the preburner was modeled computationally. The simulated 'preburner' had three two-dimensional jets entering a cavity adjacent to a baffle. The computational model employed the Patankar Spalding algorithm with upwind differencing. The turbulence model was a standard k-epsilon model with wall functions. The effects of incoming boundary conditions on turbulent kinetic energy and dissipation, k and epsilon, was studied. The results indicate a very strong sensitivity to these boundary conditions over certain ranges of values.

Campbell, W.↗
NASA NTRS1986

The flowfield structure of the 3-D shock wave - Boundary layer interaction generated by a 20 deg sharp fin at Mach 3

The three-dimensional oblique shock wave-turbulent boundary layer interaction generated by a sharp fin attached to a flat plate is investigated experimentally and theoretically for Mach 3 and Reynolds number 9 x 10 to the 5th using two different models. Both models employ the three-dimensional compressible Navier-Stokes equations in mass-averaged variables; one model utilizes the algebraic turbulent eddy viscosity model of Baldwin and Lomax (1978), while the other model employs the two-equation turbulence model of Jones and Launder (1972) coupled with the wall function model of Viegas and Rubesin (1985). The computed surface pressure, surface streamlines, pitot pressure, and yaw angle profiles are found to be in good agreement with experimental data. The three-dimensional velocity fields computed by both models are in close agreement, although the eddy viscosity profiles differ significantly within the three-dimensional interaction. This result indicates that the overall structure of this three-dimensional sharp fin interaction is insensitive to the turbulence model.

Knight, D. D.↗
NASA NTRS1986

Experiments on supersonic turbulent flow development in a square duct

The nature of supersonic, turbulent, adiabatic-wall flow in a square duct is investigated experimentally over a development length of x/D between 0 and 20 for a uniform flow, Mach 3.9 condition at the duct inlet. Initial discussion centers on the duct configuration itself, which was designed specifically to minimize wave effects and nozzle-induced distortion in the flow. Total pressure contours and local skin friction coefficient distributions are presented which show that the flow develops in a manner similar to that observed for the incompressible case. In particular, undulations exist in total pressure contours within the cross plane and in transverse skin friction coefficient distributions, which are indicative of the presence of a well-defined secondary flow superimposed upon the primary flow. The results are analyzed to show that local law-of-the-wall behavior extends well into the corner region, which implies that wall functions conventionally applied in two-equation type turbulence models, when suitably defined for compressible flow, may also be applied to supersonic streamwise corner flows.

Gessner, F. B.↗
NASA NTRS1987

Solid rocket motor aft field joint flow field analysis

An efficient Navier-Stokes analysis was successfully applied to simulate the complex flow field in the vicinity of a slot in a solid rocket motor with segment joints. The capability of the computer code to resolve the flow near solid surfaces without using a wall function assumption was demonstrated. In view of the complex nature of the flow field in the vicinity of the slot, this approach is considered essential. The results obtained from these calculations provide valuable design information, which would otherwise be extremely difficult to obtain. The results of the axisymmetric calculations indicate the presence of a region of reversed axial flow at the aft-edge of the slot and show the over-pressure in the slot to be only about 10 psi. The results of the asymmetric calculations indicate that a pressure asymmetry more than two diameters downstream of the slot has no noticeable effect on the flow field in the slot. They also indicate that the circumferential pressure differential caused in the slot due to failure of a 15 deg section of the castable inhibitor will be approximately 1 psi.

Sabnis, Jayant S.↗
NASA NTRS1987

Structure of supersonic turbulent flow past a sharp fin

The three-dimensional shock-wave/turbulent-boundary-layer interaction generated by a sharp fin is examined both experimentally and theoretically at Mach 3 for a fin angle 20 deg and Reynolds number 900,000. This study represents an extension of previous research for the sharp fin configuration to stronger interactions. The experimental data include surface pressure profiles, surface streamline patterns, and boundary-layer profiles of pitot pressure and yaw angle. Two separate theoretical approaches or 'models' were employed. Both models utilize the three-dimensional compressible Navier-Stokes equations in mass-averaged variables. The theoretical approach of Knight (1984 and 1985) employs the algebraic turbulent eddy-viscosity model of Baldwin and Lomax (1978), and the theoretical model of Horstman (1984) employs the two-equation turbulence model of Jones and Launder (1972) coupled with the wall function model of Viegas et al. (1985).

Knight, Doyle D.↗
NASA NTRS1988

Computation of two-phase shear-layer flow using an Eulerian-Lagrangian analysis

A combined Eulerian-Lagrangian analysis which combines a linearized block implicit Navier-Stokes analysis for the continuous phase with a Lagrangian analysis for the discrete phase has been utilized to simulate a two-phase shear-layer flow in a two-dimensional channel flow configuration. The motion of the particles is tracked in computational coordinate space, resulting in improved computational efficiency, and interphase coupling terms for the Eulerian analysis are computed from the instantaneous distribution of the particles rather than the trajectory information. The use of the implicit Navier-Stokes analysis for the continuous phase has made it possible to efficiently use a highly stretched mesh. A low Reynolds number form of the k-epsilon turbulence model was used with near-wall resolution, thus eliminating the need to use the 'wall-function' approach. The streamwise velocity profiles for the continuous and the discrete phases have been compared with experimental data for two test cases.

Sabnis, Jayant S.↗
NASA NTRS1988

3-D stator-rotor interaction of the SSME

The time-dependent stator-rotor interaction pattern of the Space Shuttle Main Engine's High-Pressure Fuel-Side Turbopump is presently studied by a multiple-zone solution algorithm, for both a two-dimensional case with 2:3 stator-rotor blade ratio and a three-dimensional case with a 1:1 ratio. A time-centered differencing scheme is used for temporal discretization, together with a sliding-grid multiple-zone solution method and an extended two-equation turbulence model with a wall function approach. The algorithm's accuracy is tested against experimental data.

Chen, Yen-Sen↗
NASA NTRS1989

Turbulence and modeling in transonic flow

A review is made of the performance of a variety of turbulence models in the evaluation of a particular well documented transonic flow. This is done to supplement a previous attempt to calibrate and verify transonic airfoil codes by including many more turbulence models than used in the earlier work and applying the calculations to an experiment that did not suffer from uncertainties in angle of attack and was free of wind tunnel interference. It is found from this work, as well as in the earlier study, that the Johnson-King turbulence model is superior for transonic flows over simple aerodynamic surfaces, including moderate separation. It is also shown that some field equation models with wall function boundary conditions can be competitive with it.

Rubesin, Morris W.↗
NASA NTRS1988

A study of turbulent flow computations in an angled duct with a step

This paper presents a method for effective computation of turbulent flows in an angled duct which has a step at the angled section by using a grid generation technique. An effect of flow curvature was analyzed on the separating and reattaching flows. The computations are based on the Reynolds averaged transport equations which are transformed for use on a generalized curvilinear coordinate system. A special care was taken on the treatment of pressure corrections. Wall function boundary conditions for k-epsilon model were developed.

Maruszewski, J. P.↗
NASA NTRS1990

An abbreviated Reynolds stress turbulence model for airfoil flows

An abbreviated Reynolds stress turbulence model is presented for solving turbulent flow over airfoils. The model consists of two partial differential equations, one for the Reynolds shear stress and the other for the turbulent kinetic energy. The normal stresses and the dissipation rate of turbulent kinetic energy are computed from algebraic relationships having the correct asymptotic near wall behavior. This allows the model to be integrated all the way to the wall without the use of wall functions. Results for a flat plate at zero angle of attack, a NACA 0012 airfoil and a RAE 2822 airfoil are presented.

Gaffney, R. L., Jr.↗
NASA NTRS1990

Experimental investigation of turbulent flow through a circular-to-rectangular transition duct

An incompressible, turbulent, swirl-free flow through a circular-to-rectangular transition duct was studied experimentally. The cross-sectional geometry all along the duct was defined using the equation of a superellipse. The three mean velocity components and the six Reynolds stress components were measured at two axial stations downstream from the transition. It is shown that a secondary flow vortex pair which develops along the duct sidewalls significantly distorts the mean and turbulence fields. At the duct exit, the flow is not in local equilibrium, but recovers to local equilibrium conditions in the rectangular extension duct. Analysis demonstrates that conventional wall functions are not applicable at all streamwise locations in the duct.

Davis, D. O.↗