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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Fast Solution of Fully Implicit Runge--Kutta and Discontinuous Galerkin in Time for Numerical PDEs, Part II: Nonlinearities and DAEs

Fully implicit Runge--Kutta (IRK) methods have many desirable accuracy and stability properties as time integration schemes, but high-order IRK methods are not commonly used in practice with large-scale numerical PDEs because of the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic framework for solving the nonlinear equations that arise from IRK methods (and discontinuous Galerkin discretizations in time) applied to nonlinear numerical PDEs, including PDEs with algebraic constraints. Several new linearizations of the nonlinear IRK equations are developed, offering faster and more robust convergence than the often-considered simplified Newton, as well as an effective preconditioner for the true Jacobian if exact Newton iterations are desired. Inverting these linearizations requires solving a set of block 2 x 2 systems. Under quite general assumptions, it is proven that the preconditioned 2 x 2 operator's condition number is bounded by a small constant close to one, independent of the spatial discretization, spatial mesh, and time step, and with only weak dependence on the number of stages or integration accuracy. Moreover, the new method is built using the same preconditioners needed for backward Euler-type time stepping schemes, so can be readily added to existing codes. The new methods are applied to several challenging fluid flow problems, including the compressible Euler and Navier--Stokes equations, and the vorticity-streamfunction formulation of the incompressible Euler and Navier--Stokes equations. Up to 10th-order accuracy is demonstrated using Gauss IRK, while in all cases fourth-order Gauss IRK requires roughly half the number of preconditioner applications as required by standard Singly diagonally implicit Runge--Kutta methods.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1983

A finite difference solution for the propagation of sound in near sonic flows

An explicit time/space finite difference procedure is used to model the propagation of sound in a quasi one-dimensional duct containing high Mach number subsonic flow. Nonlinear acoustic equations are derived by perturbing the time-dependent Euler equations about a steady, compressible mean flow. The governing difference relations are based on a fourth-order, two-step (predictor-corrector) MacCormack scheme. The solution algorithm functions by switching on a time harmonic source and allowing the difference equations to iterate to a steady state. The principal effect of the non-linearities was to shift acoustical energy to higher harmonics. With increased source strengths, wave steepening was observed. This phenomenon suggests that the acoustical response may approach a shock behavior at at higher sound pressure level as the throat Mach number aproaches unity. On a peak level basis, good agreement between the nonlinear finite difference and linear finite element solutions was observed, even through a peak sound pressure level of about 150 dB occurred in the throat region. Nonlinear steady state waveform solutions are shown to be in excellent agreement with a nonlinear asymptotic theory.

Hariharan, S. I.↗
NASA NTRS1984

A finite difference solution for the propagation of sound in near sonic flows

An explicit time/space finite difference procedure is used to model the propagation of sound in a quasi one-dimensional duct containing high Mach number subsonic flow. Nonlinear acoustic equations are derived by perturbing the time-dependent Euler equations about a steady, compressible mean flow. The governing difference relations are based on a fourth-order, two-step (predictor-corrector) MacCormack scheme. The solution algorithm functions by switching on a time harmonic source and allowing the difference equations to iterate to a steady state. The principal effect of the non-linearities was to shift aocustical energy to higher harmonics. With increased source strength, wave steepening was observed. This phenomenon suggests that the acoustical response may approach a shock behavior at higher sound pressure level as the throat Mach number approaches unity. On a peak level basis, good agreement between the nonlinear finite difference and linear finite element solutions was observed, even through a peak sound pressure level of about 150 dB occurred in the throat region. Nonlinear steady state waveform solutions are shown to be in excellent agreement with a nonlinear asymptotic theory. Previously announced in STAR as N83-30167

Hariharan, S. I.↗
NASA NTRS1990

Three-dimensional Euler time accurate simulations of fan rotor-stator interactions

A numerical method useful to describe unsteady 3-D flow fields within turbomachinery stages is presented. The method solves the compressible, time dependent, Euler conservation equations with a finite volume, flux splitting, total variation diminishing, approximately factored, implicit scheme. Multiblock composite gridding is used to partition the flow field into a specified arrangement of blocks with static and dynamic interfaces. The code is optimized to take full advantage of the processing power and speed of the Cray Y/MP supercomputer. The method is applied to the computation of the flow field within a single stage, axial flow fan, thus reproducing the unsteady 3-D rotor-stator interaction.

Boretti, A. A.↗
NASA NTRS2003

A Reduced-Order Model for Efficient Simulation of Synthetic Jet Actuators

