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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

A new flux-conserving numerical scheme for the steady, incompressible Navier-Stokes equations

This paper is concerned with the continued development of a new numerical method, the space-time solution element (STS) method, for solving conservation laws. The present work focuses on the two-dimensional, steady, incompressible Navier-Stokes equations. Using first an integral approach, and then a differential approach, the discrete flux conservation equations presented in a recent paper are rederived. Here a simpler method for determining the flux expressions at cell interfaces is given; a systematic and rigorous derivation of the conditions used to simulate the differential form of the governing conservation law(s) is provided; necessary and sufficient conditions for a discrete approximation to satisfy a conservation law in E2 are derived; and an estimate of the local truncation error is given. A specific scheme is then constructed for the solution of the thin airfoil boundary layer problem. Numerical results are presented which demonstrate the ability of the scheme to accurately resolve the developing boundary layer and wake regions using grids which are much coarser than those employed by other numerical methods. It is shown that ten cells in the cross-stream direction are sufficient to accurately resolve the developing airfoil boundary layer.

Scott, James R.↗
DOE OSTI2014

Variational formulation of macro-particle plasma simulation algorithms

A variation formulation of macro-particle kinetic plasma models is discussed. In the electrostatic case, the use of symplectic integrators is investigated and found to offer advantages over typical generic methods. For the electromagnetic case, gauge invariance and momentum conservation are considered in detail. It is shown that, while the symmetries responsible for these conservation laws are broken in the presence of a spatial grid, the conservation laws hold in an average sense. The requirements for exact invariance are explored and it is shown that one viable option is to represent the potentials with a truncated Fourier basis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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 NTRS1983

Convergence of generalized MUSCL schemes

Semi-discrete generalizations of the second order extension of Godunov's scheme, known as the MUSCL scheme, are constructed, starting with any three point ""E" scheme. They are used to approximate scalar conservation laws in one space dimension. For convex conservation laws, each member of a wide class is proven to be a convergent approximation to the correct physical solution. Comparison with another class of high resolution convergent schemes is made.

Osher, S.↗
NASA NTRS1984

Convergence of generalized MUSCL schemes

Semi-discrete generalizations of the second order extension of Godunov's scheme, known as the MUSCL scheme, are constructed, starting with any three point E scheme. They are used to approximate scalar conservation laws in one space dimension. For convex conservation laws, each member of a wide class is proven to be a convergent approximation to the correct physical solution. Comparison with another class of high resolution convergent schemes is made.

Osher, S.↗
DOE OSTI2024

Non-conservation and conservation for different formulations of moist potential vorticity

Potential vorticity (PV) is one of the most important quantities in atmospheric science. In the absence of dissipative processes, the PV of each fluid parcel is known to be conserved, for a dry atmosphere. However, a parcel's PV is not conserved if clouds or phase changes of water occur. Recently, PV conservation laws were derived for a cloudy atmosphere, where each parcel's PV is not conserved but parcel-integrated PV is conserved, for integrals over certain volumes that move with the flow. Hence a variety of different statements are now possible for moist PV conservation and non-conservation, and in comparison to the case of a dry atmosphere, the situation for moist PV is more complex. Here, in light of this complexity, several different definitions of moist PV are compared for a cloudy atmosphere. Numerical simulations are shown for a rising thermal, both before and after the formation of a cloud. These simulations include the first computational illustration of the parcel-integrated, moist PV conservation laws. The comparisons, both theoretical and numerical, serve to clarify and highlight the different statements of conservation and non-conservation that arise for different definitions of moist PV.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Variational integrators for reduced magnetohydrodynamics

Reduced magnetohydrodynamics is a simplified set of magnetohydrodynamics equations with applications to both fusion and astrophysical plasmas, possessing a noncanonical Hamiltonian structure and consequently a number of conserved functionals. We propose a new discretisation strategy for these equations based on a discrete variational principle applied to a formal Lagrangian. The resulting integrator preserves important quantities like the total energy, magnetic helicity and cross helicity exactly (up to machine precision). As the integrator is free of numerical resistivity, spurious reconnection along current sheets is absent in the ideal case. If effects of electron inertia are added, reconnection of magnetic field lines is allowed, although the resulting model still possesses a noncanonical Hamiltonian structure. After reviewing the conservation laws of the model equations, the adopted variational principle with the related conservation laws is described both at the continuous and discrete level. We verify the favourable properties of the variational integrator in particular with respect to the preservation of the invariants of the models under consideration and compare with results from the literature and those of a pseudo-spectral code.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Sparse invariant domain preserving discontinuous Galerkin methods with subcell convex limiting

