Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “compressible flows”

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 901 records · Page 50

Calculation of three-dimensional boundary layers on rotating turbine blades

An assessment has been made of the applicability of a three dimensional boundary layer analysis to the calculation of heat transfer and streamline flow patterns on the surfaces of both stationary and rotating turbine passages. In support of this effort, an analysis has been developed to calculate a general nonorthogonal surface coordinate system for arbitrary three dimensional surfaces and also to calculate the boundary layer edge conditions for compressible flow using the surface Euler equations and experimental pressure distributions. Using available experimental data to calibrate the method, calculations are presented for the endwall, and suction surfaces of a stationary cascade and for the pressure surface of a rotating turbine blade. The results strongly indicate that the three dimensional boundary layer analysis can give good predictions of the flow field and heat transfer on the pressure, suction, and endwall surfaces in a gas turbine passage.

Anderson, O. L.↗

User's manual for three dimensional boundary layer (BL3-D) code

An assessment has been made of the applicability of a 3-D boundary layer analysis to the calculation of heat transfer, total pressure losses, and streamline flow patterns on the surface of both stationary and rotating turbine passages. In support of this effort, an analysis has been developed to calculate a general nonorthogonal surface coordinate system for arbitrary 3-D surfaces and also to calculate the boundary layer edge conditions for compressible flow using the surface Euler equations and experimental data to calibrate the method, calculations are presented for the pressure endwall, and suction surfaces of a stationary cascade and for the pressure surface of a rotating turbine blade. The results strongly indicate that the 3-D boundary layer analysis can give good predictions of the flow field, loss, and heat transfer on the pressure, suction, and endwall surface of a gas turbine passage.

Anderson, O. L.↗

CFD analyses for advanced pump design

As one of the activities of the NASA/MSFC Pump Stage Technology Team, the present effort was focused on using CFD in the design and analysis of high performance rocket engine pumps. Under this effort, a three-dimensional Navier-Stokes code was used for various inducer and impeller flow field calculations. An existing algebraic grid generation procedure was-extended to allow for nonzero blade thickness, splitter blades, and hub/shroud cavities upstream or downstream of the (main) blades. This resulted in a fast, robust inducer/impeller geometry/grid generation package. Problems associated with running a compressible flow code to simulate an incompressible flow were resolved; related aspects of the numerical algorithm (viz., the matrix preconditioning, the artificial dissipation, and the treatment of low Mach number flows) were addressed. As shown by the calculations performed under the present effort, the resulting code, in conjunction with the grid generation package, is an effective tool for the rapid solution of three-dimensional viscous inducer and impeller flows.

Dejong, F. J.↗

Design and Test of Mixed-flow Impellers IV : Experimental Results for Impeller Models MFI-1 and MFI-2 with Changes in Blade Height

Modifications A and B of impeller model MFI-1 and A, B, and C of impeller model MFI-2 were investigated experimentally in an attempt to determine what allowance in blade height should be made for boundary layer and viscous losses in an impeller designed for isentropic compressible flow. A gradual increase in blade height was arbitrarily made from inlet to outlet in anticipation of a gradual build-up of boundary layer. Apparently there was a rapid build-up of boundary layer near the inlet in the experimental case rather than a gradual one. Therefore, the proper allowance for boundary layer cannot be described from the data obtained. Decreasing the pressure gradient along the shroud by reducing the blade height allowance apparently did little to increase the overall efficiency. At the design speed of 1400 feet per second, the overall adiabatic efficiency was increased from 0.83 for the MFI-1A to 0.85 for the MFI-1B with reduction in height; however, it is indicated from the theoretical velocity distribution and outlet surveys that the increase was due to a change from decelerating to accelerating flow along the hub rather than from any change along the shroud. It is further indicated that the consequences of a thickened or separated boundary layer depend not only on the design velocity gradients but also on the shape of the passage.

Hamrick, Joseph T↗

Vortex methods for separated flows

