Search NASA⌕ Search

SEARCH · Search NASA

Results for “Compressible Euler equations”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Adaptive remeshing for compressible flow computations

The present, quality-enhancing adaptive-mesh procedure for two-dimensional Euler equation steady state solutions is implemented by means of linear triangular elements and an explicit time-stepping scheme, in conjunction with a finite element solution algorithm. The meshes thus generated typically take the form of stretched elements in the vicinity of one-dimensional flow features; a considerable variation in element size may thereby emerge which allows the desired high-quality solutions to be obtained with commensurately high computational efficiency.

Peraire, J.↗
NASA NTRS1986

Coupling between fluid dynamics and energy addition in arcjet and microwave thrusters

A new approach to numerically solving the problem of the constricted electric arcjet is presented. An Euler Implicit finite difference scheme is used to solve the full compressible Navier Stokes equations in two dimensions. The boundary and initial conditions represent the constrictor section of the arcjet, and hydrogen is used as a propellant. The arc is modeled as a Gaussian distribution across the centerline of the constrictor. Temperature, pressure and velocity profiles for steady state converged solutions show both axial and radial changes in distributions resulting from their interaction with the arc energy source for specific input conditions. The temperature rise is largest at the centerline where there is a the greatest concentration arc energy. The solution does not converge for all initial inputs and the limitations in the range of obtainable solutions are discussed.

Micci, M. M.↗
NASA NTRS1992

Unsteady blade pressures on a propfan: Predicted and measured compressibility effects

The effect of compressibility on unsteady blade pressures is studied by solving the three-dimensional Euler equations. The operation of the eight-bladed SR7L propfan at a 4.75 deg angle of attack was considered. Euler solutions were obtained for three Mach numbers, 0.6, 0.7 and 0.8, and the predicted blade pressure waveforms were compared with flight data. The comparisons show that in general, the effect of Mach number on pressure waveforms are correctly predicted. The change in pressure waveforms are minimal when the Mach number is increased from 0.6 to 0.7. Increasing the Mach number from 0.7 to 0.8 produces significant changes in predicted pressure levels. The predicted amplitudes, however, differ from measurements at some transducer locations. At all the three Mach numbers, the measured (installed propfan) pressure waveforms show a relative phase lag compared to the computed (propfan along) waveforms due to installation effects. Measured waveforms in the blade tip region show nonlinear variations which are not captured by the present numerical procedure.

Nallasamy, M.↗
DOE OSTI2020

Low-dissipation centred schemes for hyperbolic equations in conservative and non-conservative form

We propose a one-parameter family of low-dissipation centred numerical schemes for solving hyperbolic equations in conservative or non-conservative form, using finite volume or discontinuous Galerkin finite element methods. The new schemes spring out from the multi-dimensional FORCE method and are determined by a single parameter α ≥ 1. Given an increasing sequence of real numbers 1 ≤ α 1 < α 2 < . . . < α K , there corresponds a sequence of numerical schemes with stability restriction associated to a decreasing sequence of Courant numbers 1 > c 1 > c 2 > . . . > c K > 0 and a decreasing sequence of corresponding numerical viscosity functions d 1 > d 2 > . . . > d K . For a given Courant number c k ≤ 1 there is a real number α k > 0 and a corresponding stable scheme with minimal numerical viscosity d k . The proposed schemes suit very well the family of high-order discontinuous Galerkin finite element methods and the recently proposed class of ADER-TR schemes, whose orders of accuracy define decreasing sequences of Courant numbers, as the order of accuracy increases. The centred methods of this paper are stable in 2D and 3D in the frame of simultaneous updating formulae, unlike other centred methods, such as 1D FORCE, which are not. Furthermore, the schemes are highly accurate for slowly-moving waves, which is precisely the kind of waves that traditional centred methods smear disastrously. Additional features of the proposed schemes include ease of implementation and applicability to any hyperbolic system either in conservative or non-conservative form. Here, the proposed schemes are analysed and computationally assessed through a suite of test problems for a linear model system, for the Euler equations in one and two space dimensions, and for the Baer-Nunziato equations for compressible two-phase flow.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 1: Analysis description

A new computer code was developed to solve the two-dimensional or axisymmetric, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 1 is the Analysis Description, and describes in detail the governing equations, the turbulence model, the linearization of the equations and boundary conditions, the time and space differencing formulas, the ADI solution procedure, and the artificial viscosity models.

Towne, Charles E.↗
NASA NTRS1990

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 2: User's guide

A new computer code was developed to solve the two-dimensional or axisymmetric, Reynolds averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 2 is the User's Guide, and describes the program's general features, the input and output, the procedure for setting up initial conditions, the computer resource requirements, the diagnostic messages that may be generated, the job control language used to run the program, and several test cases.

Towne, Charles E.↗
NASA NTRS1990

