Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

On a Non-Reflecting Boundary Condition for Hyperbolic Conservation Laws

A non-reflecting boundary condition (NRBC) for practical computations in fluid dynamics and aeroacoustics is presented. The technique is based on the first principle of non-reflecting, plane wave propagation and the hyperbolicity of the Euler equation system. The NRBC is simple and effective, provided the numerical scheme maintains locally a C(sup 1) continuous solution at the boundary. Several numerical examples in 1D, 2D, and 3D space are illustrated to demonstrate its robustness in practical computations.

Loh, Ching Y.↗
NASA NTRS2014

On Formulations of Discontinuous Galerkin and Related Methods for Conservation Laws

A formulation for the discontinuous Galerkin (DG) method that leads to solutions using the differential form of the equation (as opposed to the standard integral form) is presented. The formulation includes (a) a derivative calculation that involves only data within each cell with no data interaction among cells, and (b) for each cell, corrections to this derivative that deal with the jumps in fluxes at the cell boundaries and allow data across cells to interact. The derivative with no interaction is obtained by a projection, but for nodal-type methods, evaluating this derivative by interpolation at the nodal points is more economical. The corrections are derived using the approximate (Dirac) delta functions. The formulation results in a family of schemes: different approximate delta functions give rise to different methods. It is shown that the current formulation is essentially equivalent to the flux reconstruction (FR) formulation. Due to the use of approximate delta functions, an energy stability proof simpler than that of Vincent, Castonguay, and Jameson (2011) for a family of schemes is derived. Accuracy and stability of resulting schemes are discussed via Fourier analyses. Similar to FR, the current formulation provides a unifying framework for high-order methods by recovering the DG, spectral difference (SD), and spectral volume (SV) schemes. It also yields stable, accurate, and economical methods.

Huynh, H. T.↗
NASA NTRS2003

On a Non-Reflecting Boundary Condition for Hyperbolic Conservation Laws

A non-reflecting boundary condition (NRBC) for practical computations in fluid dynamics and aeroacoustics is presented. The technique is based on the hyperbolicity of the Euler equation system and the first principle of plane (simple) wave propagation. The NRBC is simple and effective, provided the numerical scheme maintains locally a C(sup 1) continuous solution at the boundary. Several numerical examples in ID, 2D and 3D space are illustrated to demonstrate its robustness in practical computations.

Loh, Ching Y.↗
NASA NTRS1984

Kinematics of the CS method for the treatment of molecular collisions

The problem of the quantum treatment of nonreactive collisions of two simple molecules or of an atom and a molecule is considered mathematically. Kinematical structure theorems analogous to the Wigner-Eckart theorem for scalar operators are derived for the transition operator (T-operator) of a two-molecule system, where, in addition to the usual conservation laws, certain additional conservation laws are presumed satisfied in a collision. The additional conservation laws were examined in order to develop a testable approximation scheme (CS method) for the specification of the kinematics of a system of two rigid diatomic molecules. The rand p-helicity basis states for a two-rigid rotor system are defined and unitary transformations established between these and conventional bases. It is found that only questionable criteria exist for testing the on-the-energy-shell matrix of a T(energy) operator for the presence of properties of two of the additional conservation laws.

Hahne, G. E.↗
NASA NTRS1991

Choice of velocity variables for complex flow computation

The issue of adopting the velocity components as dependent velocity variables for the Navier-Stokes flow computations is investigated. The viewpoint advocated is that a numerical algorithm should preferably honor both the physical conservation law in differential form and the geometric conservation law in discrete form. With the use of Cartesian velocity vector, the momentum equations in curvilinear coordinates can retain the full conservation-law form and satisfy the physical conservation laws. With the curvilinear velocity components, source terms appear in differential equations and hence the full conservation law form can not be retained. In discrete expressions, algorithms based on the Cartesian components can satisfy the geometric conservation-law form for convection terms but not for viscous terms; those based on the curvilinear components, on the other hand, cannot satisfy the geometric conservation-law form for either convection or viscous terms. Several flow solutions for domain with 90 and 360 degree turnings are presented to illustrate the issues of using the Cartesian velocity components and the staggered grid arrangement.

Shyy, W.↗
NASA NTRS1998

Discretization and Preconditioning Algorithms for the Euler and Navier-Stokes Equations on Unstructured Meshes

Chapter 1 briefly reviews several related topics associated with the symmetrization of systems of conservation laws and quasi-conservation laws: (1) Basic Entropy Symmetrization Theory; (2) Symmetrization and eigenvector scaling; (3) Symmetrization of the compressible Navier-Stokes equations; and (4) Symmetrization of the quasi-conservative form of the magnetohydrodynamic (MHD) equations. Chapter 2 describes one of the best known tools employed in the study of differential equations, the maximum principle: any function f(x) which satisfies the inequality f(double prime)>0 on the interval [a,b] attains its maximum value at one of the endpoints on the interval. Chapter three examines the upwind finite volume schemes for scalar and system conservation laws. The basic tasks in the upwind finite volume approach have already been presented: reconstruction, flux evaluation, and evolution. By far, the most difficult task in this process is the reconstruction step.

