SEARCH · Engineering Papers

Results for “RIEMANN PROBLEM”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Some notes on shock resolving flux functions. Part 1: Stationary characteristics

Numerical flux functions for solving the Euler equations using exact and/or approximate solutions of the Riemann problem of gasdynamics are discussed. Under certain restrictive conditions, schemes using these flux functions produce systems of equations which can exhibit a single degree of freedom. In some instances, the solutions represented by this degree of freedom are unstable to perturbations. This local instability can seriously degrade the temporal convergence of numerical schemes. This point is demonstrated by numerical example.

Barth, Timothy J.↗

NASA NTRS1990

Hyperbolic conservation laws and numerical methods

The mathematical structure of hyperbolic systems and the scalar equation case of conservation laws are discussed. Linear, nonlinear systems and the Riemann problem for the Euler equations are also studied. The numerical methods for conservation laws are presented in a nonstandard manner which leads to large time steps generalizations and computations on irregular grids. The solution of conservation laws with stiff source terms is examined.

Leveque, Randall J.↗

NASA NTRS1990

Numerical methods for systems of conservation laws of mixed type using flux splitting

The essentially non-oscillatory (ENO) finite difference scheme is applied to systems of conservation laws of mixed hyperbolic-elliptic type. A flux splitting, with the corresponding Jacobi matrices having real and positive/negative eigenvalues, is used. The hyperbolic ENO operator is applied separately. The scheme is numerically tested on the van der Waals equation in fluid dynamics. Convergence was observed with good resolution to weak solutions for various Riemann problems, which are then numerically checked to be admissible as the viscosity-capillarity limits. The interesting phenomena of the shrinking of elliptic regions if they are present in the initial conditions were also observed.

Shu, Chi-Wang↗

NASA NTRS1990

Structure and evolution of small-amplitude intermediate shock waves

A simplified set of equations is derived that approximates the magnetohydrodynamic (MHD) Navier-Stokes equations for weakly nonlinear disturbances whose speeds are close to the MHD intermediate speed. Its shock structure solutions are then examined. The fast and slow shock solutions are uniquely specified by their Rankine-Hugoniot relations. However, the intermediate shock solutions, which are not unique, are characterized by the integral through the shock of the noncoplanar component of the magnetic field. For situations in which this integral is conserved, the Riemann problem is well defined and predicts the evolution of intermediate shocks. This analysis is substantiated by numerical computations.

Kennel, C. F.↗

NASA NTRS1990

Non-oscillatory central differencing for hyperbolic conservation laws

Many of the recently developed high resolution schemes for hyperbolic conservation laws are based on upwind differencing. The building block for these schemes is the averaging of an appropriate Godunov solver; its time consuming part involves the field-by-field decomposition which is required in order to identify the direction of the wind. Instead, the use of the more robust Lax-Friedrichs (LxF) solver is proposed. The main advantage is simplicity: no Riemann problems are solved and hence field-by-field decompositions are avoided. The main disadvantage is the excessive numerical viscosity typical to the LxF solver. This is compensated for by using high-resolution MUSCL-type interpolants. Numerical experiments show that the quality of results obtained by such convenient central differencing is comparable with those of the upwind schemes.

Nessyahu, Haim↗

NASA NTRS1992

The role of intermediate shocks in magnetic reconnection

The present study examines the structure of discontinuity layers associated with magnetic reconnection by numerically solving the Riemann problem for the evolution of an initial current sheet which separates two plasma regions with antiparallel magnetic field components in the z direction and a common guide magnetic field in the y direction. In the presence of a nonzero normal component of the magnetic field, the initial current sheet evolves into a system of MHD discontinuities. For the initial current sheet with a zero guide field, steady intermediate shocks, slow shocks, slow expansion waves, or contact discontinuity are observed to develop. For the current sheet with a nonzero guide field, time-dependent intermediate shocks, instead of steady intermediate shocks, are observed to bound the reconnection layer.

Lin, Y.↗

NASA NTRS1992

A staggered mesh finite difference scheme for the computation of compressible flows

A simple high resolution finite difference technique is presented to approximate weak solutions to hyperbolic systems of conservation laws. The method does not rely on Riemann problem solvers and is therefore easy to extend to a wide variety of problems. The overall performance (resolution and CPU requirements) is competitive, with other state-of-the-art techniques offering sharp nonoscillatory shocks and contacts. Theoretical results confirm the reliability of the approach for linear systems and nonlinear scalar equations.

