Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “viscous effects”

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 325 records · Page 18

The computation and validation of hovering rotor performance

Recent experience with the HELIX-I code is presented, and its ability to predict the flow and performance of both conventional rotors and the unconventional anhedral parabolic tip rotor utilized on the Super Puma MK2 is described. HELIX-I is a standard full-potential rotor code having the ability to efficiently predict the detailed flow on a rotor blade, including 3D, transonic, and weak viscous effects (using appropriate boundary layer analyses). The resulting code is the first full-potential CFD code with the ability to model free wake convection and the first CFD code of any type to predict hover performance.

Mba, M. N.↗

Navier-Stokes simulations of unsteady transonic flow phenomena

Numerical simulations of two classes of unsteady flows are obtained via the Navier-Stokes equations: a blast-wave/target interaction problem class and a transonic cavity flow problem class. The method developed for the viscous blast-wave/target interaction problem assumes a laminar, perfect gas implemented in a structured finite-volume framework. The approximately factored implicit scheme uses Newton subiterations to obtain the spatially and temporally second-order accurate time history of the blast-waves with stationary targets. The inviscid flux is evaluated using either of two upwind techniques, while the full viscous terms are computed by central differencing. Comparisons of unsteady numerical, analytical, and experimental results are made in two- and three-dimensions for Couette flows, a starting shock-tunnel, and a shock-tube blockage study. The results show accurate wave speed resolution and nonoscillatory discontinuity capturing of the predominantly inviscid flows. Viscous effects were increasingly significant at large post-interaction times. While the blast-wave/target interaction problem benefits from high-resolution methods applied to the Euler terms, the transonic cavity flow problem requires the use of an efficient scheme implemented in a geometrically flexible overset mesh environment. Hence, the Reynolds averaged Navier-Stokes equations implemented in a diagonal form are applied to the cavity flow class of problems. Comparisons between numerical and experimental results are made in two-dimensions for free shear layers and both rectangular and quieted cavities, and in three-dimensions for Stratospheric Observatory For Infrared Astronomy (SOFIA) geometries. The acoustic behavior of the rectangular and three-dimensional cavity flows compare well with experiment in terms of frequency, magnitude, and quieting trends. However, there is a more rapid decrease in computed acoustic energy with frequency than observed experimentally owing to numerical dissipation. In addition, optical phase distortion due to the time-varying density field is modelled using geometrical constructs. The computed optical distortion trends compare with the experimentally inferred result, but underpredicts the fluctuating phase difference magnitude.

Atwood, C. A.↗

On the structure of gaseous confined laminar diffusion flames: Numerical investigation

The structure and characteristics of gaseous confined laminar diffusion flames are investigated by numerically solving the time-dependent two-dimensional axisymmetric conservation equations. The numerical model accounts for the important chemical and physical processes involved, including axial diffusion, viscous effects, radial convection, and finite-rate chemistry. The numerical results clearly show that the flame has a finite thickness and leakage of fuel vapor into the flame zone is possible. The effect of heat release is found to induce some radial flow. Predicted flame shape and dimensions are compared to the classical Burke-Schumann flame. The numerically calculated flame is observed to be about 15 percent taller and 5 percent narrower than that of the Burke-Schumann solution under the same conditions.

Mawid, M. A.↗

Prediction of forces and moments for hypersonic flight vehicle control effectors

This research project includes three distinct phases. For completeness, all three phases of the work are briefly described in this report. The goal was to develop methods of predicting flight control forces and moments for hypersonic vehicles which could be used in a preliminary design environment. The first phase included a preliminary assessment of subsonic/supersonic panel methods and hypersonic local flow inclination methods for such predictions. While these findings clearly indicated the usefulness of such methods for conceptual design activities, deficiencies exist in some areas. Thus, a second phase of research was conducted in which a better understanding was sought for the reasons behind the successes and failures of the methods considered, particularly for the cases at hypersonic Mach numbers. This second phase involved using computational fluid dynamics methods to examine the flow fields in detail. Through these detailed predictions, the deficiencies in the simple surface inclination methods were determined. In the third phase of this work, an improvement to the surface inclination methods was developed. This used a novel method for including viscous effects by modifying the geometry to include the viscous/shock layer.

Maughmer, Mark D.↗

A robust quasi-simultaneous interaction method for a full potential flow with a boundary layer with application to wing/body configurations