The numerical solution of the Euler or Navier-Stokes equations by Lagrangian vortex methods is discussed. The mathematical background is presented in an elementary fashion and includes the relationship with traditional point-vortex studies, the convergence to smooth solutions of the Euler equations, and the essential differences between two- and three-dimensional cases. The difficulties in extending the method to viscous or compressible flows are explained. The overlap with the excellent review articles available is kept to a minimum and more emphasis is placed on the area of expertise, namely two-dimensional flows around bluff bodies. When solid walls are present, complete mathematical models are not available and a more heuristic attitude must be adopted. The imposition of inviscid and viscous boundary conditions without conformal mappings or image vortices and the creation of vorticity along solid walls are examined in detail. Methods for boundary-layer treatment and the question of the Kutta condition are discussed. Practical aspects and tips helpful in creating a method that really works are explained. The topics include the robustness of the method and the assessment of accuracy, vortex-core profiles, timemarching schemes, numerical dissipation, and efficient programming. Calculations of flows past streamlined or bluff bodies are used as examples when appropriate.

Spalart, Philippe R.↗

Surface boundary conditions for the numerical solution of the Euler equations

We consider the implementation of boundary conditions at solid walls in inviscid Euler solutions by upwind, finite-volume methods. We review some current methods for the implementation of surface boundary conditions and examine their behavior for the problem of an oblique shock reflecting off a planar surface. We show the importance of characteristic boundary conditions for this problem and introduce a method of applying the classical flux-difference splitting of Roe as a characteristic boundary condition. Consideration of the equivalent problem of the intersection of two (equal and opposite) oblique shocks was very illuminating on the role of surface boundary conditions for an inviscid flow and led to the introduction of two new boundary-condition procedures, denoted as the symmetry technique and the curvature-corrected symmetry technique. Examples of the effects of the various surface boundary conditions considered are presented for the supersonic blunt body problem and the subcritical compressible flow over a circular cylinder. Dramatic advantages of the curvature-corrected symmetry technique over the other methods are shown, with regard to numerical entropy generation, total pressure loss, drag and grid convergence.

Dadone, A.↗

Real-Gas Effects on Binary Mixing Layers

This paper presents a computational study of real-gas effects on the mean flow and temporal stability of heptane/nitrogen and oxygen/hydrogen mixing layers at supercritical pressures. These layers consist of two counterflowing free streams of different composition, temperature, and density. As in related prior studies reported in NASA Tech Briefs, the governing conservation equations were the Navier-Stokes equations of compressible flow plus equations for the conservation of total energy and of chemical- species masses. In these equations, the expressions for heat fluxes and chemical-species mass fluxes were derived from fluctuation-dissipation theory and incorporate Soret and Dufour effects. Similarity equations for the streamwise velocity, temperature, and mass fractions were derived as approximations to the governing equations. Similarity profiles showed important real-gas, non-ideal-mixture effects, particularly for temperature, in departing from the error-function profile, which is the similarity solution for incompressible flow. The temperature behavior was attributed to real-gas thermodynamics and variations in Schmidt and Prandtl numbers. Temporal linear inviscid stability analyses were performed using the similarity and error-function profiles as the mean flow. For the similarity profiles, the growth rates were found to be larger and the wavelengths of highest instability shorter, relative to those of the errorfunction profiles and to those obtained from incompressible-flow stability analysis. The range of unstable wavelengths was found to be larger for the similarity profiles than for the error-function profiles

Okong'o, Nora↗

Development of the Glenn Heat-Transfer (Glenn-HT) Computer Code to Enable Time-Filtered Navier-Stokes (TFNS) Simulations and Application to Film Cooling on a Flat Plate Through Long Cooling Tubes

Computational fluid dynamics (CFD) analysis using Reynolds-averaged Navier-Stokes (RANS) formulation for turbomachinery-related flows has enabled improved engine component designs. RANS methodology has limitations that are related to its inability to accurately describe the spectrum of flow phenomena encountered in engines. Examples of flows that are difficult to compute accurately with RANS include phenomena such as laminar/turbulent transition, turbulent mixing due to mixing of streams, and separated flows. Large eddy simulation (LES) can improve accuracy but at a considerably higher cost. In recent years, hybrid schemes that take advantage of both unsteady RANS and LES have been proposed. This study investigated an alternative scheme, the time-filtered Navier-Stokes (TFNS) method applied to compressible flows. The method developed by Shih and Liu was implemented in the Glenn-Heat-Transfer (Glenn-HT) code and applied to film-cooling flows. In this report the method and its implementation is briefly described. The film effectiveness results obtained for film cooling from a row of 30deg holes with a pitch of 3.0 diameters emitting air at a nominal density ratio of unity and two blowing ratios of 0.5 and 1.0 are shown. Flow features under those conditions are also described.