Sanders, Richard↗

NASA NTRS1992

A numerical method for systems of conservation laws of mixed type admitting hyperbolic flux splitting

The present treatment of elliptic regions via hyperbolic flux-splitting and high order methods proposes a flux splitting in which the corresponding Jacobians have real and positive/negative eigenvalues. While resembling the flux splitting used in hyperbolic systems, the present generalization of such splitting to elliptic regions allows the handling of mixed-type systems in a unified and heuristically stable fashion. The van der Waals fluid-dynamics equation is used. Convergence with good resolution to weak solutions for various Riemann problems are observed.

Shu, Chi-Wang↗

NASA NTRS1992

The Osher scheme for non-equilibrium reacting flows

An extension of the Osher upwind scheme to nonequilibrium reacting flows is presented. Owing to the presence of source terms, the Riemann problem is no longer self-similar and therefore its approximate solution becomes tedious. With simplicity in mind, a linearized approach which avoids an iterative solution is used to define the intermediate states and sonic points. The source terms are treated explicitly. Numerical computations are presented to demonstrate the feasibility, efficiency and accuracy of the proposed method. The test problems include a ZND (Zeldovich-Neumann-Doring) detonation problem for which spurious numerical solutions which propagate at mesh speed have been observed on coarse grids. With the present method, a change of limiter causes the solution to change from the physically correct CJ detonation solution to the spurious weak detonation solution.

Suresh, Ambady↗

NASA NTRS1992

Structural relations for time-dependent intermediate shocks

A quantitative description of time-dependent (MHD) intermediate shocks is formulated. In non-coplanar Riemann problems, time-dependent 2 yields 3 intermediate shocks evolve in time as a localized self-similar structure whose strength decreases as 1/sq rt t, and whose width expands as sq rt t. We derive structural relations, similar to Rankine-Hugoniot relations, between the plasma properties and the magnitude of the transverse magnetic field, which we hope can help identify time-dependent 2 yields 3 intermediate shocks at the magnetopause, in the solar wind, and elsewhere in space. The analytic forms of the structural relations agree with our numerical results.

Wu, C. C.↗

NASA NTRS1991

A staggered mesh finite difference scheme for the computation of hypersonic Euler flows

A shock capturing finite difference method for systems of hyperbolic conservation laws is presented which avoids the need to solve Riemann problems while being competitive in performance with other current methods. A staggered spatial mesh is employed, so that complicated nonlinear waves generated at cell interfaces are averaged over cell interiors at the next time level. The full method combines to form a conservative version of the modified method of characteristics. The advantages of the method are discussed, and numerical results are presented for the two-dimensional double ellipse problem.

Sanders, Richard↗

NASA NTRS1993

Ellipsoidal figures of equilibrium - Compressible models

The results of Chandrasekhar (1969) are generalized to polytropes, using a formalism based on ellipsoidal energy variational principle to construct approximate stellar equilibrium solutions and study their stability. After reviewing the energy variational method and describing the approach, several equivalent stability conditions are established and secular vs. dynamical instabilities are discussed. Then, the equilibrium structure equations are derived for isolated, rotating polytropes, and axisymmetric configurations (compressible Maclaurin spheroids) are considered. Particular attention is given to triaxial configurations, either in a state of uniform rotation (generalizing the classical Jacobi ellipsoids) or with internal fluid motions of uniform vorticity (the compressible analogues of Riemann-S ellipsoids) and to the stability of these single star configurations. The compressible generalizations of the Roche and Roche-Riemann problems for a polytrope in orbit about a point-mass companion are solved, and the generalized Darwin problem for two identical polytropes in a binary is considered.

Lai, Dong↗

NASA NTRS1994

A flux splitting scheme with high-resolution and robustness for discontinuities

