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.
In program coupled, iterative implicit finite-difference procedure solves system of equations for laminar transitional or turbulent boundary-layer flows. Program written in FORTRAN IV.
Turbulence measurements were obtained in the nearly separated flow in a 13 deg total angle of divergence conical diffuser coupled to a constant area tailpipe. Air at 207 newtons per square centimeter and 308 K provided an inlet velocity of about 51 meters per second at an inlet unit Reynolds number of 63.7 million per meter. Very high longitudinal turbulence intensities accompanied the diffusion process with peak values approaching 40 percent when normalized by the local centerline velocity. Predictions of the pressure recovery coefficient using a mixing length concept were good in the early stages of diffusion. In the latter stages of diffusion satisfactory predictions of the pressure recovery were obtained with an empirical method.
The time dependent surface shear stress has been evaluated using surface heat transfer measurements. For fully developed turbulent pipe and open channel water flows, and incompressible and compressible turbulent boundary layer air flows the measurements indicate the absolute magnitude of the surface shear stress fluctuations will be greater than two times the mean values. The root-mean-square shear stress fluctuations were of the order of 0.2 to 0.4 times the mean surface shear values. Due to these large surface shear stress fluctuations and the nonlinear relation between heat transfer and shear stress, a special technique has been developed to evaluate the measurements. It was found that the non-linear averaging errors for a hot film-surface shear stress gauge in a fully developed pipe flow was of the order of 10 percent at low velocities. A hot wire-surface shear stress gauge was employed for measurements of turbulent boundary layers in air.
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.
A performance evaluation of an active sidewall boundary-layer removal system for the Langley 0.3-m Transonic Cryogenic Tunnel (TCT) was evaluated in 1988. This system uses a compressor and two throttling digital valves to control the boundary-layer mass flow removal from the tunnel. The compressor operates near the maximum pressure ratio for all conditions. The system uses a surge prevention and flow recirculation scheme. A microprocessor based controller is used to provide the necessary mass flow and compressor pressure ratio control. Initial tests on the system indicated problems in realizing smooth mass flow control while running the compressor at high speed and high pressure ratios. An alternate method has been conceived to realize boundary-layer mass flow control which avoids the recirculation of the compressor mass flow and operation near the compressor surge point. This scheme is based on varying the speed of the compressor for a sufficient pressure ratio to provide needed mass flow removal. The system has a mass flow removal capability of about 10 percent of test section flow at M = 0.3 and 4 percent at M = 0.8. The system performance has been evaluated in the form of the compressor map, and compressor tunnel interface characteristics covering most of the 0.3-m TCT operational envelope.
The Reynolds stress equation, pressure strain correlation, and dissipative terms and diffusion are discussed in relation to turbulence modelling using the Reynolds stress model. Algebraic modeling of Reynolds stresses and calculation of the boundary layer over an axial cylinder are examined with regards to the kinetic energy model for turbulence modelling. The numerical analysis of blade and hub wall boundary layers, and an experimental study of rotor blade boundary layer in an axial flow compressor rotor are discussed. The Patankar-Spalding numerical method for two dimensional boundary layers is included.
The numerical simulation of incompressible spatially-developing shear flows poses a special challenge to computational fluid dynamicists. The Navier-Stokes equations are elliptic and boundary equations need to be specified at the inflow and outflow boundaries in order to compute the fluid properties within the region of interest. It is, however, difficult to choose inflow and outflow conditions corresponding to a given experimental situation. Furthermore the effects that changes in the boundary conditions or in the size of the computational domain may induce on the global dynamics of the flow are presently unknown. These issues are examined in light of recent developments in hydrodynamic stability theory. The particular flow considered is the spatial mixing layer but it was expected that similar phenomena were bound to occur in other cases such as channel flow, the boundary layer, etc. A short summary of local/global and absolute/convective instability concepts is given. The results of numerical simulations are presented which strongly suggest that global resonances may be triggered in domains of finite streamwise extent although the evolution of the perturbation vorticity field is everywhere locally convective. A relationship between finite domains and pressure sources which might help in devising a scheme to eliminate these difficulties is discussed.
A simple and computationally efficient algorithm for solving the unsteady three-dimensional boundary-layer equations in the time-accurate or relaxation mode is presented. Results of the new algorithm are shown to be in quantitative agreement with detailed experimental data for flow over a swept infinite wing. The separated flow over a 6:1 ellipsoid at angle of attack, and the transonic flow over a finite-wing with shock-induced 'mushroom' separation are also computed and compared with available experimental data. It is concluded that complex, separated, three-dimensional viscous layers can be economically and routinely computed using a time-relaxation boundary-layer algorithm.
An experimental investigation at Mach 4 of shock-induced turbulent boundary layer separation at the walls of axially symmetric flow passages is discussed, with particular emphasis placed on determining the shock strengths required for incipient separation. The shock waves were produced by interchangeable sting-mounted cones placed on the axes of the flow passages and aligned with the freestream flow. The interactions under study simulate those encountered in axially symmetric engine inlets of supersonic aircraft. Knowledges of the shock strengths required for boundary layer separation in inlets is important since for shocks of somewhat greater strength rather drastic alterations in the inlet flow field may occur.
A 3-D laminar and turbulent boundary-layer method is developed for compressible flow over swept wings. The governing equations and curvature terms are derived in detail for a nonorthogonal, curvilinear coordinate system. Reynolds shear-stress terms are modeled by the Cebeci-Smith eddy-viscosity formulation. The governing equations are descretized using the second-order accurate, predictor-corrector finite-difference technique of Matsuno, which has the advantage that the crossflow difference formulas are formed independent of the sign of the crossflow velocity component. The method is coupled with a full potential wing/body inviscid code (FLO-30) and the inviscid-viscous interaction is performed by updating the original wing surface with the viscous displacement surface calculated by the boundary-layer code. The number of these global iterations ranged from five to twelve depending on Mach number, sweep angle, and angle of attack. Several test cases are computed by this method and the results are compared with another inviscid-viscous interaction method (TAWFIVE) and with experimental data.
A numerical study is made of the temporal eigenvalue spectrum of the Orr-Sommerfeld equation for the Blasius boundary layer. Unlike channel flows, there is no mathematical proof that this flow has an infinite spectrum of discrete eigenvalues. The Orr-Sommerfeld equation is integrated numerically, and the eigenvalues located by tracing out the contour lines in the complex wave velocity plane on which the real and imaginary parts of the secular determinant are zero. The spectrum of plane Poiseuille flow is used as a guide to study the spectrum of an artificial two-wall flow which consists of two Blasius boundary layers. As the upper boundary of this flow moves to infinity, it is found that the portion of the spectrum with an infinite number of eigenvalues moves towards phase velocity equal to unity and the spacing between eigenvalues goes to zero. The original few eigenvalues found are the only discrete eigenvalues that exist for Blasius flow.
The prediction of boundary layer deceleration and separation in regions of adverse pressure gradients are discussed as an unsolved problem of fundamental practical importance in aerodynamics. A three layered theoretical flow model of boundary layer separation including viscous-inviscid interaction was conceived for the case of subsonic two dimensional steady laminar flow.
The viscosity-dominated dynamic flow in a row of small transverse square cavities lying submerged in a turbulent boundary layer is considered. Experiments show that the cavity flow can be excited by free-stream disturbances in a narrow frequency band which is independent of the flow speed. The turbulent boundary layer in which the cavities are submerged remains transparent to the disturbances. The cavity flow resonates when the depths of the cavity and the Stokes layer match.
This paper will investigate the validation of the NASA developed, Reynolds-averaged Navier-Stokes (RANS) flow solver, OVERFLOW, for a boundary-layer-ingesting (BLI) offset (S-shaped) inlet in transonic flow with passive and active flow control devices as well as a baseline case. Numerical simulations are compared to wind tunnel results of a BLI inlet experiment conducted at the NASA Langley 0.3-Meter Transonic Cryogenic Tunnel. Comparisons of inlet flow distortion, pressure recovery, and inlet wall pressures are performed. The numerical simulations are compared to the BLI inlet data at a free-stream Mach number of 0.85 and a Reynolds number of approximately 2 million based on the fanface diameter. The numerical simulations with and without tunnel walls are performed, quantifying tunnel wall effects on the BLI inlet flow. A comparison is made between the numerical simulations and the BLI inlet experiment for the baseline and VG vane cases at various inlet mass flow rates. A comparison is also made to a BLI inlet jet configuration for varying actuator mass flow rates at a fixed inlet mass flow rate. Overall, the numerical simulations were able to predict the baseline circumferential flow distortion, DPCP avg, very well within the designed operating range of the BLI inlet. A comparison of the average total pressure recovery showed that the simulations were able to predict the trends but had a negative 0.01 offset when compared to the experimental levels. Numerical simulations of the baseline inlet flow also showed good agreement with the experimental inlet centerline surface pressures. The vane case showed that the CFD predicted the correct trends in the circumferential distortion levels for varying inlet mass flow but had a distortion level that was nearly twice as large as the experiment. Comparison to circumferential distortion measurements for a 15 deg clocked 40 probe rake indicated that the circumferential distortion levels are very sensitive to the symmetry of the flow and that a misalignment of the vanes in the experiment could have resulted in this difference. The numerical simulations of the BLI inlet with jets showed good agreement with the circumferential inlet distortion levels for a range of jet actuator mass flow ratios at a fixed inlet mass flow rate. The CFD simulations for the jet case also predicted an average total pressure recovery offset that was 0.01 lower than the experiment as was seen in the baseline. Comparisons of the flow features for the jet cases revealed that the CFD predicted a much larger vortex at the engine fan-face when compare to the experiment.
This paper will investigate the validation of a NASA developed, Reynolds-averaged Navier-Stokes (RANS) flow solver, OVERFLOW, for a boundary-layer-ingesting (BLI) offset (S-shaped) inlet in transonic flow with passive and active flow control devices as well as the baseline case. Numerical simulations are compared to wind tunnel results of a BLI inlet conducted at the NASA Langley 0.3-Meter Transonic Cryogenic Tunnel. Comparisons of inlet flow distortion, pressure recovery, and inlet wall pressures are performed. The numerical simulations are compared to the BLI inlet data at a freestream Mach number of 0.85 and a Reynolds number of approximately 2 million based on the length of the fan-face diameter. The numerical simulations with and without wind tunnel walls are performed, quantifying effects of the tunnel walls on the BLI inlet flow measurements. The wind tunnel test evaluated several different combinations of jet locations and mass flow rates as well as a vortex generator (VG) vane case. The numerical simulations will be performed on a single jet configuration for varying actuator mass flow rates at a fix inlet mass flow condition. Validation of the numerical simulations for the VG vane case will also be performed for varying inlet mass flow rates. Overall, the numerical simulations were able to predict the baseline circumferential flow distortion, DPCPavg, very well for comparisons made within the designed operating range of the BLI inlet. However the CFD simulations did predict a total pressure recovery that was 0.01 lower than the experiment. Numerical simulations of the baseline inlet flow also showed good agreement with the experimental inlet centerline surface pressures. The vane case showed that the CFD predicted the correct trends in the circumferential distortion for varying inlet mass flow but had a distortion level that was nearly twice as large as the experiment. Comparison to circumferential distortion measurements for a 15 deg clocked 40 probe rake indicated that the circumferential distortion levels are very sensitive to the symmetry of the flow and that a miss alignment of the vanes in the experiment could have resulted in this difference. The numerical simulations of the BLI inlet with jets showed good agreement with the circumferential inlet distortion levels for a range of jet actuator mass flow ratios at a fixed inlet mass flow rate. The CFD simulations for the jet case also predicted an average total pressure recovery that was 0.01 lower than the experiment as was seen in the baseline. Comparison of the flow features the jet case revealed that the CFD predicted a much larger vortex at the engine fan-face when compare to the experiment.
Experimental data for two three-dimensional intersecting shock-wave/turbulent boundary-layer interaction flows at Mach 8.3 are presented. The test bodies, composed of two sharp fins fastened to a flat plate test bed, were designed to generate flows with varying degrees of pressure gradient, boundary-layer separation, and turning angle. The data include surface pressure and heat transfer distributions as well as mean flow field surveys both in the undisturbed and interaction regimes. The data are presented in a convenient form to be used to validate existing or future computational models of these hypersonic flows.
This research will investigate the use of Design-of-Experiments (DOE) in the development of an optimal passive flow control vane design for a boundary-layer-ingesting (BLI) offset inlet in transonic flow. This inlet flow control is designed to minimize the engine fan-face distortion levels and first five Fourier harmonic half amplitudes while maximizing the inlet pressure recovery. Numerical simulations of the BLI inlet are computed using the Reynolds-averaged Navier-Stokes (RANS) flow solver, OVERFLOW, developed at NASA. These simulations are used to generate the numerical experiments for the DOE response surface model. In this investigation, two DOE optimizations were performed using a D-Optimal Response Surface model. The first DOE optimization was performed using four design factors which were vane height and angles-of-attack for two groups of vanes. One group of vanes was placed at the bottom of the inlet and a second group symmetrically on the sides. The DOE design was performed for a BLI inlet with a free-stream Mach number of 0.85 and a Reynolds number of 2 million, based on the length of the fan-face diameter, matching an experimental wind tunnel BLI inlet test. The first DOE optimization required a fifth order model having 173 numerical simulation experiments and was able to reduce the DC60 baseline distortion from 64% down to 4.4%, while holding the pressure recovery constant. A second DOE optimization was performed holding the vanes heights at a constant value from the first DOE optimization with the two vane angles-of-attack as design factors. This DOE only required a second order model fit with 15 numerical simulation experiments and reduced DC60 to 3.5% with small decreases in the fourth and fifth harmonic amplitudes. The second optimal vane design was tested at the NASA Langley 0.3- Meter Transonic Cryogenic Tunnel in a BLI inlet experiment. The experimental results showed a 80% reduction of DPCP(sub avg), the circumferential distortion level at the engine fan-face.
This research will investigate the use of Design-of-Experiments (DOE) in the development of an optimal passive flow control vane design for a boundary-layer-ingesting (BLI) offset inlet in transonic flow. This inlet flow control is designed to minimize the engine fan face distortion levels and first five Fourier harmonic half amplitudes while maximizing the inlet pressure recovery. Numerical simulations of the BLI inlet are computed using the Reynolds-averaged Navier-Stokes (RANS) flow solver, OVERFLOW, developed at NASA. These simulations are used to generate the numerical experiments for the DOE response surface model. In this investigation, two DOE optimizations were performed using a D-Optimal Response Surface model. The first DOE optimization was performed using four design factors which were vane height and angles-of-attack for two groups of vanes. One group of vanes was placed at the bottom of the inlet and a second group symmetrically on the sides. The DOE design was performed for a BLI inlet with a free-stream Mach number of 0.85 and a Reynolds number of 2 million, based on the length of the fan face diameter, matching an experimental wind tunnel BLI inlet test. The first DOE optimization required a fifth order model having 173 numerical simulation experiments and was able to reduce the DC60 baseline distortion from 64% down to 4.4%, while holding the pressure recovery constant. A second DOE optimization was performed holding the vanes heights at a constant value from the first DOE optimization with the two vane angles-of-attack as design factors. This DOE only required a second order model fit with 15 numerical simulation experiments and reduced DC60 to 3.5% with small decreases in the fourth and fifth harmonic amplitudes. The second optimal vane design was tested at the NASA Langley 0.3-Meter Transonic Cryogenic Tunnel in a BLI inlet experiment. The experimental results showed a 80% reduction of DPCPavg, the circumferential distortion level at the engine fan face.