SEARCH · Engineering Papers

Results for “conservative scheme”

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 163 records · Page 9

All sources NASA NTRS DOE OSTI

NASA NTRS1991

A microcomputer program for estimating low altitude wind and turbulence fields

Past efforts to develop methods for objective wind analysis and to provide turbulence estimates to pilots are reviewed. The present approach involves a wind module and a turbulence module. The wind module extends the critical dividing streamline concept and past developments in mass-conserving wind interpolation schemes. The turbulence module is based on recent efforts to develop practical atmospheric turbulence parameterization schemes based on lapse rate, wind shear, surface heating, and surface roughness. The lapse rate and wind shear are readily obtained from the wind module. Surface roughness and heating come from land use information and net radiation estimates derived from cloud cover, terrain slope, local time, and latitude. This system is argued to have considerable potential for providing useful online information for many kinds of aircraft operations.

Ludwig, F. L.↗

NASA NTRS1997

Numerical Schemes for the Hamilton-Jacobi and Level Set Equations on Triangulated Domains

Borrowing from techniques developed for conservation law equations, numerical schemes which discretize the Hamilton-Jacobi (H-J), level set, and Eikonal equations on triangulated domains are presented. The first scheme is a provably monotone discretization for certain forms of the H-J equations. Unfortunately, the basic scheme lacks proper Lipschitz continuity of the numerical Hamiltonian. By employing a virtual edge flipping technique, Lipschitz continuity of the numerical flux is restored on acute triangulations. Next, schemes are introduced and developed based on the weaker concept of positive coefficient approximations for homogeneous Hamiltonians. These schemes possess a discrete maximum principle on arbitrary triangulations and naturally exhibit proper Lipschitz continuity of the numerical Hamiltonian. Finally, a class of Petrov-Galerkin approximations are considered. These schemes are stabilized via a least-squares bilinear form. The Petrov-Galerkin schemes do not possess a discrete maximum principle but generalize to high order accuracy.

Barth, Timothy J.↗

NASA NTRS1986

On symmetric and upwind TVD schemes

The performance of the upwind and symmetric total variation diminishing (TVD) schemes in viscous and inviscid airfoil steady-state calculations is considered, and the extension of the implicit second-order-accurate TVD scheme for hyperbolic systems of conservative laws in curvilinear coordinates is discussed. For two-dimensional steady-state applications, schemes are implemented in a conservative noniterative alternating direction implicit form, and results illustrate that the algorithm produces a fairly good solution for an RAE2822 airfoil calculation. The study demonstrates that the symmetric TVD scheme is as accurate as the upwind TVD scheme, while requiring less computational effort than it.

Yee, H. C.↗

NASA NTRS1992

Progress in multi-dimensional upwind differencing

Multi-dimensional upwind-differencing schemes for the Euler equations are reviewed. On the basis of the first-order upwind scheme for a one-dimensional convection equation, the two approaches to upwind differencing are discussed: the fluctuation approach and the finite-volume approach. The usual extension of the finite-volume method to the multi-dimensional Euler equations is not entirely satisfactory, because the direction of wave propagation is always assumed to be normal to the cell faces. This leads to smearing of shock and shear waves when these are not grid-aligned. Multi-directional methods, in which upwind-biased fluxes are computed in a frame aligned with a dominant wave, overcome this problem, but at the expense of robustness. The same is true for the schemes incorporating a multi-dimensional wave model not based on multi-dimensional data but on an 'educated guess' of what they could be. The fluctuation approach offers the best possibilities for the development of genuinely multi-dimensional upwind schemes. Three building blocks are needed for such schemes: a wave model, a way to achieve conservation, and a compact convection scheme. Recent advances in each of these components are discussed; putting them all together is the present focus of a worldwide research effort. Some numerical results are presented, illustrating the potential of the new multi-dimensional schemes.

Vanleer, Bram↗

NASA NTRS1992

Triangle based adaptive stencils for the solution of hyperbolic conservation laws

A triangle based total variation diminishing (TVD) scheme for the numerical approximation of hyperbolic conservation laws in two space dimensions is constructed. The novelty of the scheme lies in the nature of the preprocessing of the cell averaged data, which is accomplished via a nearest neighbor linear interpolation followed by a slope limiting procedures. Two such limiting procedures are suggested. The resulting method is considerably more simple than other triangle based non-oscillatory approximations which, like this scheme, approximate the flux up to second order accuracy. Numerical results for linear advection and Burgers' equation are presented.

Durlofsky, Louis J.↗

NASA NTRS2007

The alpha(3) Scheme - A Fourth-Order Neutrally Stable CESE Solver