The MATRICS flow solver calculates the inviscid transonic potential flow about a wing/body semi-configuration. At present, work is in progress to extend MATRICS to take viscous effects into account through coupling with a boundary layer solver. This solver, MATRICS-V, is based on robust calculation methods for the boundary layer, the outer wing flow and their interaction. MATRICS-V is intended for (inverse) design purposes. The boundary layer and wake are based on an integral formulation of the unsteady first order boundary layer equations, the inviscid method is the existing MATRICS potential flow solver, and the interaction algorithm is of the quasi-simultaneous type. The paper gives a progress report on the coupled potential-flow boundary-layer method for transonic wing/body configurations.

Vanderwees, A. J.↗

Current status of computational methods for transonic unsteady aerodynamics and aeroelastic applications

The current status of computational methods for unsteady aerodynamics and aeroelasticity is reviewed. The key features of challenging aeroelastic applications are discussed in terms of the flowfield state: low-angle high speed flows and high-angle vortex-dominated flows. The critical role played by viscous effects in determining aeroelastic stability for conditions of incipient flow separation is stressed. The need for a variety of flow modeling tools, from linear formulations to implementations of the Navier-Stokes equations, is emphasized. Estimates of computer run times for flutter calculations using several computational methods are given. Application of these methods for unsteady aerodynamic and transonic flutter calculations for airfoils, wings, and configurations are summarized. Finally, recommendations are made concerning future research directions.

Edwards, J. W.↗

Airfoil shape optimization using sensitivity analysis on viscous flow equations

An aerodynamic shape optimization method has previously been developed by the authors using the Euler equations and has been applied to supersonic-hypersonic nozzle designs. This method has also included a flowfield extrapolation (or flow prediction) method based on the Taylor series expansion of an existing CFD solution. The present paper reports on the extension of this method to the thin-layer Navier-Stokes equations in order to account for the viscous effects. Also, to test the method under highly nonlinear conditions, it has been applied to the transonic flows. Initially, the success of the flow prediction method is tested. Then, the overall method is demonstrated by optimizing the shapes of two supercritical transonic airfoils at zero angle of attack. The first one is shape optimized to achieve a minimum drag while obtaining a lift above a specified value. Whereas, the second one is shape optimized for a maximum lift while attaining a drag below a specified value. The results of these two cases indicate that the present method can produce successfully optimized aerodynamic shapes.

Eleshaky, Mohamed E.↗

Implicit multigrid techniques for compressible flows

Recent advances in the development of the diagonalized alternating direction implicit multigrid method for compressible aerodynamic problems are reviewed. These include the extension of the method originally developed for the Euler equations to include viscous effects, the computation of turbulent flows and the implementation on parallel computers of the scheme on multiblock rids.

Caughey, David A.↗

Line relaxation methods for the solution of 2D and 3D compressible flows

An implicit finite element based algorithm for the compressible Navier-Stokes equations is outlined, and the solution of the resulting equation by a line relaxation on general meshes of triangles or tetrahedra is described. The problem of generating and adapting unstructured meshes for viscous flows is reexamined, and an approach for both 2D and 3D simulations is proposed. An efficient approach appears to be the use of an implicit/explicit procedure, with the implicit treatment being restricted to those regions of the mesh where viscous effects are known to be dominant. Numerical examples demonstrating the computational performance of the proposed techniques are given.

Hassan, O.↗

A device for measuring sonic velocity and compressor Mach number

A device has been developed which measures the velocity of sound in fluids at stagnation and is especially adaptable to turbine and compressor testing for which the constituency of the working fluid may be in doubt. By utilizing the shaft frequency of a rotary compressor, the instrument can also be used to provide a direct measurement of the compressor Mach number (ratio of blade-tip velocity to inlet velocity of sound at stagnation). A Helmholtz resonator is employed in the measurement of the sound velocity. Viscous effects in the orifice of the Helmholtz resonator are shown to be important and can be taken into account with the help of a parameter obtained from Stokes solution of the flow near an oscillating wall. This parameter includes the kinematic viscosity of the fluid and the frequency of sound in the resonator. When these effects are recognized, the resonator can be calibrated to measure velocity of sound or compressor Mach number to an accuracy of better than 0.5 percent.

Huber, Paul W↗

Design and Test of Mixed-flow Impellers III : Design and Experimental Results for Impeller Model MFI-2A and Comparison with Impeller Model MFI-1A