Ameri, Ali↗

Development of the Glenn-HT Computer Code to Enable Time-Filtered Navier-Stokes (TFNS) Simulations and Application to Film Cooling on a Flat Plate Through Long Cooling Tubes

Computational fluid dynamics (CFD) analysis using Reynolds-averaged Navier-Stokes (RANS) formulation for turbomachinery-related flows has enabled improved engine component designs. RANS methodology has limitations which are related to its inability to accurately describe the spectrum of flow phenomena encountered in engines. Examples of flows that are difficult to compute accurately with RANS include phenomena such as laminarturbulent transition, turbulent mixing due to mixing of streams, and separated flows. Large eddy simulation (LES) can improve accuracy but at a considerably higher cost. In recent years, hybrid schemes which take advantage of both unsteady RANS and LES have been proposed. This study investigated an alternative scheme, the time-filtered Navier-Stokes (TFNS) method applied to compressible flows. The method developed by Shih and Liu was implemented in the Glenn-HT code and applied to film cooling flows. In this report the method and its implementation is briefly described. The film effectiveness results obtained for film cooling from a row of 30 holes with a pitch of 3.0 diameters emitting air at a nominal density ratio of unity and four blowing ratios of 0.5, 1.0, 1.5 and 2.0 are shown. Flow features under those conditions are also described.

infrared↗

Development of the Glenn-Heat-Transfer (Glenn-HT) Computer Code to Enable Time-Filtered Navier Stokes (TFNS) Simulations and Application to Film Cooling on a Flat Plate Through Long Cooling Tubes

Computational fluid dynamics (CFD) analysis using Reynolds-averaged Navier-Stokes (RANS) formulation for turbomachinery-related flows has enabled improved engine component designs. RANS methodology has limitations that are related to its inability to accurately describe the spectrum of flow phenomena encountered in engines. Examples of flows that are difficult to compute accurately with RANS include phenomena such as laminar/turbulent transition, turbulent mixing due to mixing of streams, and separated flows. Large eddy simulation (LES) can improve accuracy but at a considerably higher cost. In recent years, hybrid schemes that take advantage of both unsteady RANS and LES have been proposed. This study investigated an alternative scheme, the time-filtered Navier-Stokes (TFNS) method applied to compressible flows. The method developed by Shih and Liu was implemented in the Glenn-Heat-Transfer (Glenn-HT) code and applied to film-cooling flows. In this report the method and its implementation is briefly described. The film effectiveness results obtained for film cooling from a row of 30deg holes with a pitch of 3.0 diameters emitting air at a nominal density ratio of unity and two blowing ratios of 0.5 and 1.0 are shown. Flow features under those conditions are also described.

heat transfer↗

Calculation of free turbulent mixing by interaction approach.

The applicability of Bradshaw's interaction hypothesis to two-dimensional free shear flows was investigated. According to it, flows with velocity extrema may be considered to consist of several interacting layers. The hypothesis leads to a new expression for the shear stress which removes the usual restriction that shear stress vanishes at the velocity extremum. The approach is based on kinetic energy and the length scale equations. The compressible flow equations are simplified by restriction to low Mach numbers, and the range of their applicability is discussed. The empirical functions of the turbulence model are found here to be correlated with the spreading rate of the shear layer. The analysis demonstrates that the interaction hypothesis is a workable concept.

Morel, T.↗

Unstructured Mesh Methods for the Simulation of Hypersonic Flows

