SEARCH · Engineering Papers
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.
Various topics telative to laminar flow aircraft certification are discussed. Boundary layer stability, flaps for laminar flow airfoils, computational wing design studies, manufacturing requirements, windtunnel tests, and flow visualization are among the topics covered.
Detailed measurements of the mean flow and turbulence were made with the use of a micrometeorological facility consisting of an instrumented 76-m tall tower located within a 100-m distance from the Atlantic Ocean at Wallops Island, Virginia. Under moderately strong wind conditions, the popular neutral boundary layer flow model fails to provide an adequate description of the actual flow. In addition to detailed flow information for all wind directions, averages of the important flow parameters used for design such as vertical distribution of mean velocity, turbulence intensities and turbulence integral scales were presented for wind direction sectors with near uniform upstream terrain. Power spectra of the three velocity components for the prevailing northwesterly and southerly winds are discussed.
Papers are presented on the heat transfer and fluid mechanics in energy systems including diesel combustion analysis using rapid sampling and the numerical simulation of thermal response of architectural structures. Experiments associated with hypersonic boundary layers, the prediction of nozzle-wall pressure distribution, and the effects of nozzle geometry on an axisymmetric thrust augmentor are outlined. Predictions are made of film thickness on a roll-metering device and the thermal characteristics of multiple-surface submerged condenser modules are described. Numerical studies are made of the interaction between natural and forced convection in a complex geometry and of pressure drop and holdup in two-phase stratified flow. A numerical solution is given for viscous flows at high Reynolds numbers and for compressible boundary layer flows near separation or reattachment. The structure and transport properties in turbulent shear flows are reviewed along with the spectra of Reynolds stresses in a developing pipe flow.
The onset and development of both dynamically and convectively forced boundary layer rolls are studied with linear and nonlinear analyses of a truncated spectral model of shallow Boussinesq flow. Emphasis is given here on the energetics of the dominant roll modes, on the magnitudes of the roll-induced modifications of the initial basic state wind and temperature profiles, and on the sensitivity of the linear stability results to the use of modified profiles as basic states. It is demonstrated that the roll circulations can produce substantial changes to the cross-roll component of the initial wind profile and that significant changes in orientation angle estimates can result from use of a roll-modified profile in the stability analysis. These results demonstrate that roll contributions must be removed from observed background wind profiles before using them to investigate the mechanisms underlying actual secondary flows in the boundary layer. The model is developed quite generally to accept arbitrary basic state wind profiles as dynamic forcing. An Ekman profile is chosen here merely to provide a means for easy comparison with other theoretical boundary layer studies; the ultimate application of the model is to study observed boundary layer profiles. Results of the analytic stability analysis are validated by comparing them with results from a larger linear model. For an appropriate Ekman depth, a complete set of transition curves is given in forcing parameter space for roll modes driven both thermally and dynamically. Preferred orientation angles, horizontal wavelengths and propagation frequencies, as well as energetics and wind profile modifications, are all shown to agree rather well with results from studies on Ekman layers as well as with studies on near-neutral and convective atmospheric boundary layers.
This paper considers viscous flows with unseparated turbulent boundary layers over two-dimensional airfoils at transonic speeds. Conventional theoretical methods are based on boundary layer formulations which do not account for the effect of the curved wake and static pressure variations across the boundary layer in the trailing edge region. In this investigation an extended viscous theory is developed that accounts for both effects. The theory is based on a rational analysis of the strong turbulent interaction at airfoil trailing edges. The method of matched asymptotic expansions is employed to develop formal series solutions of the full Reynolds equations in the limit of Reynolds numbers tending to infinity. Procedures are developed for combining the local trailing edge solution with numerical methods for solving the full potential flow and boundary layer equations. Theoretical results indicate that conventional boundary layer methods account for only about 50% of the viscous effect on lift, the remaining contribution arising from wake curvature and normal pressure gradient effects.
The stability of viscous, incompressible flow along a streamwise corner, often called the corner boundary layer problem is examined. The semi-infinite boundary value problem satisfied by small amplitude disturbances in the "bending boundary layer' region is obtained. The mean secondary flow induced by the corner exhibits a flow reversal in this region. Uniformly valid "first approximations' to solutions of the governing differential equations are derived. Uniformity at infinity is achieved by a suitable choice of the large parameter and use of an approximate Langer variable. Approximations to solutions of balanced type have a phase shift across the critical layer which is associated with instabilities in the case of two dimensional boundary layer profiles.
The stability of viscous, incompressible flow along a streamwise corner, often called the corner boundary layer problem is examined. The semi-infinite boundary value problem satisfied by small amplitude disturbances in the 'bending boundary layer' region is obtained. The mean secondary flow induced by the corner exhibits a flow reversal in this region. Uniformly valid 'first approximations' to solutions of the governing differential equations are derived. Uniformity at infinity is achieved by a suitable choice of the large parameter and use of an approximate Langer variable. Approximations to solutions of balanced type have a phase shift across the critical layer which is associated with instabilities in the case of two dimensional boundary layer profiles. Previously announced in STAR as N84-17532
Radiative heating computations are performed for high speed lunar return experiments conducted in the Electric Arc Shock Tube (EAST) facility at NASA Ames Research Center. The nonequilibrium radiative transport equations are solved via NASA's in-house radiation code NEQAIR using flow field input from US3D flow solver. The post-shock flow properties for the 10 km/s Earth entry conditions are computed using the stagnation line of a blunt-body and a full facility CFD (Computational Fluid Dynamics) simulation of the EAST shock tube. The shocked gas in the blunt-body flow achieves a thermochemical equilibrium away from the shock front whereas EAST flow exhibits a nonequilibrium behavior due to strong viscous dissipation of the shock by boundary layer. The full-tube flow calculations capture the influence of the boundary layer on the shocked gas state and provide a realistic fluid dynamic input for the radiative predictions. The integrated radiance behind the shock is calculated in NEQAIR for wavelength regimes from Vacuum-UltraViolet (VUV) to InfraRed (IR), which are pertinent to the emission characteristics of high enthalpy shock waves in air. These radiance profiles are validated against corresponding EAST shots. The full-tube simulations successfully predict a sharp radiance peak at the shock front which gets smeared in the test data due to the spatial resolution in the measurements. The full facility based radiance behind the shock shows a slightly better match with the test data in the VUV and Red spectral regions, as compared to that from a blunt-body based predictions. The UV radiance is very similar for both geometries and under-predicts the test behavior. The IR test data matches better with the blunt-body based predictions where the full-tube simulations show a significant over-prediction.
The current investigation has the objective to provide data which will make it possible to obtain a better estimate regarding the roughness drag for surface waviness. The data employed for this investigation were acquired in connection with a wavy wall study which was conducted as part of an overall program to reduce the skin friction of turbulent boundary layers in external flows. The results of the present investigation show that the low-speed roughness drag of small-amplitude sinusoidal wave trains having wavelengths of the order of the boundary-layer thickness is not only a function of h/lambda (h = total wave height, lambda = wavelength), but, in addition, is also a function of the roughness Reynolds number.
A multi-point focused laser differential interferometer (FLDI) has been developed to measure density fluctuations at 16 points along a line. A pair of cylindrical lenses on the transmitter side of a conventional single-point FLDI instrument form two closely spaced (≤200 μm), orthogonally polarized, parallel laser lines at the instrument’s focus. On the receiver side of the instrument, the interference of the beams on a 16-element photodiode array results in a single line of measurements. The further addition of a Nomarski prism creates two separate measurement lines, and the addition of a second photodiode array to the instrument enables simultaneous measurements of density fluctuations along the two lines separated by several millimeters. These two lines of measurement can be conveniently oriented at any azimuthal angle relative to the instrument’s optical axis on the measurement plane, coinciding with the instrument’s focus. Two experiments were performed to demonstrate the capabilities of the instrument. In the first experiment, a laser-induced breakdown spark generated a traveling spherical shock wave, and measurements of the resulting density disturbance and wave velocity were obtained. These results were compared to high-speed schlieren images of the shock wave acquired at 400 kHz. In the second experiment, the multi-point FLDI instrument was used to measure density disturbances in the boundary layer of a flat plate in a Mach 6 freestream flow. The measurements were made along two lines, both approximately 6 mm in length, extending from the surface of the plate through the boundary layer. High-speed schlieren images were acquired at 100 kHz during separate wind tunnel runs at matching unit Reynolds numbers to visualize the unsteady boundary layer flow and compare to the FLDI measurements.
A time marching 3-D full Navier-Stokes code, called PARC3D, is validated for an experimental Mach 5 inlet configuration using the data obtained in the 10 x 10 ft supersonic wind tunnel at the NASA Lewis Research Center. For the first time, a solution is obtained for this configuration with the actual geometry, the tunnel conditions, and all the bleed zones modeled in the computation. Pitot pressure profiles and static pressures at various locations in the inlet are compared with the corresponding experimental data. The effect of bleed zones, located in different places on the inlet walls, in eliminating the low energy vortical flow generated from the 3-D shock-boundary layer interaction is simulated very well even though some approximations are used in applying the bleed boundary conditions and in the turbulence model. A further detailed study of the effect of individual bleed ports is needed to understand fully the actual mechanism of efficiently eliminating the vortical flow from the inlet. A better turbulence model would help to improve the accuracy even further in predicting the corner flow boundary layer profiles.
Density disturbance ahead of sphere in rarefied nitrogen or argon supersonic flow, and stability of viscous three-dimensional disturbances in laminar compressible boundary layer
An optimization method has been developed to design axisymmetric body shapes such as fuselages, nacelles, and external fuel tanks with increased transition Reynolds numbers in subsonic compressible flow. The new design method involves a constraint minimization procedure coupled with analysis of the inviscid and viscous flow regions, and linear stability analysis of the compressible boundary layer. Boundary-layer transition is predicted by a 'hybrid' transition criterion based on Granville's transition criterion and a criterion using linear stability theory coupled with the e(n)-method. A tiptank of a business-jet is used as an example to illustrate that the method can be utilized to design an axisymmetric body shape with extensive natural laminar flow. Boundary layer transition is predicted to occur at a transition Reynolds number of 6.04 x 10 to the 6th on the original tiptank. On the designed body shape, a transition Reynolds number of 7.22 x 10 to the 6th is predicted using compressible linear stability theory coupled with e(n)-method.
We have further developed our capabilities to analyze transition in turbine boundary layers from first principles by integrating the nonlinear parabolized stability equations (PSE) with improved initial and boundary conditions. With modified iteration schemes, we are able to proceed deeper into the transition region where skin friction coefficient and heat transfer coefficient significantly increase. Initial and boundary conditions at elevated turbulence levels can be derived by receptivity analysis. Test runs for ERCOFTAC test case T3A at 2.4\% turbulence level provide results in good agreement with the experimental data. The sharper minimum of the skin coefficient also shown by DNS results is likely due to the missing intermittency. The method has been applied to various experimentally studied turbine blades (UTRC, VKI, Zierke, Langston, Hippensteele, and others). The PSE results, though physically reasonable, do not agree as well as expected with the experimental findings. We have, therefore, performed an extensive search for the reasons of the seemingly systematic deviations. A first source of uncertainty has been found in the often insufficient documentation of the experiments (e.g. on blockage by end-wall boundary layers). However, variation of the relevant parameters does not lead to more satisfactory agreement. A second reason has been found in the "standard procedure" which considers a 2D flow at midspan and uses a panel code and subsequent boundary-layer code to obtain the laminar basic flow for the transition analysis. Comparison with the pressure distribution obtained with a 3D design code (RVC3D) shows significant three-dimensionality of the flow (e.g. in the UTRC experiments). The spanwise variation has been neglected in our original PSE code. To overcome this problem, we have developed the PSE/3D for fully 3D boundary layers to account for streamwise and spanwise variations. Since the design code does not provide the boundary-layer flow with sufficient resolution, we have generated the Euler solution and employed a 3D boundary-layer code to obtain the viscous basic flow. Although only the linear stability level of PSE/3D has been implemented so far, the discrepancies with the experiments change but do not disappear. We still find deviations between the computed and experimental variations of C(sub f), and St along the blade for laminar flow. The main reason can be seen by comparing the solution of the boundary-layer code with the viscous results of the design code. The conventional boundary-layer solution exhibits an asymptotic behavior appropriate in external aerodynamics but does not match the steep gradients of the inviscid flow through the passage and consequently provides biased results for C(sub f), and St. An attempt is currently being made to correct this deficiency. Before attempting to perform the transition analysis for the viscous flow provided by the design code, we have analyzed the implementation and "best possible" results. Code and results exhibit flaws that may negatively affect the design and are intolerable for transition analysis. Therefore, we have decided to develop a new code to obtain a reliable basis for stability and transition studies. We expect to report improved results by the time of the meeting.
The purpose of the present note is to show that on a flat plate where both the wall temperature and mean wall pressure are constant, neither of the limitations of parallel flow or of unity for the turbulent Prandtl number are required in order for the Crocco solution to apply to the turbulent boundary-layer flow. It is shown herein that this result is subject to restrictions on the magnitude of pressure fluctuations. The same analysis is generalized to show that the compressible turbulent boundary layer on an isothermal swept flat plate is independent of the spanwise flow if the molecular Prandtl number is unity.
A recently developed K-tau model for near wall turbulent flows is applied to a variety of test cases. The turbulent flows considered include the incompressible flat plate boundary layer with adverse pressure gradients, incompressible flow past a backward facing step, and the supersonic flat plate boundary layer at zero pressure gradient. Calculations are performed for this two-equation model using an anisotropic as well as isotropic eddy-viscosity. The model predictions are shown to compare quite favorably with experimental data.
Turbulence modeling continues to be one of the most difficult problems in fluid mechanics. Existing prediction methods are well developed for certain classes of simple equilibrium flows, but are still not entirely satisfactory for a large category of complex non-equilibrium flows found in engineering practice. Direct and large-eddy simulation (LES) approaches have long been believed to have great potential for the accurate prediction of difficult turbulent flows, but the associated computational cost has been prohibitive for practical problems. This remains true for direct simulation but is no longer clear for large-eddy simulation. Advances in computer hardware, numerical methods, and subgrid-scale modeling have made it possible to conduct LES for flows or practical interest at Reynolds numbers in the range of laboratory experiments. The objective of this work is to apply ES and the dynamic subgrid-scale model to the flow of a boundary layer over a concave surface.
A newly developed spectral compressible linear stability code (SPECLS) (staggered pressure mesh) is presented for analysis of shear flow stability, and applied to high speed boundary layers and free shear flows. The formulation utilizes the first application of a staggered mesh for a compressible flow analysis by a spectral technique. An order of magnitude less number of points is needed for equivalent accuracy of growth rates compared to those calculated by a finite difference formulation. Supersonic disturbances which are found to have oscillatory structures were resolved by a spectral multi-domain discretization, which requires a factor of three fewer points than the single domain spectral stability code. It is indicated, as expected, that stability of mixing layers is enhanced by viscosity and increasing Mach number. The mean flow involves a jet being injected into a quiescent gas. Higher temperatures of the injected gas is also found to enhance stability characteristics of the free shear layer.