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 631 records · Page 35

A Structured Grid Based Solution-Adaptive Technique for Complex Separated Flows

The objective of this work was to enhance the predictive capability of widely used computational fluid dynamic (CFD) codes through the use of solution adaptive gridding. Most problems of engineering interest involve multi-block grids and widely disparate length scales. Hence, it is desirable that the adaptive grid feature detection algorithm be developed to recognize flow structures of different type as well as differing intensity, and adequately address scaling and normalization across blocks. In order to study the accuracy and efficiency improvements due to the grid adaptation, it is necessary to quantify grid size and distribution requirements as well as computational times of non-adapted solutions. Flow fields about launch vehicles of practical interest often involve supersonic freestream conditions at angle of attack exhibiting large scale separate vortical flow, vortex-vortex and vortex-surface interactions, separated shear layers and multiple shocks of different intensity. In this work, a weight function and an associated mesh redistribution procedure is presented which detects and resolves these features without user intervention. Particular emphasis has been placed upon accurate resolution of expansion regions and boundary layers. Flow past a wedge at Mach=2.0 is used to illustrate the enhanced detection capabilities of this newly developed weight function.

Thornburg, Hugh↗

A block-corrected subdomain solution procedure for recirculating flow calculations

This paper describes a robust and efficient subdomain solution procedure for two-dimensional recirculating flows. The solution domain is divided into a number of overlapping subdomains, and a direct fully coupled solution is obtained for each subdomain using a sparse matrix form of LU decomposition. An effective parabolic block correction procedure, which calculates global corrections to the tentative solution by a marching technique similar to that used for boundary layer flows, is used to accelerate the convergence of the basic procedure. The use of effective block correction is found to be essential for the success of the subdomain approach on strongly recirculating flows. In a number of laminar two-dimensional flows, the new block-corrected method performed extremely well, rivaling the best direct methods in execution time, while requiring substantially less computer storage. The new method proved to be from two to ten times faster than conventional iterative methods, while requiring only a moderate increase in storage.

Braaten, M. E.↗

Comparison of predicted and experimental external heat transfer around a film cooled cylinder in crossflow

Calculations were made of the film cooling provided by rows of holes around the circumference of a cylinder in crossflow and the results were compared to experimental data. The calculations and experimental data were for conditions that simulate most of those that are typical of air cooled turbine vane leading edges. Injection was from a single and multiple rows of holes located at different angular locations from the stagnation line. The holes in the rows were angled normal to the flow direction and at a 25 degree angle to the cylinder wall. The calculations and experimental data were for several constant values of blowing ratios for all rows and for different blowing ratios for each row, representing a simulation of a common coolant plenum supply to multiple rows of holes. The calculations were made using a finite difference boundary layer code, STAN5. Contrary to initial expectations that injection would trip the boundary layer flow into the turbulent regime, the results indicated that the high free stream acceleration apparently kept the flow laminar for holes in the first 45 degrees past stagnation. The trend in Stanton number reduction due to coolant injection was predicted with generally good agreement at the lower blowing rates, but for multile rows of holes, agreement was poor beyond the first row.

Stepka, F. S.↗

Flow Observations with Tufts and Lampblack of the Stalling of Four Typical Airfoil Sections in the NACA Variable-density Tunnel

A preliminary investigation of the stalling processes of four typical airfoil sections was made over the critical range of the Reynolds Number. Motion pictures were taken of the movements of small silk tufts on the airfoil surface as the angle of attack increased through a range of angles including the stall. The boundary-layer flow also at certain angles of attack was indicated by the patterns formed by a suspension of lampblack in oil brushed onto the airfoil surface. These observations were analyzed together with corresponding force-test measurements to derive a picture of the stalling processes of airfoils.

Abbott, Ira H↗

Development of a Low-Reynolds Number, Nonlinear kappa-epsilon Model for the Reduced Navier-Stokes Equations