A new reduced-order model of multidimensional synthetic jet actuators that combines the accuracy and conservation properties of full numerical simulation methods with the efficiency of simplified zero-order models is proposed. The multidimensional actuator is simulated by solving the time-dependent compressible quasi-1-D Euler equations, while the diaphragm is modeled as a moving boundary. The governing equations are approximated with a fourth-order finite difference scheme on a moving mesh such that one of the mesh boundaries coincides with the diaphragm. The reduced-order model of the actuator has several advantages. In contrast to the 3-D models, this approach provides conservation of mass, momentum, and energy. Furthermore, the new method is computationally much more efficient than the multidimensional Navier-Stokes simulation of the actuator cavity flow, while providing practically the same accuracy in the exterior flowfield. The most distinctive feature of the present model is its ability to predict the resonance characteristics of synthetic jet actuators; this is not practical when using the 3-D models because of the computational cost involved. Numerical results demonstrating the accuracy of the new reduced-order model and its limitations are presented.

Yamaleev, Nail K.↗
DOE OSTI2024

Georgia Tech Accelerated, Compressed, and Regularized Compute of Kinetic-based PDEs (Final Report)

This report summarizes the collaborative effort between Lawrence Livermore National Laboratory and Georgia Tech to enhance the BoBa library for tensor train computation in PDE solvers, with a target on kinetic equations and their continuum limits. We aimed to reduce computational cost and memory usage by replacing traditional array-based computations with tensor trains. We examined the compressibility of time-evolving solutions to the Euler equations with discontinuities. We also explored using the first invsicid and linear regularization of the compressible flow equations via the information geometric regularization (IGR). We explored this in a tensor train formulation. To identify that inverse terms in the IGR equations pose problems for tensor train formulations and investigate efficient methods for batched inversion of tensor trains.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1992

Practical design optimization of wing/body configurations using the Euler equations

The development of a practical method for the aerodynamic design of isolated wing and wing/body configurations is achieved through the coupling of existing computational fluid dynamics (CFD) analysis codes and a quasi-Newton numerical optimization method. The direct design method is generalized to treat three-dimensional aerodynamic optimization problems subject to inviscid, rotational, compressible flow conditions imposed by the Euler equations. The method couples either the FLO57 or the TEAM flow solver with a modified version of the QNMDIF numerical optimization algorithm. The method is applied, but is not limited, to supersonic design problems. A case study is presented illustrating the method's effectiveness in maximizing the lift-to-drag ratio, subject to a variety of constraints, of selected supersonic configurations at cruise conditions.

Reuther, J.↗
NASA NTRS2000

Progress in the Development of a Class of Efficient Low Dissipative High Order Shock-capturing Methods

In a series of papers, Olsson (1994, 1995), Olsson & Oliger (1994), Strand (1994), Gerritsen Olsson (1996), Yee et al. (1999a,b, 2000) and Sandham & Yee (2000), the issue of nonlinear stability of the compressible Euler and Navier-Stokes Equations, including physical boundaries, and the corresponding development of the discrete analogue of nonlinear stable high order schemes, including boundary schemes, were developed, extended and evaluated for various fluid flows. High order here refers to spatial schemes that are essentially fourth-order or higher away from shock and shear regions. The objective of this paper is to give an overview of the progress of the low dissipative high order shock-capturing schemes proposed by Yee et al. (1999a,b, 2000). This class of schemes consists of simple non-dissipative high order compact or non-compact central spatial differencings and adaptive nonlinear numerical dissipation operators to minimize the use of numerical dissipation. The amount of numerical dissipation is further minimized by applying the scheme to the entropy splitting form of the inviscid flux derivatives, and by rewriting the viscous terms to minimize odd-even decoupling before the application of the central scheme (Sandham & Yee). The efficiency and accuracy of these scheme are compared with spectral, TVD and fifth- order WENO schemes. A new approach of Sjogreen & Yee (2000) utilizing non-orthogonal multi-resolution wavelet basis functions as sensors to dynamically determine the appropriate amount of numerical dissipation to be added to the non-dissipative high order spatial scheme at each grid point will be discussed. Numerical experiments of long time integration of smooth flows, shock-turbulence interactions, direct numerical simulations of a 3-D compressible turbulent plane channel flow, and various mixing layer problems indicate that these schemes are especially suitable for practical complex problems in nonlinear aeroacoustics, rotorcraft dynamics, direct numerical simulation or large eddy simulation of compressible turbulent flows at various speeds including high-speed shock-turbulence interactions, and general long time wave propagation problems. These schemes, including entropy splitting, have also been extended to freestream preserving schemes on curvilinear moving grids for a thermally perfect gas (Vinokur & Yee 2000).

Yee, H. C.↗
NASA NTRS1976

Diffraction of a shock wave by a compression corner; regular and single Mach reflection

The two dimensional, time dependent Euler equations which govern the flow field resulting from the injection of a planar shock with a compression corner are solved with initial conditions that result in either regular reflection or single Mach reflection of the incident planar shock. The Euler equations which are hyperbolic are transformed to include the self similarity of the problem. A normalization procedure is employed to align the reflected shock and the Mach stem as computational boundaries to implement the shock fitting procedure. A special floating fitting scheme is developed in conjunction with the method of characteristics to fit the slip surface. The reflected shock, the Mach stem, and the slip surface are all treated as harp discontinuities, thus, resulting in a more accurate description of the inviscid flow field. The resulting numerical solutions are compared with available experimental data and existing first-order, shock-capturing numerical solutions.

