SEARCH · Search NASA

Results for “general relativity equations & solutions”

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 91 records · Page 5

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Non-existence of dynamical perturbations of Schwarzschild with vanishing self-dual part

A recent proposal of Ashtekar (1985) for the canonical quantisation of gravity involves a new choice of a (complex) polarisation on the (real) phase space of general relativity. In order for the 'four dimensional' version of Ashtekar's proposal to be viable, it is necessary that for an arbitrary (real) background solution there exist a sufficient number of (complex) solutions of the linearised Einstein equation with perturbed Weyl spinor having vanishing self-dual part. It is shown that this is not the case by explicitly demonstrating that for the Schwarzschild spacetime, aside from stationary perturbations, there are no such linearised solutions which can be expressed as superpositions of modes of real frequencies. The viability of the '3+1 version' of Ashtekar's proposal is not affected by these considerations.

Wald, R. M.↗

NASA NTRS1980

Variations on a spherically symmetrical accretion flow

Spherically symmetrical accretion onto a compact gravitating object is considered from a purely hydrodynamical viewpoint while disregarding heating and cooling processes. The study derives the radial perturbation equation, expresses it in a Eulerian form which is simpler than the Lagrangian form considered, for example, by Balazs (1972), and discusses the relation between the two formulations. This perturbation equation is used for the stability analysis of standing waves as well as for the study of short wavelength traveling waves which are a generalization of Parker's (1966) traveling wave solutions in the sense that they are not restricted to isothermal flows. These traveling waves can cause fluctuations in the mass accretion rate to grow large, but will not drive the unperturbed solution away from its assumed stationary profile. The waves conserve the total energy contained in the perturbation.

Petterson, J. A.↗

NASA NTRS1982

Particle trajectory computer program for icing analysis of axisymmetric bodies

General aviation aircraft and helicopters exposed to an icing environment can accumulate ice resulting in a sharp increase in drag and reduction of maximum lift causing hazardous flight conditions. NASA Lewis Research Center (LeRC) is conducting a program to examine, with the aid of high-speed computer facilities, how the trajectories of particles contribute to the ice accumulation on airfoils and engine inlets. This study, as part of the NASA/LeRC research program, develops a computer program for the calculation of icing particle trajectories and impingement limits relative to axisymmetric bodies in the leeward-windward symmetry plane. The methodology employed in the current particle trajectory calculation is to integrate the governing equations of particle motion in a flow field computed by the Douglas axisymmetric potential flow program. The three-degrees-of-freedom (horizontal, vertical, and pitch) motion of the particle is considered. The particle is assumed to be acted upon by aerodynamic lift and drag forces, gravitational forces, and for nonspherical particles, aerodynamic moments. The particle momentum equation is integrated to determine the particle trajectory. Derivation of the governing equations and the method of their solution are described in Section 2.0. General features, as well as input/output instructions for the particle trajectory computer program, are described in Section 3.0. The details of the computer program are described in Section 4.0. Examples of the calculation of particle trajectories demonstrating application of the trajectory program to given axisymmetric inlet test cases are presented in Section 5.0. For the examples presented, the particles are treated as spherical water droplets. In Section 6.0, limitations of the program relative to excessive computer time and recommendations in this regard are discussed.

Frost, Walter↗

NASA NTRS1947

On the one-dimensional theory of steady compressible fluid flow in ducts with friction and heat addition

Steady, diabatic (nonadiabatic), frictional, variable-area flow of a compressible fluid is treated in differential form on the basis of the one-dimensional approximation. The basic equations are first stated in terms of pressure, temperature, density, and velocity of the fluid. Considerable simplification and unification of the equations are then achieved by choosing the square of the local Mach number as one of the variables to describe the flow. The transformed system of equations thus obtained is first examined with regard to the existence of a solution. It is shown that, in general, a solution exists whose calculation requires knowledge only of the variation with position of any three of the dependent variables of the system. The direction of change of the flow variables can be obtained directly from the transformed equations without integration. As examples of this application of the equations, the direction of change of the flow variables is determined for two special flows. In the particular case when the local Mach number m = 1, a special condition must be satisfied by the flow if a solution is to exist. This condition restricts the joint rate of variation of heating, friction, and area at m = 1. Further analysis indicates that when a solution exists at this point it is not necessarily unique. Finally it is shown that the physical phenomenon of choking, which is known to occur in certain simple flow situations, is related to restrictions imposed on the variables by the form of the transformed equations. The phenomenon of choking is thus given a more general significance in that the transformed equations apply to a more general type of flow than has hitherto been treated. (author)

