Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “law of the wall”

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 37 records · Page 2

Scaling and modeling of three-dimensional, end-wall, turbulent boundary layers

The method of matched asymptotic expansion was employed to identify the various subregions in three dimensional, turbomachinery end wall turbulent boundary layers, and to determine the proper scaling of these regions. The two parts of the b.l. investigated are the 3D pressure driven part over the endwall, and the 3D part located at the blade/end wall juncture. Models are proposed for the 3d law of the wall and law of the wake. These models and the data of van den Berg and Elsenaar and of Mueller are compared and show good agreement between models and experiments.

Goldberg, U. C.↗

A Generalized Wall Function

The asymptotic solutions, described by Tennekes and Lumley (1972), for surface flows in a channel, pipe or boundary layer at large Reynolds numbers are revisited. These solutions can be extended to more complex flows such as the flows with various pressure gradients, zero wall stress and rough surfaces, etc. In computational fluid dynamics (CFD), these solutions can be used as the boundary conditions to bridge the near-wall region of turbulent flows so that there is no need to have the fine grids near the wall unless the near-wall flow structures are required to resolve. These solutions are referred to as the wall functions. Furthermore, a generalized and unified law of the wall which is valid for whole surface layer (including viscous sublayer, buffer layer and inertial sublayer) is analytically constructed. The generalized law of the wall shows that the effect of both adverse and favorable pressure gradients on the surface flow is very significant. Such as unified wall function will be useful not only in deriving analytic expressions for surface flow properties but also bringing a great convenience for CFD methods to place accurate boundary conditions at any location away from the wall. The extended wall functions introduced in this paper can be used for complex flows with acceleration, deceleration, separation, recirculation and rough surfaces.

Shih, Tsan-Hsing↗

A three-dimensional turbulent separated flow and related mesurements

The applicability of and the limits on the applicability of 11 near wall similarity laws characterizing three-dimensional turbulent boundary layer flows were determined. A direct force sensing local wall shear stress meter was used in both pressure-driven and shear-driven three-dimensional turbulent boundary layers, together with extensive mean velocity field and wall pressure field data. This resulted in a relatively large number of graphical comparisons of the predictive ability of 10 of these 11 similarity models relative to measured data over a wide range of flow conditions. Documentation of a complex, separated three-dimensional turbulent flow as a standard test case for evaluating the predictive ability of numerical codes solving such flows is presented.

Pierce, F. J.↗

New scaling law for turbulent boundary layers with high surface mass transfer

A fundamental understanding of the phenomena occurring in the turbulent boundary layer in the presence of surface blowing is limited, and considerable disagreements persist even in describing primary quantities, such as the boundary layer profile. The theories based on the linear boundary layer equations show that the thickness of the sublayer increases in the presence of surface blowing; therefore, the viscous sublayer and law of the wall modify. In this study, direct numerical simulations (DNS) of turbulent boundary layers with uniform surface mass transfers are carried out in order to scale the velocity profile. Emphasis is placed on moderate to high mass transfer rates, which are relevant to the most common hybrid rockets configuration. DNS data are used to establish a functional law of the wall and a law of wake by means of the relation between the wall shear stress and surface mass transfer. Analysis of the mean kinetic energy budget shows that the magnitude of turbulent kinetic energy increases by surface mass transfer, and the production rate extends significantly in the inner layer as the injection rate increases. DNS data of various surface blowing are used to complete the closure of turbulence kinetic energy equation and develop an eddy viscosity model. The predicted turbulent kinetic energy and eddy viscosity agree with DNS data for moderate to high blowing rates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Heat Transfer Through Turbulent Friction Layers

The "general Prandtl number" Pr(exp 1) - A(sub q)/A Pr, aside from the Reynolds number determines the ratio of turbulent to molecular heat transfer, and the temperature distribution in turbulent friction layers. A(sub q) = exchange coefficient for heat; A = exchange coefficient for momentum transfer. A formula is derived from the equation defining the general Prandtl number which describes the temperature as a function of the velocity. For fully developed thermal boundary layers all questions relating to heat transfer to and from incompressible fluids can be treated in a simple manner if the ratio of the turbulent shear stress to the total stress T(sub t)/T in the layers near the wall is known, and if the A(sub q)/A can be regarded as independent of the distance from the wall. The velocity distribution across a flat smooth channel and deep into the laminar sublayer was measured for isothermal flow to establish the shear stress ratio T(sub t)/T and to extend the universal wall friction law. The values of T(sub t)/T which resulted from these measurements can be approximately represented by a linear function of the velocity in the laminar-turbulent transition zone. The effect of the temperature relationship of the material values on the flow near the wall is briefly analyzed. It was found that the velocity at the laminar boundary (in contrast to the thickness of the laminar layer) is approximately independent of the temperature distribution. The temperature gradient at the wall and the distribution of temperature and heat flow in the turbulent friction layers were calculated on the basis of the data under two equations. The derived formulas and the figures reveal the effects of the Prandtl number, the Reynolds number, the exchange quantities and the temperature relationship of the material values.