A mixed-flow impeller was designed to give a prescribed blade-surface velocity distribution at mean blade height for a given hub-shroud profile. The blade shape at mean blade height, which was produced by the prescribed velocity distribution, was extended by means of radial lines to form the composite blade shape from hub to shroud. The resulting blade was relatively thick; therefore, it was necessary to retain the inverse blade taper which resulted from extension of the radial lines in order to prevent merging or near merging of the separate blades near the hub. For the first test version of the impeller, designated the MFI-2A, the blade height was arbitrarily made greater than that for the basic impeller (the MFI-2) to allow for viscous effects. At design equivalent speed of 1400 feet per second the peak pressure ratio and maximum adiabatic efficiency were 3.95 and 79 percent, respectively. The adiabatic efficiency of the MFI-2A is four points lower than that for impeller model MFI-1A, but because of the higher slip factor for the MFI-2A, the pressure ratios are approximately equal. The procedures followed in the design of the MFI-1A and MFI-2A were, in general, the same; and, although the prescribed initial condition resulted in geometrical configurations that were quite dissimilar, the resulting performance characteristics compare favorably with designs for which considerable development work has been necessary.

COMPRESSION AND COMPRESSORS↗

A CFD-based aerodynamic design procedure for hypersonic wind-tunnel nozzles

A new procedure which unifies the best of current classical design practices, computational fluid dynamics (CFD), and optimization procedures is demonstrated for designing the aerodynamic lines of hypersonic wind-tunnel nozzles. The new procedure can be used to design hypersonic wind tunnel nozzles with thick boundary layers where the classical design procedure has been shown to break down. An efficient CFD code, which solves the parabolized Navier-Stokes (PNS) equations using an explicit upwind algorithm, is coupled to a least-squares (LS) optimization procedure. A LS problem is formulated to minimize the difference between the computed flow field and the objective function, consisting of the centerline Mach number distribution and the exit Mach number and flow angle profiles. The aerodynamic lines of the nozzle are defined using a cubic spline, the slopes of which are optimized with the design procedure. The advantages of the new procedure are that it allows full use of powerful CFD codes in the design process, solves an optimization problem to determine the new contour, can be used to design new nozzles or improve sections of existing nozzles, and automatically compensates the nozzle contour for viscous effects as part of the unified design procedure. The new procedure is demonstrated by designing two Mach 15, a Mach 12, and a Mach 18 helium nozzles. The flexibility of the procedure is demonstrated by designing the two Mach 15 nozzles using different constraints, the first nozzle for a fixed length and exit diameter and the second nozzle for a fixed length and throat diameter. The computed flow field for the Mach 15 least squares parabolized Navier-Stokes (LS/PNS) designed nozzle is compared with the classically designed nozzle and demonstrates a significant improvement in the flow expansion process and uniform core region.

Korte, John J.↗

Numerical investigations in three-dimensional internal flows

In the present reporting period, the 3D version of the OVERFLOW code was used to solve the flow within the internal portion of the supersonic inlet. The internal portion of this inlet is bounded by an inflow plane containing the leading edge of the sidewalls, the sidewalls, the ramp and cowl surfaces and an outflow plane just downstream of the minimum geometric area of the inlet. Boundary layer bleed was used in the two-dimensional calculations discussed in the previous progress report and that same bleed was applied in the present study. For reference, this bleed corresponds to locations designated as R2 and R3 in the Mach 5 inlet model test. Using the GRIDGEN code, a three dimensional grid was generated that accounted for the viscous effects expected to occur on the sidewall, as well as those known to occur on the ramp and cowl surfaces. The internal flow grid size was 141 streamwise by 101 cross stream by 71 in the lateral direction between sidewalls. Since the flow entering the inlet was not symmetrical, the inlet was solved from sidewall to sidewall (without using a symmetry plane). In addition to the short sidewalls proposed in the Langley geometry database, a set of shorter sidewalls was also investigated in the present study and was shown to have beneficial effects with respect to the flow distortion exiting the supersonic inlet. In addition to these calculations, additional 3D solutions using the OVERFLOW code were obtained for the flow downstream of the throat of the supersonic inlet, including a terminal shock wave system produced by a backpressured subsonic diffuser.

Rose, William C.↗

In-flight pressure distributions and skin-friction measurements on a subsonic transport high-lift wing section

