Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 235 records · Page 13

Generalized fractional approach to solving partial differential equations with arbitrary dispersion relations

The dynamics of physical systems are typically modeled by partial differential equations (PDEs). Because of the complexity in physical systems, PDE models rely on approximations or limiting cases. To model the full complexity of physical systems, it is necessary to use multiscale approaches where appropriate models are used at each level. Alternatively, complexity can be directly tackled by fractional integrodifferential equations that must be derived for each system. Here, we argue that complexity can be restored in PDEs by describing them from a fractional calculus perspective. Rather than deriving a fractional integrodifferential equation, we reinterpret the dispersion relation of the system by use of the Riesz definition, which contains the required information relating the energy and momentum space of the system and thus fully describes their dynamics. The approach is demonstrated by two examples: the Landau–Lifshitz equation in a 1D ferromagnetic chain and a modified KdV equation supporting surface gravity waves or Euler dispersion. The presented approach is applicable to fluids, soft matter, and solid-state matter and can be readily generalized to higher dimensions and more complex systems. While numerical calculations are needed to determine the fractional operator, the approach is analytical and can be utilized to determine analytical solutions and investigate nonlinear problems.

97 MATHEMATICS AND COMPUTING↗

Revisiting the gluon density from the Balitsky-Kovchegov equation with kinematical constraints and large 𝑥 terms

We perform analysis of the small 𝑥 nonlinear evolution equation formulated in momentum space supplemented by higher order terms. The equation is defined in the wide range of transverse momentum and longitudinal momentum fractions extending previous studies performed in Kutak and Kwiecinski [Screening effects in the ultrahigh-energy neutrino interactions, Eur. Phys. J. C 29, 521 (2003).] and in Kutak and Stasto [Unintegrated gluon distribution from modified BK equation, Eur. Phys. J. C 41, 343 (2005).]. The linear part of the equation is motivated by the renormalization group improved small 𝑥 approach, which accounts for resummation of higher orders, and includes collinear splitting functions and kinematical constraints. The solution to the equation is then used to perform the fit to deep inelastic scattering reduced cross section data.

Perturbative QCD↗

Fokker-Planck Equation Governing the Distribution of Walkers in Auxiliary-Field Quantum Monte Carlo

Auxiliary-field quantum Monte Carlo (AFQMC) is typically formulated as an open-ended random walk in an overcomplete space of Slater determinants, implemented through a Langevin equation. However, the explicit form of the underlying Fokker-Planck equation governing the walker population distribution has remained unknown. Here, in this Letter, we derive the Fokker-Planck equation for AFQMC and propose a novel numerical scheme to solve it. The solution of the Fokker-Planck equation reveals the wave function actually sampled by the AFQMC algorithm. Interestingly, we find that even when the exact ground state is used as a guiding wave function in constrained path AFQMC, contrary to the common assumption, the wave function sampled by AFQMC is not exact. Beyond clarifying several fundamental aspects of AFQMC, the availability of a Fokker-Planck equation formulation opens new avenues for systematically improving its accuracy, which we outline in this Letter.

Monte Carlo methods↗

Computer program for solving compressible nonsimilar-boundary-layer equations for laminar, transitional, or turbulent flows of a perfect gas

A computer program is described which solves the compressible laminar, transitional, or turbulent boundary-layer equations for planar or axisymmetric flows. Three-point implicit difference relations are used to reduce the momentum and energy equations to finite-difference form. These equations are solved simultaneously without iteration. Turbulent flow is treated by the inclusion of either a two-layer eddy-viscosity model or a mixing-length formulation. The eddy conductivity is related to the eddy viscosity through a static turbulent Prandtl number which may be an arbitrary function of the distance from the wall boundary. The transitional boundary layer is treated by the inclusion of an intermittency function which modifies the fully turbulent model. The laminar-boundary-layer equations are recovered when the intermittency is zero, and the fully turbulent equations are solved when the intermittency is unity.

Price, J. M.↗

An equation of state for oxygen and nitrogen 2

Preliminary equations of state for oxygen and nitrogen are presented which are the result of least squares fitting of the available P-density-T data, heat capacities at constant volume, and data to establish the criteria for phase equilibrium, simultaneously. The equation for nitrogen is applicable to 10,000 atmospheres at temperatures from 70 K to 1300 K, and the equation for oxygen to 350 atmospheres for temperatures from 56 K to 323 K. The development of the functional form of the equation of state, data weighting methods, and the techniques employed for least squares fitting of related data are discussed. Comparisons of the equation of state developed with selected values of experimental P-density-T data, and heat capacity values are included for both fluids.

