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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Parallel CE/SE Computations via Domain Decomposition

This paper describes the parallelization strategy and achieved parallel efficiency of an explicit time-marching algorithm for solving conservation laws. The Space-Time Conservation Element and Solution Element (CE/SE) algorithm for solving the 2D and 3D Euler equations is parallelized with the aid of domain decomposition. The parallel efficiency of the resultant algorithm on a Silicon Graphics Origin 2000 parallel computer is checked.

Himansu, Ananda↗
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↗
DOE OSTI2014

Shock dynamics of phase diagrams

A thermodynamic phase transition denotes a drastic change of state of a physical system due to a continuous change of thermodynamic variables, as for instance pressure and temperature. The classical van der Waals equation of state is the simplest model that predicts the occurrence of a critical point associated with the gas–liquid phase transition. Nevertheless, below the critical temperature theoretical predictions of the van der Waals theory significantly depart from the observed physical behaviour. We develop a novel approach to classical thermodynamics based on the solution of Maxwell relations for a generalised family of nonlocal entropy functions. This theory provides an exact mathematical description of discontinuities of the order parameter within the phase transition region, it explains the universal form of the equations of state and the occurrence of triple points in terms of the dynamics of nonlinear shock wave fronts. -- Highlights: •A new generalisation of van der Waals equation of state. •Description of phase transitions in terms of shock dynamics of state curves. •Proof of the universality of equations of state for a general class of models. •Interpretation of triple points as confluence of classical shock waves. •Correspondence table between thermodynamics and nonlinear conservation laws.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Gauge natural formulation of conformal gravity

We consider conformal gravity as a gauge natural theory. We study its conservation laws and superpotentials. We also consider the Mannheim and Kazanas spherically symmetric vacuum solution and discuss conserved quantities associated to conformal and diffeomorphism symmetries.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

On One Integrable Discrete System

In this paper, we study a system of nonlinear equations on a square graph related to the affine algebra A{sup (1)}{sub 1} . This system is the simplest representative of the class of discrete systems corresponding to affine Lie algebras. We find the Lax representation and construct hierarchies of higher symmetries. In neighborhoods of singular points ⋋ = 0 and ⋋ = ∞, we construct formal asymptotic expansions of eigenfunctions of the Lax pair and, based on these expansions, find series of local conservation laws for the system considered.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Angular momentum effects in fission

Background: The role of angular momentum in fission has long been discussed but the observable effects are difficult to quantify. Purpose: Here we discuss a variety of effects associated with angular momentum in fission and present quantitative illustrations. Methods: We employ the fission simulation model freya , which is well suited for this purpose because it obeys all conservation laws, including linear and angular momentum conservation at each step of the process. We first discuss the implementation of angular momentum in freya and then assess particular observables, including various correlated observables. We also study potential effects of neutron-induced fission of the low-lying isomeric state of U 235 relative to the ground state. Results: The fluctuations inherent in the fission process ensure that the spin of the initial compound nucleus has only a small influence on the fragment spins, which are therefore nearly uncorrelated. There is a marked correlation between the spin magnitude of the fission fragments and the photon multiplicity. We also consider the dynamical anisotropy caused by the rotation of an evaporating fragment and study especially the distribution of the projected neutron-neutron opening angles, showing that while it is dominated by the effect of the evaporation recoils, it is possible to extract the signal of the dynamical anisotropy by means of a Fourier decomposition. Finally, we note that the use of an isomeric target, U 235 m ( n th ,f), may enhance the symmetric yields and can thus result in higher neutron multiplicities for low total fragment kinetic energies. Conclusions: While the initial angular momentum of the fissioning nucleus tends to have little effect on the observables, those of the produced fragments influence the emitted neutrons and photons in a significant and correlated manner which may be exploited experimentally to elucidate the fission process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2012

Development of a Robust and Efficient Parallel Solver for Unsteady Turbomachinery Flows