HIGH-SPEED↗

NASA NTRS2023

Another Geometric Interpretation of Cramer’s Rule

We develop a geometric interpretation of Cramer’s rule as a generalization of projection onto orthogonal basis vectors using the rows of the adjugate. This interpretation makes connections between elementary linear algebra concepts like the solution to linear equations, inner products, and projections. Such connections are useful for introducing broader concepts related to Hilbert spaces and geometric algebras like Grassman algebra. Such connections were essential for the author’s mathematical education as an engineer.

linear algebra↗

NASA NTRS1991

A Fourier-Chebyshev pseudospectral method for solving steady 3-D Navier-Stokes and heat equations in cylindrical cavities

A Fourier-Chebyshev pseudospectral method for solving steady 3D Navier-Stokes equations in cylindrical cavities is presented and discussed. The general method is pseudo-unsteady and uses a semi-implicit finite difference scheme for the time discretization. The generalized ADI (Alternating Direction Implicit) procedure is then applied to reduce the problem to a successive solution of one-dimensional problems. The spatial approximation uses a Fourier-Galerkin approximation in the periodic direction and a Chebyshev-collocation approximation in the other directions. Difficulties related to the pressure are surmounted by using the artifical compressibility method. A suitable variable change was chosen to avoid the problem of singularity at the axis generated by cylindrical coordinates. The method is first tested on an advection-diffusion equation and then on the Navier-Stokes equations. Finally, the method is illustrated by a convection problem of a differentially heated fluid.

Pulicani, J. P.↗

NASA NTRS1988

Solution of an integral equation encountered in studies on radiative transfer in completely absorbing leaf canopies

The leaf normal distribution function is related to the extinction coefficient of direct solar radiation by a Fredholm integral equation of the first kind. The integral equation is solved by the method of constrained least squares. The solutions are stable and correct only in extreme cases and, in general, it is not possible to recover an arbitrary leaf normal distribution function.

Myneni, Ranga B.↗

NASA NTRS1994

A linear method for analyzing lightning field changes

A constrained, least-squares method for analyzing multiple-station measurements of lightning field changes (delta Es) is introduced. Previous methods have attempted to fit the spatial pattern of lightning delta Es using nonlinear models, such as a point charge (Q) or a point dipole (P) model. With the linear method, the delta Es are described not by models but by a general volume charge distribution that is deposited on a large (40 x 40 x 20 cu km) Cartesian grid above the measuring network. A linear system of equations is used to relate the measured delta Es to the charges that are deposited at each grid point. With this approach, the information content of the measurements can be quantified by an eigenanalysis of the covariance matrix of the linear system. Constraints can be used to reduce the infinity of possible solutions to the linear system and also to reduce systematic biases that can be introduced by the method of solution. It is shown that a Landweber iterative method, derived from the general method of steepest descent, can be used to solve the linear system and that the resulting volume charge distributions are generally consistent with computer-simulated charge sources, when these sources are over the measuring network. The Landweber iteration has also provided solutions for natural lightning events that are consistent with Q- and P-model results.

Koshak, William J.↗

NASA NTRS1994

Boundary element analysis on vector and parallel computers

Boundary element analysis (BEA) can be characterized as a numerical technique that generally shifts the computational burden in the analysis toward numerical integration and the solution of nonsymmetric and either dense or blocked sparse systems of algebraic equations. Researchers have explored the concept that the fundamental characteristics of BEA can be exploited to generate effective implementations on vector and parallel computers. In this paper, the results of some of these investigations are discussed. The performance of overall algorithms for BEA on vector supercomputers, massively data parallel single instruction multiple data (SIMD), and relatively fine grained distributed memory multiple instruction multiple data (MIMD) computer systems is described. Some general trends and conclusions are discussed, along with indications of future developments that may prove fruitful in this regard.

Kane, J. H.↗

NASA NTRS1989

The evolution of energetic particles and the emitted radiation in solar flares

