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 NTRS2000

Detection of Extraterrestrial Ecology (Exoecology)

Researchers in the Astrobiology Technology Branch at Ames Research Center have begun investigating alternate concepts for the detection of extraterrestrial life. We suggest searching for extraterrestrial ecology, exoecology, as well as for extraterrestrial biology, exobiology. Ecology describes the interactions of living things with their environment. All ecosystems are highly constrained by their environment and are constrained by well-known system design principles. Ecology could exist wherever there is an energy source and living I things have discovered some means to capture, store, and use the available energy. Terrestrial ecosystems use as energy sources, light, organic molecules, and in thermal vents and elsewhere, simple inorganic molecules. Ecosystem behavior is controlled by matter and energy conservation laws and can be described by linear and nonlinear dynamic systems theory. Typically in an ecosystem different molecules are not in chemical equilibrium and scarce material is conserved, stored, or recycled. Temporal cycles and spatial variations are often observed. These and other -eneral principles of exoecology can help guide the search for extraterrestrial life. The chemical structure observed in terrestrial biology may be highly contingent on evolutionary accidents. Oxygen was not always abundant on Earth. Primitive sulfur bacteria use hydrogen sulfide and sulfur to perform photosynthesis instead of water and oxygen. Astrobiologists have assumed, for the sake of narrowing and focusing our life detection strategies, that extraterrestrial life will have detailed chemical similarities with terrestrial life. Such assumptions appear very reasonable and they allow us to design specific and highly sensitive life detection experiments. But the fewer assumptions we make, the less chance we have of being entirely wrong The best strategy for the detection of extraterrestrial life could be a mixed strategy. We should use detailed assumptions based on terrestrial biology to guide some but not all future searches for alien life. The systems principles of exoecology seem much more fundamental and inescapable than the terrestrial biology analogies of exobiology. We should search for exoecology as well as exobiology.

Jones, Harry↗
NASA NTRS2018

High-Order Stabilized Finite Elements on Dynamic Meshes

The development of dynamic mesh capability for turbulent flow simulations using the Streamlined Upwind Petrov-Galerkin (SUPG) discretization is described. The current work extends previous research to include high-order spatial accuracy, including the satisfaction of the discrete geometric conservation law (GCL) on curved elements. Two closely-related schemes are described and the ability of these schemes to satisfy the GCL, while also maintaining temporal accuracy and conservation is assessed. Studies indicate that although one scheme discretizes the time derivative in conservative form, both schemes exhibit temporal conservation errors that decrease according to the expected design order of accuracy. The source of the temporal conservation errors is examined, and it is demonstrated that many finite-volume and finite-element schemes can also be expected to have difficulty strictly satisfying conservation in time. The effects on conservation are examined and, while present in the simulations, are seen to be negligible for the problems considered.

Anderson, W. Kyle↗
NASA NTRS2001

On the Conservation and Convergence to Weak Solutions of Global Schemes

In this paper we discuss the issue of conservation and convergence to weak solutions of several global schemes, including the commonly used compact schemes and spectral collocation schemes, for solving hyperbolic conservation laws. It is shown that such schemes, if convergent boundedly almost everywhere, will converge to weak solutions. The results are extensions of the classical Lax-Wendroff theorem concerning conservative schemes.

Carpenter, Mark H.↗
DOE OSTI2022

Rotating cosmologies: classical and quantum

Here, we revisit spatially flat, anisotropic cosmologies within the framework of mini-superspace. Putting special emphasis on the symmetries of the mini-superspace action and on the associated conservation laws, we unveil a new class of rotating cosmologies driven by solid matter. Their rotating is physical, in that it is characterized in an invariant way in terms of a conserved angular momentum. Along the way, we confirm the results of Bartolo et al. regarding the slow decay of anisotropies for solid inflation. We then use our minisuperspace approach as a laboratory to address certain puzzles of quantum cosmology — among these, how to characterize the spacetime symmetries of a quantum state at the level of the wavefunction of the universe. For the case of a solid driven cosmology, this question seems better defined than in more standard cases. Other questions remain unanswered, though; in particular, the general question of how to operate a minisuperspace-like truncation of degrees of freedom that is consistent at the quantum level.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1995

The direction of heat flow in an expanding current free stellar atmosphere

