Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conservation form”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Numerical solution of transonic full stream function equations in conservation form

The stream function equation in conservation form is solved iteratively based on the artificial compressibility method. The density is not a unique function of the mass flux. In order to avoid the ambiguity near the sonic line, the density is updated in terms of the velocity, which is obtained through a simple integration of a first order equation step by step in the flow field. Iteration algorithms and finite difference approximations are discussed and numerical results of both conservative and nonconservative calculations are presented.

Hafez, M. M.↗
NASA NTRS1984

Implicit treatment of the unsteady full potential equation in conservation form

An implicit, conservative treatment for the unsteady full potential equation in two-dimensions is presented. The method employs a local time linearization for density, and introduces flux biasing concepts based on sonic conditions for the generation of artificial viscosity to capture shocks without any overshoots. The boundary condition is treated implicitly using a splitting procedure consistent with the approximate factorization scheme. This allows for extremely large Courant numbers, even for nonorthogonal grid at the body. The method has application not only to unsteady problems, but also to generate the starting blunt body solution for a supersonic full potential marching code. Results are presented for flows over cylinders, spheres and airfoils. Comparisons are made with available Euler and full potential results, and are in excellent agreement.

Shankar, V.↗
DOE OSTI2020

Low-dissipation centred schemes for hyperbolic equations in conservative and non-conservative form

We propose a one-parameter family of low-dissipation centred numerical schemes for solving hyperbolic equations in conservative or non-conservative form, using finite volume or discontinuous Galerkin finite element methods. The new schemes spring out from the multi-dimensional FORCE method and are determined by a single parameter α ≥ 1. Given an increasing sequence of real numbers 1 ≤ α 1 < α 2 < . . . < α K , there corresponds a sequence of numerical schemes with stability restriction associated to a decreasing sequence of Courant numbers 1 > c 1 > c 2 > . . . > c K > 0 and a decreasing sequence of corresponding numerical viscosity functions d 1 > d 2 > . . . > d K . For a given Courant number c k ≤ 1 there is a real number α k > 0 and a corresponding stable scheme with minimal numerical viscosity d k . The proposed schemes suit very well the family of high-order discontinuous Galerkin finite element methods and the recently proposed class of ADER-TR schemes, whose orders of accuracy define decreasing sequences of Courant numbers, as the order of accuracy increases. The centred methods of this paper are stable in 2D and 3D in the frame of simultaneous updating formulae, unlike other centred methods, such as 1D FORCE, which are not. Furthermore, the schemes are highly accurate for slowly-moving waves, which is precisely the kind of waves that traditional centred methods smear disastrously. Additional features of the proposed schemes include ease of implementation and applicability to any hyperbolic system either in conservative or non-conservative form. Here, the proposed schemes are analysed and computationally assessed through a suite of test problems for a linear model system, for the Euler equations in one and two space dimensions, and for the Baer-Nunziato equations for compressible two-phase flow.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1975

Transonic potential flow calculations using conservation form

A method is presented for the solution of the full potential equation in conservation form. This assures that a weak solution satisfying proper isentropic jump conditions is obtained, giving an improved representation of shock waves in comparison with earlier nonconservative schemes. The method uses the concept of artificial viscosity to produce a stable difference scheme in the supersonic zone, and artificial time to generate a convergent iterative scheme for the solution of the difference equations. Results of calculations using the nonconservative and conservative schemes are compared.

Jameson, A.↗
NASA NTRS1981

Solution of the transonic full potential equation in conservative form using an implicit algorithm

The paper presents numerical solutions of the full potential equation in conservative form. The iteration scheme used is a fully implicit approximate factorization technique and provides a significant improvement in convergence speed relative to standard successive line overrelaxation algorithms. The spatial differencing algorithm is centrally differenced in both subsonic and supersonic regions to maintain stability. This effectively approximates rotated differencing, thereby greatly improving the reliability of the algorithm.

Holst, T.↗
NASA NTRS1980

An alternating-direction-implicit algorithm for the unsteady potential equation in conservation form

An implicit finite difference scheme for an efficient computation of unsteady potential flow about airfoils is presented. The formulation uses density and velocity potential as dependent variables, and is cast in conservation form to assure the theoretically correct determination of shockwave location and speed. To enable boundary conditions to be imposed directly on the airfoil surface, a time varying sheared rectilinear coordinate transformation is employed. Calculated time history solutions on a pulsating airfoil are compared with the results of another unsteady transonic code. It is concluded that the method has excellent numerical stability and gives accurate solutions with sharply resolved shocks.