Stewart, R. B.↗

An equation of state for oxygen and nitrogen

Preliminary equations of state are presented for oxygen and nitrogen which provide accurate representations of the available P-density-T data for both fluids. The equation for nitrogen is applicable for temperatures from 70 K to 1300 K at pressures to 10,000 atmospheres, and the equation for oxygen for temperatures from 70 K to 323 K at pressures to 350 atmospheres. Deviations of calculated densities from representative experimental data are included. A volume-explicit equation of state for oxygen to be used in estimating density values in the range of applicability of the equation of state is also presented.

Stewart, R. B.↗

Brownian motion from Boltzmann's equation.

Two apparently disparate lines of inquiry in kinetic theory are shown to be equivalent: (1) Brownian motion as treated by the (stochastic) Langevin equation and Fokker-Planck equation; and (2) Boltzmann's equation. The method is to derive the kinetic equation for Brownian motion from the Boltzmann equation for a two-component neutral gas by a simultaneous expansion in the density and mass ratios.

Montgomery, D.↗

Similarity and generalized finite-difference solutions of parabolic partial differential equations.

Techniques are presented for obtaining generalized finite-difference solutions to partial differential equations of the parabolic type. It is shown that the advantages of similarity in the solution of similar problems are generally not lost if the solution to the original partial differential equations is effected in the physical plane by finite-difference methods. The analysis results in a considerable saving in computational effort in the solution of both similar and nonsimilar problems. Several examples, including both the heat-conduction equation and the boundary-layer equations, are given. The analysis also provides a practical means of estimating the accuracy of finite-difference solutions to parabolic equations.

Clausing, A. M.↗

A rigorous cosmic-ray transport equation with no restrictions on particle energy.

A new transport equation for the cosmic-ray omnidirectional intensity is obtained. This equation follows exactly from the coupled pair of differential moment equations we presented earlier. It can be characterized as a nonlocal convection-diffusion equation in which the usual transport coefficients are replaced by time integral operators. The nonlocal equation is shown to reduce to the standard convection-diffusion form if the adiabatic approximation can be applied. In general, the adiabatic approximation does not apply; however, by going to the limits of infinite and zero gyroradius and, in addition, applying the adiabatic approximation, the large- and small-gyroradius transport theories due originally to Jokipii are regained. The validity of these theories as asymptotic limits and as approximate theories in the interplanetary magnetic field is discussed.

Klimas, A. J.↗

Power-spectral-density relationship for retarded differential equations

The power spectral density (PSD) relationship between input and output of a set of linear differential-difference equations of the retarded type with real constant coefficients and delays is discussed. The form of the PSD relationship is identical with that applicable to unretarded equations. Since the PSD relationship is useful if and only if the system described by the equations is stable, the stability must be determined before applying the PSD relationship. Since it is sometimes difficult to determine the stability of retarded equations, such equations are often approximated by simpler forms. It is pointed out that some common approximations can lead to erroneous conclusions regarding the stability of a system and, therefore, to the possibility of obtaining PSD results which are not valid.

Barker, L. K.↗

Classical seventh-, sixth-, and fifth-order Runge-Kutta-Nystrom formulas with stepsize control for general second-order differential equations

Runge-Kutta-Nystrom formulas of the seventh, sixth, and fifth order were derived for the general second order (vector) differential equation written as the second derivative of x = f(t, x, the first derivative of x). The formulas include a stepsize control procedure, based on a complete coverage of the leading term of the local truncation error in x, and they require no more evaluations per step than the earlier Runge-Kutta formulas for the first derivative of x = f(t, x). The developed formulas are expected to be time saving in comparison to the Runge-Kutta formulas for first-order differential equations, since it is not necessary to convert the second-order differential equations into twice as many first-order differential equations. The examples shown saved from 25 percent to 60 percent more computer time than the earlier formulas for first-order differential equations, and are comparable in accuracy.

Fehlberg, E.↗

Differential equation of exospheric lateral transport and its application to terrestrial hydrogen

The differential equation description of exospheric lateral transport of Hodges and Johnson is reformulated to extend its utility to light gases. Accuracy of the revised equation is established by applying it to terrestrial hydrogen. The resulting global distributions for several static exobase models are shown to be essentially the same as those that have been computed by Quessette using an integral equation approach. The present theory is subsequently used to elucidate the effects of nonzero lateral flow, exobase rotation, and diurnal tidal winds on the hydrogen distribution. Finally it is shown that the differential equation of exospheric transport is analogous to a diffusion equation. Hence it is practical to consider exospheric transport as a continuation of thermospheric diffusion, a concept that alleviates the need for an artificial exobase dividing thermosphere and exosphere.

