Search NASA⌕ Search

SEARCH · Search NASA

Results for “EULER EQUATION”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Finite Volume Discretization of the Euler Equations in Pronghorn

Modeling flow and heat transfer in high temperature gas reactors (HTGR) requires the ability to model a wide range of flow speeds from slow (natural convection), to intermediate (forced-flow conditions), to supersonic regimes (depressurization) for a wide range of geometries including the pebble bed, upper and lower plenum, and risers. In previous work, Pronghorn has effectively modeled low-to-medium speed flows in scenarios such as the one described in the two-dimensional PBMR-400 benchmark, using its finite-element-based streamline-upwind Petrov-Galerkin (SUPG) stabilized implementation of the Euler equations. However, limitations of this method become apparent when dealing with more complicated geometries (e.g. imposing slip boundary conditions at nodes belonging to two different boundaries) and when gas speeds are fast enough for shocks and supersonic flow to occur. For these problems, the finite-element-based solver lacks robustness and is plagued by slow iterative convergence or even divergence. In order to address these challenges, the Pronghorn code at INL has been updated with new, modified versions of its original equations. The new Pronghorn models are built on the finite volume method with a Harten-Lax-van Leer-Contact (HLLC) Riemann solver based numerical flux method, which (1) allows imposing slip boundary conditions much more robustly and (2) performs well for a wide range of flow speeds. The finite-volume-based flow solver will form the basis for a robust coarse-mesh thermal-hydraulics capability in Pronghorn.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1994

Canonical-variables multigrid method for steady-state Euler equations

A novel approach for the solution of inviscid flow problems for subsonic compressible flows is described. The approach is based on canonical forms of the equations, in which subsystems governed by hyperbolic operators are separated from those governed by elliptic ones. The discretizations used as well as the iterative techniques for the different subsystems are inherently different. Hyperbolic parts, which describe, in general, propagation phenomena, are discretized using upwind schemes and are solved by marching techniques. Elliptic parts, which are directionally unbiased, are discretized using h-elliptic central discretizations, and are solved by pointwise relaxations together with coarse grid acceleration. The resulting discretization schemes introduce artificial viscosity only for the hyperbolic parts of the system; thus a smaller total artificial viscosity is used, while the multigrid solvers used are much more efficient. Solutions of the subsonic compressible Euler equations are achieved at the same efficiency as the full potential equation.

Taasan, Shlomo↗
NASA NTRS1985

A simulation of rotor-stator interaction using the Euler equations and patched grids

An unsteady Euler code to study rotor-stator interaction problem was developed. The code uses patched grids that move relative to each other to simulate the motion of the rotor airfoils with respect to the stator airfoils. The Osher integration scheme is used in conjunction with an implicit relaxation approach. The scheme is second order accurate in space and time, and is also TVD in each spatial direction. The numerical results were found to be periodic in time, thus demonstrating the capability of the integration and zonal schemes in simulating periodic time dependent flow. The pressure contours obtained are almost oscillation free because of the TVD nature of the scheme. A new procedure was developed to simulate flows about bodies that move relative to each other. This capability should prove to be very useful in the areas of rotor-stator interaction, propeller-nacelle interaction, and helicopter rotor-fuselage interaction.

Rai, M. M.↗
NASA NTRS1998

A Preconditioning Method for Shape Optimization Governed by the Euler Equations

We consider a classical aerodynamic shape optimization problem subject to the compressible Euler flow equations. The gradient of the cost functional with respect to the shape variables is derived with the adjoint method at the continuous level. The Hessian (second order derivative of the cost functional with respect to the shape variables) is approximated also at the continuous level, as first introduced by Arian and Ta'asan (1996). The approximation of the Hessian is used to approximate the Newton step which is essential to accelerate the numerical solution of the optimization problem. The design space is discretized in the maximum dimension, i.e., the location of each point on the intersection of the computational mesh with the airfoil is taken to be an independent design variable. We give numerical examples for 86 design variables in two different flow speeds and achieve an order of magnitude reduction in the cost functional at a computational effort of a full solution of the analysis partial differential equation (PDE).

Arian, Eyal↗
NASA NTRS1985

Spectral methods for the Euler equations. II - Chebyshev methods and shock fitting

The Chebyshev spectral collocation method for the Euler gasdynamic equations is described. It is used with shock fitting to compute several two-dimensional gasdynamic flows. Examples include a shock/acoustic wave interaction, a shock/vortex interaction, and the classical blunt-body problem. With shock fitting, the spectral method has a clear advantage over second-order finite differences in that equivalent accuracy can be obtained with far fewer grid points.

Hussaini, M. Y.↗
NASA NTRS1984

Spectral methods for the Euler equations: Chebyshev methods and shock-fitting