Previous work at NASA LeRC has shown that flow distortions in aircraft engine inlet ducts can be significantly reduced by mounting vortex generators, or small wing sections, on the inside surface of the engine inlet. The placement of the vortex generators is an important factor in obtaining the optimal effect over a wide operating envelope. In this regard, the only alternative to a long and expensive test program which would search out this optimal configuration is a good prediction procedure which could narrow the field of search. Such a procedure has been developed in collaboration with NASA LeRC, and results obtained by NASA personnel indicate that it shows considerable promise for predicting the viscous turbulent flow in engine inlet ducts in the presence of vortex generators. The prediction tool is a computer code which numerically solves the reduced Navier-Stokes equations and so is commonly referred to as RNS3D. Obvious deficiencies in RNS3D have been addressed in previous work. Primarily, it is known that the predictions of the mean velocity field of a turbulent boundary layer flow approaching separation are not in good agreement with data. It was suggested that the use of an algebraic mixing-length turbulence model in RNS3D is at least partly to blame for this. Additionally, the current turbulence model includes an assumption of isotropy which will ultimately fail to capture turbulence-driven secondary flow known to exist in noncircular ducts.

Boger, David A.↗

Nonlinear Evolution and Breakdown of Azimuthally Compact Crossflow Vortex Pattern over a Yawed Cone

Hypersonic boundary-layer flows over a circular cone at moderate angle of incidence can support strong crossflow instability in between the windward and leeward rays on the plane of symmetry. Due to the more efficient excitation of stationary crossflow vortices by surface roughness, a possible path to transition in such flows corresponds to rapid amplification of the high-frequency secondary instabilities of finite amplitude stationary crossflow vortices. In the present paper, the previous analyses of crossflow instability over a 7- degree half-angle, yawed circular cone in a Mach 6 free stream have been extended to the nonlinear evolution of azimuthally localized crossflow vortex packets and the amplification characteristics and nonlinear breakdown of high-frequency secondary instabilities associated with those packets. A comparison between plane marching PSE and direct Navier-Stokes simulations (DNS) reveals favorable agreement in regard to mode shapes, most amplified disturbance frequencies, and N-factor evolution. In contrast, the quasi-parallel predictions are found to result in severe underprediction of the N-factors. The direct numerical simulations also indicate that the breakdown of secondary instabilities in a 3D hypersonic boundary layer shares certain common features with the previous computations of crossflow transition over subsonic swept wings.

Vortex structure↗

Turbulence processes and simple closure schemes

The problem of closure in turbulence in the case of two-point correlations resides in the existence of two unknowns E and W, the energy spectrum function and the transfer function, respectively, in the spectrum equation. In the case of weak turbulence, W is negligible. In case of higher correlations, closure can be effective by neglecting the inertia term in the highest order term used. Specifying a certain number of spectra at an initial time is also a way of getting around the closure problem. A simple case of turbulent shear flow is then considered, where two-point correlation equations are used and the velocity is broken into mean and fluctuating components. This yields a differential equation for the energy spectrum, the three terms of which are the energy spectrum, production term and dissipation term. They are plotted for a particular time. Similar analyses and comparisons with experiment are made for pipe and boundary layer flows.

Deissler, R. G.↗

Solution of 3-dimensional time-dependent viscous flows. Part 3: Application to turbulent and unsteady flows

A numerical scheme is developed for solving the time dependent, three dimensional compressible viscous flow equations to be used as an aid in the design of helicopter rotors. In order to further investigate the numerical procedure, the computer code developed to solve an approximate form of the three dimensional unsteady Navier-Stokes equations employing a linearized block implicit technique in conjunction with a QR operator scheme is tested. Results of calculations are presented for several two dimensional boundary layer flows including steady turbulent and unsteady laminar cases. A comparison of fourth order and second order solutions indicate that increased accuracy can be obtained without any significant increases in cost (run time). The results of the computations also indicate that the computer code can be applied to more complex flows such as those encountered on rotating airfoils. The geometry of a symmetric NACA four digit airfoil is considered and the appropriate geometrical properties are computed.

Weinberg, B. C.↗

The effect of a leading-edge slat on the dynamic stall of an oscillating airfoil

The dynamic stall characteristics of a slatted airfoil were investigated experimentally on a 2-ft-chord airfoil oscillating in pitch at M = 0.2 for a range of reduced frequency and mean angle of oscillation. The slat produced a flow that remained attached to the airfoil for angles well above those normally attained by the retreating blade of a helicopter during high speed flight. The dynamic stall vortex usually associated with these flight conditions was completely eliminated for all angles under 30 deg. Instantaneous surface pressure, lift, and pitching moment data are presented as a function of incidence throughout the oscillation cycle; a detailed analysis of instantaneous boundary-layer flow behavior for the various test conditions is included.