This report describes the research work undertaken at the Massachusetts Institute of Technology. The aim of this research is to identify effective algorithms and methodologies for the efficient and routine solution of hypersonic viscous flows about re-entry vehicles. For over ten years we have received support from NASA to develop unstructured mesh methods for Computational Fluid Dynamics. As a result of this effort a methodology based on the use, of unstructured adapted meshes of tetrahedra and finite volume flow solvers has been developed. A number of gridding algorithms flow solvers, and adaptive strategies have been proposed. The most successful algorithms developed from the basis of the unstructured mesh system FELISA. The FELISA system has been extensively for the analysis of transonic and hypersonic flows about complete vehicle configurations. The system is highly automatic and allows for the routine aerodynamic analysis of complex configurations starting from CAD data. The code has been parallelized and utilizes efficient solution algorithms. For hypersonic flows, a version of the, code which incorporates real gas effects, has been produced. One of the latest developments before the start of this grant was to extend the system to include viscous effects. This required the development of viscous generators, capable of generating the anisotropic grids required to represent boundary layers, and viscous flow solvers. In figures I and 2, we show some sample hypersonic viscous computations using the developed viscous generators and solvers. Although these initial results were encouraging, it became apparent that in order to develop a fully functional capability for viscous flows, several advances in gridding, solution accuracy, robustness and efficiency were required. As part of this research we have developed: 1) automatic meshing techniques and the corresponding computer codes have been delivered to NASA and implemented into the GridEx system, 2) a finite element algorithm for the solution of the viscous compressible flow equations which can solve flows all the way down to the incompressible limit and that can use higher order (quadratic) approximations leading to highly accurate answers, and 3) and iterative algebraic multigrid solution techniques.

Peraire, Jaime↗

Interferometric investigations of compressible dynamic stall over a transiently pitching airfoil

The compressible dynamic stall flowfield over a NACA 0012 airfoil transiently pitching from 0 to 60 deg at a constant rate under compressible flow conditions has been studied using real-time interferometry. A quantitative description of the overall flowfield, including the finer details of dynamic stall vortex formation, growth, and the concomitant changes in the airfoil pressure distribution, has been provided by analyzing the interferograms. For Mach numbers above 0.4, small multiple shocks appear near the leading edge and are present through the initial stages of dynamic stall. Dynamic stall was found to occur coincidentally with the bursting of the separation bubble over the airfoil. Compressibility was found to confine the dynamic stall vortical structure closer to the airfoil surface. The measurements show that the peak suction pressure coefficient drops with increasing freestream Mach number, and also it lags the steady flow values at any given angle of attack. As the dynamic stall vortex is shed, an anti-clockwise vortex is induced near the trailing edge, which actively interacts with the post-stall flow.

Chandrasekhara, M. S.↗

LES of Temporally Evolving Mixing Layers by an Eighth-Order Filter Scheme

An eighth-order filter method for a wide range of compressible flow speeds (H.C. Yee and B. Sjogreen, Proceedings of ICOSAHOM09, June 22-26, 2009, Trondheim, Norway) are employed for large eddy simulations (LES) of temporally evolving mixing layers (TML) for different convective Mach numbers (Mc) and Reynolds numbers. The high order filter method is designed for accurate and efficient simulations of shock-free compressible turbulence, turbulence with shocklets and turbulence with strong shocks with minimum tuning of scheme parameters. The value of Mc considered is for the TML range from the quasi-incompressible regime to the highly compressible supersonic regime. The three main characteristics of compressible TML (the self similarity property, compressibility effects and the presence of large-scale structure with shocklets for high Mc) are considered for the LES study. The LES results using the same scheme parameters for all studied cases agree well with experimental results of Barone et al. (2006), and published direct numerical simulations (DNS) work of Rogers & Moser (1994) and Pantano & Sarkar (2002).

Hadjadj, A↗

Instantaneous velocity field imaging instrument for supersonic reacting flows

