SEARCH · Engineering Papers

Results for “Conservation equations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,351 records · Page 75

All sources NASA NTRS DOE OSTI

NASA NTRS1972

Interplanetary shock normals

A scheme has been developed to improve the estimate of an interplanetary shock normal by using both magnetic field and plasma data from a single spacecraft. Calculation of the basic shock model employs a subset of the eight magnetohydrodynamic conservation relations for a shock in an isotropic medium. This subset consists of six equations that are devoid of pressure and temperature terms. A sigma-weighted least squares loss function technique is used to best fit the overdetermination equations with respect to the eleven parameters of the system. This procedure usually yields a normal much more accurate than the one obtained using magnetic field values alone. An example is given of the use of the technique on the data for the Pioneer 7 shock of 29 August 1966.

Lepping, R. P.↗

NASA NTRS1975

Numerical solution of the Navier-Stokes equations with topography

A finite difference scheme for solving the equations of fluid motion in a generalized coordinate system has been constructed. The scheme conserves mass and all the first integral moments of the motion. The scheme also advectively 'almost conserves' second moments, in that the magnitude of implicit numerical smoothing is typically about an order smaller than explicit viscosity and diffusion. Calculations with the model support the theoretical conjecture that the difference scheme is stable whenever the analogous Cartesian scheme is stable. The scheme has been used to calculate dry atmospheric convection due to differential heating between top and bottom of mountainous terrain. The general small-scale characteristics of mountain up-slope winds have been simulated. In addition, the results have demonstrated the crucial role played by the eddy diffusivities and the environmental stability, in determining both the quantitative and the qualitative features of the circulation.

Gal-Chen, T.↗

NASA NTRS1975

Gravitational theories with nonzero divergence of the energy-momentum tensor

Assuming that the divergence of the energy-momentum tensor is nonzero leads to a class of theories with consistent field equations and gauge conditions as well as compatibility with the Newtonian limit of the conservation laws. Both the Einstein and the Brans-Dicke theories are used as models, but the extension to other viable theories such as vector-metric and two-metric theories is possible. One particularly interesting theory emerges that agrees with the ordinary Brans-Dicke theory except for the post-Newtonian parameter zeta sub 2, which predicts nonconservation of total momentum. Unfortunately, no accurate experimental limits for this parameter are known. It thus remains for future experiments in lunar-laser ranging to test this theory.

Smalley, L. L.↗

NASA NTRS1981

An explicit mixed numerical method for mesoscale model

A mixed numerical method has been developed for mesoscale models. The technique consists of a forward difference scheme for time tendency terms, an upstream scheme for advective terms, and a central scheme for the other terms in a physical system. It is shown that the mixed method is conditionally stable and highly accurate for approximating the system of either shallow-water equations in one dimension or primitive equations in three dimensions. Since the technique is explicit and two time level, it conserves computer and programming resources.

Hsu, H.-M.↗

NASA NTRS1987

The terminal area simulation system. Volume 1: Theoretical formulation

A three-dimensional numerical cloud model was developed for the general purpose of studying convective phenomena. The model utilizes a time splitting integration procedure in the numerical solution of the compressible nonhydrostatic primitive equations. Turbulence closure is achieved by a conventional first-order diagnostic approximation. Open lateral boundaries are incorporated which minimize wave reflection and which do not induce domain-wide mass trends. Microphysical processes are governed by prognostic equations for potential temperature water vapor, cloud droplets, ice crystals, rain, snow, and hail. Microphysical interactions are computed by numerous Orville-type parameterizations. A diagnostic surface boundary layer is parameterized assuming Monin-Obukhov similarity theory. The governing equation set is approximated on a staggered three-dimensional grid with quadratic-conservative central space differencing. Time differencing is approximated by the second-order Adams-Bashforth method. The vertical grid spacing may be either linear or stretched. The model domain may translate along with a convective cell, even at variable speeds.

Proctor, F. H.↗

NASA NTRS1988

Analysis of the spectral vanishing viscosity method for periodic conservation laws

The convergence of the spectral vanishing method for both the spectral and pseudospectral discretizations of the inviscid Burgers' equation is analyzed. It is proven that this kind of vanishing viscosity is responsible for a spectral decay of those Fourier coefficients located toward the end of the computed spectrum; consequently, the discretization error is shown to be spectrally small independent of whether the underlying solution is smooth or not. This in turn implies that the numerical solution remains uniformly bounded and convergence follows by compensated compactness arguments.

Maday, Yvon↗

NASA NTRS1990

Symmetrization of conservation laws with entropy for high-temperature hypersonic computations