Carr, L. W.↗

Effects of Cross-Sectional Shape, Solidity, and Distribution of Heat-Transfer Coefficient on the Torsional Stiffness of Thin Wings Subjected to Aerodynamic Heating

A study has been made of the effects of varying the shape, solidity, and heat-transfer coefficient of thin wings with regard to their influence on the torsional-stiffness reduction induced by aerodynamic heating. The variations in airfoil shape include blunting, flattening, and combined blunting and flattening of a solid wing of symmetrical double-wedge cross section. Hollow double-wedge wings of constant skin thickness with and without internal webs also are considered. The effects of heat-transfer coefficients appropriate for laminar and turbulent flow are investigated in addition to a step transition along the chord from a lower to a higher constant value of heat-transfer coefficient. From the results given it is concluded that the flattening of a solid double wedge decreases the reduction in torsional stiffness while slight degrees of blunting increase the loss. The influence of chordwise variations in heat-transfer coefficient due to turbulent and laminar boundary-layer flow on the torsional stiffness of solid wings is negligible. The effect of a step transition in heat-transfer coefficient along the chord of a solid wing can, however, become appreciable. The torsional-stiffness reduction of multiweb and hollow double-wedge wings is substantially less than that calculated for a solid wing subjected to the same heating conditions.

Thomson, Robert G.↗

Implicit Large-Eddy Simulations of Zero-Pressure Gradient, Turbulent Boundary Layer

A set of direct simulations of zero-pressure gradient, turbulent boundary layer flows are conducted using various span widths (62-630 wall units), to document their influence on the generated turbulence. The FDL3DI code that solves compressible Navier-Stokes equations using high-order compact-difference scheme and filter, with the standard recycling/rescaling method of turbulence generation, is used. Results are analyzed at two different Re values (500 and 1,400), and compared with spectral DNS data. They show that a minimum span width is required for the mere initiation of numerical turbulence. Narrower domains ((is) less than 100 w.u.) result in relaminarization. Wider spans ((is) greater than 600 w.u.) are required for the turbulent statistics to match reference DNS. The upper-wall boundary condition for this setup spawns marginal deviations in the mean velocity and Reynolds stress profiles, particularly in the buffer region.

Turblent boundary layer↗

Experimental study of low Reynolds number nozzles

High-performance electrothermal thrusters operate in a low nozzle-throat Reynolds number regime. Under these conditions, the flow boundary layer occupies a large volume inside the nozzle, contributing to large viscous losses. Four nozzles (conical, bell, trumpet, and modified trumpet) and a sharp-edged orifice were evaluated over a Reynolds number range of 500 to 9000 with unheated nitrogen and hydrogen. The nozzles showed significant decreases in specific impulse efficiency with decreasing Reynolds number. At Reynolds numbers less than 1000, all four nozzles were probably filled with a large boundary layer. The discharge coefficient decreased with Reynolds number in the same manner as the specific impulse efficiency. The bell and modified trumpet nozzles had discharge coefficients 4 to 8 percent higher than those of the cone or trumpet nozzles. The Two-Dimensional Kinetics (TDK) nozzle analysis computer program was used to predict nozzle performance. The results were then compared to the experimental results in order to determine the accuracy of the program within this flow regime.

Grisnik, Stanley P.↗

Experimental study of low Reynolds number nozzles

High-performance electrothermal thrusters operate in a low nozzle-throat Reynolds number regime. Under these conditions, the flow boundary layer occupies a large volume inside the nozzle, contributing to large viscous losses. Four nozzles (conical, bell, trumpet, and modified trumpet) and a sharp-edged orifice were evaluated over a Reynolds number range of 500 to 9000 with unheated nitrogen and hydrogen. The nozzles showed significant decreases in specific impulse efficiency with decreasing Reynolds number. At Reynolds numbers less than 1000, all four nozzles were probably filled with a large boundary layer. The discharge coefficient decreased with Reynolds number in the same manner as the specific impulse efficiency. The bell and modified trumpet nozzles had discharge coefficients 4 to 8 percent higher than those of the cone or trumpet nozzles. The Two-Dimensional Kinetics (TDK) nozzle analysis computer program was used to predict nozzle performance. The results were then compared to the experimental results in order to determine the accuracy of the program within this flow regime.