The energy budgets of the corona are substantially influenced by losses attributed to thermal conduction losses down the transition region, towards the chromosphere. Although Fourier law for such a plasma regime is now universally criticized as inappropriate, the direction of this heat flow seems to be certain. The purpose of this paper is to point out, based on conservation laws, the direction in which the heat must flow in this regime, provided it is known that there is a time independent expansion devoid of parallel currents. It will be argued that the direction of the heat flux is from the base of the transition region up into the corona, rather than the opposite. Under the circumstances mentioned above, the direction of the heat flow should be opposite to the gradient of the electron pressure gradient. Although the transition region is modeled as isobaric, this is only approximately true. Such a likely scenario makes the heat flow almost certainly to flow from the lower temperature base of the TR up into the incipient corona if a steady state current free expansion is to be maintained.

Scudder, J. D.↗
NASA NTRS2002

The Delta(dot) B = O Constraint vs. 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. On the other hand, the MHD equations can be derived from basic principles in either conservative or non-conservative form. The non-conservative system has a better conditioned eigensystem. The Delta(dot)B = 0 constraint of the A4HD equations is only an initial condition constraint. One does not need the Delta(dot)B condition to close the MHD system. We formulate our new low dissipative high order 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 Delta(dot)B 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.↗
DOE OSTI2021

On the violation of the zeroth law of turbulence in space plasmas

The zeroth law of turbulence states that, for fixed energy input into large-scale motions, the statistical steady state of a turbulent system is independent of microphysical dissipation properties. This behaviour, which is fundamental to nearly all fluid-like systems from industrial processes to galaxies, occurs because nonlinear processes generate smaller and smaller scales in the flow, until the dissipation – no matter how small – can thermalise the energy input. Using direct numerical simulations and theoretical arguments, we show that in strongly magnetised plasma turbulence such as that recently observed by the Parker Solar Probe spacecraft, the zeroth law is routinely violated. Namely, when such turbulence is ‘imbalanced’ – when the large-scale energy input is dominated by Alfvénic perturbations propagating in one direction (the most common situation in space plasmas) – nonlinear conservation laws imply the existence of a ‘barrier’ at scales near the ion gyroradius. This causes energy to build up over time at large scales. The resulting magnetic-energy spectra bear a strong resemblance to those observed in situ , exhibiting a sharp, steep kinetic transition range above and around the ion-Larmor scale, with flattening at yet smaller scales. The effect thus offers a possible solution to the decade-long puzzle of the position and variability of ion-kinetic spectral breaks in plasma turbulence. The existence of the ‘barrier’ also suggests that, how a plasma is forced at large scales (the imbalance) may have a crucial influence on thermodynamic properties such as the ion-to-electron heating ratio.

Physics↗
DOE OSTI2022

Momentum conservation in current drive and alpha-channeling-mediated rotation drive

Alpha channeling uses waves to extract hot ash from a fusion plasma, transferring energy from the ash to the wave. It has been proposed that this process could create a radial electric field, efficiently driving E × B rotation. However, existing theories ignore the nonresonant particles, which play a critical role in enforcing momentum conservation in quasilinear theory. Because cross field charge transport and momentum conservation are fundamentally linked, this non-consistency throws the rotation drive into question. This paper has two main goals. First, we provide a pedantic and cohesive introduction to the recently developed simple, general, self-consistent quasilinear theory for electrostatic waves that explains the torques which allow for current drive parallel to the magnetic field, and charge extraction across it; a theory that has largely resolved the question of rotation drive by alpha channeling. We show how the theory reveals a fundamental difference between the reaction of nonresonant particles to plane waves that grow in time vs steady-state waves that have a nonuniform spatial structure, allowing rotation drive in the latter case while precluding it in the former, and we review the local and global conservation laws that lead to this result. Second, we provide two new results in support of the theory. First, we provide a novel two-particle Hamiltonian model that rigorously establishes the relationship between charge transport and momentum conservation. Second, we compare the new quasilinear theory to the oscillation-center theories of ponderomotive forces, showing how the latter often obscure the time-dependent nonresonant recoil, but ultimately lead to similar results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

A new flux conserving Newton's method scheme for the two-dimensional, steady Navier-Stokes equations

A new numerical method is developed for the solution of the two-dimensional, steady Navier-Stokes equations. The method that is presented differs in significant ways from the established numerical methods for solving the Navier-Stokes equations. The major differences are described. First, the focus of the present method is on satisfying flux conservation in an integral formulation, rather than on simulating conservation laws in their differential form. Second, the present approach provides a unified treatment of the dependent variables and their unknown derivatives. All are treated as unknowns together to be solved for through simulating local and global flux conservation. Third, fluxes are balanced at cell interfaces without the use of interpolation or flux limiters. Fourth, flux conservation is achieved through the use of discrete regions known as conservation elements and solution elements. These elements are not the same as the standard control volumes used in the finite volume method. Fifth, the discrete approximation obtained on each solution element is a functional solution of both the integral and differential form of the Navier-Stokes equations. Finally, the method that is presented is a highly localized approach in which the coupling to nearby cells is only in one direction for each spatial coordinate, and involves only the immediately adjacent cells. A general third-order formulation for the steady, compressible Navier-Stokes equations is presented, and then a Newton's method scheme is developed for the solution of incompressible, low Reynolds number channel flow. It is shown that the Jacobian matrix is nearly block diagonal if the nonlinear system of discrete equations is arranged approximately and a proper pivoting strategy is used. Numerical results are presented for Reynolds numbers of 100, 1000, and 2000. Finally, it is shown that the present scheme can resolve the developing channel flow boundary layer using as few as six to ten cells per channel width, depending on the Reynolds number.

