Search NASA⌕ Search

SEARCH · Search NASA

Results for “nonlinear conservation laws”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Relaxation and Preconditioning for High Order Discontinuous Galerkin Methods with Applications to Aeroacoustics and High Speed Flows

This project is about the investigation of the development of the discontinuous Galerkin finite element methods, for general geometry and triangulations, for solving convection dominated problems, with applications to aeroacoustics. Other related issues in high order WENO finite difference and finite volume methods have also been investigated. methods are two classes of high order, high resolution methods suitable for convection dominated simulations with possible discontinuous or sharp gradient solutions. In [18], we first review these two classes of methods, pointing out their similarities and differences in algorithm formulation, theoretical properties, implementation issues, applicability, and relative advantages. We then present some quantitative comparisons of the third order finite volume WENO methods and discontinuous Galerkin methods for a series of test problems to assess their relative merits in accuracy and CPU timing. In [3], we review the development of the Runge-Kutta discontinuous Galerkin (RKDG) methods for non-linear convection-dominated problems. These robust and accurate methods have made their way into the main stream of computational fluid dynamics and are quickly finding use in a wide variety of applications. They combine a special class of Runge-Kutta time discretizations, that allows the method to be non-linearly stable regardless of its accuracy, with a finite element space discretization by discontinuous approximations, that incorporates the ideas of numerical fluxes and slope limiters coined during the remarkable development of the high-resolution finite difference and finite volume schemes. The resulting RKDG methods are stable, high-order accurate, and highly parallelizable schemes that can easily handle complicated geometries and boundary conditions. We review the theoretical and algorithmic aspects of these methods and show several applications including nonlinear conservation laws, the compressible and incompressible Navier-Stokes equations, and Hamilton-Jacobi-like equations.

Shu, Chi-Wang↗
NASA NTRS2013

Entropy Stable Spectral Collocation Schemes for the Navier-Stokes Equations: Discontinuous Interfaces

Nonlinear entropy stability and a summation-by-parts framework are used to derive provably stable, polynomial-based spectral collocation methods of arbitrary order. The new methods are closely related to discontinuous Galerkin spectral collocation methods commonly known as DGFEM, but exhibit a more general entropy stability property. Although the new schemes are applicable to a broad class of linear and nonlinear conservation laws, emphasis herein is placed on the entropy stability of the compressible Navier-Stokes equations.

Carpenter, Mark H.↗
NASA NTRS1990

A kinetic equation with kinetic entropy functions for scalar conservation laws

A nonlinear kinetic equation is constructed and proved to be well-adapted to describe general multidimensional scalar conservation laws. In particular, it is proved to be well-posed uniformly in epsilon - the microscopic scale. It is also shown that the proposed kinetic equation is equipped with a family of kinetic entropy functions - analogous to Boltzmann's microscopic H-function, such that they recover Krushkov-type entropy inequality on the macroscopic scale. Finally, it is proved by both - BV compactness arguments in the one-dimensional case, that the local density of kinetic particles admits a continuum limit, as it converges strongly with epsilon below 0 to the unique entropy solution of the corresponding conservation law.

Perthame, Benoit↗
NASA NTRS2003

Finite Volume Methods: Foundation and Analysis

Finite volume methods are a class of discretization schemes that have proven highly successful in approximating the solution of a wide variety of conservation law systems. They are extensively used in fluid mechanics, porous media flow, meteorology, electromagnetics, models of biological processes, semi-conductor device simulation and many other engineering areas governed by conservative systems that can be written in integral control volume form. This article reviews elements of the foundation and analysis of modern finite volume methods. The primary advantages of these methods are numerical robustness through the obtention of discrete maximum (minimum) principles, applicability on very general unstructured meshes, and the intrinsic local conservation properties of the resulting schemes. Throughout this article, specific attention is given to scalar nonlinear hyperbolic conservation laws and the development of high order accurate schemes for discretizing them. A key tool in the design and analysis of finite volume schemes suitable for non-oscillatory discontinuity capturing is discrete maximum principle analysis. A number of building blocks used in the development of numerical schemes possessing local discrete maximum principles are reviewed in one and several space dimensions, e.g. monotone fluxes, E-fluxes, TVD discretization, non-oscillatory reconstruction, slope limiters, positive coefficient schemes, etc. When available, theoretical results concerning a priori and a posteriori error estimates are given. Further advanced topics are then considered such as high order time integration, discretization of diffusion terms and the extension to systems of nonlinear conservation laws.

Barth, Timothy↗
NASA NTRS1989

Second-order accurate nonoscillatory schemes for scalar conservation laws

Explicit finite difference schemes for the computation of weak solutions of nonlinear scalar conservation laws is presented and analyzed. These schemes are uniformly second-order accurate and nonoscillatory in the sense that the number of extrema of the discrete solution is not increasing in time.

Huynh, Hung T.↗
NASA NTRS1989

Second-order accurate nonoscillatory schemes for scalar conservation laws

Explicit finite difference schemes for the computation of weak solutions of nonlinear scalar conservation laws is presented and analyzed. These schemes are uniformly second-order accurate and nonoscillatory in the sense that the number of extrema of the discrete solution is not increasing in time.

