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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1972

A study of the liquid-vapor phase change of mercury based on irreversible thermodynamics.

The object of this work is to determine the transport coefficients which appear in linear irreversible-thermodynamic rate equations of a phase change. An experiment which involves the steady-state evaporation of mercury was performed to measure the principal transport coefficient appearing in the mass-rate equation and the coupling transport coefficient appearing in both the mass-rate equation and the energy-rate equation. The principal transport coefficient sigma, usually termed the 'condensation' or 'evaporation' coefficient, is found to be approximately 0.9, which is higher than that measured previously in condensation-of-mercury experiments. The experimental value of the coupling coefficient K does not agree with the value predicted from Schrage's kinetic analysis of the phase change. A modified kinetic analysis in which the Onsager reciprocal law and the conservation laws are invoked is presented which removes this discrepancy but which shows that the use of Schrage's equation for predicting mass rates of phase change is a good approximation.

Adt, R. R., Jr.↗
NASA NTRS1991

A Runge-Kutta discontinuous finite element method for high speed flows

A Runge-Kutta discontinuous finite element method is developed for hyperbolic systems of conservation laws in two space variables. The discontinuous Galerkin spatial approximation to the conservation laws results in a system of ordinary differential equations which are marched in time using Runge-Kutta methods. Numerical results for the two-dimensional Burger's equation show that the method is (p+1)-order accurate in time and space, where p is the degree of the polynomial approximation of the solution within an element and is capable of capturing shocks over a single element without oscillations. Results for this problem also show that the accuracy of the solution in smooth regions is unaffected by the local projection and that the accuracy in smooth regions increases as p increases. Numerical results for the Euler equations show that the method captures shocks without oscillations and with higher resolution than a first-order scheme.

Bey, Kim S.↗
NASA NTRS1999

Bidirectional Reflectance of a Macroscopically Flat, High-Albedo Particulate Surface: An Efficient Radiative Transfer Solution and Applications to Regoliths

Many remote sensing applications rely on accurate knowledge of the bidirectional reflection function (BRF) of surfaces composed of discrete, randomly positioned scattering particles. Theoretical computations of BRFs for plane-parallel particulate layers are usually reduced to solving the radiative transfer equation (RTE) using one of existing exact or approximate techniques. Since semi-empirical approximate approaches are notorious for their low accuracy, violation of the energy conservation law, and ability to produce unphysical results, the use of numerically exact solutions of RTE has gained justified popularity. For example, the computation of BRFs for macroscopically flat particulate surfaces in many geophysical publications is based on the adding-doubling (AD) and discrete ordinate (DO) methods. A further saving of computer resources can be achieved by using a more efficient technique to solve the plane-parallel RTE than the AD and DO methods. Since many natural particulate surfaces can be well represented by the model of an optically semi-infinite, homogeneous scattering layer, one can find the BRF directly by solving the Ambartsumian's nonlinear integral equation using a simple iterative technique. In this way, the computation of the internal radiation field is avoided and the computer code becomes highly efficient and very accurate and compact. Furthermore, the BRF thus obtained fully obeys the fundamental physical laws of energy conservation and reciprocity. In this paper, we discuss numerical aspects and the computer implementation of this technique, examine the applicability of the Henyey-Greenstein phase function and the sigma-Eddington approximation in BRF and flux calculations, and describe sample applications demonstrating the potential effect of particle shape on the bidirectional reflectance of flat regolith surfaces. Although the effects of packing density and coherent backscattering are currently neglected, they can also be incorporated. The FORTRAN implementation of the technique is available on the World Wide Web, and can be applied to a wide range of remote sensing problems. BRF computations for undulated (macroscopically rough) surfaces are more complicated and often rely on time consuming Monte Carlo procedures. This approach is especially inefficient for optically thick, weakly absorbing media (e.g., snow and desert surfaces at visible wavelengths since a photon may undergo many internal scattering events before it exists the medium or is absorbed. However, undulated surfaces can often be represented as collections of locally flat tilted facets characterized by the BRF found from the traditional plane parallel RTE. In this way the MOnte Carlo procedure could be used only to evaluate the effects of surface shadowing and multiple surface reflections, thereby bypassing the time-consuming ray tracing inside the medium and providing a great savings of CPU time.