Reichardt, H.↗

Streamline curvature effects on turbulent boundary layers

A theoretical tool has been developed for predicting, in a nonempirical manner, effects of streamline curvature and coordinate-system rotation on turbulent boundary layers. The second-order closure scheme developed by Wilcox and Traci has been generalized for curved streamline flow and for flow in a rotating coordinate system. A physically based straightforward argument shows that curvature/rotation primarily affects the turbulent mixing energy; the argument yields suitable curvature/rotation terms which are added to the mixing-energy equation. Singular-perturbation solutions valid in the wall layer of a curved-wall boundary layer and a fully developed rotating channel flow demonstrate that, with the curvature/rotation terms, the model predicts the curved-wall and the rotating coordinate system laws of the wall. Results of numerical computations of curved-wall boundary layers and of rotating channel flow show that curvature/rotation effects can be computed accurately with second-order closure.

Wilcox, D. C.↗

Convex curvature effects on the heated turbulent boundary layer

A convexly curved and isothermally heated wall with a 45-cm radius of curvature is subjected to turbulent boundary layer flow measurements in order to determine wall heat transfer rates and mean velocity and temperature profiles. Significant curvature effects are noted, with Stanton number and skin friction coefficient reductions of 35-40 percent by comparison with flat plate values for the same momentum or enthalpy thickness Reynolds numbers. Profiles of mean velocity and temperature show a more rapid growth of the wake regions, and a shortening of the log-linear region, as a result of curvature. Turbulent Prandtl numbers deduced from the mean temperature profiles under the assumption of a wall thermal law were found to be increased by 40-50 percent by this strong convex curvature.

Simon, T. W.↗

Improved Nonequilibrium Algebraic Model Of Turbulence

Blend of previous models predicts pressure distributions more accurately. Improved algebraic model represents some of time-averaged effects of turbulence in transonic flow of air over airfoil. Based partly on comparisons among various eddy-viscosity formulations for turbulence and partly on premise that law of wall more universally valid in immediate region of surface in presence of adverse gradient of pressure than mixing-length theory and original Johnson and King model.

Johnson, D. A.↗

An enrichment wall modeling framework for spectral element methods

In the present work, a first-of-its-kind enrichment wall-model is developed within the spectral element method (SEM) framework for large-eddy simulations (LES) of wall-bounded turbulent flows. The method augments the polynomial solution in the wall-adjacent elements with an analytical law-of-the-wall enrichment function representing the mean velocity near the wall. In the solution representation, this enrichment function captures the large gradients in the boundary layer, which allows the polynomial modes to represent the turbulent fluctuations. The enriched solution is able to resolve the shear stress at the wall without any modification to the no-slip wall boundary conditions, which allows for greater accuracy in the near-wall region compared to traditional methods. The enrichment wall modeling approach is implemented in a high-order SEM computational fluid dynamics solver, Nek5000, and its performance is assessed in turbulent channel flow wall-modeled LES for a range of Reynolds numbers. It is demonstrated that the enrichment wall-model improves solution accuracy on under-resolved near-wall grids as compared to traditional shear stress wall-models.

42 ENGINEERING↗

A modified wall wake velocity profile for turbulent compressible boundary layers.

The law of the wake and the law of the wall in incompressible turbulent boundary layers formulated by Coles (1956) and their use by Mathews et al. (1970) in the development of a wall-wake representation of the velocity profile in a form applicable for isoenergetic compressible boundary layers are extended to a modified wall-wake velocity profile for turbulent compressible boundary layers. The modified wall-wake profile is shown to provide good representations of experimental velocity distributions.

Sun, C.-C.↗

A universal velocity transformation for boundary layers with pressure gradients

The logarithmic law of the wall does not capture the mean flow when a boundary layer is subjected to a strong pressure gradient. In such a boundary layer, the mean flow is affected by the spatio-temporal history of the imposed pressure gradient; and accounting for history effects remains a challenge. This work aims to develop a universal mean flow scaling for boundary layers subjected to arbitrary adverse or/and favourable pressure gradients. Here, we derive from the Navier–Stokes equation a velocity transformation that accounts for the history effects and maps the mean flow to the canonical law of the wall. The transformation is tested against channel flows with a suddenly imposed adverse or favourable pressure gradient, boundary layer flows subjected to an adverse pressure gradient, and Couette–Poiseuille flows with a streamwise pressure gradient. It is found that the transformed velocity profiles follow closely the equilibrium law of the wall.