Hodges, R. R., Jr.↗

Singular perturbation equations for flexible satellites

The dynamic model of a flexible satellite with n generalized structural position coordinates requires the solution of a set of N coupled nonlinear ordinary or partial differential equations. For single composite body satellites, N is equal to (n + 3). For dual-spin systems, N is equal to (n + 9). These equations usually involve time derivatives up to the second order. For large values of n, linearization of the system has so far been the only practicable way of solution. The present study shows a method of avoiding this linearization by reducing the N equations to three or nine nonlinear, coupled, first-order ordinary differential equations involving only the angular velocities of the composite bodies. The solutions for these angular velocities lead to linear equations in the n generalized structural position coordinates, which can then be solved by known methods.

Huang, T. C.↗

An equation of state for oxygen and nitrogen

Recent measurements of the thermodynamic properties of oxygen and nitrogen have been used to develop a single equation of state for both fluids; with appropriate coefficients, the equation describes the pressure-density-temperature surface for each fluid within the estimated experimental uncertainty of the data. The formulation uses least squares curve fitting of the equation of state to the available experimental data (including pressure-density-temperature data, saturation data, and measured isochoric heat capacities). New vapor pressure equations and equations for the ideal gas heat capacities of oxygen and nitrogen are included for use in calculation of derived thermodynamic properties.

Jacobsen, R. T.↗

Nonlinear equations of motion for the elastic bending and torsion of twisted nonuniform rotor blades

The equations of motion are developed by two complementary methods, Hamilton's principle and the Newtonian method. The resulting equations are valid to second order for long, straight, slender, homogeneous, isotropic beams undergoing moderate displacements. The ordering scheme is based on the restriction that squares of the bending slopes, the torsion deformation, and the chord/radius and thickness/radius ratios are negligible with respect to unity. All remaining nonlinear terms are retained. The equations are valid for beams with mass centroid axis and area centroid (tension) axis offsets from the elastic axis, nonuniform mass and stiffness section properties, variable pretwist, and a small precone angle. The strain-displacement relations are developed from an exact transformation between the deformed and undeformed coordinate systems. These nonlinear relations form an important contribution to the final equations. Several nonlinear structural and inertial terms in the final equations are identified that can substantially influence the aeroelastic stability and response of hingeless helicopter rotor blades.

Hodges, D. H.↗

Almost periodic solutions to difference equations

The theory of Massera and Schaeffer relating the existence of unique almost periodic solutions of an inhomogeneous linear equation to an exponential dichotomy for the homogeneous equation was completely extended to discretizations by a strongly stable difference scheme. In addition it is shown that the almost periodic sequence solution will converge to the differential equation solution. The preceding theory was applied to a class of exponentially stable partial differential equations to which one can apply the Hille-Yoshida theorem. It is possible to prove the existence of unique almost periodic solutions of the inhomogeneous equation (which can be approximated by almost periodic sequences) which are the solutions to appropriate discretizations. Two methods of discretizations are discussed: the strongly stable scheme and the Lax-Wendroff scheme.

Bayliss, A.↗

Attenuation of sound in ducts with acoustic treatment: A generalized approximate equation

A generalized approximate equation for duct lining sound attenuation is presented. The specification of two parameters, the maximum possible attenuation and the optimum wall acoustic impedance is shown to completely determine the sound attenuation for any acoustic mode at any selected wall impedance. The equation is based on the nearly circular shape of the constant attenuation contours in the wall acoustic impedance plane. For impedances far from the optimum, the equation reduces to Morse's approximate expression. The equation can be used for initial acoustic liner design. Not least important is the illustrative nature of the solutions which provide an understanding of the duct propagation problem usually obscured in the exact calculations. Sample calculations using the approximate attenuation equation show that the peak and the bandwidth of the sound attenuation spectrum can be represented by quite simple functions of the ratio of actual wall acoustic resistance to optimum resistance.

Rice, E. J.↗

The method of Ritz applied to the equation of Hamilton

Without any reference to the theory of differential equations, the initial value problem of the nonlinear, nonconservative double pendulum system is solved by the application of the method of Ritz to the equation of Hamilton. Also shown is an example of the reduction of the traditional eigenvalue problem of linear, homogeneous, differential equations of motion to the solution of a set of nonhomogeneous algebraic equations. No theory of differential equations is used. Solution of the time-space path of the linear oscillator is demonstrated and compared to the exact solution.

Bailey, C. D.↗