Mishchenko, Michael I.↗
DOE OSTI2016

Dynamics of the diffusive DM-DE interaction – Dynamical system approach

We discuss dynamics of a model of an energy transfer between dark energy (DE) and dark matter (DM) . The energy transfer is determined by a non-conservation law resulting from a diffusion of dark matter in an environment of dark energy. The relativistic invariance defines the diffusion in a unique way. The system can contain baryonic matter and radiation which do not interact with the dark sector. We treat the Friedman equation and the conservation laws as a closed dynamical system. The dynamics of the model is examined using the dynamical systems methods for demonstration how solutions depend on initial conditions. We also fit the model parameters using astronomical observation: SNIa, H ( z ), BAO and Alcock-Paczynski test. We show that the model with diffuse DM-DE is consistent with the data.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1997

The Implicit and Explicit alpha-mu Schemes

Artificial numerical dissipation is an important issue in large Reynolds number computations. In such computations, the artificial dissipation inherent in traditional numerical schemes can overwhelm the physical dissipation and yield inaccurate results on meshes of practical size. In the present work, the space-time conservation element and solution element method is used to construct new and accurate numerical schemes such that artificial numerical dissipation will not overwhelm physical dissipation. Specifically, these schemes have the property that numerical dissipation vanishes when the physical viscosity goes to zero. These new schemes therefore accurately model the physical dissipation even when it is extremely small. The method of space-time conservation element and solution element, currently under development, is a nontraditional numerical method for solving conservation laws. The method is developed on the basis of local and global flux conservation in a space-time domain, in which space and time are treated in a unified manner. Explicit solvers for model and fluid dynamic conservation laws have previously been investigated. In this paper, we introduce a new concept in the design of implicit schemes, and use it to construct two highly accurate solvers for a convection-diffusion equation. The two schemes become identical in the pure convection case, and in the pure diffusion case. The implicit schemes are applicable over the whole Reynolds number range, from purely diffusive equations to purely inviscid (convective) equations. The stability and consistency of the schemes are analyzed, and some numerical results are presented. It is shown that, in the inviscid case, the new schemes become explicit and their amplification factors are identical to those of the Leapfrog scheme. On the other hand, in the pure diffusion case, their principal amplification factor becomes the amplification factor of the Crank-Nicolson scheme. We also construct an explicit solver with the treatment of diffusion being based on that in the implicit solvers. The explicit solver has only a CFL stability limitation on the Courant number, yet it retains the second-order spatial accuracy of the implicit schemes.

Chang, Sin-Chung↗
NASA NTRS2004

Balanced Central Schemes for the Shallow Water Equations on Unstructured Grids

We present a two-dimensional, well-balanced, central-upwind scheme for approximating solutions of the shallow water equations in the presence of a stationary bottom topography on triangular meshes. Our starting point is the recent central scheme of Kurganov and Petrova (KP) for approximating solutions of conservation laws on triangular meshes. In order to extend this scheme from systems of conservation laws to systems of balance laws one has to find an appropriate discretization of the source terms. We first show that for general triangulations there is no discretization of the source terms that corresponds to a well-balanced form of the KP scheme. We then derive a new variant of a central scheme that can be balanced on triangular meshes. We note in passing that it is straightforward to extend the KP scheme to general unstructured conformal meshes. This extension allows us to recover our previous well-balanced scheme on Cartesian grids. We conclude with several simulations, verifying the second-order accuracy of our scheme as well as its well-balanced properties.

Bryson, Steve↗
DOE OSTI2016

Variational formulations of guiding-center Vlasov-Maxwell theory

