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 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS1977

Pure field theories and MACSYMA algorithms

A pure field theory attempts to describe physical phenomena through singularity-free solutions of field equations resulting from an action principle. The physics goes into forming the action principle and interpreting specific results. Algorithms for the intervening mathematical steps are sketched. Vacuum general relativity is a pure field theory, serving as model and providing checks for generalizations. The fields of general relativity are the 10 components of a symmetric Riemannian metric tensor; those of the Einstein-Straus generalization are the 16 components of a nonsymmetric. Algebraic properties are exploited in top level MACSYMA commands toward performing some of the algorithms of that generalization. The light cone for the theory as left by Einstein and Straus is found and simplifications of that theory are discussed.

Ament, W. S.↗

NASA NTRS1982

Matter in general relativity

Two theories of matter in general relativity, the fluid theory and the kinetic theory, were studied. Results include: (1) a discussion of various methods of completing the fluid equations; (2) a method of constructing charged general relativistic solutions in kinetic theory; and (3) a proof and discussion of the incompatibility of perfect fluid solutions in anisotropic cosmologies. Interpretations of NASA gravitational experiments using the above mentioned results were started. Two papers were prepared for publications based on this work.

Ray, J. R.↗

NASA NTRS1986

Radiative transfer in a sphere illuminated by a parallel beam - An integral equation approach

The problem of multiple scattering of nonpolarized light in a planetary body of arbitrary shape illuminated by a parallel beam is formulated using the integral equation approach. There exists a simple functional whose stationarity condition is equivalent to solving the equation of radiative transfer and whose value at the stationary point is proportional to the differential cross section. The analysis reveals a direct relation between the microscopic symmetry of the phase function for each scattering event and the macroscopic symmetry of the differential cross section for the entire planetary body, and the interconnection of these symmetry relations and the variational principle. The case of a homogeneous sphere containing isotropic scatterers is investigated in detail. It is shown that the solution can be expanded in a multipole series such that the general spherical problem is reduced to solving a set of decoupled integral equations in one dimension. Computations have been performed for a range of parameters of interest, and illustrative examples of applications to planetary problems as provided.

Shia, R.-L.↗

NASA NTRS1982

General curvilinear coordinate systems

The basic ideas of the construction and use of numerically-generated boundary-fitted coordinate systems for the numerical solution of partial differential equations are discussed. With such coordinate systems, all computation can be done on a fixed square grid in the rectangular transformed region regardless of the shape or movement of the physical boundaries. A number of different types of configurations for the transformed region and the basic transformation relations from a cartesian system to a general curvilinear system are given. The material of this paper is applicable to all types of coordinate system generation.

Thompson, J. P.↗

NASA NTRS1991

Analysis of lightning field changes produced by Florida thunderstorms