Scott, James R.↗
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 NTRS1989

Algorithm For Hypersonic Flow In Chemical Equilibrium

Implicit, finite-difference, shock-capturing algorithm calculates inviscid, hypersonic flows in chemical equilibrium. Implicit formulation chosen because overcomes limitation on mathematical stability encountered in explicit formulations. For dynamical portion of problem, Euler equations written in conservation-law form in Cartesian coordinate system for two-dimensional or axisymmetric flow. For chemical portion of problem, equilibrium state of gas at each point in computational grid determined by minimizing local Gibbs free energy, subject to local conservation of molecules, atoms, ions, and total enthalpy. Major advantage: resulting algorithm naturally stable and captures strong shocks without help of artificial-dissipation terms to damp out spurious numerical oscillations.

Palmer, Grant↗
NASA NTRS2005

Multi-Dimensional Spectral Difference Method for Unstructured Grids

A new, high-order, conservative, and efficient method for conservation laws on unstructured grids is developed. It combines the best features of structured and unstructured grid methods to attain computational efficiency and geometric flexibility; it utilizes the concept of discontinuous and high-order local representations to achieve conservation and high accuracy; and it is based on the finite-difference formulation for simplicity. Universal reconstructions are obtained by distributing unknowns in a geometrically similar manner for all unstructured cells. Placements of the unknown and flux points with various order of accuracy are given for the line, triangular and tetrahedral elements. The data structure of the new method permits an optimum use of cache memory, resulting in further computational efficiency on modern computers. A new pointer system is developed that reduces memory requirements and simplifies programming for any order of accuracy. Numerical solutions are presented and compared with the exact solutions for wave propagation problems in both two and three dimensions to demonstrate the capability of the method. Excellent agreement has been found. The method is simpler and more efficient than previous discontinuous Galerkin and spectral volume methods for unstructured grids.

Liu, Yen↗
NASA NTRS1994

Green's formula and variational principles for cosmic-ray transport with application to rotating and shearing flows

Green's theorem and Green's formula for the diffusive cosmic-ray transport equation in relativistic flows are derived. Green's formula gives the solution of the transport equation in terms of the Green's function of the adjoint transport equation, and in terms of distributed sources throughout the region R of interest, plus terms involving the particle intensity and streaming on the boundary. The adjoint transport equation describes the time-reversed particle transport. An Euler-Lagrange variational principle is then obtained for both the mean scattering frame distribution function f, and its adjoint f(dagger). Variations of the variational functional with respect to f(dagger) yield the transport equation, whereas variations of f yield the adjoint transport equation. The variational principle, when combined with Noether's theorem, yields the conservation law associated with Green's theorem. An investigation of the transport equation for steady, azimuthal, rotating flows suggests the introduction of a new independent variable H to replace the comoving frame momentum variable p'. For the case of rigid rotating flows, H is conserved and is shown to be analogous to the Hamiltonian for a bead on a rigidly rotating wire. The variable H corresponds to a balance between the centrifugal force and the particle inertia in the rotating frame. The physical interpretation of H includes a discussion of nonrelativistic and special relativistic rotating flows as well as the cases of aziuthal, differentially rotating flows about Schwarzs-child and Kerr black holes. Green's formula is then applied to the problem of the acceleration of ultra-high-energy cosmic rays by galactic rotation. The model for galactic rotation assumes an angular velocity law Omega = Omega(sub 0)(omega(sub 0)/omega), where omega denotes radial distance from the axis of rotation. Green's functions for the galactic rotation problem are used to investigate the spectrum of accelerated particles arising from monoenergetic and truncated power-law sources. We conclude that it is possible to accelerate particles beyond the knee by galactic rotation, but not in sufficient number to adequately explain the observed spectrum.

Webb, G. M.↗
DOE OSTI2025

Moment-preserving Monte-Carlo Coulomb collision method for particle codes