Huynh, Hung T.↗
NASA NTRS1986

Semi-implicit and fully implicit shock-capturing methods for hyperbolic conservation laws with stiff source terms

Some numerical aspects of finite-difference algorithms for nonlinear multidimensional hyperbolic conservation laws with stiff nonhomogenous (source) terms are discussed. If the stiffness is entirely dominated by the source term, a semi-implicit shock-capturing method is proposed provided that the Jacobian of the soruce terms possesses certain properties. The proposed semi-implicit method can be viewed as a variant of the Bussing and Murman point-implicit scheme with a more appropriate numerical dissipation for the computation of strong shock waves. However, if the stiffness is not solely dominated by the source terms, a fully implicit method would be a better choice. The situation is complicated by problems that are higher than one dimension, and the presence of stiff source terms further complicates the solution procedures for alternating direction implicit (ADI) methods. Several alternatives are discussed. The primary motivation for constructing these schemes was to address thermally and chemically nonequilibrium flows in the hypersonic regime. Due to the unique structure of the eigenvalues and eigenvectors for fluid flows of this type, the computation can be simplified, thus providing a more efficient solution procedure than one might have anticipated.

Yee, H. C.↗
NASA NTRS1987

Semi-implicit and fully implicit shock-capturing methods for hyperbolic conservation laws with stiff source terms

Some numerical aspects of finite-difference algorithms for nonlinear multidimensional hyperbolic conservation laws with stiff nonhomogeneous (source) terms are discussed. If the stiffness is entirely dominated by the source term, a semi-implicit shock-capturing method is proposed provided that the Jacobian of the source terms possesses certain properties. The proposed semi-implicit method can be viewed as a variant of the Bussing and Murman point-implicit scheme with a more appropriate numerical dissipation for the computation of strong shock waves. However, if the stiffness is not solely dominated by the source terms, a fully implicit method would be a better choice. The situation is complicated by problems that are higher than one dimension, and the presence of stiff source terms further complicates the solution procedures for alternating direction implicit (ADI) methods. Several alternatives are discussed. The primary motivation for constructing these schemes was to address thermally and chemically nonequilibrium flows in the hypersonic regime. Due to the unique structure of the eigenvalues and eigenvectors for fluid flows of this type, the computation can be simplified, thus providing a more efficient solution procedure than one might have anticipated.

Yee, H. C.↗
NASA NTRS1989

A class of high resolution explicit and implicit shock-capturing methods

An attempt is made to give a unified and generalized formulation of a class of high resolution, explicit and implicit shock capturing methods, and to illustrate their versatility in various steady and unsteady complex shock wave computations. Included is a systematic review of the basic design principle of the various related numerical methods. Special emphasis is on the construction of the basis nonlinear, spatially second and third order schemes for nonlinear scalar hyperbolic conservation laws and the methods of extending these nonlinear scalar schemes to nonlinear systems via the approximate Riemann solvers and the flux vector splitting approaches. Generalization of these methods to efficiently include equilibrium real gases and large systems of nonequilibrium flows are discussed. Some issues concerning the applicability of these methods that were designed for homogeneous hyperbolic conservation laws to problems containing stiff source terms and shock waves are also included. The performance of some of these schemes is illustrated by numerical examples for 1-, 2- and 3-dimensional gas dynamics problems.

Yee, H. C.↗
NASA NTRS1989

Local error estimates for discontinuous solutions of nonlinear hyperbolic equations

Let u(x,t) be the possibly discontinuous entropy solution of a nonlinear scalar conservation law with smooth initial data. Suppose u sub epsilon(x,t) is the solution of an approximate viscosity regularization, where epsilon greater than 0 is the small viscosity amplitude. It is shown that by post-processing the small viscosity approximation u sub epsilon, pointwise values of u and its derivatives can be recovered with an error as close to epsilon as desired. The analysis relies on the adjoint problem of the forward error equation, which in this case amounts to a backward linear transport with discontinuous coefficients. The novelty of this approach is to use a (generalized) E-condition of the forward problem in order to deduce a W(exp 1,infinity) energy estimate for the discontinuous backward transport equation; this, in turn, leads one to an epsilon-uniform estimate on moments of the error u(sub epsilon) - u. This approach does not follow the characteristics and, therefore, applies mutatis mutandis to other approximate solutions such as E-difference schemes.

Tadmor, Eitan↗
NASA NTRS1995

Multiresolution Representation Using Biorthogonal Multiwavelets

We generalize Harten's multiresolution representation to biorthogonal multiwavelets. Several variants are considered. For example, a given array of discrete point values is transformed to point values and derivatives or point 'values and cell averages'. Compact Hermite interpolation is used in the decomposition and reconstruction algorithm. The resulting basis functions that are symmetric or skewsymmetric, compact, and smooth with optimal order accuracy. Harten's approach has several advantages: the multiresolution scheme is inherently discrete, non-periodic boundary conditions are easy to implement, and the representation can be extended to unstructured grids in bounded domains. We demonstrate the compression features of the new mutliwavelets by application to variable scale piecewise smooth functions with jump discontinuities typical of numerical solutions of nonlinear hyperbolic conservation laws.