Flight experiments are being conducted as part of a multiphased subsonic transport high-lift research program for correlation with wind-tunnel and computational results. The NASA Langley Transport Systems Research Vehicle (B737-100 aircraft) is used to obtain in-flight flow characteristics at full-scale Reynolds numbers to contribute to the understanding of 3-D high-lift, multi-element flows including attachment-line transition and relaminarization, confluent boundary-layer development, and flow separation characteristics. Flight test results of pressure distributions and skin friction measurements were obtained for a full-chord wing section including the slat, main-wing, and triple-slotted, Fowler flap elements. Test conditions included a range of flap deflections, chord Reynolds numbers (10 to 21 million), and Mach numbers (0.16 to 0.40). Pressure distributions were obtained at 144 chordwise locations of a wing section (53-percent wing span) using thin pressure belts over the slat, main-wing, and flap elements. Flow characteristics observed in the chordwise pressure distributions included leading-edge regions of high subsonic flows, leading-edge attachment-line locations, slat and main-wing cove-flow separation and reattachment, and trailing-edge flap separation. In addition to the pressure distributions, limited skin-friction measurements were made using Preston-tube probes. Preston-tube measurements on the slat upper surface suggested relaminarization of the turbulent flow introduced by the pressure belt on the slat leading-edge surface when the slat attachment line was laminar. Computational analysis of the in-flight pressure measurements using two-dimensional, viscous multielement methods modified with simple-sweep theory showed reasonable agreement. However, overprediction of the pressures on the flap elements suggests a need for better detailed measurements and improved modeling of confluent boundary layers as well as inclusion of three-dimensional viscous effects in the analysis.

Yip, Long P.↗

Experimental investigation of nozzle/plume aerodynamics at hypersonic speeds

Work continued on the improvement of 16-Inch Shock Tunnel. This comprised studies of ways of improving driver gas ignition, an improved driver gas mixing system, an axial translation system for the driver tube, improved diaphragm materials (carbon steel vs. stainless steel), a copper liner for the part of the driven tube near the nozzle, the use of a buffer gas between the driver and driven gases, the use of N2O in the driven tube, the use of a converging driven tube, operation of the facility as a non-reflected shock tunnel and expansion tube, operation with heated hydrogen or helium driver gas, the use of detonations in the driver and the construction of an enlarged test section. Maintenance and developmental work continued on the scramjet combustor continued. New software which greatly speeds up data analysis has been written and brought on line. In particular, software which provides very rapid generation of model surface heat flux profiles has been brought on line. A considerable amount of theoretical work was performed in connection with upgrading the 16 Inch Shock Tunnel Facility. A one-dimensional Godunov code for very high velocities and any equation of state is intended to add viscous effects in studying the operation of the Shock Tunnel and also of two-stage light gas guns.

Bogdanoff, David W.↗

Transient Ejector Analysis (TEA) code user's guide

A FORTRAN computer program for the semi analytic prediction of unsteady thrust augmenting ejector performance has been developed, based on a theoretical analysis for ejectors. That analysis blends classic self-similar turbulent jet descriptions with control-volume mixing region elements. Division of the ejector into an inlet, diffuser, and mixing region allowed flexibility in the modeling of the physics for each region. In particular, the inlet and diffuser analyses are simplified by a quasi-steady-analysis, justified by the assumption that pressure is the forcing function in those regions. Only the mixing region is assumed to be dominated by viscous effects. The present work provides an overview of the code structure, a description of the required input and output data file formats, and the results for a test case. Since there are limitations to the code for applications outside the bounds of the test case, the user should consider TEA as a research code (not as a production code), designed specifically as an implementation of the proposed ejector theory. Program error flags are discussed, and some diagnostic routines are presented.

Drummond, Colin K.↗

Design of a Low Aspect Ratio Transonic Compressor Stage Using CFD Techniques

A transonic compressor stage has been designed for the Naval Postgraduate School Turbopropulsion Laboratory. The design relied heavily on CFD techniques while minimizing conventional empirical design methods. The low aspect ratio (1.2) rotor has been designed for a specific head ratio of .25 and a tip relative inlet Mach number of 1.3. Overall stage pressure ratio is 1.56. The rotor was designed using an Euler code augmented by a distributed body force model to account for viscous effects. This provided a relatively quick-running design tool, and was used for both rotor and stator calculations. The initial stator sections were sized using a compressible, cascade panel code. In addition to being used as a case study for teaching purposes, the compressor stage will be used as a research stage. Detailed measurements, including non-intrusive LDV, will be compared with the design computations, and with the results of other CFD codes, as a means of assessing and improving the computational codes as design tools.

Sanger, Nelson L.↗

Computational fluid dynamics studies of nuclear rocket performance

A CFD analysis of a low pressure nuclear rocket concept is presented with the use of an advanced chemical kinetics, Navier-Stokes code. The computations describe the flow field in detail, including gas dynamic, thermodynamic and chemical properties, as well as global performance quantities such as specific impulse. Computational studies of several rocket nozzle shapes are conducted in an attempt to maximize hydrogen recombination. These Navier-Stokes calculations, which include real gas and viscous effects, predict lower performance values than have been reported heretofore.

Stubbs, Robert M.↗