SEARCH · Engineering Papers

Results for “general relativity equations and 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 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Space-time dependent thermal conductivity in nonlocal thermal transport

Nonlocal thermal transport is generally described by the Peierls-Boltzmann transport equation (PBE). However, solving the PBE for a general space-time dependent problem remains a challenging task due to the high dimensionality of the integro-differential equation. In this work, we present a direct solution to the space-time dependent PBE with a linearized collision matrix using an eigendecomposition method. We show that there exists a generalized Fourier-type relation that links heat flux to the local temperature, and this constitutive relation defines a thermal conductivity that depends on both time and space. Combining this approach with ab initio calculations of phonon properties, we demonstrate that the space-time dependent thermal conductivity gives rise to an oscillatory response in temperature in a transient grating geometry in high thermal conductivity materials. The present solution method allows us to extend the reach of our computational capability for heat conduction to space-time dependent nondiffusive transport regimes. Here, this capability will not only enable a more accurate interpretation of thermal measurements that observe nonlocal thermal transport, but also enhance our physical understanding of nonlocal thermal transport in high thermal conductivity materials that are promising candidates for nanoscale thermal management applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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↗

DOE OSTI2021

Growth of accretion driven scalar hair around Kerr black holes

We report that Scalar fields around compact objects are of interest for scalar-tensor theories of gravity and dark matter models consisting of a massive scalar, e.g., axions. We study the behavior of a scalar field around a Kerr black hole with nontrivial asymptotic boundary conditions - both nonzero density and nonzero angular momentum. Starting from an initial radially homogeneous configuration, a scalar cloud is accreted, which asymptotes to known stationary configurations over time. We study the cloud growth for different parameters including black hole spin, scalar field mass, and the scalar field density and angular momentum far from the black hole. We characterize the transient growth of the mass and angular momentum in the cloud, and the spatial profile of the scalar around the black hole, and relate the results of fully nonlinear simulations to an analytic perturbative expansion. We also highlight the potential for these accreted clouds to create monochromatic gravitational wave signals - similar to the signals from superradiant clouds, although significantly weaker in amplitude.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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

DOE OSTI2024

Quasinormal modes of naked singularities in presence of nonlinear scalar fields

We study linear perturbations against static spherically symmetric background configurations of general relativity with a real scalar field (SF), which is minimally coupled with gravity; it is nonlinear due to the presence of the self-action potential. The background solutions have a naked singularity at the center of the configuration. The focus is on the stability of the background and fundamental frequencies of the quasinormal modes (QNM) of the axial perturbations in the Regge-Wheeler gauge. The problem is reduced to one hyperbolic master equation with an effective potential 𝑊eff, which turns out to be positive for a general non-negative SF potential; this ensures the linear stability with respect to this kind of perturbations. For numerical simulations, the SF potential was chosen in the power-law form 𝑉⁡(𝜙) ∼ 𝜙 2⁢𝑛 with 2 < 𝑛 ≤40. We extracted the fundamental frequencies of QNM for different 𝑛 and various sets of the background configuration parameters. The results show that even for a small background SF, there is a significant difference between the fundamental frequencies and ones in case of the Schwarzschild background. Furthermore, the results are also compared with the case of the Fisher-Janis-Newman-Winicour background dealing with a massless linear scalar field.

General relativity↗

DOE OSTI2020

Extended Lagrangian Born–Oppenheimer molecular dynamics using a Krylov subspace approximation

It is shown how the electronic equations of motion in extended Lagrangian Born–Oppenheimer molecular dynamics simulations can be integrated using low-rank approximations of the inverse Jacobian kernel. This kernel determines the metric tensor in the harmonic oscillator extension of the Lagrangian that drives the evolution of the electronic degrees of freedom. The proposed kernel approximation is derived from a pseudoinverse of a low-rank estimate of the Jacobian, which is expressed in terms of a generalized set of directional derivatives with directions that are given from a Krylov subspace approximation. The approach allows a tunable and adaptive approximation that can take advantage of efficient preconditioning techniques. The proposed kernel approximation for the integration of the electronic equations of motion makes it possible to apply extended Lagrangian first-principles molecular dynamics simulations to a broader range of problems, including reactive chemical systems with numerically sensitive and unsteady charge solutions. This can be achieved without requiring exact full calculations of the inverse Jacobian kernel in each time step or relying on iterative non-linear self-consistent field optimization of the electronic ground state prior to the force evaluations as in regular direct Born–Oppenheimer molecular dynamics. We note the low-rank approximation of the Jacobian is directly related to Broyden’s class of quasi-Newton algorithms and Jacobian-free Newton–Krylov methods and provides a complementary formulation for the solution of nonlinear systems of equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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

DOE OSTI2020

Lorentz Symmetry Group, Retardation, Intergalactic Mass Depletion and Mechanisms Leading to Galactic Rotation Curves

The general theory of relativity (GR) is symmetric under smooth coordinate transformations, also known as diffeomorphisms. The general coordinate transformation group has a linear subgroup denoted as the Lorentz group of symmetry, which is also maintained in the weak field approximation to GR. The dominant operator in the weak field equation of GR is thus the d’Alembert (wave) operator, which has a retarded potential solution. Galaxies are huge physical systems with dimensions of many tens of thousands of light years. Thus, any change at the galactic center will be noticed at the rim only tens of thousands of years later. Those retardation effects are neglected in the present day galactic modelling used to calculate rotational velocities of matter in the rims of the galaxy and surrounding gas. The significant differences between the predictions of Newtonian instantaneous action at a distance and observed velocities are usually explained by either assuming dark matter or by modifying the laws of gravity (MOND). In this paper, we will show that, by taking general relativity seriously without neglecting retardation effects, one can explain the radial velocities of galactic matter in the M33 galaxy without postulating dark matter. It should be stressed that the current approach does not require that velocities v are high; in fact, the vast majority of galactic bodies (stars, gas) are substantially subluminal—in other words, the ratio of $\frac{v}{c} \ll 1$. Typical velocities in galaxies are 100 km/s, which makes this ratio 0.001 or smaller. However, one should consider the fact that every gravitational system, even if it is made of subluminal bodies, has a retardation distance, beyond which the retardation effect cannot be neglected. Every natural system, such as stars and galaxies and even galactic clusters, exchanges mass with its environment, for example, the sun loses mass through solar wind and galaxies accrete gas from the intergalactic medium. This means that all natural gravitational systems have a finite retardation distance. The question is thus quantitative: how large is the retardation distance? For the M33 galaxy, the velocity curve indicates that the retardation effects cannot be neglected beyond a certain distance, which was calculated to be roughly 14,000 light years; similar analysis for other galaxies of different types has shown similar results. We demonstrate, using a detailed model, that this does not require a high velocity of gas or stars in or out of the galaxy and is perfectly consistent with the current observational knowledge of galactic and extra galactic material content and dynamics

79 ASTRONOMY AND ASTROPHYSICS↗

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↗