97 MATHEMATICS AND COMPUTING↗

Nonadiabatic and three-dimensional effects in compressible turbulent boundary layers

A defect stream function formulation for nonadiabatic flow with small crossflow is developed. The first-integral property of this formulation provides for two removal of the streamline curvature term in the governing equation so that the form of the reduced equation for small crossflow is the same as that for two dimensional flow. The combined law of the wall and wake is used in place of the no-slip boundary condition. The tangential velocity equation for law-of-the-wall flow is shown to be the same for three-dimensions as for two when the Boussinesq approximation applies, and a closed form solution for the crossflow angle in the inner region is obtained. Analytic solutions for nonadiabatic, compressible, equilibrium flow with a Clauser outer-region eddy-viscosity model are obtained, and excellent agreement with experimental skin friction and velocity profile data for nonadiabatic, compressible flat-plate flow is achieved. An analytic solution for a linear inner-region eddy-viscosity model is also obtained; the wake function part of this solution is found to be inconsistent with the empirically established law of the wake.

Barnwell, Richard W.↗

Investigation of turbulent flow in highly curved ducts with application to turbomachinery components

Numerical predictions were performed using a semi-elliptic calculation procedure for the case of turbulent flow in passage through a 90 deg bend of square cross section. Two versions of the isotropic turbulent viscosity two equation k-epsilon model were used. The first, the wall function model (WFM), employs the logarithmic law-of-the-wall relation and the notion of equilibrium flow to set all the necessary boundary conditions at the first grid point adjacent to a solid wall. The second model, the Van Driest model (VDM), employs Prandtl's original mixing length formulation, in conjunction with Van Driest's semi-empirical relation for the mixing length, to calculate the turbulent viscosity in the near wall regions of the flow. In this case, boundary conditions for k and epsilon, required to calculate these quantities in the core of the flow, are obtained by matching the mixing length and Reynolds number model formulations in an overlapping region of the flow near the walls. In both cases the results obtained show an improvement over earlier calculations using an elliptic numerical procedure. This is attributed to the finer grids possible in the present work. Of the two models, the VDM formulation shows better overall conformity with the mean flow measurements. Neither model reproduces well the details of the stress distribution as a result of the implied isotropic turbulent viscosity.

Humphrey, J. A. C.↗

Geometrically Flexible and Efficient Flow Analysis of High Speed Vehicles Via Domain Decomposition, Part 1: Unstructured-Grid Solver for High Speed Flows

The ability to solve the equations governing the hypersonic turbulent flow of a real gas on unstructured grids using a spatially-elliptic, 2nd-order accurate, cell-centered, finite-volume method has been recently implemented in the VULCAN-CFD code. This paper describes the key numerical methods and techniques that were found to be required to robustly obtain accurate solutions to hypersonic flows on non-hex-dominant unstructured grids. The methods and techniques described include: an augmented stencil, weighted linear least squares, cell-average gradient method, a robust multidimensional cell-average gradient-limiter process that is consistent with the augmented stencil of the cell-average gradient method and a cell-face gradient method that contains a cell skewness sensitive damping term derived using hyperbolic diffusion based concepts. A data-parallel matrix-based symmetric Gauss-Seidel point-implicit scheme, used to solve the governing equations, is described and shown to be more robust and efficient than a matrix-free alternative. In addition, a y+ adaptive turbulent wall boundary condition methodology is presented. This boundary condition methodology is deigned to automatically switch between a solve-to-the-wall and a wall-matching-function boundary condition based on the local y+ of the 1st cell center off the wall. The aforementioned methods and techniques are then applied to a series of hypersonic and supersonic turbulent flat plate unit tests to examine the efficiency, robustness and convergence behavior of the implicit scheme and to determine the ability of the solve-to-the-wall and y+ adaptive turbulent wall boundary conditions to reproduce the turbulent law-of-the-wall. Finally, the thermally perfect, chemically frozen, Mach 7.8 turbulent flow of air through a scramjet flow-path is computed and compared with experimental data to demonstrate the robustness, accuracy and convergence behavior of the unstructured-grid solver for a realistic 3-D geometry on a non-hex-dominant grid.

White, Jeffery A.↗

The multi-zone calculation of turbomachinery flows. II - The multi-zone calculation of the turbulent, two-specie flow through the SSME HPFTP first and second stage cavities