The variational formulations of guiding-center Vlasov-Maxwell theory based on Lagrange, Euler, and Euler-Poincaré variational principles are presented. Each variational principle yields a different approach to deriving guiding-center polarization and magnetization effects into the guiding-center Maxwell equations. The conservation laws of energy, momentum, and angular momentum are also derived by Noether method, where the guiding-center stress tensor is now shown to be explicitly symmetric.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The tensor-train stochastic finite volume method for uncertainty quantification

The stochastic finite volume method offers an efficient one-pass approach for assessing uncertainty in hyperbolic conservation laws. Still, it struggles with the curse of dimensionality when dealing with multiple stochastic variables. Here, we introduce the stochastic finite volume method within the tensor-train framework to counteract this limitation. This integration, however, comes with its own set of difficulties, mainly due to the propensity for shock formation in hyperbolic systems. To overcome these issues, we have developed a tensor-train-adapted stochastic finite volume method that employs a global WENO reconstruction, making it suitable for such complex systems. This approach represents the first step in designing tensor-train techniques for hyperbolic systems and conservation laws involving shocks.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Dissipation-preserving spectral element method for damped seismic wave equations

This article describes the extension of the conformal symplectic method to solve the damped acoustic wave equation and the elastic wave equations in the framework of the spectral element method. The conformal symplectic method is a variation of conventional symplectic methods to treat non-conservative time evolution problems, which has superior behaviors in long-time stability and dissipation preservation. To reveal the intrinsic dissipative properties of the model equations, we first reformulate the original systems in their equivalent conformal multi-symplectic structures and derive the corresponding conformal symplectic conservation laws. We thereafter separate each system into a conservative Hamiltonian system and a purely dissipative ordinary differential equation system. Based on the splitting methodology, we solve the two subsystems respectively. The dissipative one is cheaply solved by its analytic solution. While for the conservative system, we combine a fourth-order symplectic Nyström method in time and the spectral element method in space to cover the circumstances in realistic geological structures involving complex free-surface topography. The Strang composition method is adopted thereby to concatenate the corresponding two parts of solutions and generate the completed conformal symplectic method. A relative larger Courant number than that of the traditional Newmark scheme is found in the numerical experiments in conjunction with a spatial sampling of approximately 5 points per wavelength. A benchmark test for the damped acoustic wave equation validates the effectiveness of our proposed method in precisely capturing dissipation rate. The classical Lamb problem is used to demonstrate the ability of modeling Rayleigh wave in elastic wave propagation. More comprehensive numerical experiments are presented to investigate the long-time simulation, low dispersion and energy conservation properties of the conformal symplectic methods in both the attenuating homogeneous and heterogeneous media.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1994

Steady Secondary Momentum and Enthalpy Streaming in the Pulse Tube Refrigerator

Our study investigates the steady secondary streaming that occurs in the pulse tube refrigerator through application of the conservation laws and support by flow visualization experiments. A series expansion of the mass, momentum and energy conservation laws for an axisymmetric system yields a set of linear differential equations. The series solution is expanded in the small parameter epsilon = 1/Str where Str is the Strouhal number. The O(1) basic flow problem is solved for the case of a compressible ideal gas oscillating within a cylindrical tube. The phase and amplitude of the boundary conditions on axial velocity are treated as independent parameters. The basic solution for the gas domain is solved for the case of temperature and heat flux coupling to the tube wall. The temperature coupled solution shows how the Fourier number of the tube wall affects the gas temperature, phase angle, and time-averaged enthalpy flow. The basic solution is also examined for the time-averaged Reynolds stresses which are shown to be the driving force that produces steady secondary momentum streaming. Calculated solutions to the steady momentum streaming are shown to be in qualitative agreement with flow visualization experiments.

Lee, J. M.↗
DOE OSTI2016

On the quantum Landau collision operator and electron collisions in dense plasmas