Bart, Timothy J.↗
NASA NTRS1990

Shock-layer bounds for a singularly perturbed equation

The size of the shock-layer governed by a conservation law is studied. The conservation law is a parabolic reaction-convection-diffusion equation with a small parameter multiplying the diffusion term and convex flux. Rigorous upper and lower bounding functions for the solution of the conservation law are established based on maximum-principle arguments. The bounding functions demonstrate that the size of the shock-layer is proportional to the parameter multiplying the diffusion term.

Scroggs, Jeffrey S.↗
DOE OSTI2020

The geometric theory of charge conservation in particle-in-cell simulations

In recent years, several gauge-symmetric particle-in-cell (PIC) methods have been developed whose simulations of particles and electromagnetic fields exactly conserve charge. While it is rightly observed that these methods’ gauge symmetry gives rise to their charge conservation, this causal relationship has generally been asserted via ad hoc derivations of the associated conservation laws. In this work, we develop a comprehensive theoretical grounding for charge conservation in gauge-symmetric Lagrangian and Hamiltonian PIC algorithms. For Lagrangian variational PIC methods, we apply Noether’s second theorem to demonstrate that gauge symmetry gives rise to a local charge conservation law as an off-shell identity. For Hamiltonian splitting methods, we show that the momentum map establishes their charge conservation laws. We define a new class of algorithms – gauge-compatible splitting methods – that exactly preserve the momentum map associated with a Hamiltonian system’s gauge symmetry – even after time discretization. This class of algorithms affords splitting schemes a decided advantage over alternative Hamiltonian integrators. We apply this general technique to design a novel, explicit, symplectic, gauge-compatible splitting PIC method, whose momentum map yields an exact local charge conservation law. Finally, our study clarifies the appropriate initial conditions for such schemes and examines their symplectic reduction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1946

Some Effects of Compressibility on the Flow Through Fans and Turbines

The laws of conservation of mass, momentum, and energy are applied to the compressible flow through a two-dimensional cascade of airfoils. A fundamental relation between the ultimate upstream and downstream flow angles, the inlet Mach number, and the pressure ratio across the cascade is derived. Comparison with the corresponding relation for incompressible flow shows large differences. The fundamental relation reveals two ranges of flow angles and inlet Mach numbers, for which no ideal pressure ratio exists. One of these nonideal operating ranges is analogous to a similar type in incompressible flow. The other is characteristic only of compressible flow. The effect of variable axial-flow area is treated. Some implications of the basic conservation laws in the case of nonideal flow through cascades are discussed.

Perl, W.↗
DOE OSTI2021

A moving discontinuous Galerkin finite element method with interface condition enforcement for compressible flows

A variation of moving discontinuous Galerkin finite element method with interface condition enforcement (MDG-ICE) is developed for solving the compressible Euler equations. The MDG-ICE method, originating from the work of Corrigan et al. [1], [2], [3], [4], is based on the space-time DG formulation, where both flow field and grid geometry are considered as independent variables and the conservation laws are enforced both on discrete elements and element interfaces. The element conservation laws are solved in the standard discontinuous solution space to determine conservative quantities, while the interface conservation is enforced using a variational formulation in a continuous space to determine discrete grid geometry. The resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation can then be solved in a least-squares sense, leading to an unconstrained nonlinear least-squares problem that is regularized and solved by Levenberg-Marquardt method. A number of numerical experiments for both 1D unsteady and 2D steady state compressible flow problems are conducted to assess the accuracy and robustness of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to implicitly detect and track all types of discontinuities via interface conservation enforcement and satisfy the conservation law on both elements and interfaces via grid movement and grid management, demonstrating that an exponential rate of convergence for Sod and Lax-Harden shock tube problems can be achieved and highly accurate solutions without overheating to both double-rarefaction wave and Noh problems can be obtained.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Hidden Killing fields, geometric symmetries and black hole mergers

Highlights: • Construction of geometric models with a variable degree of symmetry. • Formulation of Energy and Momentum Balance Laws. • Description of backreactions on spacetime geometries caused by dissipative effects. • Analytic description of accretion and/or merger processes in GR. • Exact model for the merger of two extremal Reissner–Nordström black holes. In the present work, using the recently introduced framework of local geometric deformations, special types of vector fields – so-called hidden Killing vector fields – are constructed, which solve the Killing equation not globally, but only locally, i.e. in local subregions of spacetime. Taking advantage of the fact that the vector fields coincide locally with Killing fields and therefore allow the consideration of integral laws that convert into exact physical conservation laws on local scales, balance laws in dynamical systems without global Killing symmetries are derived that mimic as closely as possible the conservation laws for energy and angular momentum of highly symmetric models. The utility of said balance laws is demonstrated by a concrete geometric example, namely a toy model for the binary merger of two extremal Reissner–Nordström black holes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1987

Lectures series in computational fluid dynamics