A flux splitting scheme is proposed for the general nonequilibrium flow equations with an aim at removing numerical dissipation of Van-Leer-type flux-vector splittings on a contact discontinuity. The scheme obtained is also recognized as an improved Advection Upwind Splitting Method (AUSM) where a slight numerical overshoot immediately behind the shock is eliminated. The proposed scheme has favorable properties: high-resolution for contact discontinuities; conservation of enthalpy for steady flows; numerical efficiency; applicability to chemically reacting flows. In fact, for a single contact discontinuity, even if it is moving, this scheme gives the numerical flux of the exact solution of the Riemann problem. Various numerical experiments including that of a thermo-chemical nonequilibrium flow were performed, which indicate no oscillation and robustness of the scheme for shock/expansion waves. A cure for carbuncle phenomenon is discussed as well.

Wada, Yasuhiro↗

NASA NTRS1994

High-order centered difference methods with sharp shock resolution

In this paper we consider high-order centered finite difference approximations of hyperbolic conservation laws. We propose different ways of adding artificial viscosity to obtain sharp shock resolution. For the Riemann problem we give simple explicit formulas for obtaining stationary one and two-point shocks. This can be done for any order of accuracy. It is shown that the addition of artificial viscosity is equivalent to ensuring the Lax k-shock condition. We also show numerical experiments that verify the theoretical results.

Gustafsson, Bertil↗

NASA NTRS2003

A Stochastic Collocation Algorithm for Uncertainty Analysis

This report describes a stochastic collocation method to adequately handle a physically intrinsic uncertainty in the variables of a numerical simulation. For instance, while the standard Galerkin approach to Polynomial Chaos requires multi-dimensional summations over the stochastic basis functions, the stochastic collocation method enables to collapse those summations to a one-dimensional summation only. This report furnishes the essential algorithmic details of the new stochastic collocation method and provides as a numerical example the solution of the Riemann problem with the stochastic collocation method used for the discretization of the stochastic parameters.

Mathelin, Lionel↗

NASA NTRS2008

Magnetohydrodynamic Effects in Propagating Relativistic Ejecta: Reverse Shock and Magnetic Acceleration

We solve the Riemann problem for the deceleration of arbitrarily magnetized relativistic ejecta injected into a static unmagnetized medium. We find that for the same initial Lorentz factor, the reverse shock becomes progressively weaker with increasing magnetization s (the Poynting-to-kinetic energy flux ratio), and the shock becomes a rarefaction wave when s exceeds a critical value, sc, defined by the balance between the magnetic pressure in the ejecta and the thermal pressure in the forward shock. In the rarefaction wave regime, we find that the rarefied region is accelerated to a Lorentz factor that is significantly larger than the initial value. This acceleration mechanism is due to the strong magnetic pressure in the ejecta.

Mizuno, Y.↗

NASA NTRS2010

Evolution of Advection Upstream Splitting Method Schemes

This paper focuses on the evolution of advection upstream splitting method(AUSM) schemes. The main ingredients that have led to the development of modern computational fluid dynamics (CFD) methods have been reviewed, thus the ideas behind AUSM. First and foremost is the concept of upwinding. Second, the use of Riemann problem in constructing the numerical flux in the finite-volume setting. Third, the necessity of including all physical processes, as characterised by the linear (convection) and nonlinear (acoustic) fields. Fourth, the realisation of separating the flux into convection and pressure fluxes. The rest of this review briefly outlines the technical evolution of AUSM and more details can be found in the cited references. Keywords: Computational fluid dynamics methods, hyperbolic systems, advection upstream splitting method, conservation laws, upwinding, CFD

Liou, Meng-Sing↗

NASA NTRS1987

Comparative study of high-resolution shock-capturing schemes for a real gas

Recently developed second-order explicit shock-capturing methods, in conjunction with generalized flux-vector splittings, and a generalized approximate Riemann solver for a real gas are studied. The comparisons are made on different one-dimensional Riemann (shock-tube) problems for equilibrium air with various ranges of Mach numbers, densities and pressures. Six different Riemann problems are considered. These tests provide a check on the validity of the generalized formulas, since theoretical prediction of their properties appears to be difficult because of the non-analytical form of the state equation. The numerical results in the supersonic and low-hypersonic regimes indicate that these produce good shock-capturing capability and that the shock resolution is only slightly affected by the state equation of equilibrium air. The difference in shock resolution between the various methods varies slightly from one Riemann problem to the other, but the overall accuracy is very similar. For the one-dimensional case, the relative efficiency in terms of operation count for the different methods is within 30%. The main difference between the methods lies in their versatility in being extended to multidimensional problems with efficient implicit solution procedures.

Montagne, J.-L.↗