A multi-zone Navier-Stokes methodology to calculate the two-specie flow through the first and second stage cooling cavities of the Space Shuttle Main Engine (SSME) high pressure fuel turbopump (HPFTP) is developed. A simplified two-component fluid formulation is used to model the interaction of coolant and hot gas. Johnston's secant approximation is used to define an appropriate near wall velocity for use in a three-dimensional law of the wall. The basic Navier-Stokes algorithm used is a finite-volume, predictor-corrector algorithm which uses a pressure correction technique. A multi-zone method is used to partition each cavity into easily handled subdomains. The results show that coolant flow is pumped up the turbine wheel for both cavities, creating a region of large temperature gradients on the turbine shank.

Williams, M.↗

In-Flight Capability for Evaluating Skin-Friction Gages and Other Near-Wall Flow Sensors

An 8-in.-square boundary-layer sensor panel has been developed for in-flight evaluation of skin-friction gages and other near-wall flow sensors on the NASA Dryden Flight Research Center F-15B/Flight Test Fixture (FTF). Instrumentation on the sensor panel includes a boundary-layer rake, temperature sensors, static pressure taps, and a Preston tube. Space is also available for skin-friction gages or other near-wall flow sensors. Pretest analysis of previous F-15B/FTF flight data has identified flight conditions suitable for evaluating skin-friction gages. At subsonic Mach numbers, the boundary layer over the sensor panel closely approximates the two-dimensional (2D), law-of-the-wall turbulent boundary layer, and skin-friction estimates from the Preston tube and the rake (using the Clauser plot method) can be used to evaluate skin-friction gages. At supersonic Mach numbers, the boundary layer over the sensor panel becomes complex, and other means of measuring skin friction are needed to evaluate the accuracy of new skin-friction gages. Results from the flight test of a new rubber-damped skin-friction gage confirm that at subsonic Mach numbers, nearly 2D, law-of-the-wall turbulent boundary layers exist over the sensor panel. Sensor panel data also show that this new skin-friction gage prototype does not work in flight.

Bui, Trong T.↗

Mean-flow measurements in a supersonic three-dimensional turbulent boundary layer

Measurements are presented of the mean flow in a supersonic turbulent boundary layer subjected to a constant weak transverse pressure gradient. The temperature and longitudinal velocity and the axial shear stress were only slightly affected by the three-dimensionality but a clearly defined crossflow component appeared, in the manner suggested by theory and confirmed in earlier experiments. The flow deflection angle and transverse-velocity component achieved a self-preserving form so long as the transverse-pressure gradient remained constant, and both achieved a maximum at the sublayer edge; the deflection angle seemed to decrease again between the latter and the surface. An empirical relation was found between the pressure-gradient strength and the maximum in the crossflow, and the dependence of the latter on distance from the surface was used to test analytical crossflow predictions. The data are in general agreement with Van Den Berg's law of the wall. The data also support the so-called parabolic law following a relaxation distance, especially if inviscid gradients are accounted for, and if the normal coordinate is contracted by the compressibility transformation.

Demetriades, A.↗

Development of a two-equation turbulence model for hypersonic flows: Evaluation of a low Reynolds number correction to the Kappa - epsilon two equation compressible turbulence model - Volume 1

The objective of the current research is the development of an improved k-epsilon two-equation compressible turbulence model for turbulent boundary layer flows experiencing strong viscous-inviscid interactions. The development of an improved model is important in the design of hypersonic vehicles such as the National Aerospace Plane (NASP) and the High Speed Civil Transport (HSCT). Improvements have been made to the low Reynolds number functions in the eddy viscosity and dissipation of solenoidal dissipation of the k-epsilon turbulence mode. These corrections offer easily applicable modifications that may be utilized for more complex geometries. The low Reynolds number corrections are functions of the turbulent Reynolds number and are therefore independent of the coordinate system. The proposed model offers advantages over some current models which are based upon the physical distance from the wall, that modify the constants of the standard model, or that make more corrections than are necessary to the governing equations. The code has been developed to solve the Favre averaged, boundary layer equations for mass, momentum, energy, turbulence kinetic energy, and dissipation of solenoidal dissipation using Keller's box scheme and the Newton spatial marching method. The code has been validated by removing the turbulent terms and comparing the solution with the Blasius solution, and by comparing the turbulent solution with an existing k-epsilon model code using wall function boundary conditions. Excellent agreement is seen between the computed solution and the Blasius solution, and between the two codes. The model has been tested for both subsonic and supersonic flat-plate turbulent boundary layer flow by comparing the computed skin friction with the Van Driest II theory and the experimental data of Weighardt; by comparing the transformed velocity profile with the data of Weighardt, and the Law of the Wall and the Law of the Wake; and by comparing the computed results of an adverse pressure gradient with the experimental data of Fernando and Smits. Good agreement is obtained with the experimental correlations for all flow conditions.

Knight, Doyle D.↗