The lecture notes cover the basic principles of computational fluid dynamics (CFD). They are oriented more toward practical applications than theory, and are intended to serve as a unified source for basic material in the CFD field as well as an introduction to more specialized topics in artificial viscosity and boundary conditions. Each chapter in the test is associated with a videotaped lecture. The basic properties of conservation laws, wave equations, and shock waves are described. The duality of the conservation law and wave representations is investigated, and shock waves are examined in some detail. Finite difference techniques are introduced for the solution of wave equations and conservation laws. Stability analysis for finite difference approximations are presented. A consistent description of artificial viscosity methods are provided. Finally, the problem of nonreflecting boundary conditions are treated.

Thompson, Kevin W.↗
NASA NTRS1978

The 3 K blackbody radiation, Dirac's Large Numbers Hypothesis, and scale-covariant cosmology

A program is described which is intended to derive a generalized system of gravitational equations that allow (but do not require) G to vary, to use the 3-K blackbody radiation to fix the relation between G and the gauge function, and to employ Dirac's (1937) Large Numbers Hypothesis to derive the geometry of the universe. Einstein's equations are retained in their total integrity, but the specification is made that they are valid only when gravitational units are used. A scale-invariant form of Einstein's equations is obtained, and from this are derived the energy conservation law, the baryon-number conservation law, and the appropriate cosmological equations. Dirac's proposals of 1937 and 1973 are incorporated into the formalism, and a gauge based on consolidation of the 3-K blackbody radiation is presented. A unique solution for the geometry of the universe is determined for zero curvature solely from the 3-K radiation and the Large Numbers Hypothesis; this solution predicts a deceleration parameter exactly equal to unity.

Canuto, V.↗
DOE OSTI2020

Energy and momentum conservation in the Euler–Poincaré formulation of local Vlasov–Maxwell-type systems

The action principle by Low for the classic Vlasov–Maxwell system contains a mix of Eulerian and Lagrangian variables. This renders the Noether analysis of reparametrization symmetries inconvenient, especially since the well-known energy- and momentum-conservation laws for the system are expressed in terms of Eulerian variables only. While an Euler–Poincaré formulation of Vlasov–Maxwell-type systems, effectively starting with Low's action and using constrained variations for the Eulerian description of particle motion, has been known for a while Cendra et al, it is hard to come by a documented derivation of the related energy- and momentum-conservation laws in the spirit of the Euler–Poincaré machinery. To our knowledge only one such derivation exists in the literature so far, dealing with the so-called guiding-center Vlasov–Darwin system Sugama et al. The present exposition discusses a generic class of local Vlasov–Maxwell-type systems, with a conscious choice of adopting the language of differential geometry to exploit the Euler–Poincaré framework to its full extent. After reviewing the transition from a Lagrangian picture to an Eulerian one, we demonstrate how symmetries generated by isometries in space lead to conservation laws for linear- and angular-momentum density and how symmetry by time translation produces a conservation law for energy density. Here, we also discuss what happens if no symmetries exist. Finally, two explicit examples will be given—the classic Vlasov–Maxwell and the drift-kinetic Vlasov–Maxwell—and the results expressed in the language of regular vector calculus for familiarity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1999

Acoustic Energy Estimates in Inhomogeneous Moving Media

In ducted fan engine noise research, there is a need for defining a simple and easy to use acoustic energy conservation law to help in quantification of noise control techniques. There is a well known conservation law relating acoustic energy and acoustic energy flux in the case of an isentropic irrotational flow. Several different approaches have been taken to generalize this conservation law. For example, Morfey finds an identity by separating out the irrotational part of the perturbed flow. Myers is able to find a series of indentities by observing an algebraic relationship between the basic conservation of energy equation for a background flow and the underlying equations of motion. In an approximate sense, this algebraic relationship is preserved under perturbation. A third approach which seems to have not been pursued in the literature is a result known as Noether's theorem. There is a Lagrangian formulation for the Euler equation of fluid mechanics. Noether's theorem says that any group action that leaves the Lagrangian action invariant leads to a conserved quantity. This presentation will include a survey of current results regarding acoustic energy and preliminary results on the symmetries of the Lagrangian.

Farassat, F.↗
NASA NTRS1983

A numerical method based on the Fourier-Fourier transform approach for modeling 1-D electron plasma evolution

A numerical method is presented for studying one-dimensional electron plasma evolution under typical interplanetary conditions. The method applies the Fourier-Fourier transform approach to a plasma model that is a generalization of the electrostatic Vlasov-Poisson system of equations. Conservation laws that are modified to include the plasma model generalization and also the boundary effects of nonperiodic solutions are given. A new conservation law for entropy in the transformed space is then introduced. These conservation laws are used to verify the numerical solutions. A discretization error analysis is presented. Two numerical instabilities and the methods used for their suppression are treated. It is shown that in interplanetary plasma conditions, the bump-on-tail instability produces significant excitation of plasma oscillations at the Bohm-Gross frequency and its second harmonic. An explanation of the second harmonic excitation is given in terms of wave-wave coupling during the growth phase of the instability.

Klimas, A. J.↗