Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conservation equations”

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 847 records · Page 47
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Application of a nonisentropic full potential method to AGARD standard airfoils

An entropy-correction method for the unsteady full potential equation is presented. The unsteady potential equation is modified to model the entropy jumps across shock waves. The conservative form of the modified equation is solved in generalized coordinates using an implicit, approximate factorization method. A flux-biasing differencing method, which generates the proper amounts of artificial viscosity in supersonic regions, is used to discretize the flow equations in space. Calculated results are presented for the NLR 7301, NACA 0012, and NACA 64A010A airfoils. Comparisons of the present method and solutions of the Euler equations are presented for the NLR 7301 airfoil, and comparisons of the present method and experimental data are presented for all three airfoils. The comparisons show that the present method more accurately models solutions of the Euler equations and experiment than does the isentropic potential formulation. In addition, it is shown that modeling shock-generated entropy extends the range of validity of the full potential method.

Whitlow, Woodrow, Jr.↗
DOE OSTI2026

A sweeping positivity-preserving high-order finite difference WENO scheme for Euler equations

We develop a simple, high-order, conservative and robust positivity-preserving sweeping procedure for the density and the nonlinear pressure function in the compressible Euler equations. Using the scaling limiter in Zhang and Shu (J Comput Phys 229:3091–3120, 2010), we obtain a non-trivial extension of the scalar sweeping technique in Liu et al. (J Sci Comput 73:1028–1071, 2017) for the positivity of pressure. The sweeping procedure developed in this paper is a post-processing technique, which can be applied to any concave functions of the conserved variables in hyperbolic conservation law systems. Thus, it has applications beyond the Euler equations. This procedure preserves positivity and conservation of physical quantities without destroying the accuracy of the underlying scheme. The algorithm works for general schemes including finite difference, finite volume and discontinuous Galerkin methods; however, in this paper we focus on finite difference weighted essentially non-oscillatory (WENO) methods. As a result, we provide numerical tests of the fifth-order finite difference WENO scheme to demonstrate the accuracy and robustness of the technique.

Compressible Euler equations↗
NASA NTRS1993

A semi-Lagrangian approach to the shallow water equation

We present a formulation of the shallow water equations that emphasizes the conservation of potential vorticity. A locally conservative semi-Lagrangian time-stepping scheme is developed, which leads to a system of three coupled PDE's to be solved at each time level. We describe a smoothing analysis of these equations, on which an effective multigrid solver is constructed. Some results from applying this solver to the static version of these equations are presented.

Bates, J. R.↗
NASA NTRS1988

Properties of finite difference models of non-linear conservative oscillators

Finite-difference (FD) approaches to the numerical solution of the differential equations describing the motion of a nonlinear conservative oscillator are investigated analytically. A generalized formulation of the Duffing and modified Duffing equations is derived and analyzed using several FD techniques, and it is concluded that, although it is always possible to contstruct FD models of conservative oscillators which are themselves conservative, caution is required to avoid numerical solutions which do not accurately reflect the properties of the original equation.

Mickens, R. E.↗
DOE OSTI2021

Fluctuations and Conservation Laws

We report on recent progress concerning effects of global conservation laws on cumulants of conserved quantities. Specifically, we will relate for an arbitrary equation of state cumulants of a conserved charge measured in a subvolume of a thermal system with the corresponding grandcanonical susceptibilities, taking into account exact global conservation of that charge. Applications to actual measurement at the RHIC and LHC as well as extensions to multiple conserved charges will be discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1980

An alternative approach to the numerical simulation of steady inviscid flow

A numerical procedure for the efficient simulation of steady inviscid flow is described and its utility is demonstrated. The method is uniformly valid for application in the subsonic, transonic and supersonic flow regimes. It does not rely on the introduction of additional assumptions beyond those necessary to obtain the Euler equations from the Navier-Stokes equations, nor does it make use of a time-asymptotic solution of the unsteady equations of motion. Application of the herein-defined surrogate equation technique allows the formulation of stable, fully-conservative, type-dependent finite difference equations for use in obtaining numerical solutions to systems of first-order partial differential equations, such as the steady-state Euler equations or their various approximations. Computational results are presented for the full Euler equations used to simulate rotational subsonic flow and for the transonic small disturbance equations. For the latter case, a computational efficiency greater than that obtained by means of the standard perturbation potential approach is indicated.

Johnson, G. M.↗
NASA NTRS1987

Conservative finite volume solutions of a linear hyperbolic transport equation in two and three dimensions using multiple grids

The feasibility of the multiple grid technique is investigated by solving linear hyperbolic equations for simple two- and three-dimensional cases. The results are compared with exact solutions and those obtained from the single grid calculations. It is demonstrated that the technique works reasonably well when two grid systems contain grid cells of comparative sizes. The study indicates that use of the multiple grid does not introduce any significant error and that it can be used to attack more complex problems.

Tiwari, Surendra N.↗
NASA NTRS1992