Chipman, R. R.↗
NASA NTRS1980

Implicit computations of unsteady transonic flow governed by the full-potential equation in conservation form

An alternating-direction implicit algorithm is presented for solving the conservative, full-potential equation for unsteady, transonic flow. A new development is the time-linearization of the density function. This linearization reduces the solution process from one of solving a system of two equations at each mesh point to one of solving a single equation. Two sample cases are computed. First, a one-dimensional traveling shock wave is computed and compared with the analytic solution. Second, a two-dimensional case is computed of a flow field that results from a thickening and subsequently thinning airfoil. The resulting flow field, which includes a traveling shock wave, is compared to the flow field obtained from the low-frequency, small-disturbance, transonic equation.

Goorjian, P. M.↗
NASA NTRS1980

An implicit algorithm for the transonic full-potential equation in conservative form

A fast, implicit approximate factorization algorithm for the solution of the conservative full-potential equation for transonic flow in two and three dimensions is presented. Stability in supersonic regions is maintained by the use of an upwind evaluation of the density coefficient along all coordinate directions, providing an effective upwind difference of the streamwise terms for any orientation of the velocity vector and thereby enhancing the reliability of the algorithm. The algorithm is shown to provide rapid convergence for the computation of certain difficult two-dimensional test cases, including cases with fishtail shock patterns, demonstrating the reliability and efficiency of the procedure. Surface pressure coefficient distributions obtained by the present method are also found to be in good agreement with those computed by successive-line overrelaxation and a hybrid direct-solver/successive-line overrelaxation scheme, with significant reductions in CPU time required. A three-dimensional solution for a swept wing mounted between parallel walls is also presented which demonstrates the high convergence rate of the algorithm in three dimensions as well as two.

Holst, T.↗
NASA NTRS1986

Numerical modeling of rotor flows with a conservative form of the full-potential equations

A computer program has been developed to solve a three-dimensional conservative formulation of the full-potential equation. Its ability to solve transonic, unsteady rotor flows is demonstrated by comparison to forward flight non-lifting pressure data at low to moderate advance ratios. A 'split potential' formulation has been added to the code which incorporates known vorticity fields into the full-potential calculation. Using this methodology, rotor wake contributions have been incorporated into the computer code. Pressure results are presented for lifting rotors in hover. These results are compared to experimental data as well as to other predictions.

Strawn, R. C.↗
DOE OSTI2026

Neural entropy-stable conservative flux form neural networks for learning hyperbolic conservation laws

We propose a neural entropy-stable conservative flux form neural network (NESCFN) for learning hyperbolic conservation laws and their associated entropy functions directly from solution trajectories, without requiring any predefined numerical discretization. While recent neural network architectures have successfully integrated classical numerical principles into learned models, most rely on prior knowledge of the governing equations or assume a fixed discretization. Our approach removes this dependency by embedding entropy-stable design principles into the learning process itself, enabling the discovery of physically consistent dynamics in a fully data-driven setting. By jointly learning both the flux function and a corresponding entropy, NESCFN promotes conservation and entropy dissipation, which is critical for long-term stability and fidelity in the system of hyperbolic conservation laws. Furthermore, numerical results demonstrate that the method achieves stability and conservation over extended time horizons and accurately captures shock propagation speeds, even without oracle access to future-time solution profiles in the training data.

Conservative flux form↗
NASA NTRS1983

On the symmetric form of systems of conservation laws with entropy

The present investigation has the objective to review the general structure of system of conservation laws with entropy, giving particular attention to a presentation of symmetric formulations of the equations of gas dynamics. Systems of hyperbolic conservation laws are considered, taking into account the case in which a system of equations is symmetric hyperbolic, questions of symmetrization, and entropy functions. The Euler equations for a polytropic gas in conservation form are considered. The symmetric form retains the conservation properties of the equations. Weak solutions remain, therefore, unchanged. Of particular computational interest is the possibility of using the symmetric form.

Harten, A.↗
DOE OSTI2026

Entropy Stable Conservative Flux Form Neural Networks