The conservation element and solution element (CESE) development is driven by a belief that a solver should (i) enforce conservation laws in both space and time, and (ii) be built from a non-dissipative (i.e., neutrally stable) core scheme so that the numerical dissipation can be controlled effectively. To provide a solid foundation for a systematic CESE development of high order schemes, in this paper we describe a new 4th-order neutrally stable CESE solver of the advection equation Theta u/Theta + alpha Theta u/Theta x = 0. The space-time stencil of this two-level explicit scheme is formed by one point at the upper time level and three points at the lower time level. Because it is associated with three independent mesh variables u(sup n) (sub j), (u(sub x))(sup n) (sub j) , and (uxz)(sup n) (sub j) (the numerical analogues of u, Theta u/Theta x, and Theta(exp 2)u/Theta x(exp 2), respectively) and four equations per mesh point, the new scheme is referred to as the alpha(3) scheme. As in the case of other similar CESE neutrally stable solvers, the alpha(3) scheme enforces conservation laws in space-time locally and globally, and it has the basic, forward marching, and backward marching forms. These forms are equivalent and satisfy a space-time inversion (STI) invariant property which is shared by the advection equation. Based on the concept of STI invariance, a set of algebraic relations is developed and used to prove that the alpha(3) scheme must be neutrally stable when it is stable. Moreover it is proved rigorously that all three amplification factors of the alpha(3) scheme are of unit magnitude for all phase angles if |v| <= 1/2 (v = alpha delta t/delta x). This theoretical result is consistent with the numerical stability condition |v| <= 1/2. Through numerical experiments, it is established that the alpha(3) scheme generally is (i) 4th-order accurate for the mesh variables u(sup n) (sub j) and (ux)(sup n) (sub j); and 2nd-order accurate for (uxx)(sup n) (sub j). However, in some exceptional cases, the scheme can achieve perfect accuracy aside from round-off errors.

Chang, Sin-Chung↗

NASA NTRS2008

The a(3) Scheme--A Fourth-Order Space-Time Flux-Conserving and Neutrally Stable CESE Solver

The CESE development is driven by a belief that a solver should (i) enforce conservation laws in both space and time, and (ii) be built from a non-dissipative (i.e., neutrally stable) core scheme so that the numerical dissipation can be controlled effectively. To initiate a systematic CESE development of high order schemes, in this paper we provide a thorough discussion on the structure, consistency, stability, phase error, and accuracy of a new 4th-order space-time flux-conserving and neutrally stable CESE solver of an 1D scalar advection equation. The space-time stencil of this two-level explicit scheme is formed by one point at the upper time level and three points at the lower time level. Because it is associated with three independent mesh variables (the numerical analogues of the dependent variable and its 1st-order and 2ndorder spatial derivatives, respectively) and three equations per mesh point, the new scheme is referred to as the a(3) scheme. Through the von Neumann analysis, it is shown that the a(3) scheme is stable if and only if the Courant number is less than 0.5. Moreover, it is established numerically that the a(3) scheme is 4th-order accurate.

Chang, Sin-Chung↗

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

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 NTRS1991

High order filtering methods for approximating hyperbolic systems of conservation laws

The essentially nonoscillatory (ENO) schemes, while potentially useful in the computation of discontinuous solutions of hyperbolic conservation-law systems, are computationally costly relative to simple central-difference methods. A filtering technique is presented which employs central differencing of arbitrarily high-order accuracy except where a local test detects the presence of spurious oscillations and calls upon the full ENO apparatus to remove them. A factor-of-three speedup is thus obtained over the full-ENO method for a wide range of problems, with high-order accuracy in regions of smooth flow.

Lafon, F.↗

NASA NTRS1976

Numerical solution of nonlinear partial differential equations of mixed type

A review is presented of some recently developed numerical methods for the solution of nonlinear equations of mixed type. The methods considered use finite difference approximations to the differential equation. Central difference formulas are employed in the subsonic zone and upwind difference formulas are used in the supersonic zone. The relaxation method for the small disturbance equation is discussed and a description is given of difference schemes for the potential flow equation in quasi-linear form. Attention is also given to difference schemes for the potential flow equation in conservation form, the analysis of relaxation schemes by the time dependent analogy, the accelerated iterative method, and three-dimensional calculations.

Jameson, A.↗

NASA NTRS1986

Higher-order flux difference splitting schemes for the Euler equations using upstream interpolations

A class of explicit two time-level, 2p + 1 space-point, (2p 1)-th order, upwind-biased flux difference splitting schemes are proposed for the numerical advection based on Lagrange's interpolation, and the method is an accord with the physical domain of dependence. A normalized Jacobian coefficient matrix is introduced to convert the schemes to hyperbolic systems of conservation laws, and approaches to make the higher-order schemes total variation stable are discussed. Accuracy and stability of the present schemes are examined, and implicit total variation diminishing schemes are developed for steady-state calculations.Application to gasdynamic problems for both steady and unsteady flows covering a wide range of Mach numbers is considered, and results for a blast wave passing a cylinder, and head-on collision of two blast waves over a circular arc, are presented. The flow patterns were found to be symmetric, and good resolution of flow structures was obtained.

