Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary layer flow”

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 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Cross correlation and length scales in turbulent flows near surfaces

Two kinds of length scales are used in turbulent flows; 'functional length scales' such as mixing length, dissipation length L(sub epsilon), etc., and 'flow-field length scales' derived from cross correlations of velocity, pressure, etc. in the flow. Some connection between these scales are derived here. We first consider the cross correlation R(sub vv)(y,y(sub 1)) of the normal components u at two heights y, y(sub 1) above a rigid surface, normalized by the velocity y(sub 1) (greater than y). For shear-free boundary layers it is found theoretically, and in field and numerical experiments that R(sub vv) approximately equals y/y(sub 1). For shear layers it is also found that R(sub vv) approximately equals f(y/y(sub 1)) less than or equal to y,y(sub 1). This function f differs slightly between low Reynolds number numerical simulations and field experiments. The lateral structure defined by R(sub vv)(y,r(sub 3); y(sub 1),0) is also self similar and shows that the eddies centered at about y(sub 1) appear to have constant lateral width a(sub 3) above and below y(sub 1), where a(sub 3, sup +) approximately equals 7+1/(1.4dU(sup +)/dy(sup +)), when normalized on u(sub *) and v, where U is the mean velocity. Results for L(sub epsilon, sup -1) from direct numerical simulation are found to compare well with the formula L(sub epsilon, sup -1) = A(sub B)/y + A(sub S)dU/dy/v, for unidirectional and reversing turbulent boundary layers and channel flow, except near where dU/dy approximately equals 0. The conclusion is that the large-scale eddy structure and length scales in these flows are determined by a combination of shear and blocking, and that the vertical component of turbulence has a self-similar structure in both kinds of boundary layer.

Hunt, J. C. R.↗
DOE OSTI2023

Demonstration of RANS models with wall functions in the spectral element code Nek5000

The spectral element based computational fluid dynamics (CFD) code Nek5000 has been traditionally used for high-fidelity applications, such as direct numerical simulation (DNS) and large eddy simulation (LES). These techniques require very fine numerical resolution to accurately capture turbulent fluctuations which can be prohibitively expensive for users without access to leadership class computing facilities. For broader application and adoption, significant effort has been invested to develop Reynolds-averaged Navier–Stokes (RANS) capabilities in Nek5000. Here, this work presents details of the implementation and demonstration of the standard wall functions for the κ–τ model in Nek5000. Results using the wall-modeled approach are compared to a wall-resolved approach for cases with negligible pressure gradient, viz., channel flow, pipe flow and flow in a reactor subchannel. Results show reasonably good agreement between the two approaches for friction factor and Nusselt number. Some expected differences are identified near the wall. These cases demonstrate the potential for significant computational savings by using much coarser meshes for the wall-modeled approach, with only minor differences between the predicted result. Additionally, several Reynolds numbers up to 1,000,000 are demonstrated for pipe flow and predicted friction factors and Nusselt numbers compared well to available correlations, with the worst below 10%. As the Reynolds number is increased, better agreement is observed between the correlations and the wall-modeled approach. In addition, flow in a molten salt fast reactor (MSFR) core is considered which features an adverse pressure gradient and flow separation. It showcases the inability of standard wall functions to accurately predict flows with adverse pressure gradients. The results, however, match reasonably well in trend in regions of the flow where the boundary layer is attached. Ongoing research is dedicated to include a pressure gradient correction to wall functions to improve the accuracy of flows with separation or reattachment and adverse or favorable pressure gradients.

42 ENGINEERING↗
NASA NTRS1982

Three-dimensional flow field in the tip region of a compressor rotor passage. II - Turbulence properties

The turbulence properties in the annulus wall region of an axial flow compressor rotor was measured using a triaxial, hot-wire probe rotating with the rotor. The flow was surveyed across the entire passage at five axial locations (leading edge, 1/4 chord, 1/2 chord, 3/4 chord, and the trailing edge location) and at six radial locations in a low-speed compressor rotor. The data derived include all three components of turbulence intensity and three components of turbulence stress. A comprehensive interpretation of the data with emphasis on features related to rotation, leakage flow, annulus wall boundary layer, and blade boundary layer interactions is included. All the components of turbulent intensities and stresses are found to be high in the leakage-flow mixing region. The radial component of intensities and stresses is found to be much higher than the corresponding streamwise components. The turbulent spectra clearly reveal the decay process of the inlet-guide-vane wake within the rotor passage.

Lakshminarayana, B.↗
NASA NTRS2017

Computational Design and Analysis of a Transonic Natural Laminar Flow Wing for a Wind Tunnel Model