We propose an entropy-stable conservative flux form neural network (CFN) to predict the dynamics of unknown governing conservation laws. The design of the network is based on the entropy-stable, second-order, and non-oscillatory Kurganov-Tadmor (KT) scheme. The proposed entropy-stable CFN, hereafter referred to as ESCFN, uses slope limiting as a denoising mechanism, ensuring accurate predictions in both noisy and sparse observation environments, as well as in both smooth and discontinuous regions. Importantly, our method is designed to predict long term dynamics of the unknown conservation law exclusively from a short temporal window of observed data, that is, without oracle knowledge of the PDE or later-time solution profiles. Numerical experiments demonstrate that the ESCFN achieves both stability and conservation while maintaining accuracy over extended time domains, and successfully predicts shock propagation speeds in long-term simulations. Furthermore, it is also robust to both noisy and sparse data environments.

Hyperbolic conservation laws↗
NASA NTRS1991

Flux-vector splitting algorithm for chain-rule conservation-law form

A flux-vector splitting algorithm with Newton-Raphson iteration was developed for the 'full compressible' Navier-Stokes equations cast in chain-rule conservation-law form. The algorithm is intended for problems with deforming spatial domains and for problems whose governing equations cannot be cast in strong conservation-law form. The usefulness of the algorithm for such problems was demonstrated by applying it to analyze the unsteady, two- and three-dimensional flows inside one combustion chamber of a Wankel engine under nonfiring conditions. Solutions were obtained to examine the algorithm in terms of conservation error, robustness, and ability to handle complex flows on time-dependent grid systems.

Shih, T. I.-P.↗
NASA NTRS1990

Viscous shock profiles and primitive formulations

Weak solutions of hyperbolic systems in primitive (non-conservation) form for which a consistent conservation form exists are considered. It is shown that primitive formulations, shock relations are not uniquely defined by the states to either side of the shock but also depend on the viscous path connecting the two. Scheme-dependent high order correction terms are derived that enforce consistent viscous shock profiles. The resulting primitive algorithm is conservative to the order of approximation. One dimensional Euler calculations of flows containing strong shocks clearly show that conservation errors in primitive flow calculations are of comparable quality.

Karni, S.↗
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 NTRS2002

The Delta x B = 0 Constraint Versus Minimization of Numerical Errors in MHD Simulations

The MHD equations are a system of non-strictly hyperbolic conservation laws. The non-convexity of the inviscid flux vector resulted in corresponding Jacobian matrices with undesirable properties. It has previously been shown by Powell et al. (1995) that an 'almost' equivalent MHD system in non-conservative form can be derived. This non-conservative system has a better conditioned eigensystem. Aside from Powell et al., the MHD equations can be derived from basic principles in either conservative or non-conservative form. The Delta x B = 0 constraint of the MHD equations is only an initial condition constraint, it is very different from the incompressible Navier-Stokes equations in which the divergence condition is needed to close the system (i.e., to have the same number of equations and the same number of unknown). In the MHD formulations, if Delta x B = 0 initially, all one needs is to construct appropriate numerical schemes that preserve this constraint at later time evolutions. In other words, one does not need the Delta x B condition to close the MHD system. We formulate our new scheme together with the Cargo & Gallice (1997) form of the MHD approximate Riemann solver in curvilinear grids for both versions of the MHD equations. A novel feature of our new method is that the well-conditioned eigen-decomposition of the non-conservative MHD equations is used to solve the conservative equations. This new feature of the method provides well-conditioned eigenvectors for the conservative formulation, so that correct wave speeds for discontinuities are assured. The justification for using the non-conservative eigen-decomposition to solve the conservative equations is that our scheme has a better control of the numerical error associated with the divergence of the magnetic condition. Consequently, computing both forms of the equations with the same eigen-decomposition is almost equivalent. It will be shown that this approach, using the non-conservative eigensystem when solving the conservative equations, also works well in the context of standard shock-capturing schemes.

Yee, H. C.↗
NASA NTRS1995

Nonlinear (time domain) and linearized (time and frequency domain) solutions to the compressible Euler equations in conservation law form

Two linearized solvers (time and frequency domain) based on a high resolution numerical scheme are presented. The basic approach is to linearize the flux vector by expressing it as a sum of a mean and a perturbation. This allows the governing equations to be maintained in conservation law form. A key difference between the time and frequency domain computations is that the frequency domain computations require only one grid block irrespective of the interblade phase angle for which the flow is being computed. As a result of this and due to the fact that the governing equations for this case are steady, frequency domain computations are substantially faster than the corresponding time domain computations. The linearized equations are used to compute flows in turbomachinery blade rows (cascades) arising due to blade vibrations. Numerical solutions are compared to linear theory (where available) and to numerical solutions of the nonlinear Euler equations.

Sreenivas, Kidambi↗