The evolution of accelerated particle distributions in a magnetized plasma and the resulting radiation are calculated, and the results are applied to solar flares. To study the radiation on timescales of order the particle lifetimes, the evolution of the particle distribution is determined by the use of the Fokker-Planck equation including Coulomb collisions and magnetic mirroring. Analytic solution to the equations are obtained for limiting cases such as homogeneous injection in a homogeneous plasma, and for small pitch angle. These analytic solutions are then used to place constraints on flare parameters such as density, loop length, and the injection timescale for very short implusive solar flares. For general particle distributions in arbitrary magnetic field and background density, the equation is solved numerically. The relative timing of microwaves and X-rays during individual flares is investigated. A number of possible sources for excessive microwave flux are discussed including a flattening in the electron spectrum above hard X-ray energies, thermal synchrotron emission, and trapping of electron by converging magnetic fields. Over shorter timescales, the Fokker-Planck equation is solved numerically to calculate the temporal evolution of microwaves and X-rays from nonthermal thick target models. It is shown that magnetic trapping will not account for the observed correlation of microwaves of approximately 0.15 seconds behind X-rays in flares with rapid time variation, and thus higher energy electrons must be accelerated later than lower energy electrons.

Lu, Edward Tsang↗

NASA NTRS1976

A new unified approach for analyzing wing-body-tail configurations with control surfaces

A general theory for steady and unsteady, subsonic and supersonic potential aerodynamics for complex configurations is presented. Special attention is given to the theoretical formulation and the corresponding numerical implementation for coplanar interfering surfaces. Applying the Green's function method to the equation of the velocity potential and discretizing the spatial problem by using the finite-element technique, yields a set of differential-delay equations in time relating the potential to the normal wash. For fully unsteady flow, the motion is assumed to consist of constant subsonic or supersonic speed for time t less than or equal to 0 (steady state) and of small perturbations around the steady state for time t greater than 0; the solution is obtained in the Laplace domain. From the potential, the aerodynamic pressure and the generalized forces are evaluated. The program SOUSSA (Steady, Oscillatory and Unsteady Subsonic and Supersonic Aerodynamics) is briefly described. Numerical results obtained with SOUSSA are presented.

Tseng, K.↗

NASA NTRS1981

Mesh generation by conformal and quasiconformal mappings

It is pointed out that many recent advances in the finite-difference solution of elliptic equations have been limited to regions whose boundary contours coincide with coordinate lines of the Cartesian coordinate system. The reason for this is related to the fact that in the case of an arbitrary curvilinear coordinate system the original equation becomes much more complex. However, there is no added complexity if an orthogonal coordinate system is generated from a conformal mapping. In the present investigation, a finite difference method developed for the construction of conformal mappings has been generalized to construct quasi-conformal mappings. It is expected that the use of more sophisticated numerical algorithms could lead to improvements in both speed and accuracy. Quasi-conformal mappings have applications not only in the solution of elliptic equations but also in other areas such as orthogonal mesh generation on surfaces and the solution of certain fluid flow problems.

Mastin, C. W.↗

NASA NTRS1980

Collisional redistribution of radiation. II - The effects of degeneracy on the equations of motion for the density matrix. III - The equation of motion for the correlation function and the scattered spectrum

The effect of correlations between an absorber atom and perturbers in the binary-collision approximation are applied to degenerate atomic systems. A generalized absorption profile which specifies the final state of the atom after an absorption event is related to the total intensities of Rayleigh scattering and fluorescence from the atom. It is suggested that additional dynamical information to that obtainable from ordinary absorption experiments is required in order to describe redistributed atomic radiation. The scattering of monochromatic radiation by a degenerate atom is computed in a binary-collision approximation; an equation of motion is derived for the correlation function which is valid outside the quantum-regression regime. Solutions are given for the weak-field conditions in terms of generalized absorption and emission profiles that depend on the indices of the atomic multipoles.

Burnett, K.↗

NASA NTRS1972

Magnetic latitude effects in the solar wind

The Weber-Davis model of the solar wind is generalized to include the effects of latitude. The principal assumptions of high electrical conductivity, rotational symmetry, the polytropic relation between pressure and density, and a flow-alined field in a system rotating with the sun, are retained. An approximate solution to the resulting equations for spherical boundary conditions at the base of the corona indicates a small component of latitudinal flow toward the solar poles at large distances from the sun as result of latitudinal magnetic forces.