PROTEUS two-dimensional Navier-Stokes computer code, version 1.0. Volume 3: Programmer's reference

A new computer code was developed to solve the 2-D or axisymmetric, Reynolds-averaged, unsteady compressible Navier-Stokes equations in strong conservation law form. The thin-layer or Euler equations may also be solved. Turbulence is modeled using an algebraic eddy viscosity model. The objective was to develop a code for aerospace applications that is easy to use and easy to modify. Code readability, modularity, and documentation were emphasized. The equations are written in nonorthogonal body-fitted coordinates, and solved by marching in time using a fully-coupled alternating-direction-implicit procedure with generalized first- or second-order time differencing. All terms are linearized using second-order Taylor series. The boundary conditions are treated implicitly, and may be steady, unsteady, or spatially periodic. Simple Cartesian or polar grids may be generated internally by the program. More complex geometries require an externally generated computational coordinate system. The documentation is divided into three volumes. Volume 3 is the Programmer's Reference, and describes the program structure, the FORTRAN variables stored in common blocks, and the details of each subprogram.

Towne, Charles E.↗
NASA NTRS1979

Calculated global conditions for the unsteady Euler's flow equations

This paper considers the numerical solution of a fluid dynamical system with time-dependent oscillating boundaries. The problem detailed, a compressible and inviscid flow in a one-dimensional pipe with vibrating ends, is examined to develop techniques needed to solve the more complex unsteady flows in turbomachine cascades. The problem is shown to be mathematically ill-posed but by complementing the governing equations with energy conservation, a unique solution is obtained. An artificial damping term is introduced into the governing equations to improve the performance of the time-marching integration.

Kuzo, D. M.↗
NASA NTRS1990

A review of reduced Navier-Stokes computations for compressible viscous flows

A reduced form of the Navier-Stokes equations, defined by a single composite of the Euler, boundary layer, and triple deck approximations, is considered for the computation of viscous interacting flows. Global pressure or pseudopotential relaxation methods are integrated with coupled sparse matrix direct solvers or coupled strongly implicit ILU algorithms to efficiently capture sharp shocks and regions of recirculation. Solutions are obtained for a variety of 2D and 3D geometries and for Mach numbers (M) spanning the range from incompressible (M = 0) to supersonic (M = 6) flow.

Rubin, S. G.↗
NASA NTRS1985

Numerical solution of Euler's equation by perturbed functionals

A perturbed functional iteration has been developed to solve nonlinear systems. It adds at each iteration level, unique perturbation parameters to nonlinear Gauss-Seidel iterates which enhances its convergence properties. As convergence is approached these parameters are damped out. Local linearization along the diagonal has been used to compute these parameters. The method requires no computation of Jacobian or factorization of matrices. Analysis of convergence depends on properties of certain contraction-type mappings, known as D-mappings. In this article, application of this method to solve an implicit finite difference approximation of Euler's equation is studied. Some representative results for the well known shock tube problem and compressible flows in a nozzle are given.

Dey, S. K.↗
NASA NTRS1985

An adaptive finite element procedure for compressible high speed flows

A practical finite element-based solution procedure for high-speed inviscid compressible flow problems is described. The method provides time-accurate solutions to the coompressible Euler equations, and is computationally more efficient than the one-step Taylor-Galerkin approach and better suited for implementation on the modern generation of vector computers. The method is coupled to an adaptive mesh refinement process that enables steady state solutions of improved quality to be obtained.

Loehner, R.↗
NASA NTRS1993

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.↗
NASA NTRS2000

A Finite-Element Approach for Modeling Inviscid and Viscous Compressible Flows using Prismatic Grids

The Galerkin finite-element method is used to solve the Euler and Navier-Stokes equations on prismatic meshes. It is shown that the prismatic grid is advantageous for correctly and efficiently capturing the boundary layers in high Reynolds number flows. It can be captured accurately because of the ability to cluster grid points normal to the body. The efficiency derives from the implicit treatment of the normal direction. To treat the normal direction implicitly, a semi-implicit Runge-Kutta time stepping scheme is developed. The semi-implicit algorithm is validated on simple geometries for inviscid and viscous flows and its convergence history is compared to that of the explicit Runge-Kutta scheme. The semi-implicit scheme is shown to be a factor of 3 to 4 faster in terms of CPU time to convergence.

Pandya, S. A.↗
NASA NTRS2003

A Study of the Impact of Variations on Aerodynamic Flow in Gas Turbine Engines via Monte-Carlo Simulations