In this paper, we develop high-order nodal discontinuous Galerkin (DG) methods for hyperbolic conservation laws that satisfy invariant domain preserving properties using subcell flux corrections and convex limiting. These methods are based on a subcell flux corrected transport (FCT) methodology that involves blending a high-order target scheme with a robust, low-order invariant domain preserving method that is obtained using a graph viscosity technique. Furthermore, the new low-order discretizations are based on sparse stencils which do not increase with the polynomial degree of the high-order DG method. As a result, the accuracy of the low-order method does not degrade when used with high-order target methods. The method is applied to both scalar conservation laws, for which the discrete maximum principle is naturally enforced, and to systems of conservation laws such as the Euler equations, for which positivity of density and a minimum principle for specific entropy are enforced. Numerical results are presented on a number of benchmark test cases.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Implicit shock tracking for unsteady flows by the method of lines

Here, a recently developed high-order implicit shock tracking (HOIST) framework for resolving discontinuous solutions of inviscid, steady conservation laws is extended to the unsteady case. Central to the framework is an optimization problem which simultaneously computes a discontinuity-aligned mesh and the corresponding high-order approximation to the flow, which provides nonlinear stabilization and a high-order approximation to the solution. This work extends the implicit shock tracking framework to the case of unsteady conservation laws using a method of lines discretization via a diagonally implicit Runge-Kutta method by “solving a steady problem at each timestep”. We formulate and solve an optimization problem that produces a feature-aligned mesh and solution at each Runge-Kutta stage of each timestep, and advance this solution in time by standard Runge-Kutta update formulas. A Rankine-Hugoniot based prediction of the shock location together with a high-order, untangling mesh smoothing procedure provides a high-quality initial guess for the optimization problem at each time, which results in rapid convergence of the sequential quadratic programing (SQP) optimization solver. This method is shown to deliver highly accurate solutions on coarse, high-order discretizations without nonlinear stabilization and recover the design accuracy of the Runge-Kutta scheme. We demonstrate this framework on a series of inviscid, unsteady conservation laws in both one- and two- dimensions. We also verify that our method is able to recover the design order of accuracy of our time integrator in the presence of a strong discontinuity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Re-appraisal and extension of the Gratton-Vargas two-dimensional analytical snowplow model of plasma focus. II. Looking at the singularity

The Gratton-Vargas snowplow model, recently revisited and expanded [S. K. H. Auluck, Phys. Plasmas 20, 112501 (2013)], has given rise to significant new insights into some aspects of the Dense Plasma Focus (DPF), in spite of being a purely kinematic description having no reference to plasma phenomena. It is able to provide a good fit to the experimental current waveforms in at least 4 large facilities. It has been used for construction of a local curvilinear frame of reference, in which conservation laws for mass, momentum, and energy can be reduced to effectively-one-dimensional hyperbolic conservation law equations. Its utility in global parameter optimization of device parameters has been demonstrated. These features suggest that the Gratton-Vargas model deserves a closer look at its supposed limitations near the singular phase of the DPF. This paper presents a discussion of its development near the device axis, based on the original work of Gratton and Vargas, with some differences. It is shown that the Gratton-Vargas partial differential equation has solutions for times after the current singularity, which exhibit an expanding bounded volume (which can serve as model of an expanding plasma column) and decreasing dynamic inductance of the discharge, in spite of having no built-in hydrodynamics. This enables the model to qualitatively reproduce the characteristic shape of the current derivative in DPF experiments without reference to any plasma phenomena, such as instabilities, anomalous resistance, or reflection of hydrodynamic shock wave from the axis. The axial propagation of the solution exhibits a power-law dependence on the dimensionless time starting from the time of singularity, which is similar to the power-law relations predicted by theory of point explosions in ideal gases and which has also been observed experimentally.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

A lattice Maxwell system with discrete space–time symmetry and local energy–momentum conservation

A lattice Maxwell system is developed here with gauge-symmetry, symplectic structure and discrete space–time symmetry. Noether's theorem for Lie group symmetries is generalized to discrete group symmetries for the lattice Maxwell system. As a result, the lattice Maxwell system is shown to admit a discrete local energy–momentum conservation law corresponding to the discrete space–time symmetry. A lattice model that respects all local conservation laws and geometric structures is as good as and probably more preferable than standard models on continuous space–time. It can also be viewed as an effective algorithm for the governing differential equations on continuous space–time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1985