The technical tasks conducted to develop and demonstrate a new gas velocity measurement technique for high enthalpy reacting flows is described. The technique is based on Doppler-shifted Planar Laser-induced Fluorescence (PLIF) imaging of the OH radical. The imaging approach permits, in principle, single-shot measurements of the 2-D distribution of a single velocity component in the measurement plane, and is thus a technique of choice for applications in high enthalpy transient flow facilities. In contrast to previous work in this area, the present program demonstrated an approach which modified the diagnostic technique to function under the constraints of practical flow conditions of engineering interest, rather than vice-versa. In order to accomplish the experimental demonstrations, the state-of-the-art in PLIF diagnostic techniques was advanced in several ways. Each of these tasks is described in detail and is intended to serve as a reference in supporting the transition of this new capability to the fielded PLIF instruments now installed at several national test facilities. Among the new results of general interest in LlF-based flow diagnostics, a detailed set of the first measurements of the collisional broadening and shifting behavior of OH (1,0) band transitions in H7-air combustion environments is included. Such measurements are critical in the design of a successful strategy for PLIF velocity imaging; they also relate to accurate concentration and temperature measurements, particularly in compressible flow regimes. Furthermore, the results shed new light on the fundamental relationship between broadening and energy transfer collisions in OH A(sup 2)Sigma(+)v(sup ') = 1. The first single-pulse, spectrally-resolved measurements of the output of common pulsed dye lasers were also produced during the course of this effort. As with the OH broadening measurements, these data are a significant aspect of a successful velocity imaging strategy, and also have potential implications for many other LIF measurement techniques. Our results indicated the need to modify the commercially available laser cavity in order to accommodate the constraints imposed by typical SCRAMJET combustion characteristics as well as to increase the instrument's velocity dynamic range to span an intra-image range in excess of 2 km/s. The various technical efforts were brought together in a series of experiments demonstrating the applicability of the technique in a high pressure, high temperature H2-air combustion system. The resultant images were compared with 2-D flow simulations in order to determine the accuracy of the instrument. Mean velocity imaging in flows with an axis of symmetry was demonstrated with an accuracy of +/- 50 m/s out of an intra-image dynamic range of 1600 m/s, including reversed flow. A more complex configuration amenable to single-shot imaging in flows without an axis of symmetry was also demonstrated. Limitations imposed by available equipment resulted in an accuracy of about +/- 200 m/s out of 1750 m/s in these demonstrations. Minor modifications to the present configuration were suggested to improve this performance. Each technical task is described in detail, along with significance of the results for the overall imaging velocimeter configuration. This report should allow the user community to integrate this new measurement capability in their existing instrumentation platforms.

Allen, M. G.↗

Simulation of Inviscid Compressible Multi-Phase Flow with Condensation

Condensation of vapours in rapid expansions of compressible gases is investigated. In the case of high temperature gradients the condensation will start at conditions well away from thermodynamic equilibrium of the fluid. In those cases homogeneous condensation is dominant over heterogeneous condensation. The present work is concerned with development of a simulation tool for computation of high speed compressible flows with homogeneous condensation. The resulting ow solver should preferably be accurate and robust to be used for simulation of industrial flows in general geometries.

Kelleners, Philip↗

Data-parallel lower-upper relaxation method for reacting flows

The implicit lower-upper symmetric Gauss-Seidel (LU-SGS) method of Yoon and Jameson is modified for use on massively parallel computers. The method has been implemented on the Thinking Machines CM-5 and the MasPar MP-1 and MP-2, where large percentages of the theoretical peak floating point performance are obtained. It is shown that the new data-parallel LU relaxation method has better convergence properties than the original method for two different inviscid compressible flow simulations. The convergence is also improved for five-species reacting air computations. The performance of the method on various partitions of the CM-5 and on the MasPar computers is discussed. The new method shows promise for the efficient simulation of very large perfect gas and reacting flows.

Candler, Graham V.↗

PeleLMeX: an AMR Low Mach Number Reactive Flow Simulation Code without level sub-cycling

PeleLMeX simulates chemically reacting low Mach number flows with block-structured adaptive mesh refinement (AMR). The code is built upon the AMReX library, which provides the underlying data structures and tools to manage and operate on them across massively parallel computing architectures. PeleLMeX algorithmic features are inherited from its predecessor PeleLM but key improvements allow representation of more complex physical processes. Together with its compressible flow counterpart PeleC, the thermo-chemistry library PelePhysics and the multi-physics library PeleMP, it forms the Pele suite of open-source reactive flow simulation codes.

97 MATHEMATICS AND COMPUTING↗