Binary-pairing Monte-Carlo methods are widely used in particle-in-cell codes to capture effects of small angle Coulomb collisions. These methods preserve momentum and energy exactly when the simulation particles have equal weights. However, when the interacting particles are of varying weight, these physical conservation laws are only preserved on average. Here, we 1) extend these methods to weighted particles such that the scattering physics is correct on average, and 2) describe a new method for adjusting the particle velocities post scatter to restore exact conservation of momentum and energy. In conclusion, the efficacy of the model is illustrated with various test problems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

Entropy conservative finite element schemes

The question of entropy stability for discrete approximations to hyperbolic systems of conservation laws is studied. The amount of numerical viscosity present in such schemes is quantified and related to their entropy stability by means of comparison. To this end, two main ingredients are used: entropy variables and the construction of certain entropy conservative schemes in terms of piecewise-linear finite element approximations. It is then shown that conservative schemes are entropy stable, if and (for three-point schemes) only if, they contain more numerical viscosity than the abovementioned entropy conservation ones.

Tadmor, E.↗
DOE OSTI2024

A Local Macroscopic Conservative (LoMaC) Low Rank Tensor Method for the Vlasov Dynamics

Abstract In this paper, we propose a novel Local Macroscopic Conservative (LoMaC) low rank tensor method for simulating the Vlasov-Poisson (VP) system. The LoMaC property refers to the exact local conservation of macroscopic mass, momentum and energy at the discrete level. This is a follow-up work of our previous development of a conservative low rank tensor approach for Vlasov dynamics ( arXiv:2201.10397 ). In that work, we applied a low rank tensor method with a conservative singular value decomposition to the high dimensional VP system to mitigate the curse of dimensionality, while maintaining the local conservation of mass and momentum. However, energy conservation is not guaranteed, which is a critical property to avoid unphysical plasma self-heating or cooling. The new ingredient in the LoMaC low rank tensor algorithm is that we simultaneously evolve the macroscopic conservation laws of mass, momentum and energy using a flux-difference form with kinetic flux vector splitting; then the LoMaC property is realized by projecting the low rank kinetic solution onto a subspace that shares the same macroscopic observables by a conservative orthogonal projection. The algorithm is extended to the high dimensional problems by hierarchical Tuck decomposition of solution tensors and a corresponding conservative projection algorithm. Extensive numerical tests on the VP system are showcased for the algorithm’s efficacy.

Guo, Wei↗
DOE OSTI2023

Collective dynamics of polarized spin-half fermions in relativistic heavy-ion collisions

Standard relativistic hydrodynamics has been successful in describing the properties of the strongly interacting matter produced in the heavy-ion collision experiments. Recently, there has been a significant theoretical advancement in this field to explain spin polarization of hadrons emitted in these processes. Although current models have successfully explained some of the experimental data based on the coupling between spin polarization and vorticity of the medium, they still lack a clear understanding of the differential measurements. This is commonly interpreted as an indication that the spin needs to be treated as an independent degree of freedom whose dynamics is not entirely bound to flow circulation. In particular, if the spin is a macroscopic property of the system, in equilibrium its dynamics should follow hydrodynamic laws. Here, we develop a framework of relativistic hydrodynamics which includes spin degrees of freedom from the quantum kinetic theory for Dirac fermions and use it for modeling the dynamics of matter. Following experimental observations, we assume that the polarization effects are small and derive conservation laws for the net baryon current, the energy–momentum tensor and the spin tensor based on the de Groot–van Leeuwen–van Weert definitions of these currents. We present various properties of the spin polarization tensor and its components, analyze the propagation properties of the spin polarization components, and derive the spin-wave velocity for arbitrary statistics. We find that only the transverse spin components propagate, analogously to the electromagnetic waves. Finally, using our framework, we study the space–time evolution of the spin polarization for the systems respecting certain space–time symmetries and calculate the mean spin polarization per particle, which can be compared to the experimental data. We find that, for some observables, our spin polarization results agree qualitatively with the experimental findings and other model calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1981

Numerical approximation of boundary conditions with applications to inviscid equations of gas dynamics

A comprehensive overview of the state of the art of well-posedness and stability analysis of difference approximations for initial boundary value problems of the hyperbolic type is presented. The applicability of recent theoretical development to practical calculations for nonlinear gas dynamics is examined. The one dimensional inviscid gas dynamics equations in conservation law form are selected for numerical experiments. The class of implicit schemes developed from linear multistep methods in ordinary differential equations is chosen and the use of linear extrapolation as an explicit or implicit boundary scheme is emphasized. Specification of boundary data in the primitive variables and computation in terms of the conservative variables in the interior is discussed. Some numerical examples for the quasi-one-dimensional nozzle are given.

Yee, H. C.↗