The traditional design and analysis practice for advanced propulsion systems relies heavily on expensive full-scale prototype development and testing. Over the past decade, use of high-fidelity analysis and design tools such as CFD early in the product development cycle has been identified as one way to alleviate testing costs and to develop these devices better, faster and cheaper. In the design of advanced propulsion systems, CFD plays a major role in defining the required performance over the entire flight regime, as well as in testing the sensitivity of the design to the different modes of operation. Increased emphasis is being placed on developing and applying CFD models to simulate the flow field environments and performance of advanced propulsion systems. This necessitates the development of next generation computational tools which can be used effectively and reliably in a design environment. The turbomachinery simulation capability presented here is being developed in a computational tool called Loci-STREAM [1]. It integrates proven numerical methods for generalized grids and state-of-the-art physical models in a novel rule-based programming framework called Loci [2] which allows: (a) seamless integration of multidisciplinary physics in a unified manner, and (b) automatic handling of massively parallel computing. The objective is to be able to routinely simulate problems involving complex geometries requiring large unstructured grids and complex multidisciplinary physics. An immediate application of interest is simulation of unsteady flows in rocket turbopumps, particularly in cryogenic liquid rocket engines. The key components of the overall methodology presented in this paper are the following: (a) high fidelity unsteady simulation capability based on Detached Eddy Simulation (DES) in conjunction with second-order temporal discretization, (b) compliance with Geometric Conservation Law (GCL) in order to maintain conservative property on moving meshes for second-order time-stepping scheme, (c) a novel cloud-of-points interpolation method (based on a fast parallel kd-tree search algorithm) for interfaces between turbomachinery components in relative motion which is demonstrated to be highly scalable, and (d) demonstrated accuracy and parallel scalability on large grids (approx 250 million cells) in full turbomachinery geometries.

West, Jeff↗
NASA NTRS2019

Entropy Stable Method for the Euler Equations Revisited: Central Differencing via Entropy Splitting and SBP

The two decades old high order central differencing via entropy splitting and summation-by-parts (SBP) difference boundary closure of Ols- son & Oliger (1994), Gerritsen & Olsson (1996), and Yee et al. (2000) is revisited. The entropy splitting is a form of skew-symmetric splitting of the nonlinear Euler flux derivatives. Central differencing applied to the entropy splitting form of the Euler flux derivatives together with SBP difference operators will, hereafter, be referred to as entropy split schemes. This study is prompted by the recent growing interest in numerical methods for which a discrete entropy conservation law holds, a discrete global entropy conservation can be proved and/or the numerical method possesses a stable entropy in the framework of SBP difference operators and L2-energy norm estimate. The objective of this paper is to recast the entropy split scheme as the re- cent definition of an entropy stable method for central differencing with SBP operators for both periodic and non-periodic boundary conditions for non- linear Euler equations. Standard high order spatial central differencing as well as high order central spatial DRP (dispersion relation preserving) spatial differencing is part of the entropy stable methodology framework. Long time integration of 2D and 3D test cases is included to show the comparison of this efficient entropy stable method with the Tadmor-type of entropy conservative methods. Studies also include the comparison among the three skew-symmetric splittings on their nonlinear stability and accuracy performance without added numerical dissipations for smooth flows. These are, namely, entropy splitting, Ducros et al. splitting and the Kennedy & Grub- ber splitting.

Sjogreen, Bjorn↗
DOE OSTI2016

Fokker-Planck equation in the presence of a uniform magnetic field

The Fokker-Planck equation in the presence of a uniform magnetic field is derived which has the same form as the case of no magnetic field but with different Fokker-Planck coefficients. The coefficients are calculated explicitly within the binary collision model, which are free from infinite sums of Bessel functions. They can be used to investigate relaxation and transport phenomena conveniently. The kinetic equation is also manipulated into the Landau form from which it is straightforward to compare with previous results and prove the conservation laws.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Deviations from the Gutenberg–Richter law on account of a random distribution of block sizes

This paper studies properties of a continuum with structure. The characteristic size of the structure governs the fact that difference relations are nonautomatically transformed into differential ones. It is impossible to consider an infinitesimal volume of a body, to which the major conservation laws could be applied, because the minimum representative volume of the body must contain at least a few elementary microstructures. The corresponding equations of motion are equations of infinite order, solutions of which include, along with usual sound waves, unusual waves with abnormally low velocities without a lower limit. It is shown that in such media weak perturbations can increase or decrease outside the limits. The number of complex roots of the corresponding dispersion equation, which can be interpreted as the number of unstable solutions, depends on the specific surface of cracks and is an almost linear dependence on a logarithmic scale, as in the seismological Gutenberg–Richter law. If the distance between one pore (crack) to another one is a random value with some distribution, we must write another dispersion equation and examine different scenarios depending on the statistical characteristics of the random distribution. In this case, there are sufficient deviations from the Gutenberg–Richter law and this theoretical result corresponds to some field and laboratory observations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Strongly intensive quantities

Analysis of fluctuations of hadron production properties in collisions of relativistic particles profits from use of measurable intensive quantities which are independent of system size variations. The first family of such quantities was proposed in 1992; another is introduced in this paper. Furthermore we present a proof of independence of volume fluctuations for quantities from both families within the framework of the grand canonical ensemble. These quantities are referred to as strongly intensive ones. Influence of conservation laws and resonance decays is also discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Reduced MHD in Astrophysical Applications: Two-dimensional or Three-dimensional?