A time-accurate algorithm for chemical non-equilibrium viscous flows at all speeds

A time-accurate, coupled solution procedure is described for the chemical nonequilibrium Navier-Stokes equations over a wide range of Mach numbers. This method employs the strong conservation form of the governing equations, but uses primitive variables as unknowns. Real gas properties and equilibrium chemistry are considered. Numerical tests include steady convergent-divergent nozzle flows with air dissociation/recombination chemistry, dump combustor flows with n-pentane-air chemistry, nonreacting flow in a model double annular combustor, and nonreacting unsteady driven cavity flows. Numerical results for both the steady and unsteady flows demonstrate the efficiency and robustness of the present algorithm for Mach numbers ranging from the incompressible limit to supersonic speeds.

Shuen, J.-S.↗
NASA NTRS1993

A coupled implicit method for chemical non-equilibrium flows at all speeds

The present time-accurate coupled-solution procedure addresses the chemical nonequilibrium Navier-Stokes equations over a wide Mach-number range uses, in conjunction with the strong conservation form of the governing equations, five unknown primitive variables. The numerical tests undertaken address steady convergent-divergent nozzle flows with air dissociation/recombination, dump combustor flows with n-pentane/air chemistry, and unsteady nonreacting cavity flows.

Shuen, Jian-Shun↗
NASA NTRS1982

Numerical Solution for Navier-Stokes Equations

Carefully selected blend of computational techniques solves complete set of equations for viscous, unsteady, hypersonic flow in general curvilinear coordinates. New algorithm has tested computation of axially directed flow about blunt body having shape similar to that of such practical bodies as wide-body aircraft or artillery shells. Method offers significant computational advantages because of conservation-law form of equations and because it reduces amount of metric data required.

Warsi, Z. U. A.↗
DOE OSTI2026

Parallel-in-Time Solution of Allen-Cahn Equations by Integrating Operator Learning into the Parareal Method

While recent advances in deep learning have shown promising efficiency gains in solving time-dependent partial differential equations (PDEs), matching the accuracy of conventional numerical solvers still remains a challenge. One strategy to improve the accuracy of deep learning-based solutions for time-dependent PDEs is to use the learned model as the coarse propagator in the Parareal method and a traditional numerical method as the fine solver. However, successful integration of deep learning into the Parareal method requires consistency between the coarse and fine solvers, particularly for PDEs exhibiting rapid changes such as sharp transitions. Here, to ensure this consistency, we propose using convolutional neural networks (CNNs) to learn the fully discrete time-stepping operator defined by the same numerical scheme employed as the fine solver. We demonstrate the effectiveness of the proposed method in solving the classical and mass-conservative Allen–Cahn (AC) equations. Through iterative updates in the Parareal algorithm, our approach achieves a significant computational speedup compared to traditional fine solvers while converging to high-accuracy solutions. Our results highlight that the proposed hybrid Parareal algorithm effectively accelerates simulations, particularly when implemented on multiple GPUs, and converges to the desired accuracy in only a few iterations. Another advantage of our method is that the CNN model is trained on trajectory-based data generated from random initial conditions, such that the trained model can be used to solve the AC equations with various initial conditions without retraining. This work demonstrates the potential of integrating neural network methods into parallel-in-time frameworks for efficient and accurate simulations of time-dependent PDEs.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1981

On various treatments of potential equations at shocks

The potential equation that governs transonic inviscid flowfields is elliptic in subsonic regions and hyperbolic in supersonic regions. The transition from supersonic to subsonic flow may occur either continuously or across a surface of discontinuity known as a shock. Because of the changing type of governing equation and the presence of the surfaces of discontinuity, special treatments of the governing equations are needed at shocks to preserve mass conservation and to ensure satisfaction of the governing equation. In the so-called shock-capturing techniques, shocks are captured by smearing the discontinuity over several mesh cells in the computational domain. Most existing conservative schemes conserve mass flux throughout the flowfield, including shock boundaries; however, they also introduce zero-order errors in the approximation to the governing equation. It is shown that the zero-order errors, which do not diminish as the mesh spacing approaches zero, can cause discrepancies in the prediction of shock strength and location. These parameters are extremely sensitive features of transonic inviscid flowfield calculations. It is also shown that the zero-order errors can be avoided by developing higher-order schemes that properly model the governing equation at shocks.

Chen, L. T.↗
NASA NTRS2001

High Order Difference Method for Low Mach Number Aeroacoustics