The quantum Landau collision operator, which extends the widely used Landau/Fokker-Planck collision operator to include quantum statistical effects, is discussed. The quantum extension can serve as a reference model for including electron collisions in non-equilibrium dense plasmas, in which the quantum nature of electrons cannot be neglected. In this paper, the properties of the Landau collision operator that have been useful in traditional plasma kinetic theory and plasma transport theory are extended to the quantum case. We outline basic properties in connection with the conservation laws, the H-theorem, and the global and local equilibrium distributions. We discuss the Fokker-Planck form of the operator in terms of three potentials that extend the usual two Rosenbluth potentials. We establish practical closed-form expressions for these potentials under local thermal equilibrium conditions in terms of Fermi-Dirac and Bose-Einstein integrals. We study the properties of linearized quantum Landau operator, and extend two popular approximations used in plasma physics to include collisions in kinetic simulations. We apply the quantum Landau operator to the classic test-particle problem to illustrate the physical effects embodied in the quantum extension. We present useful closed-form expressions for the electron-ion momentum and energy transfer rates. Throughout the paper, similarities and differences between the quantum and classical Landau collision operators are emphasized.

Daligault, Jérôme↗
DOE OSTI2021

The importance of Noether's theorem for the way we think today

In the summer of 1918, German mathematician Emmy Noether published the two-part theorem that now bears her name, establishing profound connections between symmetries and conservation laws and between symmetries and interactions. These insights exert a pervasive influence in physics; they give meaning to conservation laws that elevates them beyond useful empirical rules, and they underlie all our theories of the fundamental interactions. In Noether's later career, her papers, lectures, and personal interactions with students and colleagues drove the development of abstract algebra, establishing her in the pantheon of twentieth-century mathematicians.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1989

An Entropy-Based Approach to Nonlinear Stability

Many numerical methods used in computational fluid dynamics (CFD) incorporate an artificial dissipation term to suppress spurious oscillations and control nonlinear instabilities. The same effect can be accomplished by using upwind techniques, sometimes augmented with limiters to form Total Variation Diminishing (TVD) schemes. An analysis based on numerical satisfaction of the second law of thermodynamics allows many such methods to be compared and improved upon. A nonlinear stability proof is given for discrete scalar equations arising from a conservation law. Solutions to such equations are bounded in the L sub 2 norm if the second law of thermodynamics is satisfied in a global sense over a periodic domain. It is conjectured that an analogous statement is true for discrete equations arising from systems of conservation laws. Analysis and numerical experiments suggest that a more restrictive condition, a positive entropy production rate in each cell, is sufficient to exclude unphysical phenomena such as oscillations and expansion shocks. Construction of schemes which satisfy this condition is demonstrated for linear and nonlinear wave equations and for the one-dimensional Euler equations.

Merriam, Marshal L.↗
DOE OSTI2024

Cycle equivalence classes, orthogonal Weingarten calculus, and the mean field theory of memristive systems

Abstract It has been recently noted that for a class of dynamical systems with explicit conservation laws represented via projector operators, the dynamics can be understood in terms of lower dimensional equations. This is the case, for instance, of memristive circuits. Memristive systems are important classes of devices with wide-ranging applications in electronic circuits, artificial neural networks, and memory storage. We show that such mean-field theories can emerge from averages over the group of orthogonal matrices, interpreted as cycle-preserving transformations applied to the projector operator describing Kirchhoff’s laws. Our results provide insights into the fundamental principles underlying the behavior of resistive and memristive circuits and highlight the importance of conservation laws for their mean-field theories. In addition, we argue that our results shed light on the nature of the critical avalanches observed in quasi-two-dimensional nanowires as boundary phenomena.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2002

Interchange Method in Compressible Magnetized Couette Flow: Magnetorotational and Magnetoconvective Instabilities