Foundations for the numerical solution of the Euler equations

The Navier-Stokes equations represent an extremely good model of the physical phenomena encountered in most aeronautical problems. However, the computational resource needed to solve the Navier-Stokes equations are so large that even with today's supercomputers, it is necessary to make use of simpler models. A large number of external aerodynamic problems can be accurately described by a simpler model. This model consists of an outer inviscid flow plus a boundary-layer thickness correction for the vehicle shape. The outer inviscid model may be represented by the potential equation or by the Euler equation. The present paper provides the foundations for the numerical solution of the Euler equations. The governing equations are considered, taking into account conservation laws, the medium, the differential form of the conservation laws, generalized solutions, shock-fitting, and characteristics. Attention is also given to initial and boundary conditions, existence and uniqueness, and rotational phenomena.

Salas, M. D.↗
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.↗
DOE OSTI2013

Conjugate flow action functionals

We present a new general framework to construct an action functional for a non-potential field theory. The key idea relies on representing the governing equations relative to a diffeomorphic flow of curvilinear coordinates which is assumed to be functionally dependent on the solution field. Such flow, which will be called the conjugate flow, evolves in space and time similarly to a physical fluid flow of classical mechanics and it can be selected in order to symmetrize the Gâteaux derivative of the field equations with respect to suitable local bilinear forms. This is equivalent to requiring that the governing equations of the field theory can be derived from a principle of stationary action on a Lie group manifold. By using a general operator framework, we obtain the determining equations of such manifold and the corresponding conjugate flow action functional. In particular, we study scalar and vector field theories governed by second-order nonlinear partial differential equations. The identification of transformation groups leaving the conjugate flow action functional invariant could lead to the discovery of new conservation laws in fluid dynamics and other disciplines.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Effects of the guide field on electron distribution functions in the diffusion region of asymmetric reconnection

Crescent electron distribution functions in the electron diffusion region in asymmetric reconnection with a guide field (strength less than 50% of the upstream magnetic field) are investigated by theory and fully kinetic simulations. Electron motion in the electron diffusion region is a combination of meandering across the current sheet and gyration around the guide field. Combining the two canonical momentum conservation laws in a one-dimensional model with the energy conservation law, we derive the shape of crescent electron distribution functions in both velocity planes: perpendicular to the magnetic field and coplanner with the magnetic field. There are two major effects of the guide field on crescent distribution functions: (1) widening the opening angle of the crescent in a reduced distribution function and (2) the linear cut-off of a distribution function nonparallel to the magnetic field. We also discuss a translational mapping (perpendicular to the linear cut-off) of distribution functions between the X-line and another position. In addition, considering the gyration due to the normal magnetic field, we predict the shape of distribution functions in general positions in the electron diffusion region. Finally, the predictions will be useful to identify features of magnetic reconnection in a guide field by satellite observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A mathematical model of asynchronous data flow in parallel computers *

Abstract We present a simplified model of data flow on processors in a high-performance computing framework involving computations necessitating inter-processor communications. From this ordinary differential model, we take its asymptotic limit, resulting in a model which treats the computer as a continuum of processors and data flow as an Eulerian fluid governed by a conservation law. We derive a Hamilton–Jacobi equation associated with this conservation law for which the existence and uniqueness of solutions can be proven. We then present the results of numerical experiments for both discrete and continuum models; these show a qualitative agreement between the two and the effect of variations in the computing environment’s processing capabilities on the progress of the modelled computation.

Barnard, Richard C.↗
DOE OSTI2018

Perturbative variational formulation of the Vlasov-Maxwell equations

The perturbative variational formulation of the Vlasov-Maxwell equations is presented up to the third order in the perturbation analysis. From the second and third-order Lagrangian densities, the first-order and second-order Vlasov-Maxwell equations are expressed in gauge-invariant and gauge-independent forms, respectively. Upon deriving the reduced second-order Vlasov-Maxwell Lagrangian for the linear nonadiabatic gyrokinetic Vlasov-Maxwell equations, the reduced Lagrangian densities for the linear drift-wave equation and the linear hybrid kinetic-magnetohydrodynamic (MHD) equations are derived, with their associated wave-action conservation laws obtained by the Noether method. The exact wave-action conservation law for the linear hybrid kinetic-MHD equations is written explicitly. To conclude, a new form of the third-order Vlasov-Maxwell Lagrangian is derived in which ponderomotive effects play a crucial role.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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↗