Originally proposed as an efficient approach to computation of nonlinear dynamics in tokamak fusion research devices, reduced magnetohydrodynamics (RMHD) has subsequently found application in studies of coronal heating, flux tube dynamics, charged particle transport, and, in general, as an approximation to describe plasma turbulence in space physics and astrophysics. Given the diverse set of derivations available in the literature, there has emerged some level of discussion and a lack of consensus regarding the completeness of RMHD as a turbulence model, and its applicability in contexts such as the solar wind. Some of the key issues in this discussion are examined here, emphasizing that RMHD is properly neither 2D nor fully 3D, being rather an incomplete representation that enforces at least one family of extraneous conservation laws.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Glucans monomer-exchange dynamics as an open chemical network

We describe the oligosaccharides-exchange dynamics performed by the so-called D-enzymes on polysaccharides. To mimic physiological conditions, we treat this process as an open chemical network by assuming some of the polymer concentrations fixed (chemostatting). We show that three different long-time behaviors may ensue: equilibrium states, nonequilibrium steady states, and continuous growth states. We dynamically and thermodynamically characterize these states and emphasize the crucial role of conservation laws in identifying the chemostatting conditions inducing them.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2013

Functional entropy variables: A new methodology for deriving thermodynamically consistent algorithms for complex fluids, with particular reference to the isothermal Navier–Stokes–Korteweg equations

We propose a new methodology for the numerical solution of the isothermal Navier–Stokes–Korteweg equations. Our methodology is based on a semi-discrete Galerkin method invoking functional entropy variables, a generalization of classical entropy variables, and a new time integration scheme. We show that the resulting fully discrete scheme is unconditionally stable-in-energy, second-order time-accurate, and mass-conservative. We utilize isogeometric analysis for spatial discretization and verify the aforementioned properties by adopting the method of manufactured solutions and comparing coarse mesh solutions with overkill solutions. Various problems are simulated to show the capability of the method. Our methodology provides a means of constructing unconditionally stable numerical schemes for nonlinear non-convex hyperbolic systems of conservation laws.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

A Freestream-Preserving High-Order Finite-Volume Method for Mapped Grids with Adaptive-Mesh Refinement

A fourth-order accurate finite-volume method is presented for solving time-dependent hyperbolic systems of conservation laws on mapped grids that are adaptively refined in space and time. Novel considerations for formulating the semi-discrete system of equations in computational space combined with detailed mechanisms for accommodating the adapting grids ensure that conservation is maintained and that the divergence of a constant vector field is always zero (freestream-preservation property). Advancement in time is achieved with a fourth-order Runge-Kutta method.

42 ENGINEERING↗
DOE OSTI2015

Energy and momentum preserving Coulomb collision model for kinetic Monte Carlo simulations of plasma steady states in toroidal fusion devices

A kinetic Monte Carlo model suited for self-consistent transport studies is proposed and tested. The Monte Carlo collision operator is based on a widely used model of Coulomb scattering by a drifting Maxwellian and a new algorithm enforcing the momentum and energy conservation laws. The difference to other approaches consists in a specific procedure of calculating the background Maxwellian parameters, which does not require ensemble averaging and, therefore, allows for the use of single-particle algorithms. This possibility is useful in transport balance (steady state) problems with a phenomenological diffusive ansatz for the turbulent transport, because it allows a direct use of variance reduction methods well suited for single particle algorithms. In addition, a method for the self-consistent calculation of the electric field is discussed. Results of testing of the new collision operator using a set of 1D examples, and preliminary results of 2D modelling in realistic tokamak geometry, are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Equilibration in fermionic systems

Highlights: •Local equilibration in finite fermion systems is based on a nonlinear diffusion equation. •The fermionic diffusion equation is solved analytically in closed form. •The exact solutions agree with the numerical results. •Particle-number and energy conservation are shown to be fulfilled. •The model is applied to hadrons at low energies, and relativistic quarks. -- Abstract: The time evolution of a finite fermion system towards local statistical equilibrium is investigated using analytical solutions of a nonlinear partial differential equation that had been derived earlier from the Boltzmann collision term. The solutions of this fermionic diffusion equation are rederived in closed form, evaluated exactly for simplified initial conditions, and applied to hadron systems at low energies in the MeV-range, as well as to quark systems at relativistic energies in the TeV-range where antiparticle production is abundant. Conservation laws for particle number including created antiparticles, and for the energy are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