We obtain the general forms of the axisymmetric stability criteria in a magnetized compressible Couette flow using an energy variational principle, the so-called interchange or Chandrasekhar s met hod, which we applied successfully in the incompressible case. This formulation accounts for the simultaneous presence of gravity, rotation, a toroidal magnetic field, a weak axial magnetic field, entropy gradients, and density gradients in the initial equilibrium state. The power of the method lies in its simplicity which allows us to derive extremely compact and physically clear expressions for the relevant stability criteria despite the inclusion of so many physical effects. In the implementation of the method, all the applicable conservation laws are explicitly taken into account during the variations of a quantity with dimensions of energy which we call the free energy function. As in the incompressible case, the presence of an axial field invalidates the conservation laws of angular momentum and azimuthal magnetic flux and introduces instead isorotation and axial current conservation along field lines. Our results are therefore markedly different depending on whether an axial magnetic field is present, and generalize in two simple expressions all previously known, partial stability criteria for the appearance of magnetorotational instability. Furthermore, the coupling between magnetic tension and buoyancy and its influence to the dynamics of nonhomoentropic magnetized flows becomes quite clear from our results. In the limits of plane-parallel atmospheres and homoentropic flows, our formulation easily recovers the stability criteria for suppression of convective and Parker instabilities, as well as some related special cases studied over 40 years ago by Newcomb and Tserkovnikov via laborious variational techniques.

Christodoulou, Dimitris M.↗
NASA NTRS2020

Entropy Stable Nonconforming Discretizations with the Summation-By-Parts Property for Curvilinear Coordinates

The entropy conservative/stable algorithm of Friedrichet al.(2018) for hyperbolic conservation laws on nonconforming p-refined/coarsened Cartesian grids, is extended to curvilinear grids for the compressible Euler equations. The primary focus is on constructing appropriate coupling procedures across the curvilinear nonconforming interfaces. A simple and flexible approach is proposed that uses interpolation operators from one element to the other. On the element faces,the analytic metrics are used to construct coupling terms, while metric terms in the volume are approximated to satisfy a discretization of the geometric conservation laws. The resulting scheme is entropy conservative/stable, elementwise conservative, and freestream preserving. The accuracy and stability properties of the resulting numerical algorithm are shown to be comparable to those ofthe original conforming scheme (∼p+ 1 convergence) in the context of the isentropic Euler vortex and the inviscid Taylor–Green vortex problems on manufactured high order grids.

David Del Rey Fernandez↗
DOE OSTI2019

Colloquium: A Century of Noether's Theorem

In the summer of 1918, Emmy Noether published the theorem that now bears her name, establishing a profound two-way connection between symmetries and conservation laws. The influence of this insight is pervasive in physics; it underlies all of our theories of the fundamental interactions and gives meaning to conservation laws that elevates them beyond useful empirical rules. Noether's papers, lectures, and personal interactions with students and colleagues drove the development of abstract algebra, establishing her in the pantheon of twentieth-century mathematicians. This essay traces her path from Erlangen through G\"ottingen to a brief but happy exile at Bryn Mawr College in Pennsylvania, illustrating the importance of "Noether's Theorem" for the way we think today. The text draws on a colloquium presented at Fermilab on 15 August 2018.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2016

Hydromechanical Modeling of Unsaturated Flow in Double Porosity Media

Geomaterials with aggregated structure or containing fissures often exhibit a bimodal pore size distribution that can be viewed as two coexisting pore regions of different scales. The double-porosity concept enables continuum modeling of such materials by considering two interacting pore scales satisfying relevant conservation laws. This paper develops a thermodynamically consistent framework for hydromechanical modeling of unsaturated flow in double-porosity media. With an explicit treatment of the two pore scales, conservation laws are formulated incorporating an effective stress tensor that is energy-conjugate to the rate of deformation tensor of the solid matrix. A constitutive framework is developed on the basis of energy-conjugate pairs identified in the first law of thermodynamics, which is then incorporated into a three-field mixed finite-element formulation for double-porosity media. In conclusion, numerical simulations of laboratory- and field-scale problems are presented to demonstrate the impact of double porosity on the resulting hydromechanical responses.

58 GEOSCIENCES↗