Grisnik, Stanley P.↗

Turbulent flow between a rotating disk and a parallel wall

A scalar effective viscosity method is used to calculate the incompressible flow in the gap between two infinite, parallel smooth disks, one of them rotating and the other stationary. The flow is calculated for two values of the gap Reynolds number omega(s squared)/nu, namely 296 and 2852, where omega is the disk angular velocity, s the gap width, and nu the kinematic viscosity of the fluid; and over a range of radial Reynolds number omega(r squared)/nu from 0 to 10 million, where r is the radius from the axis of rotation. A numerical shooting technique is used to match the two boundary-layer flows arising on the disk and wall. Gradual transition is assumed for radial Reynolds numbers from 160,000 to 250,000. The results are in good agreement with Daily and Nece's measurements of velocity profiles and disk friction drag for enclosed rotating disks, and at sufficiently large radius the local skin friction approaches that of turbulent Couette flow.

Cooper, P.↗

Feasibility study of optical boundary layer transition detection method

A high sensitivity differential interferometer was developed to locate the region where the boundary layer flow undergoes transition from laminar to turbulent. Two laboratory experimental configurations were used to evaluate the performance of the interferometer: open shear layer, and low speed wind tunnel turbulent spot configuration. In each experiment, small temperature fluctuations were introduced as the signal source. Simultaneous cold wire measurements were compared with the interferometer data. The comparison shows that the interferometer is sensitive to very weak phase variations in the order of 0.001 the laser wavelength. An attempt to detect boundary layer transition over a flat plate at NASA-Langley Unitary Supersonic Wind Tunnel using the interferometer system was performed. The phase variations during boundary layer transition in the supersonic wind tunnel were beyond the minimum signal-to-noise level of the instrument.

Azzazy, M.↗

Phase Synchronization and Desynchronization of Structural Response Induced by Turbulent and External Sound

Acoustic and turbulent boundary layer flow loadings over a flexible structure are used to study the spatial-temporal dynamics of the response of the structure. The stability of the spatial synchronization and desynchronization by an active external force is investigated with an array of coupled transducers on the structure. In the synchronous state, the structural phase is locked, which leads to the formation of spatial patterns while the amplitude peaks exhibit chaotic behaviors. Large amplitude, spatially symmetric loading is superimposed on broadband, but in the desynchronized state, the spectrum broadens and the phase space is lost. The resulting pattern bears a striking resemblance to phase turbulence. The transition is achieved by using a low power external actuator to trigger broadband behaviors from the knowledge of the external acoustic load inducing synchronization. The changes are made favorably and efficiently to alter the frequency distribution of power, not the total power level. Before synchronization effects are seen, the panel response to the turbulent boundary layer loading is discontinuously spatio-temporally correlated. The stability develops from different competing wavelengths; the spatial scale is significantly shorter than when forced with the superimposed external sound. When the external sound level decreases and the synchronized phases are lost, changes in the character of the spectra can be linked to the occurrence of spatial phase transition. These changes can develop broadband response. Synchronized responses of fuselage structure panels have been observed in subsonic and supersonic aircraft; results from two flights tests are discussed.

Maestrello, Lucio↗

On a finite element CFD algorithm for compressible, viscous and turbulent aerodynamic flows

This paper develops and analyses individual construction aspects of an efficient and accurate finite element algorithm for prediction of viscous and turbulent flow fields of impact in aerodynamics. The theoretical construction employs a Taylor weak statement (TWS) for coincident embedding of stability mechanisms within a classic Galerkin finite element formulation of semidiscrete approximation error orthogonalization. A wide variety of the stabilizing mechanisms of independently derived CFD algorithms are contained within the TWS theory. An implicit construction that meets the requirement of efficient convergence to steady state is developed. The theoretical asymptotic error estimates of the TWS finite element algorithm for supersonic and viscous boundary layer flows are verified. Application to a three-dimensional turbulent flow is cited.

Baker, A. J.↗