SEARCH · Engineering Papers

Results for “Equation”

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 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Differential equations for cosmological correlators

Cosmological fluctuations retain a memory of the physics that generated them in their spatial correlations. The strength of correlations varies smoothly as a function of external kinematics, which is encoded in differential equations satisfied by cosmological correlation functions. In this work, we provide a broader perspective on the origin and structure of these differential equations. As a concrete example, we study conformally coupled scalar fields in a power-law cosmology. The wavefunction coefficients in this model have integral representations, with the integrands being the product of the corresponding flat-space results and “twist factors” that depend on the cosmological evolution. Similar twisted integrals arise for loop amplitudes in dimensional regularization, and their recent study has led to the discovery of rich mathematical structures and powerful new tools for computing multi-loop Feynman integrals in quantum field theory. The integrals of interest in cosmology are also part of a finite-dimensional basis of master integrals, which satisfy a system of first-order differential equations. We develop a formalism to derive these differential equations for arbitrary tree graphs. The results can be represented in graphical form by associating the singularities of the differential equations with a set of graph tubings. Upon differentiation, these tubings grow in a local and predictive fashion. In fact, a few remarkably simple rules allow us to predict — by hand — the equations for all tree graphs. While the rules of this “kinematic flow” are defined purely in terms of data on the boundary of the spacetime, they reflect the physics of bulk time evolution. We also study the analogous structures in tr ϕ 3 theory, and see some glimpses of hidden structure in the sum over planar graphs. This suggests that there is an autonomous combinatorial or geometric construction from which cosmological correlations, and the associated spacetime, emerge.

Cosmological models↗

DOE OSTI2025

Validating corrosion models: A comparison of governing equations

Experimental validation of Finite Element Method (FEM) models varying electrochemical governing equations, inclusion of chemical reactions, and time on the resultant damage profile for two galvanic couples is explored. Two anode materials (Magnesium AZ31 and Carbon Steel) in contact with a cathode (Stainless Steel 304 L) were modeled in/exposed to NaCl (1 and 0.1 M respectively for the anode materials) for up to one week. The physics approach, inclusion of chemical reactions, and the boundary conditions required to accurately represent the damage profile in FEM models depended on the galvanic couple materials and, ultimately, the corrosion rate. For high rates of corrosion (i.e., magnesium anode), the Nernst-Planck equation with Electroneutrality was sufficient to describe the damage, while, for low rates of corrosion (i.e., carbon steel anode), the Laplace equation was sufficient. In all cases, the most complete governing equation (Nernst-Planck-Poisson Equation) was not necessary to accurately describe the damage. Precipitation reactions in solution also played a critical role in the predicted damage profile, especially for high corrosion rate systems. Finally, for short time periods (< 6 h), the choice of governing equations does not significantly influence damage profile results. Overall, the choice of physics to reduce error in simulations relies on the boundary conditions, geometry, conductivity of the solution, electrochemical potential differences, and time of exposure. The above results are discussed with regard to accuracy and computational savings.

Carbon steel↗

DOE OSTI2025

Euler equations and the Sod shock tube problem

The Euler equations are a subset of the magnetohydrodynamic (MHD) equations in the infinitely collisional, unmagnetized limit. MHD modeling is central to many areas of plasma physics, ranging from low-temperature glow discharges to inertial confinement fusion. An important aspect of the Euler equations is their ability to describe states with discontinuities, such as shock waves. A standard benchmark test for numerical implementation of the Euler equations is the Sod shock tube. In this test, the system is initialized at rest with a pressure and density discontinuity, which results in a shock wave traveling into the low-pressure region and a rarefaction wave traveling into the high-pressure region. Starting with the presentation of the Euler equations, a numerical algorithm is presented here to solve these equations in one dimension. This is followed by an overview of the Sod shock tube problem that includes the precise initial setup and the analytic solution. Finally, the analytic solution is compared with results from numerical simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1971

Non-partial wave treatment of reactive and non-reactive scattering Coupled integral equation formalism.

Coupled integral equations are derived for the full scattering amplitudes for both reactive and nonreactive channels. The equations do not involve any partial wave expansion and are obtained using channel operators for reactive and nonreactive collisions. These coupled integral equations are similar in nature to equations derived for purely nonreactive collisions of structureless particles. Using numerical quadrature techniques, these equations may be reduced to simultaneous algebraic equations which may then be solved.

Hayes, E. F.↗

NASA NTRS1976

A quasi-linear kinetic equation for cosmic rays in the interplanetary medium

A kinetic equation for interplanetary cosmic rays is set up with the aid of weak-plasma-turbulence theory for an idealized radially symmetric model of the interplanetary magnetic field. As a starting point, this treatment invokes the Vlasov equation instead of the traditional Fokker-Planck equation. Quasi-linear theory is applied to obtain a momentum diffusion equation for the heliocentric frame of reference which describes the interaction of cosmic rays with convecting magnetic irregularities in the solar-wind plasma. Under restricted conditions, the well-known equation of solar modulation can be obtained from this kinetic equation.

Luhmann, J. G.↗

NASA NTRS1976