A natural laminar flow (NLF) wind tunnel model has been designed and analyzed for a wind tunnel test in the National Transonic Facility (NTF) at the NASA Langley Research Center. The NLF design method is built into the CDISC design module and uses a Navier-Stokes flow solver, a boundary layer profile solver, and stability analysis and transition prediction software. The NLF design method alters the pressure distribution to support laminar flow on the upper surface of wings with high sweep and flight Reynolds numbers. The method addresses transition due to attachment line contamination/transition, Gortler vortices, and crossflow and Tollmien-Schlichting modal instabilities. The design method is applied to the wing of the Common Research Model (CRM) at transonic flight conditions. Computational analysis predicts significant extents of laminar flow on the wing upper surface, which results in drag savings. A 5.2 percent scale semispan model of the CRM NLF wing will be built and tested in the NTF. This test will aim to validate the NLF design method, as well as characterize the laminar flow testing capabilities in the wind tunnel facility.

Lynde, Michelle N.↗
NASA NTRSDate not supplied

Adaptive mesh refinement and turbulence modeling

Adaptive mesh refinement is shown to be essential in turbulence modeling. Adaptive Mesh Refinement (AMR) is a more computationally efficient means to obtain a discrete approximation to a continuous boundary value problem of a specified accuracy than classic isotropic grid refinement. Previous application of this methodology in the assessment of turbulence models suggested that the field variable solution on the interior of the domain was more sensitive to grid refinement than would be suggested by the sensitivity to grid refinement of surface quantities. Revisiting challenging high speed flow field revealed that issues believed to be fundamental turbulence model issues appear to be related to an under-resolved boundary layer edge – a flow detail previously considered unimportant.

Hypersonics↗
NASA NTRS2013

Numerical Investigation of a Fuselage Boundary Layer Ingestion Propulsion Concept

In the present study, a numerical assessment of the performance of fuselage boundary layer ingestion (BLI) propulsion techniques was conducted. This study is an initial investigation into coupling the aerodynamics of the fuselage with a BLI propulsion system to determine if there is sufficient potential to warrant further investigation of this concept. Numerical simulations of flow around baseline, Boundary Layer Controlled (BLC), and propelled boundary layer controlled airships were performed. Computed results showed good agreement with wind tunnel data and previous numerical studies. Numerical simulations and sensitivity analysis were then conducted on four BLI configurations. The two design variables selected for the parametric study of the new configurations were the inlet area and the inlet to exit area ratio. Current results show that BLI propulsors may offer power savings of up to 85% over the baseline configuration. These interim results include the simplifying assumption that inlet ram drag is negligible and therefore likely overstate the reduction in power. It has been found that inlet ram drag is not negligible and should be included in future analysis.

Elmiligui, Alaa A.↗
NASA NTRS1974

Computation of steady axisymmetric flow using a one-dimensional time-dependent method

An iterative numerical method for computing steady, three dimensional, viscous, compressible flow fields, about aerodynamic bodies was studied. In order to develop the iterative method economically, the primary emphasis was directed towards supersonic, axisymmetric flow. However, the technique readily extends to three spatial dimensions. The viscous flow field about a cone-cylinder-flare body was calculated and compared to existing experimental data. Numerical predictions of the cone boundary layer and the flow field shock wave structure agreed with corresponding measurements. A separation was calculated at the cylinder-flare junction in six iterations; however, the size of the vortex did not correspond to the measured size. It was estimated that fifty iterations would be required to properly define the vortex. It was concluded that the iteration technique is of limited value for plane two dimensional and axisymmetrix flows, but of great value for three-dimensional flows.

Walitt, L.↗
NASA NTRS1993

On the instability of flow in a streamwise corner

The linear stability of an incompressible laminar flow in the blending boundary layer between the boundary layer in a 90 deg streamwise corner and a Blasius boundary layer well away from the corner is examined using a locally parallel flow approximation. It is shown that the influence of the outer boundary conditions associated with oblique modes of disturbances which are anti-symmetric about the bisector plane have a profound effect on the stability of the flow. As a result, in good agreement with observation, the critical streamwise Reynolds number, associated with a spanwise location is significantly reduced as the corner is approached, being R(sub Cr) = 60 approximately for spanwise distance of z* = 6x*R(sup -1) from the corner compared with R(sub Cr) = 322 approximately for z* = 20x*R(sup -1), where x* measures downstream distance from the leading edges. At R = 600, the growth rate of the most amplified mode of disturbance at the former location is over six times greater than that at the latter; the corresponding wave angle at the two locations is respectively 44 deg and 5 deg, approximately.

Dhanak, Manhar R.↗
NASA NTRS1989

Coherent structures in a boundary layer and shear layer of a turbulent backward-facing step flow