Warming, Robert F.↗
NASA NTRS1985

Implicit Total Variation Diminishing (TVD) schemes for steady-state calculations

The novel implicit and unconditionally stable, high resolution Total Variation Diminishing (TVD) scheme whose application to steady state calculations is presently examined is a member of a one-parameter family of implicit, second-order accurate systems developed by Harten (1983) for the computation of weak solutions for one-dimensional hyperbolic conservation laws. The scheme will not generate spurious oscillations for a nonlinear scalar equation and a constant coefficient system. Numerical experiments for a quasi-one-dimensional nozzle problem show that the experimentally determined stability limit correlates exactly with the theoretical stability limit for the nonlinear scalar hyberbolic conservation laws.

Yee, H. C.↗
NASA NTRS1991

The P1-RKDG method for two-dimensional Euler equations of gas dynamics

A class of nonlinearly stable Runge-Kutta local projection discontinuous Galerkin (RKDG) finite element methods for conservation laws is investigated. Two dimensional Euler equations for gas dynamics are solved using P1 elements. The generalization of the local projections, which for scalar nonlinear conservation laws was designed to satisfy a local maximum principle, to systems of conservation laws such as the Euler equations of gas dynamics using local characteristic decompositions is discussed. Numerical examples include the standard regular shock reflection problem, the forward facing step problem, and the double Mach reflection problem. These preliminary numerical examples are chosen to show the capacity of the approach to obtain nonlinearly stable results comparable with the modern nonoscillatory finite difference methods.

Cockburn, Bernardo↗
NASA NTRS1985

Implicit TVD schemes for hyperbolic conservation laws in curvilinear coordinates

The Harten (1983, 1984) total variation-diminishing (TVD) schemes, constituting a one-parameter explicit and implicit, second-order-accurate family, have the property of not generating spurious oscillations when applied to one-dimensional, nonlinear scalar hyperbolic conservation laws and constant coefficient hyperbolic systems. These methods are presently extended to the multidimensional hyperbolic conservation laws in curvilinear coordinates. Means by which to linearize the implicit operator and solution strategies, in order to improve the computation efficiency of the implicit algorithm, are discussed. Numerical experiments with steady state airfoil calculations indicate that the proposed linearized implicit TVD schemes are accurate and robust.

Yee, H. C.↗
NASA NTRS1994

A Class of High-Resolution Explicit and Implicit Shock-Capturing Methods

The development of shock-capturing finite difference methods for hyperbolic conservation laws has been a rapidly growing area for the last decade. Many of the fundamental concepts, state-of-the-art developments and applications to fluid dynamics problems can only be found in meeting proceedings, scientific journals and internal reports. This paper attempts to give a unified and generalized formulation of a class of high-resolution, explicit and implicit shock capturing methods, and to illustrate their versatility in various steady and unsteady complex shock waves, perfect gases, equilibrium real gases and nonequilibrium flow computations. These numerical methods are formulated for the purpose of ease and efficient implementation into a practical computer code. The various constructions of high-resolution shock-capturing methods fall nicely into the present framework and a computer code can be implemented with the various methods as separate modules. Included is a systematic overview of the basic design principle of the various related numerical methods. Special emphasis will be on the construction of the basic nonlinear, spatially second and third-order schemes for nonlinear scalar hyperbolic conservation laws and the methods of extending these nonlinear scalar schemes to nonlinear systems via the approximate Riemann solvers and flux-vector splitting approaches. Generalization of these methods to efficiently include real gases and large systems of nonequilibrium flows will be discussed. Some perbolic conservation laws to problems containing stiff source terms and terms and shock waves are also included. The performance of some of these schemes is illustrated by numerical examples for one-, two- and three-dimensional gas-dynamics problems. The use of the Lax-Friedrichs numerical flux to obtain high-resolution shock-capturing schemes is generalized. This method can be extended to nonlinear systems of equations without the use of Riemann solvers or flux-vector splitting approaches and thus provides a large savings for multidimensional, equilibrium real gases and nonequilibrium flow computations.

Yee, H. C.↗
NASA NTRS1976

An implicit finite-difference algorithm for hyperbolic systems in conservation-law form

An implicit finite-difference scheme is developed for the efficient numerical solution of nonlinear hyperbolic systems in conservation-law form. The algorithm is second-order time-accurate, noniterative, and in a spatially factored form. Second- or fourth-order central and second-order one-sided spatial differencing are accommodated within the solution of a block tridiagonal system of equations. Significant conceptual and computational simplifications are made for systems whose flux vectors are homogeneous functions (of degree one), e.g., the Eulerian gasdynamic equations. Conservative hybrid schemes, which switch from central to one-sided spatial differencing whenever the local characteristic speeds are of the same sign, are constructed to improve the resolution of weak solutions. Numerical solutions are presented for a nonlinear scalar model equation and the two-dimensional Eulerian gasdynamic equations.

Beam, R. M.↗
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.↗