Vijayashankar, V. S.↗
NASA NTRS1977

Diffraction of a shock wave by a compression corner. II - Single Mach reflection

The two-dimensional time-dependent Euler equations which govern the flow field resulting from the interaction of a planar shock with a compression corner are solved for initial conditions which result in single Mach reflection of the incident planar shock. The Euler equations are first transformed to include the self-similarity of the flow field. A second transformation is employed to normalize the distances between the ramp and the reflected shock and between the wall and the Mach stem. The resulting equations in strong conservation-law form are solved using a second-order discontinuity-fitting finite-difference approach. The results are compared with experimental interferograms and existing first-order shock-capturing numerical solutions.

Shankar, V.↗
NASA NTRS1979

The lambda-scheme

A method for integrating the Euler equations of gas dynamics for compressible flows in any hyperbolic case is presented. This method is applied to the Mach number distribution over a stretch of an infinite duct having a variable cross section, and to the distribution in a channel opening into a vacuum with the Mach number equalling 1.04. An example of the ability of this method to handle two-dimensional unsteady flows is shown using the steady shock-and-isobars pattern reached asymptotically about an ablated blunt body with a free stream Mach number equalling 12. A final example is presented where the technique is applied to a three-dimensional steady supersonic flow, with a Mach number of 2 and an angle of attack of 5 deg.

Moretti, G.↗
NASA NTRS1982

A natural conservative flux difference splitting for the hyperbolic systems of gasdynamics

Gasdynamic hyperbolic systems are treated by a novel conservative flux difference splitting upwind method which is applicable to explicit or implicit and iterative or direct schemes, for marching in time or space with the Euler or compressible Navier-Stokes equations. The method is able to capture sharp shocks correctly by maintaining global conservation. The embedded characteristics formulation is cast in the volumetric internal energy primitive variable, which is effective for real as well as perfect gases. Higher order upwind formulae are constructed from the distribution of pieces of the simple internal flux difference splitting, and the method reverts to first order at the appropriate points in order not to violate the domain of dependence by differencing across discontinuities. At insignificant computational overhead, switching is done by simple algebraic extensions of the truth functions that determine upwind direction for the first order scheme.

Lombard, C. K.↗
NASA NTRS1985

On symmetric and upwind TVD schemes

A class of explicit and implicit total variation diminishing (TVD) schemes for the compressible Euler and Navier-Stokes equations was developed. They do not generate spurious oscillations across shocks and contact discontinuities. In general, shocks can be captured within 1 to 2 grid points. For the inviscid case, these schemes are divided into upwind TVD schemes and symmetric (nonupwind) TVD schemes. The upwind TVD scheme is based on the second-order TVD scheme. The symmetric TVD scheme is a generalization of Roe's and Davis' TVD Lax-Wendroff scheme. The performance of these schemes on some viscous and inviscid airfoil steady-state calculations is investigated. The symmetric and upwind TVD schemes are compared.

Yee, H. C.↗
NASA NTRS1985

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

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

Assessment of a 3-D boundary layer analysis to predict heat transfer and flow field in a turbine passage

An assessment was made of the applicability of a three dimensional boundary layer analysis of heat transfer, total pressure losses, and streamline flow patterns on the surfaces of both stationary and rotating turbine passages. In support of this effort, an analysis was 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. 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 three dimensional 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.↗
NASA NTRS1988

The use of direct solvers for compressible flow computations

The feasibility of using a vectorized direct solver for the compressible Euler and Navier-Stokes equations is examined for both single-grid and multi-grid strategies. A procedure is developed for comparing the computational effort required for the direct method with that of the vertical line Gauss-Seidel iteration scheme in order to provide a criteria for choosing between the two techniques. The direct method is shown to have a relatively wide range of application on a vector processor with large memory. Results for both inviscid and viscous test problems in two-dimensions are examined.

Riggins, David↗
NASA NTRS1988

Finite element Euler computations in three dimensions

A two-step explicit FEM solution algorithm for the three-dimensional compressible Euler and Navier-Stokes equations based on unstructured triangular and tetrahedral grids is described and demonstrated. The method represents an extension and refinement of the algorithms presented by Loehner et al. (1984 and 1985), Peraire et al. (1987), and Morgan et al. (1987). The formulation and numerical implementation are outlined; the mesh generation, data structures, and adaptive remeshing are explained; and results for a two-dimensional airfoil, a three-dimensional engine air intake, a B747 in landing configuration, and a generic fighter aircraft are presented in extensive graphics and discussed in detail.

Peraire, Jaime↗