Yang, J. Y.↗

NASA NTRS1989

Thickening oscillation of a delta wing using Navier-Stokes and Navier-displacement equations

The problem of unsteady, supersonic, locally conical, vortical flow around a delta wing undergoing thickening oscillation is solved using the unsteady, thin-layer Navier-Stokes equations and the unsteady linearized Navier-displacement equations. The unsteady, thin-layer Navier-Stokes equations are solved using the implicit approximate-factorization finite-volume scheme to compute the conservative components of the flow vector field. With the conservative components known at any time step, the linearized, Navier-displacement equations are solved using the alternating, direction-implicit scheme to obtain the grid points displacements due to known displacement boundary conditions. A grid-displacements limiter, in the form of a low mesh Reynolds number, is used to limit grid-folding in regions of highly reversed flow.

Chuang, Hsin-Kung A.↗

NASA NTRS1986

Finite-volume scheme for transonic potential flow about airfoils and bodies in an arbitrarily shaped channel

A conservative finite-volume difference scheme is developed for the potential equation to solve transonic flow about airfoils and bodies in an arbitrarily shaped channel. The scheme employs a mesh which is a nearly conformal O mesh about the airfoil and nearly orthogonal at the channel walls. The mesh extends to infinity upstream and downstream, where the mapping is singular. Special procedures are required to treat the singularities at infinity, including computation of the metrics near those points. Channels with exit areas different from inlet areas are solved; a body with a sting mount is an example of such a case.

South, Jerry C., Jr.↗

NASA NTRS1992

Low-to-high altitude predictions of three-dimensional ablative reentry flowfields

A new coupled-chemistry PNS scheme has been developed to study 3D hypersonic reentry flows over sphere cones at small angles of attack. A new quasi-steady ablation model for Teflon was used to predict ablative reentry flow fields under laminar, transitional, and turbulent conditions. This solution scheme uses a new two-step solution scheme for solving the species conservation equations under low-altitude conditions. As an example, the scheme is used to study hypersonic flow over a 6-deg sphere-cone configuration at 50 and 125 kft conditions under zero and nonzero angle of attack conditions.

Bhutta, Bilal A.↗

NASA NTRS1972

A two-spacecraft test of a single spacecraft method of estimating shock normals

By assuming the validity of a subset of the Rankine-Hugoniot conservation relations for interplanetary (IP) shocks in an isotropic medium it has been demonstrated, in principle, that improved shock normals can be calculated by using a least-squares technique on combined magnetic field and plasma data from a single spacecraft. The scheme uses those six conservation relations not involving pressure and temperature. This paper deals with a test of the scheme by examining in detail a shock across which the magnetic field changed direction by a small amount (approximately 10 deg). On January 26, 1968 at about 1430 UT this shock was observed by the plasma and magnetic field instruments in Explorers 33 and 35. The spacecraft were 76.6 and 56.9 R sub E sunward of the earth, respectively (and 43.5 R sub E from each other), and therefore well outside the earth's bow shock region, a necessary condition for a valid test.

Lepping, R. P.↗

NASA NTRS1988

High-resolution shock-capturing schemes for inviscid and viscous hypersonic flows

A class of implicit Total Variation Diminishing (TVD) type algorithms suitable for transonic and supersonic multidimensional Euler and Navier-Stokes equations was extended to hypersonic computations. The improved conservative shock-capturing schemes are spatially second- and third-order, and are fully implicit. They can be first- or second-order accurate in time and are suitable for either steady or unsteady calculations. Enhancement of stability and convergence rate for hypersonic flows is discussed. With the proper choice of the temporal discretization and suitable implicit linearization, these schemes are fairly efficient and accurate for very complex two-dimensional hypersonic inviscid and viscous shock interactions. This study is complimented by a variety of steady and unsteady viscous and inviscid hypersonic blunt-body flow computations. Due to the inherent stiffness of viscous flow problems, numerical experiments indicated that the convergence rate is in general slower for viscous flows than for inviscid steady flows.

Yee, H. C.↗

NASA NTRS1990

High-resolution shock-capturing schemes for inviscid and viscous hypersonic flows

Hypersonic computations are presently conducted with an extension of a class of high-resolution implicit TVD algorithms suited to transonic multidimensional Euler and Navier-Stokes equations. These conservative shock-capturing schemes, which are spatially second- and third-order, may be first- and second-order accurate in time and suitable for either steady or unsteady calculations. Attention is given to the enhancement of hypersonic flows' convergence rate and stability; accuracy and efficiency is achieved by these means for very complex two-dimensional hypersonic viscous and inviscid shock interactions.

Yee, H. C.↗