A high order finite difference method with improved accuracy and stability properties for computational aeroacoustics (CAA) at low Mach numbers is proposed. The Euler equations are split into a conservative and a symmetric non- conservative portion to allow the derivation of a generalized energy estimate. Since the symmetrization is based on entropy variables, that splitting of the flux derivatives is referred to as entropy splitting. Its discretization by high order central differences was found to need less numerical dissipation than conventional conservative schemes. Owing to the large disparity of acoustic and stagnation quantities in low Mach number aeroacoustics, the split Euler equations are formulated in perturbation form. The unknowns are the small changes of the conservative variables with respect to their large stagnation values. All nonlinearities and the conservation form of the conservative portion of the split flux derivatives can be retained, while cancellation errors are avoided with its discretization opposed to the conventional conservative form. The finite difference method is third-order accurate at the boundary and the conventional central sixth-order accurate stencil in the interior. The difference operator satisfies the summation by parts property analogous to the integration by parts in the continuous energy estimate. Thus, strict stability of the difference method follows automatically. Spurious high frequency oscillations are suppressed by a characteristic-based filter similar to but without limiter. The time derivative is approximated by a 4-stage low-storage second-order explicit Runge-Kutta method. The method has been applied to simulate vortex sound at low Mach numbers. We consider the Kirchhoff vortex, which is an elliptical patch of constant vorticity rotating with constant angular frequency in irrotational flow. The acoustic pressure generated by the Kirchhoff vortex is governed by the 2D Helmholtz equation, which can be solved analytically using separation of variables.

Mueller, B.↗
NASA NTRS1985

Approximate-factorization schemes for solving the transonic full-potential equation

The present paper provides a general discussion of approximate-factorization techniques applied to the transonic full-potential equation. Giving particular attention to the AF2 approximate-factorization scheme. This scheme was first introduced by Ballhaus and Steger (1975) for solving the low-frequency (unsteady), transonic small-disturbance equation. The full-potential equation algorithm is examined, taking into account the governing equations, grid generation, the artificial density scheme (spatial differencing), the alternating direction implicit scheme, the AF2 iteration scheme, temporal damping, and boundary conditions. Computed results are also presented. It is shown that fast, fully-implicit algorithms of the approximate-factorization variety are both efficient and reliable for solving the conservative full-potential equation.

Holst, T. L.↗
NASA NTRS1983

High resolution applications of the Osher upwind scheme for the Euler equations

The 'Osher' scheme was introduced by Osher (1981). It represents an upwind finite-difference method for hyperbolic systems of conservation laws, including the Euler equations. In studies conducted by Osher (1981) and Osher and Solomon (1982), the method was applied to the nonisentropic form of the Euler equations in one dimension and the isentropically restricted form in two spatial dimensions, both in Cartesian coordinates. Chakravarthy and Osher (1982) have shown an approach for extending the Osher scheme to the Euler equations written for general geometries, taking into account the use of mappings to arbitrary curvilinear coordinate systems. The present investigation is concerned with the high resolution extension of the Osher scheme to second-order accuracy. Results are presented for several example problems, giving attention to quasi-one-dimensional Laval nozzle flow, a one-dimensional shock tube problem, and supersonic flow over a cylinder.

Chakravarthy, S. R.↗
NASA NTRS2011

Orion Parachute Riser Cutter Development

This paper presents the tests and analytical approach used on the development of a steel riser cutter for the CEV Parachute Assembly System (CPAS) used on the Orion crew module. Figure 1 shows the riser cutter and the steel riser bundle which consists of six individual cables. Due to the highly compressed schedule, initial unavailability of the riser material and the Orion Forward Bay mechanical constraints, JSC primarily relied on a combination of internal ballistics analysis and LS-DYNA simulation for this project. Various one dimensional internal ballistics codes that use standard equation of state and conservation of energy have commonly used in the development of CAD devices for initial first order estimates and as an enhancement to the test program. While these codes are very accurate for propellant performance prediction, they usually lack a fully defined kinematic model for dynamic predictions. A simple piston device can easily and accurately be modeled using an equation of motion. However, the accuracy of analytical models is greatly reduced on more complicated devices with complex external loads, nonlinear trajectories or unique unlocking features. A 3D finite element model of CAD device with all critical features included can vastly improve the analytical ballistic predictions when it is used as a supplement to the ballistic code. During this project, LS-DYNA structural 3D model was used to predict the riser resisting load that was needed for the ballistic code. A Lagrangian model with eroding elements shown in Figure 2 was used for the blade, steel riser and the anvil. The riser material failure strain was fine tuned by matching the dent depth on the anvil with the actual test data. LS-DYNA model was also utilized to optimize the blade tip design for the most efficient cut. In parallel, the propellant type and the amount were determined by using CADPROG internal ballistics code. Initial test results showed a good match with LS-DYNA and CADPROG simulations. Final paper will present a detailed roadmap from initial ballistic modeling and LS-DYNA simulation to the performance testing. Blade shape optimization study will also be presented.

Oguz, Sirri↗
NASA NTRS1977

Equations and closure methods

Basic differential equations governing compressible turbulent boundary layer flow are reviewed, including conservation of mass and energy, momentum equations derived from Navier-Stokes equations, and equations of state. Closure procedures were broken down into: (1) simple or zeroth-order methods, (2) first-order or mean field closure methods, and (3) second-order or mean turbulence field methods.

Source record↗