A wind tunnel experiment has been carried out at the NASA Ames Research Center to analyze the evolution of coherent structures from a boundary layer to a shear layer in a turbulent, backward-facing, step flow. A miniature X-wire/cold-wire probe has been used in conjunction with two arrays of cold wires, one aligned in the plane of main shear and the other in the spanwise direction of the flow, to detect and characterize delta-scale organized structures in the outer regions of the flow and to provide detailed information concerning these structures. Kinematic features of the events associated with the large scale structures were analyzed and topological pictures of the evolving flow, as well as the contributions to the Reynolds shear stress components are presented.

Jovic, Srba↗
DOE OSTI2022

Adaptive Determination of the Optimal Exchange Location in Wall-Modeled Large-Eddy Simulation

Wall-modeled large-eddy simulation introduces a modeling interface (or exchange location) separating the wall-modeled layer from the rest of the domain. The current state-of-the-art is to rely on user expertise when choosing where to place this modeling interface, whether this choice is tied to the grid or not. This paper presents a postprocessing algorithm that determines the exchange location systematically. The algorithm is based on a model for the error in the predicted wall shear stress and a model for the computational cost, and then finds the exchange location that minimizes a combination of the two. Here, the algorithm is tested both a priori and a posteriori using an equilibrium wall model for the flow over a wall-mounted hump, a boundary layer in an adverse pressure gradient, and a shock/boundary-layer interaction. The algorithm produces exchange locations that mostly agree with what an experienced user would suggest, with thinner wall-modeled layers in nonequilibrium flow regions and thicker wall-modeled layers where the boundary layer is closer to equilibrium. This suggests that the algorithm should be useful in simulations of realistic and highly complex geometries.

42 ENGINEERING↗
NASA NTRS1975

Solutions for incompressible separated boundary layers including viscous-inviscid interaction

Numerical solutions are presented for the laminar and turbulent boundary-layer equations for incompressible flows with separation and reattachment. The separation angularity is avoided by using an inverse technique in which the displacement thickness is prescribed and the pressure is deduced from the resulting solution. The turbulent results appear qualitatively correct despite the use of a two-layer eddy-viscosity model which is generally assumed appropriate only for mild-pressure-gradient flows. A new viscous-inviscid interaction technique is presented in which the inviscid flow is solved inversely by prescribing the pressure from the boundary-layer solution and deducing the new displacement thickness from the solution of a Cauchy integral. Calculations are presented using this interaction procedure for a laminar flow in which separation and reattachment occur on a solid surface.

Carter, J. E.↗
NASA NTRS2020

Simulation and Modeling of Cold-Wall Hypersonic Turbulent Boundary Layers on Flat Plate

Direct numerical simulations (DNS) of flat-plate, zero-pressure-gradient turbulent boundary layers are presented for nominal freestream Mach numbers of 11 and 14 and a highly cooled wall (wall-to-recovery temperature of approximately 0.2). The flow conditions of the DNS are representative of the experimental data for a Mach 11.1 turbulent boundary layer on a flat plate that was tested at Calspan–University of Buffalo Research Center (CUBRC) and the operational conditions of the AEDC Hypervelocity Tunnel No. 9 at Mach 14. The wall shear stress and turbulent heat flux predicted by DNS show good comparisons with those measured at CUBRC and those modeled with a Baldwin-Lomax model. The generated DNS database extends the range of Reynolds numbers from those typical of previous DNS studies for hypersonic boundary-layer flows and confirms the validity of Morkovin’s hypothesis within the relatively broader range of Reynolds numbers considered in this study. A priori assessment of turbulent heat flux models with DNS data shows that the algebraic energy flux (AEF) formulation gives better predictions of the streamwise turbulent heat flux than the commonly used constant turbulent Prandtl number model while both models give similarly good predictions of the wall-normal turbulent heat flux. A posteriori assessment shows that the AEF model improves the prediction of the temperature profiles while keeping a good prediction of the velocity field.

Junji Huang↗
NASA NTRS1983

Progress toward the analysis of supersonic inlet flows

A zonal flow analysis procedure was developed to predict the flow through the supersonic diffuser of an axisymmetric mixed compression inlet at angle-of-attack. In this analysis, the inlet flow is divided into three types of regions, each with different dominant flow phenomena. These are the inviscid supersonic core, boundary layer, and shock/boundary layer interaction flows. An appropriate analysis was selected or developed for the three-dimensional flow in each type of region. Procedures were developed to interface these analyses for the overall inlet flow analysis. This analysis was applied to an inlet operating at M = 2.58 at several angle-of-attack conditions. Comparisons are presented between computed and measured flow properties for the inlet and for the component analysis flows. Extensions of the present procedure to include the terminal shock and subsonic diffuser flows are recommended. Desirable experiments for evaluation of the inlet analysis procedure or the component analyses and to support improved modeling or extension of the inlet analysis are defined and recommended.