A new method is introduced for inferring the charges deposited in a lightning flash. Lightning-caused field changes (delta E's) are described by a more general volume charge distribution than is defined on a large cartesian grid system centered above the measuring networks. It is shown that a linear system of equations can be used to relate delta E's at the ground to the values of charge on this grid. It is possible to apply more general physical constraints to the charge solutions, and it is possible to access the information content of the delta E data. Computer-simulated delta E inversions show that the location and symmetry of the charge retrievals are usually consistent with the known test sources.

Koshak, William John↗

NASA NTRS2003

Implicit Approaches for Moving Boundaries in a 3-D Cartesian Method

This work considers numerical simulation of three-dimensional flows with time-evolving boundaries. Such problems pose a variety of challenges for numerical schemes, and have received a substantial amount of attention in the recent literature. Since such simulations are unsteady, time-accurate solution of the governing equations is required. In special cases, the body motion can be treated by a uniform rigid motion of the computational domain. For the more general situation of relative-body motion, however, this simplification is unavailable and the simulations require a mechanism for ensuring that the mesh evolves with the moving boundaries. This involves a "remeshing" of the computational domain (either localized or global) at each physical timestep, and places a premium on both the speed and robustness of the remeshing algorithms. This work presents a method which includes unsteady flow simulation, rigid domain motion, and relative body motion using a time-evolving Cartesian grid system in three dimensions.

Murman, Scott M.↗

NASA NTRS2020

Solution of the Nonlinear High-Fidelity Generalized Method of Cells Micromechanics Relations via Order-Reduction Techniques

The High-Fidelity Generalized Method of Cells (HFGMC) is one technique, distinct from traditional finite element approaches, for accurately simulating nonlinear composite material behavior. In this work, the HFGMC global system of equations for doubly-periodic repeating unit cells with nonlinear constituents has been reduced in size through the novel application of a Petrov-Galerkin Proper Orthogonal Decomposition order-reduction scheme in order to improve its computational efficiency. Order-reduced models of an E-glass/Nylon 12 composite led to a 4.8 to 6.3x speedup in the equation assembly/solution runtime while maintaining model accuracy. This corresponded to a 21 to 38 percent reduction in total runtime. The significant difference in assembly/solution and total runtimes was attributed to the evaluation of integration point inelastic field quantities; this step was identical between the unreduced and order-reduced models. Nonetheless, order-reduced techniques offer the potential to significantly improve the computational efficiency of multiscale calculations.

Ricks, Trenton M.↗

NASA NTRS1974

Evaluation of approximate relations for Delta /Q/ using a numerical solution of the Boltzmann equation

Data obtained from a numerical solution of the Boltzmann equation for shock-wave structure are used to test the accuracy of accepted approximate expressions for the two moments of the collision integral Delta (Q) for general intermolecular potentials in systems with a large translational nonequilibrium. The accuracy of the numerical scheme is established by comparison of the numerical results with exact expressions in the case of Maxwell molecules. They are then used in the case of hard-sphere molecules, which are the furthest-removed inverse power potential from the Maxwell molecule; and the accuracy of the approximate expressions in this domain is gauged. A number of approximate solutions are judged in this manner, and the general advantages of the numerical approach in itself are considered.

Nathenson, M.↗

NASA NTRS1972

A new two-metric theory of gravity with prior geometry

A theory is presented of gravity which has the same post-Newtonian (PN) as that of general relativity. The field equations, and a calculation of the PN limit of the theory are discussed along with the equations of stellar structure for static spherically symmetric stars. A special exterior spherically symmetric solution, time dependent solutions, conservations, and gravitational waves are analyzed.

Lightman, A. P.↗

NASA NTRS1989

Stability and dynamic analysis of a slender column with curved longitudinal stiffeners

The results of a stability design study are presented for a slender column with curved longitudinal stiffeners for large space structure applications. Linear stability analyses are performed using a link-plate representation of the stiffeners to determine stiffener local buckling stresses. Results from a set of parametric analyses are used to determine an approximate explicit expression for stiffener local buckling in terms of its geometric parameters. This expression along with other equations governing column stability and mass are assembled into a determinate system describing minimum mass stiffened column design. An iterative solution is determined to solve this system and a computer program incorporating this routine is presented. Example design problems are presented which verify the solution accuracy and illustrate the implementation of the solution routine. Also, observations are made which lead to a greatly simplified first iteration design equation relating the percent increase in column mass to the percent increase in column buckling load. From this, generalizations are drawn as to the mass savings offered by the stiffened column concept. Finally, the percent increase in fundamental column vibration frequency due to the addition of deployable stiffeners is studied.

Lake, Mark S.↗

NASA NTRS1995

A similarity hypothesis for the two-point correlation tensor in a temporally evolving plane wake

The analysis demonstrated that the governing equations for the two-point velocity correlation tensor in the temporally evolving wake admit similarity solutions, which include the similarity solutions for the single-point moment as a special case. The resulting equations for the similarity solutions include two constants, beta and Re(sub sigma), that are ratios of three characteristic time scales of processes in the flow: a viscous time scale, a time scale characteristic of the spread rate of the flow, and a characteristic time scale of the mean strain rate. The values of these ratios depend on the initial conditions of the flow and are most likely measures of the coherent structures in the initial conditions. The occurrences of these constants in the governing equations for the similarity solutions indicates that these solutions, in general, will only be the same for two flows if these two constants are equal (and hence the coherent structures in the flows are related). The comparisons between the predictions of the similarity hypothesis and the data presented here and elsewhere indicate that the similarity solutions for the two-point correlation tensors provide a good approximation of the measures of those motions that are not significantly affected by the boundary conditions caused by the finite extent of real flows. Thus, the two-point similarity hypothesis provides a useful tool for both numerical and physical experimentalist that can be used to examine how the finite extent of real flows affect the evolution of the different scales of motion in the flow.

Ewing, D. W.↗

NASA NTRS1981

Stresses in adhesively bonded joints - A closed-form solution

The general plane strain problem of adhesively bonded structures consisting of two different, orthotropic adherends is considered, under the assumption that adherend thicknesses are constant and small in relation to the lateral dimensions of the bonded region, so that they may be treated as plates. The problem is reduced to a system of differential equations for the adhesive stresses which is solved in closed form, with a single lap joint and a stiffened plate under various loading conditions being considered as examples. It is found that the plate theory used in the analysis not only predicts the correct trend for adhesive stresses but gives surprisingly accurate results, the solution being obtained by assuming linear stress-strain relations for the adhesive.

Delale, F.↗

NASA NTRS1982

Numerical computation of exponential matrices using the Cayley-Hamilton theorem

A method for computing exponential matrices, which often arise naturally in the solution of systems of linear differential equations, is developed. An exponential matrix is generated as a linear combination of a finite number (equal to the matrix order) of matrices, the coefficients of which are scalar infinite sums. The method can be generalized to apply to any formal power series of matrices. Attention is focused upon the exponential function, and the matrix exponent is assumed tri-diagonal in form. In such cases, the terms in the coefficient infinite sums can be extracted, as recursion relations, from the characteristic polynomial of the matrix exponent. Two numerical examples are presented in some detail: (1) the three dimensional infinitesimal rotation rate matrix, which is skew symmetric, and (2) an N-dimensional tri-diagonal and symmetric finite difference matrix which arises in the numerical solution of the heat conduction partial differential equation. In the second example, the known eigenvalues and eigenvectors of the finite difference matrix permit an analytical solution for the exponential matrix, through the theory of diagonalization and similarity transformations, which is used for independent verification. The convergence properties of the scalar infinite summations are investigated for finite difference matrices of various orders up to ten, and it is found that the number of terms required for convergence increases slowly with the order of the matrix.

Walden, H.↗

NASA NTRS1989

Corrections to the thin wall approximation in general relativity

The question is considered whether the thin wall formalism of Israel applies to the gravitating domain walls of a lambda phi(exp 4) theory. The coupled Einstein-scalar equations that describe the thick gravitating wall are expanded in powers of the thickness of the wall. The solutions of the zeroth order equations reproduce the results of the usual Israel thin wall approximation for domain walls. The solutions of the first order equations provide corrections to the expressions for the stress-energy of the wall and to the Israel thin wall equations. The modified thin wall equations are then used to treat the motion of spherical and planar domain walls.

Garfinkle, David↗

NASA NTRS1972

Application of a stochastic inverse to the geophysical inverse problem

The inverse problem for gross earth data can be reduced to an undertermined linear system of integral equations of the first kind. A theory is discussed for computing particular solutions to this linear system based on the stochastic inverse theory presented by Franklin. The stochastic inverse is derived and related to the generalized inverse of Penrose and Moore. A Backus-Gilbert type tradeoff curve is constructed for the problem of estimating the solution to the linear system in the presence of noise. It is shown that the stochastic inverse represents an optimal point on this tradeoff curve. A useful form of the solution autocorrelation operator as a member of a one-parameter family of smoothing operators is derived.

Jordan, T. H.↗

NASA NTRS1992

On the instability of boundary layers on heated flat plates

The stability of a boundary layer on a heated flat plate is investigated in the linear regime. The flow is shown to be unstable to longitudinal vortex structures which in general develop in a nonparallel manner in the streamwise direction. Solutions of the nonparallel equations are obtained numerically at O(1) values of the appropriate stability parameter, i.e., the Grashof number. The particular cases investigated relate to the situations when the instability is induced by localized or distributed wall roughness or nonuniform wall heating. The case when the vortices are induced by freestream disturbances is also considered. The fastest growing mode is found to be governed by a quasi-parallel theory at high wavenumbers. The wavenumber and growth rate of the fastest growing mode are found in closed form. At low wavenumbers the vortex instability is shown to be closely related to Tollmein-Schlichting waves. The effect of wall heating or cooling on the latter type of instability is discussed.

Hall, Philip↗

NASA NTRS1982

An integration scheme to determine the dynamic response of a launch vehicle with several payloads

The coupled equations of motion of a system composed of a launch vehicle and multiple payloads are derived. In the process, it is shown how superfluous interface degrees of freedom on the booster side can be accommodated in the formulation. The discrete system equations are directly integrated, avoiding the solution of an expensive system eigenvalue problem. A modified Newmark-Chan-Beta numerical integration scheme is used to obtain the response. The unique form of the equations of motion allow for a quick and accurate solution both for the response and the internal loads. Although completely general, the presented technique is most effective when the number of interface degrees of freedom is relatively small compared to the overall number of degrees of freedom of the system. The technique is applied to the case of the S.T.S.-S.T.-OMS Kit system. Finally, an evaluation of the presented method is included.

Engels, R. C.↗

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.↗