The unsteady compressible inviscid flow is characterized by the conservations of mass, momentum, and energy; or simply the Euler equations. In this paper, a study of the subsonic one-dimensional Euler equations with local preconditioning is presented with a modal analysis approach. Specifically, this study investigates the behavior of airflow in a gas turbine engine using the specified conditions at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine, under the impact of variations in pressure, velocity, temperature, and density at low Mach numbers. Two main questions that motivate this research are: 1) Is there any aerodynamic problem with the existing gas turbine engines that could impact aircraft performance? 2) If yes, what aspect of a gas turbine engine could be improved via design to alleviate that impact and to optimize aircraft performance. This paper presents an initial attempt to the flow behavior in terms (perturbation) using simulation outputs from a customer-deck model obtained from Pratt&Whitney, (i.e., pressure, temperature, velocity, density) about their mean states at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine. Flow behavior is analyzed for the high pressure compressor and combustion chamber employing the conditions on their left and right boundaries. In the same fashion, similar analyses are carried out for the high and low-pressure turbines. In each case, the eigenfrequencies that are obtained for different boundary conditions are examined closely based on their probabilistic distributions, a result of a Monte Carlo 10,000-sample simulation. Furthermore, the characteristic waves and eave response are analyzed and contrasted among different cases, with and without preconditioners. The results reveal the existence of flow instabilities due to the combined effect of variations and excessive pressures; which are clearly the case in the combustion chamber and high-pressure turbine. Finally a discussion is presented on potential impacts of the instabilities and what can be improved via design to alleviate them for a better aircraft performance.

Ngo, Khiem Viet↗
NASA NTRS2004

Impact of Variations on 1-D Flow in Gas Turbine Engines via Monte Carlo Simulations

The unsteady compressible inviscid flow is characterized by the conservations of mass, momentum, and energy; or simply the Euler equations. In this paper, a study of the subsonic one-dimensional Euler equations with local preconditioning is presented using a modal analysis approach. Specifically, this study investigates the behavior of airflow in a gas turbine engine using the specified conditions at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine, to determine the impact of variations in pressure, velocity, temperature, and density at low Mach numbers. Two main questions motivate this research: 1) Is there any aerodynamic problem with the existing gas turbine engines that could impact aircraft performance? 2) If yes, what aspect of a gas turbine engine could be improved via design to alleviate that impact and to optimize aircraft performance? This paper presents an initial attempt to model the flow behavior in terms of their eigenfrequencies subject to the assumption of the uncertainty or variation (perturbation). The flow behavior is explored using simulation outputs from a customer-deck model obtained from Pratt & Whitney. Variations of the main variables (i.e., pressure, temperature, velocity, density) about their mean states at the inflow and outflow boundaries of the compressor, combustion chamber, and turbine are modeled. Flow behavior is analyzed for the high-pressure compressor and combustion chamber utilizing the conditions on their left and right boundaries. In the same fashion, similar analyses are carried out for the high-pressure and low-pressure turbines. In each case, the eigenfrequencies that are obtained for different boundary conditions are examined closely based on their probabilistic distributions, a result of a Monte Carlo 10,000 sample simulation. Furthermore, the characteristic waves and wave response are analyzed and contrasted among different cases, with and without preconditioners. The results reveal the existence of flow instabilities due to the combined effect of variations and excessive pressures in the case of the combustion chamber and high-pressure turbine. Finally, a discussion is presented on potential impacts of the instabilities and what can be improved via design to alleviate them for a better aircraft performance.

Ngo, Khiem Viet↗
NASA NTRS1992

A numerical resolution study of high order essentially non-oscillatory schemes applied to incompressible flow

High order essentially non-oscillatory (ENO) schemes, originally designed for compressible flow and in general for hyperbolic conservation laws, are applied to incompressible Euler and Navier-Stokes equations with periodic boundary conditions. The projection to divergence-free velocity fields is achieved by fourth order central differences through Fast Fourier Transforms (FFT) and a mild high-order filtering. The objective of this work is to assess the resolution of ENO schemes for large scale features of the flow when a coarse grid is used and small scale features of the flow, such as shears and roll-ups, are not fully resolved. It is found that high-order ENO schemes remain stable under such situations and quantities related to large-scale features, such as the total circulation around the roll-up region, are adequately resolved.

Weinan, E.↗
NASA NTRS1994

A numerical resolution study of high order essentially non-oscillatory schemes applied to incompressible flow

High order essentially non-oscillatory (ENO) schemes, originally designed for compressible flow and in general for hyperbolic conservation laws, are applied to incompressible Euler and Navier-Stokes equations with periodic boundary conditions. The projection to divergence-free velocity fields is achieved by fourth-order central differences through fast Fourier transforms (FFT) and a mild high-order filtering. The objective of this work is to assess the resolution of ENO schemes for large scale features of the flow when a coarse grid is used and small scale features of the flow, such as shears and roll-ups, are not fully resolved. It is found that high-order ENO schemes remain stable under such situations and quantities related to large scale features, such as the total circulation around the roll-up region, are adequately resolved.

Weinan, E.↗