Paynter, G. C.↗
NASA NTRS1959

Use of a Stanton Tube for Skin-Friction Measurements

A small total-pressure tube resting against a flat-plate surface was used as a Stanton tube and calibrated as a skin-friction meter at various subsonic and supersonic speeds. Laminar flow was maintained for the supersonic runs at a Mach number M(sub infinity) of 2. At speeds between M(sub infinity) = 1.33 and M(sub infinity) = 1.87, the calibrations were carried-out in a turbulent boundary layer. The subsonic flows were found to be in transition. The skin-friction readings of a floating-element type of balance served as the reference values against which the Stanton tube was calibrated. A theoretical model was developed which, for moderate values of the shear parameter tau, accurately predicts the performance of the Stanton tube in subsonic and supersonic flows. A "shear correction factor" was found to explain the deviations from the basic model when T became too large. Compressibility effects were important only in the case of turbulent supersonic flows, and they did not alter the form of the calibration curve. The test Reynolds numbers, based on the distance from the leading edge and free-stream conditions, ranged from 70,000 to 875,000. The turbulent-boundary-layer Reynolds numbers, based on momentum thickness, varied between 650 and 2,300. Both laminar and turbulent velocity profiles were taken and the effect of pressure gradient on the calibration was investigated.

Abarbanel, S. S.↗
NASA NTRS1994

A spectrally accurate boundary-layer code for infinite swept wings

This report documents the development, validation, and application of a spectrally accurate boundary-layer code, WINGBL2, which has been designed specifically for use in stability analyses of swept-wing configurations. Currently, we consider only the quasi-three-dimensional case of an infinitely long wing of constant cross section. The effects of streamwise curvature, streamwise pressure gradient, and wall suction and/or blowing are taken into account in the governing equations and boundary conditions. The boundary-layer equations are formulated both for the attachment-line flow and for the evolving boundary layer. The boundary-layer equations are solved by marching in the direction perpendicular to the leading edge, for which high-order (up to fifth) backward differencing techniques are used. In the wall-normal direction, a spectral collocation method, based upon Chebyshev polynomial approximations, is exploited. The accuracy, efficiency, and user-friendliness of WINGBL2 make it well suited for applications to linear stability theory, parabolized stability equation methodology, direct numerical simulation, and large-eddy simulation. The method is validated against existing schemes for three test cases, including incompressible swept Hiemenz flow and Mach 2.4 flow over an airfoil swept at 70 deg to the free stream.

Pruett, C. David↗
NASA NTRS1995

Turbulance boundary conditions for shear flow analysis, using the DTNS flow solver

The effects of different turbulence boundary conditions were examined for two classical flows: a turbulent plane free shear layer and a flat plate turbulent boundary layer with zero pressure gradient. The flow solver used was DTNS, an incompressible Reynolds averaged Navier-Stokes solver with k-epsilon turbulence modeling, developed at the U.S. Navy David Taylor Research Center. Six different combinations of turbulence boundary conditions at the inflow boundary were investigated: In case 1, 'exact' k and epsilon profiles were used; in case 2, the 'exact' k profile was used, and epsilon was extrapolated upstream; in case 3, both k and epsilon were extrapolated; in case 4, the turbulence intensity (I) was 1 percent, and the turbulent viscosity (mu(sub t)) was equal to the laminar viscosity; in case 5, the 'exact' k profile was used and mu(sub t) was equal to the laminar viscosity; in case 6, the I was 1 percent, and epsilon was extrapolated. Comparisons were made with experimental data, direct numerical simulation results, or theoretical predictions as applicable. Results obtained with DTNS showed that turbulence boundary conditions can have significant impacts on the solutions, especially for the free shear layer.

Mizukami, M.↗
NASA NTRS1991

Effect of non-equilibrium flow chemistry and surface catalysis on surface heating to AFE

The effect of nonequilibrium flow chemistry on the surface temperature distribution over the forebody heat shield on the Aeroassisted Flight Experiment (AFE) vehicle was investigated using a reacting boundary-layer code. Computations were performed by using boundary-layer-edge properties determined from global iterations between the boundary-layer code and flow field solutions from a viscous shock layer (VSL) and a full Navier-Stokes solution. Surface temperature distribution over the AFE heat shield was calculated for two flight conditions during a nominal AFE trajectory. This study indicates that the surface temperature distribution is sensitive to the nonequilibrium chemistry in the shock layer. Heating distributions over the AFE forebody calculated using nonequilibrium edge properties were similar to values calculated using the VSL program.

Stewart, David A.↗