Winge, C. R., Jr.↗

NASA NTRS1980

Stress intensity factors in two bonded elastic layers containing cracks perpendicular to and on the interface. Part 1: Analysis

The basic crack problem which is essential for the study of subcritical crack propagation and fracture of layered structural materials is considered. Because of the apparent analytical difficulties, the problem is idealized as one of plane strain or plane stress. An additional simplifying assumption is made by restricting the formulation of the problem to crack geometries and loading conditions which have a plane of symmetry perpendicular to the interface. The general problem is formulated in terms of a coupled system of four integral equations. For each relevant crack configuration of practical interest, the singular behavior of the solution near and at the ends and points of intersection of the cracks is investigated and the related characteristic equations are obtained. The edge crack terminating at and crossing the interface, the T-shaped crack consisting of a broken layer and a delamination crack, the cross-shaped crack which consists of a delamination crack intersecting a crack which is perpendicular to the interface, and a delamination crack initiating from a stress-free boundary of the bonded layers are some of the practical crack geometries considered.

Lu, M. C.↗

NASA NTRS1983

Stress intensity factors in two bonded elastic layers containing cracks perpendicular to and on the interface. I Analysis. II - Solution and results

The basic crack problem which is essential for the study of subcritical crack propagation and fracture of layered structural materials is considered. Because of the apparent analytical difficulties, the problem is idealized as one of plane strain or plane stress. An additional simplifying assumption is made by restricting the formulation of the problem to crack geometries and loading conditions which have a plane of symmetry perpendicular to the interface. The general problem is formulated in terms of a coupled systems of four integral equations. For each relevant crack configuration of practical interest, the singular behavior of the solution near and at the ends and points of intersection of the cracks is investigated and the related characteristic equations are obtained. The edge crack terminating at and crossing the interface, the T-shaped crack consisting of a broken layer and a delamination crack, the cross-shaped crack which consists of a delamination crack intersecting a crack which is perpendicular to the interface, and a delamination crack initiating from a stress-free boundary of the bonded layers are some of the practical crack geometries considered. Previously announced in STAR as N80-18428 and N80-18429

Lu, M.-C.↗

NASA NTRS1994

Variations in the modal characteristics of a telescopically deploying beam

The equations of motion for a two-segment deploying telescopic beam are derived through application of Lagrange's equation. The outer tube of the beam is fixed at one end and the inner tube slides freely relative to the fixed segment. The resulting nonlinear, non-autonomous set of equations is linearized and simplified to the standard Euler-Bernoulli partial differential equations for an elastic beam by freezing the deployment process at various stages of deployment, and examining the small amplitude and natural modes of vibration of the resulting configuration. Application of the natural boundary conditions and compatibility of motion relations for the two segments in their common region of overlap leads to a transcendental characteristic equation in the frequency parameter Beta(L). Numerical solution of the equation for the characteristic roots determines the modal frequencies, and the corresponding mode shapes are obtained from the general solution of the Euler-Bernoulli equation tailored to the natural boundary conditions. Sample results of modal frequencies and shapes are presented for various stages of deployment and discussed. It is shown that for all intermediate stages of deployment (between 0 and 100 percent) the spectral distribution is drastically altered by the appearance of regions of very closely spaced modal frequencies. The sources of this modal agglomeration are explored.

Amos, Anthony K.↗

NASA NTRS1979

Development of a Model and Computer Code to Describe Solar Grade Silicon Production Processes

Mathematical models and computer codes based on these models, which allow prediction of the product distribution in chemical reactors for converting gaseous silicon compounds to condensed-phase silicon were developed. The following tasks were accomplished: (1) formulation of a model for silicon vapor separation/collection from the developing turbulent flow stream within reactors of the Westinghouse (2) modification of an available general parabolic code to achieve solutions to the governing partial differential equations (boundary layer type) which describe migration of the vapor to the reactor walls, (3) a parametric study using the boundary layer code to optimize the performance characteristics of the Westinghouse reactor, (4) calculations relating to the collection efficiency of the new AeroChem reactor, and (5) final testing of the modified LAPP code for use as a method of predicting Si(1) droplet sizes in these reactors.

Srivastava, R.↗