The Chebyshev spectral collocation method for the Euler gas-dynamic equations is described. It is used with shock fitting to compute several two-dimensional, gas-dynamic flows. Examples include a shock-acoustic wave interaction, a shock/vortex interaction, and the classical blunt body problem. With shock fitting, the spectral method has a clear advantage over second order finite differences in that equivalent accuracy can be obtained with far fewer grid points.

Hussaini, M. Y.↗
NASA NTRS1985

Convergence acceleration of iterative solutions of Euler equations for transonic flow computations

In this paper, two acceleration techniques for Euler calculations are investigated. The first technique is an extrapolation procedure based on the Power Method; it is applicable when the iterative matrix has dominant eigenvalues. Both real and complex conjugate roots are allowed. The second technique is a generalization of the Minimal Residual Method, where the extrapolation step consists of a weighted combination of the corrections at different iteration levels and the weights are chosen to minimize the Lz norm of the residual. Numerical results, using Jameson's Runge-Kutta Multigrid Code, are presented. The extra computational work to apply either technique is negligible and the extra storage is not a problem on current supercomputers.

Hafez, M.↗
NASA NTRS1993

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

The block adaptive multigrid method applied to the solution of the Euler equations

In the present study, a scheme capable of solving very fast and robust complex nonlinear systems of equations is presented. The Block Adaptive Multigrid (BAM) solution method offers multigrid acceleration and adaptive grid refinement based on the prediction of the solution error. The proposed solution method was used with an implicit upwind Euler solver for the solution of complex transonic flows around airfoils. Very fast results were obtained (18-fold acceleration of the solution) using one fourth of the volumes of a global grid with the same solution accuracy for two test cases.

Pantelelis, Nikos↗
NASA NTRS2003

Robust and Simple Non-Reflecting Boundary Conditions for the Euler Equations: A New Approach Based on the Space-Time CE/SE Method

This paper reports on a significant advance in the area of non-reflecting boundary conditions (NRBCs) for unsteady flow computations. As a part of the development of the space-time conservation element and solution element (CE/SE) method, sets of NRBCs for 1D Euler problems are developed without using any characteristics-based techniques. These conditions are much simpler than those commonly reported in the literature, yet so robust that they are applicable to subsonic, transonic and supersonic flows even in the presence of discontinuities. In addition, the straightforward multidimensional extensions of the present 1D NRBCs have been shown numerically to be equally simple and robust. The paper details the theoretical underpinning of these NRBCs, and explains their unique robustness and accuracy in terms of the conservation of space-time fluxes. Some numerical results for an extended Sod's shock-tube problem, illustrating the effectiveness of the present NRBCs are included, together with an associated simple Fortran computer program. As a preliminary to the present development, a review of the basic CE/SE schemes is also included.

Chang, Sin-Chung↗
NASA NTRS2005

Robust and Simple Non-Reflecting Boundary Conditions for the Euler Equations - A New Approach based on the Space-Time CE/SE Method

This paper reports on a significant advance in the area of nonreflecting boundary conditions (NRBCs) for unsteady flow computations. As a part of t he development of t he space-time conservation element and solution element (CE/SE) method, sets of NRBCs for 1D Euler problems are developed without using any characteristics- based techniques. These conditions are much simpler than those commonly reported in the literature, yet so robust that they are applicable to subsonic, transonic and supersonic flows even in the presence of discontinuities. In addition, the straightforward multidimensional extensions of the present 1D NRBCs have been shown numerically to be equally simple and robust. The paper details the theoretical underpinning of these NRBCs, and explains t heir unique robustness and accuracy in terms of t he conservation of space-time fluxes. Some numerical results for an extended Sod's shock-tube problem, illustrating the effectiveness of the present NRBCs are included, together with an associated simple Fortran computer program. As a preliminary to the present development, a review of the basic CE/SE schemes is also included.

Chang, S.-C.↗
NASA NTRS1976

Transonic solutions of the Euler equations by the finite volume method

An investigation is conducted of the time-accurate convergence of representative transonic flows to a steady state under given constraints of time-accuracy. Factored explicit and implicit difference operators are used to accelerate the calculations. Attention is given to flow at a Mach number of 1.35 past a circular cylinder and supersonic flow past a NACA 0012 airfoil for three different Mach numbers. Questions of transonic wave behavior are considered along with the equations of motion and the characteristics of the mesh network.

Rizzi, A.↗
NASA NTRS1979

A two-dimensional unsteady Euler-equation solver for flow regions with arbitrary boundaries

A new technique is described for solving supersonic fluid dynamic problems containing multiple regions of continuous flow, each bounded by a permeable or impermeable surface. Region boundaries are, in general, arbitrarily shaped and time dependent. Discretization of such a region for solution by conventional finite difference procedures is accomplished using an elliptic solver which alleviates the dependence on a particular base coordinate system. Multiple regions are coupled together through the boundary conditions. The technique has been applied to a variety of problems including the shock diffraction and pointed wedge with detached bow shock.

Hindman, R. G.↗