Results of Hughes, France, and Mallet are generalized to conservation law systems taking into account high-temperature effects. Symmetric forms of different equation sets are derived in terms of entropy variables. First, the case of a general divariant gas is studied; it can be specialized to the usual Navier-Stokes equations, as well as to situations where the gas is vibrationally excited, and undergoes equilibrium chemical reactions. The case of gas in thermochemical nonequilibrium is considered next. Transport phenomena, and in particular mass diffusion, are examined in the framework of symmetric advective-diffusive systems.

Chalot, F.↗

NASA NTRS1999

Bidirectional Reflectance of Flat, Optically Thick Particulate Layers: An Efficient Radiative Transfer Solution and Applications to Snow and Soil Surfaces

We describe a simple and highly efficient and accurate radiative transfer technique for computing bidirectional reflectance of a macroscopically flat scattering layer composed of nonabsorbing or weakly absorbing, arbitrarily shaped, randomly oriented and randomly distributed particles. The layer is assumed to be homogeneous and optically semi-infinite, and the bidirectional reflection function (BRF) is found by a simple iterative solution of the Ambartsumian's nonlinear integral equation. As an exact Solution of the radiative transfer equation, the reflection function thus obtained fully obeys the fundamental physical laws of energy conservation and reciprocity. Since this technique bypasses the computation of the internal radiation field, it is by far the fastest numerical approach available and can be used as an ideal input for Monte Carlo procedures calculating BRFs of scattering layers with macroscopically rough 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 at http://ww,,v.giss.nasa.gov/-crmim/brf.html and can be applied to a wide range of remote sensing, engineering, and biophysical problems. We also examine the potential effect of ice crystal shape on the bidirectional reflectance of flat snow surfaces and the applicability of the Henyey-Greenstein phase function and the 6-Eddington approximation in calculations for soil surfaces.

Mishchenko, Michael I.↗

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↗

NASA NTRS2014

On Formulations of Discontinuous Galerkin and Related Methods for Conservation Laws

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

Huynh, H. T.↗

NASA NTRS1988

A diagonally inverted LU implicit multigrid scheme

A new Diagonally Inverted LU Implicit scheme is developed within the framework of the multigrid method for the three-dimensional unsteady Euler equations. The matrix systems that are to be inverted in the LU scheme are treated by local diagonalizing transformations that decouple them into systems of scalar equations. Unlike the Diagonalized ADI method, the time accuracy of the LU scheme is not reduced since the diagonalizing procedure does not destroy time conservation. Even more importantly, this diagonalization significantly reduces the computational effort required to solve the LU approximation and therefore transforms it into a more efficient method of numerically solving the three-dimensional Euler equations.

Yokota, J. W.↗

NASA NTRS1995

On the Conservative Interface Treatment for Multi-Block Viscous Flow Computations

A pressure-based multi-block computational method is developed for solving the incompressible Navier-Stokes equations in general curvilinear grid systems. The scheme is based on the semi-implicit type flow solver with the staggered grid. Issues concerning the mass and momentum flux treatments at the discontinuous grid interface are addressed. Systematic numerical experiments for different interface treatments involving (1) straightforward interpolation, (2) globally conservative scheme, and (3) locally conservative scheme have been conducted. It is demonstrated that mass conservation has to be maintained locally, at the grid interface, with accuracy compatible with that of the scheme used in interior domain. Direct interpolation or globally conservative interface treatment of mass flux can not yield solutions with desirable accuracy.

Liu, J.↗

NASA NTRS1993

Analysis of the Harrier forebody/inlet design using computational techniques

Under the support of this Cooperative Agreement, computations of transonic flow past the complex forebody/inlet configuration of the AV-8B Harrier II have been performed. The actual aircraft configuration was measured and its surface and surrounding domain were defined using computational structured grids. The thin-layer Navier-Stokes equations were used to model the flow along with the Chimera embedded multi-grid technique. A fully conservative, alternating direction implicit (ADI), approximately-factored, partially flux-split algorithm was employed to perform the computation. An existing code was altered to conform with the needs of the study, and some special engine face boundary conditions were developed. The algorithm incorporated the Chimera technique and an algebraic turbulence model in order to deal with the embedded multi-grids and viscous governing equations. Comparison with experimental data has yielded good agreement for the simplifications incorporated into the analysis. The aim of the present research was to provide a methodology for the numerical solution of complex, combined external/internal flows. This is the first time-dependent Navier-Stokes solution for a geometry in which the fuselage and inlet share a wall. The results indicate the methodology used here is a viable tool for transonic aircraft modeling.

Chow, Chuen-Yen↗

DOE OSTI2022

The kinetic origin of the fluid helicity—A symmetry in the kinetic phase space