A tensor formulation of the equation of transfer for spherically symmetric flows

A tensor formulation of the equation of radiative transfer is derived in a seven-dimensional Riemannian space such that the resulting equation constitutes a divergence in any coordinate system. After being transformed to a spherically symmetric comoving coordinate system, the transfer equation contains partial derivatives in angle and frequency, as well as optical depth due to the effects of aberration and the Doppler shift. However, by virtue of the divergence form of this equation, the divergence theorem may be applied to yield a numerical differencing scheme which is expected to be stable and to conserve luminosity. It is shown that the equation of transfer derived by this method in a Lagrangian coordinate system may be reduced to that given by Castor (1972), although it is, of course, desirable to leave the equation in divergence form.

Haisch, B. M.↗

NASA NTRS1977

Symbolic computation of recurrence equations for the Chebyshev series solution of linear ODE's

If a linear ordinary differential equation with polynomial coefficients is converted into integrated form then the formal substitution of a Chebyshev series leads to recurrence equations defining the Chebyshev coefficients of the solution function. An explicit formula is presented for the polynomial coefficients of the integrated form in terms of the polynomial coefficients of the differential form. The symmetries arising from multiplication and integration of Chebyshev polynomials are exploited in deriving a general recurrence equation from which can be derived all of the linear equations defining the Chebyshev coefficients. Procedures for deriving the general recurrence equation are specified in a precise algorithmic notation suitable for translation into any of the languages for symbolic computation. The method is algebraic and it can therefore be applied to differential equations containing indeterminates.

Geddes, K. O.↗

NASA NTRS1977

Derivation of an integral equation for transonic flows

It is noted that the nonlinear partial differential equation for the perturbation velocity potential and boundary conditions describing steady inviscid compressible transonic flow past a thin two-dimensional airfoil can be transformed into a singular integrodifferential equation and that differentiation of the latter yields an integral equation. Two forms of this integral equation currently exist: one for the singularity that is enclosed in an infinitely long strip of vanishing thickness and the other for the singularity that is enclosed in a vanishing circle. In the present article, a more general integral equation is derived by enclosing the singularity in a vanishing rectangular cavity of arbitrary aspect ratio. The two existing forms of this equation are deduced as special cases distinguished by the respective values for the aspect ratio (infinity for the first form and unity for the second).

Ogana, W.↗

NASA NTRS1976

Scaled equation of state parameters for gases in the critical region

In the light of recent theoretical developments, the paper presents an accurate characterization of anomalous thermodynamic behavior of xenon, helium 4, helium 3, carbon dioxide, steam and oxygen in the critical region. This behavior is associated with long range fluctuations in the system and the physical properties depend primarily on a single variable, namely, the correlation length. A description of the thermodynamic behavior of fluids in terms of scaling laws is formulated, and the two successfully used scaled equations of state (NBS equation and Linear Model parametric equation) are compared. Methods for fitting both equations to experimental equation of state data are developed and formulated, and the optimum fit for each of the two scaled equations of the above gases are presented and the results are compared. By extending the experimental data for the above one-component fluids to partially miscible binary liquids, superfluid liquid helium, ferromagnets and solids exhibiting order-disorder transitions, the principle of universality is concluded. Finally by using this principle, the critical regions for nine additional fluids are described.

Sengers, J. M. H. L.↗

NASA NTRS1979

Algebraic methods for the solution of some linear matrix equations

The characterization of polynomials whose zeros lie in certain algebraic domains (and the unification of the ideas of Hermite and Lyapunov) is the basis for developing finite algorithms for the solution of linear matrix equations. Particular attention is given to equations PA + A'P = Q (the Lyapunov equation) and P - A'PA = Q the (discrete Lyapunov equation). The Lyapunov equation appears in several areas of control theory such as stability theory, optimal control (evaluation of quadratic integrals), stochastic control (evaluation of covariance matrices) and in the solution of the algebraic Riccati equation using Newton's method.

Djaferis, T. E.↗

NASA NTRS1978

Stochastic partial differential equations in turbulence related problems

The theory of stochastic partial differential equations (PDEs) and problems relating to turbulence are discussed by employing the theories of Brownian motion and diffusion in infinite dimensions, functional differential equations, and functional integration. Relevant results in probablistic analysis, especially Gaussian measures in function spaces and the theory of stochastic PDEs of Ito type, are taken into account. Linear stochastic PDEs are analyzed through linearized Navier-Stokes equations with a random forcing. Stochastic equations for waves in random media as well as model equations in turbulent transport theory are considered. Markovian models in fully developed turbulence are discussed from a stochastic equation viewpoint.

Chow, P.-L.↗

NASA NTRS1981

Surrogate-equation technique for simulation of steady inviscid flow

A numerical procedure for the iterative solution of inviscid flow problems is described, and its utility for the calculation of steady subsonic and transonic flow fields is demonstrated. Application of the surrogate equation technique defined herein allows the formulation of stable, fully conservative, type dependent finite difference equations for use in obtaining numerical solutions to systems of first order partial differential equations, such as the steady state Euler equations. Steady, two dimensional solutions to the Euler equations for both subsonic, rotational flow and supersonic flow and to the small disturbance equations for transonic flow are presented.