Helicity, a topological degree that measures the winding and linking of vortex lines, is preserved by ideal (barotropic) fluid dynamics. In the context of the Hamiltonian description, the helicity is a Casimir invariant characterizing a foliation of the associated Poisson manifold. Casimir invariants are special invariants that depend on the Poisson bracket, not on the particular choice of the Hamiltonian. The total mass (or particle number) is another Casimir invariant, whose invariance guarantees the mass (particle) conservation (independent of any specific choice of the Hamiltonian). In a kinetic description (e.g., that of the Vlasov equation), the helicity is no longer an invariant (although the total mass remains a Casimir of the Vlasov’s Poisson algebra). The implication is that some “kinetic effect” can violate the constancy of the helicity. To elucidate how the helicity constraint emerges or submerges, we examine the fluid reduction of the Vlasov system; the fluid (macroscopic) system is a “sub-algebra” of the kinetic (microscopic) Vlasov system. In the Vlasov system, the helicity can be conserved if a special helicity symmetry condition holds. To put it another way, breaking helicity symmetry induces a change in the helicity. We delineate the geometrical meaning of helicity symmetry and show that for a special class of flows (the so-called epi-two-dimensional flows), the helicity symmetry is written as ∂γ = 0 for a coordinate γ of the configuration space.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1981

Optimal time splitting for two- and three-dimensional Navier-Stokes equations with mixed derivatives

A new explicit, time splitting algorithm has been developed for finite difference modelling of the full two and three-dimensional time-dependent, compressible, viscous Navier-Stokes equations of fluid mechanics. The scheme is optimal in the sense that the split operators achieve their maximum allowable time step, i.e., the corresponding Courant number. The algorithm allows a conservation-form formulation. Stability is proven analytically and verified numerically. In proving stability it was shown that all nine matrix coefficients of the Navier-Stokes equations are simultaneously symmetrizable by a similarity transformation. Two such transformations and their resulting symmetric matrix coefficients are presented explicitly.

Abarbanel, S.↗

NASA NTRS1992

Kinetic description of electron beams in the solar chromosphere

We formulate the relativistic Fokker-Plank equation for a beam of accelerated electrons interacting with a partially ionized plasma. In our derivation we conserved those terms contributing to velocity diffusion and found that this effect cannot be neglected a priori. We compute the terms accounting for elastic and inelastic collisions with neutral hydrogen and helium. Collisions with neutral hydrogen are found to be dominant throughout the chromosphere, except at the uppermost layers close to the transition region. As an application, we compute the loss of energy and momentum for a power-law beam impinging on the solar chromosphere, for a particular case in which the Fokker-Planck equation can be integrated analytically. We find that most of the beam energy is deposited in a relatively thin region of the chromosphere, a result which is largely insensitive to the theoretical method employed to compute the energy deposition rate.

Gomez, Daniel O.↗

DOE OSTI2023

Analytical closure to the spatially-filtered Euler equations for shock-dominated flows

To ensure numerical stability in the vicinity of shocks, a variety of methods have been used, including shock-capturing schemes such as weighted essentially non-oscillatory schemes, as well as the addition of artificial diffusivities to the governing equations. Centered finite difference schemes are often avoided near discontinuities due to the tendency for significant oscillations. However, such schemes have desirable conservation properties compared to many shock-capturing schemes. The objective of this work is to derive all necessary viscous/diffusion terms from first principles and then demonstrate the performance of these analytical terms within a centered differencing framework. The physical Euler equations are spatially-filtered with a Gaussian-like filter. Sub-filter scale (SFS) terms arise in the momentum and energy equations. Analytical closure is provided for each of them by leveraging the jump conditions for a shock. No SFS terms are present in the continuity or species equations. Here, this approach is tested for several problems involving shocks in one and two dimensions. Implemented within a centered difference code, the SFS terms perform well for a range of flow conditions without introducing excessive diffusion.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Electron kinetics in a high- Z plasmoid

The problem of the electron dynamics on a closed magnetic field line passing through a high-Z plasmoid is considered. The electron kinetic equation is integrated over bounce motion and pitch angle, reducing the independent variables to a single adiabatic invariant plus time. Integration of the full Landau self-collision operator is carried out exactly, resulting in a nonlinear integro-differential operator in the new invariant. Conservation laws and the H theorem of the integrated self-collision operator are proven. Numerical solutions of the integrated kinetic equation are obtained with a self-consistent quasineutral electric potential, given the initial condition of a cold plasmoid immersed in a hot ambient plasma. The fact that cold electrons are deeply trapped in a potential with a parabolic peak leads to exactly 3/4 the usual rate of collisional heating by the ambient plasma, independent of any other parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