Johnson, G. M.↗

NASA NTRS1981

On various treatments of potential equations at shocks

The potential equation that governs transonic inviscid flowfields is elliptic in subsonic regions and hyperbolic in supersonic regions. The transition from supersonic to subsonic flow may occur either continuously or across a surface of discontinuity known as a shock. Because of the changing type of governing equation and the presence of the surfaces of discontinuity, special treatments of the governing equations are needed at shocks to preserve mass conservation and to ensure satisfaction of the governing equation. In the so-called shock-capturing techniques, shocks are captured by smearing the discontinuity over several mesh cells in the computational domain. Most existing conservative schemes conserve mass flux throughout the flowfield, including shock boundaries; however, they also introduce zero-order errors in the approximation to the governing equation. It is shown that the zero-order errors, which do not diminish as the mesh spacing approaches zero, can cause discrepancies in the prediction of shock strength and location. These parameters are extremely sensitive features of transonic inviscid flowfield calculations. It is also shown that the zero-order errors can be avoided by developing higher-order schemes that properly model the governing equation at shocks.

Chen, L. T.↗

NASA NTRS1982

A composite velocity procedure for the compressible Navier-Stokes equations

A new boundary-layer relaxation procedure is presented. In the spirit of the theory of matched asymptotic expansions, a multiplicative composite of the appropriate velocity representations for the inviscid and viscous regions is prescribed. The resulting equations are structured so that far from the surface of the body the momentum equations lead to the Bernoulli relation for the pressure, while the continuity equation reduces to the familiar compressible potential equation. Close to the body surface, the governing equations and solution techniques are characteristic of those describing interacting boundary-layers; although, the full Navier-Stokes equations are considered here. Laminar flow calculations for the subsonic flow over an axisymmetric boattail simulator geometry are presented for a variety of Reynolds and Mach numbers. A strongly implicit solution method is applied for the coupled velocity components.

Khosla, P. K.↗

NASA NTRS1983

Numerical computation of transonic flow governed by the full-potential equation

Numerical solution techniques for solving transonic flow fields governed by the full potential equation are discussed. In a general sense relaxation schemes suitable for the numerical solution of elliptic partial differential equations are presented and discussed with emphasis on transonic flow applications. The presentation can be divided into two general categories: An introductory treatment of the basic concepts associated with the numerical solution of elliptic partial differential equations and a more advanced treatment of current procedures used to solve the full potential equation for transonic flow fields. The introductory material is presented for completeness and includes a brief introduction (Chapter 1), governing equations (Chapter 2), classical relaxation schemes (Chapter 3), and early concepts regarding transonic full potential equation algorithms (Chapter 4).

Holst, T. L.↗

NASA NTRS1984

Fundamental Equations of the LIF Hydroxyl Measurements, Appendix a

The equation relating OH concentration to the observed detector count rate for a given laser flux is defined. That equation, which will be referred to henceforth as the central equation for absolute sensitivity, will be cast in terms of standard spectroscopic quantities such as oscillator strengths, Hoenl-London factors, rotational and vibrational quantum numbers, etc., such that it may be tested directly under laboratory conditions. The working expressions used by each of the OH laser induced fluorescence (LIF) groups will be related to this central equation, thereby establishing a common language for the comparison of experimental results and the definition of LIF system performance. The relevant structural details of the OH radical via a spectroscopic description of the energy levels and transitions employed in the LIF method are detailed. The central equation relating observed count rate to absolute concentration will be derived, and that equation will be related to the observed quantities cited by each of the experimental groups engaged in the measurement of OH in the Earth's atmosphere.

Source record↗

NASA NTRS1985

Numerical solutions of 3-dimensional Navier-Stokes equations for closed bluff-bodies

The Navier-Stokes equations are solved numerically. These equations are unsteady, compressible, viscous, and three-dimensional without neglecting any terms. The time dependency of the governing equations allows the solution to progress naturally for an arbitrary initial guess to an asymptotic steady state, if one exists. The equations are transformed from physical coordinates to the computational coordinates, allowing the solution of the governing equations in a rectangular parallelepiped domain. The equations are solved by the MacCormack time-split technique which is vectorized and programmed to run on the CDc VPS 32 computer. The codes are written in 32-bit (half word) FORTRAN, which provides an approximate factor of two decreasing in computational time and doubles the memory size compared to the 54-bit word size.

Abolhassani, J. S.↗

NASA NTRS1985

Stable boundary conditions and difference schemes for Navier-Stokes equations

The Navier-Stokes equations can be viewed as an incompletely elliptic perturbation of the Euler equations. By using the entropy function for the Euler equations as a measure of energy for the Navier-Stokes equations, it was possible to obtain nonlinear energy estimates for the mixed initial boundary value problem. These estimates are used to derive boundary conditions which guarantee L2 boundedness even when the Reynolds number tends to infinity. Finally, a new difference scheme for modelling the Navier-Stokes equations in multidimensions for which it is possible to obtain discrete energy estimates exactly analogous to those we obtained for the differential equation was proposed.

Dutt, P.↗