Engineering Papers⌕ Search

SEARCH · Engineering Papers

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.

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

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Theory of biaxial graded-index optical fiber

A biaxial graded-index fiber with a homogeneous cladding is studied. Two methods, wave equation and matrix differential equation, of formulating the problem and their respective solutions are discussed. For the wave equation formulation of the problem it is shown that for the case of a diagonal permittivity tensor the longitudinal electric and magnetic fields satisfy a pair of coupled second-order differential equations. Also, a generalized dispersion relation is derived in terms of the solutions for the longitudinal electric and magnetic fields. For the case of a step-index fiber, either isotropic or uniaxial, these differential equations can be solved exactly in terms of Bessel functions. For the cases of an istropic graded-index and a uniaxial graded-index fiber, a solution using the Wentzel, Krammers and Brillouin (WKB) approximation technique is shown. Results for some particular permittivity profiles are presented. Also the WKB solutions is compared with the vector solution found by Kurtz and Streifer. For the matrix formulation it is shown that the tangential components of the electric and magnetic fields satisfy a system of four first-order differential equations which can be conveniently written in matrix form. For the special case of meridional modes, the system of equations splits into two systems of two equations. A general iterative technique, asymptotic partitioning of systems of equations, for solving systems of differential equations is presented. As a simple example, Bessel's differential equation is written in matrix form and is solved using this asymptotic technique. Low order solutions for particular examples of a biaxial and uniaxial graded-index fiber are presented. Finally numerical results obtained using the asymptotic technique are presented for particular examples of isotropic and uniaxial step-index fibers and isotropic, uniaxial and biaxial graded-index fibers.

Kawalko, Stephen F.↗
NASA NTRS1987

A coupled mode formulation by reciprocity and a variational principle

A coupled mode formulation for parallel dielectric waveguides is presented via two methods: a reciprocity theorem and a variational principle. In the first method, a generalized reciprocity relation for two sets of field solutions satisfying Maxwell's equations and the boundary conditions in two different media, respectively, is derived. Based on the generalized reciprocity theorem, the coupled mode equations can then be formulated. The second method using a variational principle is also presented for a general waveguide system which can be lossy. The results of the variational principle can also be shown to be identical to those from the reciprocity theorem. The exact relations governing the 'conventional' and the new coupling coefficients are derived. It is shown analytically that the present formulation satisfies the reciprocity theorem and power conservation exactly, while the conventional theory violates the power conservation and reciprocity theorem by as much as 55 percent and the Hardy-Streifer (1985, 1986) theory by 0.033 percent, for example.

Chuang, Shun-Lien↗
DOE OSTI2025

Radiation GRMHD Models of Accretion onto Stellar-mass Black Holes. I. Survey of Eddington Ratios

We summarize results from a survey of radiation-dominated black hole accretion flows across a wide range of mass accretion rates, as well as two values of black hole spin and initial magnetic field geometry. These models apply an algorithm targeting direct solutions to the radiation transport equation in full general relativity and have been enabled by access to modern exascale computing systems. Super-Eddington accretion flows form geometrically thick radiation-pressure-supported disks that drive powerful equatorial outflows. A narrow funnel-shaped photosphere in the inner region results in very low radiative efficiencies in this regime. The structure of near- and sub-Eddington accretion depends on whether there is net vertical magnetic flux at the midplane of the disk. With net flux, the disk forms a thin, dense layer at the midplane surrounded by a magnetically dominated corona, whereas without net flux, the disk remains magnetically dominated everywhere. Although none of our models achieve the magnetically arrested disk regime, those with net vertical flux and a rapidly spinning black hole still produce powerful relativistic jets. Our calculations adopt simple opacity models (with scalings appropriate to stellar-mass black hole accretion). We discuss the application of our results to observations of X-ray binaries and ultraluminous X-ray sources such as Cyg X-3 and SS433. We also speculate on the application of our super-Eddington models to the interpretation of little red dots recently discovered by JWST.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1984

The effect of closed boundary conditions on a stationary dynamo

One of two boundary conditions generally assumed in solutions of the dynamo equation is related to the disappearance of the azimuthal field at the boundary. Parker (1984) points out that for the realization of this condition the field must escape freely through the surface. Escape requires that the field be detached from the gas in which it is embedded. In the case of the sun, this can be accomplished only through reconnection in the tenuous gas above the visible surface. Parker concludes that the observed magnetic activity on the solar surface permits at most three percent of the emerging flux to escape. He arrives at the conclusion that, instead of B(phi) = 0, the partial derivative of B(phi) to r is equal to zero. The present investigation is concerned with the effect of changing the boundary condition according to Parker's conclusion. Implications for the solar convection zone are discussed.

Choudhuri, A. R.↗
NASA NTRS1977

Unsteady supersonic aerodynamic theory for interfering surfaces by the method of potential gradient

A generalized solution of the hyperbolic wave equation was further developed to relate the velocity components at a field point to the potential gradient distribution in the dependence domain. Singular integrals were evaluated in closed form, with numerical integration methods for more complex but analytic functions. Idealization of the lifting surfaces by trapezoidal elements with two sides parallel to the streamlines is computationally efficient. Streamwise integrals were performed analytically, and spanwise integrals were neccessary only on element leading and trailing sides. All integrands vanish on the Mach cone. Pressure distribution on a double delta wing and generalized aerodynamic coefficients for three AGARD planforms were calculated and compared with available results.

Jones, W. P.↗
NASA NTRS1994

(abstract) Spacecraft Doppler Tracking with the Deep Space Network in the Search for Gravitational Waves

The Deep Space Network generates accurate radio science data observables for investigators who use radio links between spacecraft and the Earth to examine small changes in the phase and/or amplitude of the signal to study a wide variety of structures and phenomena in space. Several such studies are directed at aspects of the theory of general relativity such as gravitational redshift and gravitational waves. A gravitational wave is a propagating, polarized gravitational field, a ripple in the curvature of space-time. In Einstein's theory of general relativity, the waves are propagating solutions of the Einstein field equations. Their amplitudes are dimensionless strain amplitudes that change the fractional difference in distance between test masses and the rates at which separated clocks keep time. Predicted by all relativistic theories of gravity, they are extremely weak (the ratio of gravitational forces to electrical forces is about 10(sup -40)) and are generated at detectable levels only by astrophysical sources - very massive sources under violent dynamical conditions. The waves have never been detected but searches in the low-frequency band using Doppler tracking of many spacecraft have been conducted and others are being planned. Upper limits have been placed on the gravitational wave strength with the best sensitivities to date are for periodic waves being 7 x 10(sup -15).

gravitational waves theory of general relativity D↗
NASA NTRS1987

Finite difference solution for a generalized Reynolds equation with homogeneous two-phase flow

An attempt is made to relate elements of two-phase flow and kinetic theory to the modified generalized Reynolds equation and to the energy equation, in order to arrive at a unified model simulating the pressure and flows in journal bearings, hydrostatic journal bearings, or squeeze film dampers when a two-phase situation occurs due to sudden fluid depressurization and heat generation. The numerical examples presented furnish a test of the algorithm for constant properties, and give insight into the effect of the shaft fluid heat transfer coefficient on the temperature profiles. The different level of pressures achievable for a given angular velocity depends on whether the bearing is thermal or nonisothermal; upwind differencing is noted to be essential for the derivation of a realistic profile.

Braun, M. J.↗
NASA NTRS1977

Integration-free interval doubling for Riccati equation solutions

Starting with certain identities obtained by Reid (1972) and Redheffer (1962) for general matrix Riccati equations (RE's), we give various algorithms for the case of constant coefficients. The algorithms are based on two ideas - first, relate the RE solution with general initial conditions to anchored RE solutions; and second, when the coefficients are constant, the anchored solutions have a basic shift-invariance property. These ideas are used to construct an integration-free, superlinearly convergent iterative solution to the algebraic RE. Preliminary numerical experiments show that our algorithms, arranged in square-root form, provide a method that is numerically stable and appears to be competitive with other methods of solving the algebraic RE.

Sidhu, G. S.↗
NASA NTRS1995

Existence of frozen-in coordinate systems

The 'frozen-in' coordinate systems were first introduced in the works on 'reconnection' and 'magnetic barrier' theories (see review by M.l.Pudovkin and V.S.Semenov, Space Sci. Rev. 41,1 1985). The idea was to utilize the mathematical apparatus developed for 'general relativity' theory to simplify obtaining solutions to the ideal MHD equations set. Magnetic field (B), plasma velocity (v), and their vector product were used as coordinate vectors. But there exist no stationary solutions of ideal MHD set that satisfies the required boundary conditions at infinity (A.D.Chertkov, Solar Wind Seven Conf.,Pergamon Press,1992,165) having non-zero vector product of v and B where v and B originate from the same sphere. The existence of a solution is the hidden mine of the mentioned theories. The solution is constructed in the coordinate system, which is unknown and indeterminate before obtaining this solution. A substitution of the final solution must be done directly into the initial MHD set in order to check the method. One can demonstrate that 'solutions' of Petschek's problem, obtained by 'frozen-in' coordinate systems, does not satisfy just the 'frozen-in' equation, i.e. induction equation. It stems from the fact that Petschek's 're-connection' model, treated as a boundary problem, is over determined. This problem was incorrectly formulated.

Chertkov, A. D.↗
DOE OSTI2024

Amplitudes, supersymmetric black hole scattering at $\mathcal{O}\left({G}^5\right)$, and loop integration

We compute the potential-graviton contribution to the scattering amplitude, the radial action, and the scattering angle of two extremal black holes in $\mathcal{N}$ = 8 supergravity at the fifth post-Minkowskian order and to next-to-leading order in a large mass expansion (first self-force order). Properties of classical unitarity cuts allow us to focus on the integration-by-parts reduction of planar integrals, while nonplanar integrals at this order are obtained from the planar ones by straightforward manipulations. We present the solution to the differential equations for all master integrals necessary to evaluate the classical scattering amplitudes of massive scalar particles at this order in all gravitational theories, in particular in $\mathcal{N}$ = 8 supergravity, and in general relativity. Despite the appearance of higher-weight generalized polylogarithms and elliptic functions in the solution to the differential equation for master integrals, the final supergravity answer is remarkably simple and contains only (harmonic) polylogarithmic functions up to weight 2. The systematic analysis of elliptic integrals discussed here, as well as the particular organization of boundary integrals in $\mathcal{N}$ = 8 observables are independent of supersymmetry and may have wider applications, including to aspects of collider physics.

Black Holes↗
DOE OSTI2022

Numerical Solution of the Steady-State Network Flow Equations for a Non-Ideal Gas

Herein we formulate a steady-state network flow problem for non-ideal gas that relates injection rates and nodal pressures in the network to flows in pipes. For this problem, we present and prove a theorem on uniqueness of generalized solution for a broad class of non-ideal pressure-density relations that satisfy a monotonicity property. Further, we develop a Newton-Raphson algorithm for numerical solution of the steady-state problem, which is made possible by a systematic non-dimensionalization of the equations. The developed algorithm has been extensively tested on benchmark instances and shown to converge robustly to a generalized solution. Previous results [1]-[4], indicate that the steady-state network flow equations for an ideal gas are difficult to solve by the Newton-Raphson method because of its extreme sensitivity to the initial guess. In contrast, we find that non-dimensionalization of the steady-state problem is key to robust convergence of the Newton-Raphson method. We identify criteria based on the uniqueness of solutions under which the existence of a non-physical generalized solution found by a non-linear solver implies non-existence of a physical solution, i.e., infeasibility of the problem. Finally, we compare pressure and flow solutions based on ideal and non-ideal equations of state to demonstrate the need to apply the latter in practice. The solver developed in this article is open-source and is made available for both the academic and research communities as well as the industry.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Magnetic frame-dragging correction to the electromagnetic solution of a compact neutron star

ABSTRACT Neutron stars are usually modelled as spherical, rotating perfect conductors with a predominant intrinsic dipolar magnetic field anchored to their stellar crust. Due to their compactness, General Relativity corrections must be accounted for in Maxwell’s equations, leading to modified interior and exterior electromagnetic solutions. We present analytical solutions for slowly rotating magnetized neutron stars, taking into account the magnetic frame-dragging correction. For typical compactness values, i.e. Rs ∼ 0.5 [R*], we show that the new terms lead to a per cent order correction in the magnetic field orientation and strength compared to the case, with no magnetic frame-dragging correction. Also, we obtain a self-consistent redistribution of the surface azimuthal current. We verify the validity of the derived solution through two-dimensional particle-in-cell simulations of an isolated neutron star. Defining the azimuthal electric and magnetic field amplitudes during the transient phase as observables, we prove that the magnetic frame-dragging correction reduces the transient wave amplitude, as expected from the analytical solution. We show that simulations are more accurate and stable, when we include all first-order terms. The increased accuracy at lower spatiotemporal resolutions translates into a reduction in simulation runtimes.

Torres, R. (ORCID:0000000291820228)↗
DOE OSTI2025

Lectures on statistical mechanics

Presented here is a transcription of the lecture notes from Professor Allan N. Kaufman’s graduate statistical mechanics course Physics 212A and 212B at the University of California Berkeley from the 1972–1973 academic year. 212A addressed equilibrium statistical mechanics with topics: fundamentals (micro-canonical and sub-canonical ensembles, adiabatic law and action conservation, fluctuations, pressure, and virial theorem), classical fluids and other systems (equation of state, deviations from ideality, virial coefficients and van der Waals potential, canonical ensemble and partition function, quasistatic evolution, grand-canonical ensemble and partition function, chemical potential, simple model of a phase transition, quantum virial expansion, numerical simulation of equations of state, and phase transition), chemical equilibrium (systems with multiple species and chemical reactions, law of mass action, Saha equation, chemical equilibrium including ionization and excited states), and long-range interactions (including Coulomb, dipole, and gravitational interactions, Debye–Hückel theory, and shielding). 212B addressed nonequilibrium statistical mechanics with topics: fundamentals (definitions: realizations, moments, characteristic function, and discrete variables), Brownian motion (Langevin equation, fluctuation–dissipation theorem, spatial diffusion, Boltzmann’s H-theorem), Liouville and Klimontovich equations, Landau equation (derivation, elaboration, and H-theorem, and irreversibility), Markov processes and Fokker–Planck equation (derivations of the Fokker–Planck equation and a master equation), linear response and transport theory (linear Boltzmann equation, linear response theory of Kubo and Mori, relation of entropy production to electrical conductivity, transport relations and coefficients, normal mode solutions of the transport equations, sketch of a generalized Langevin equation method for transport theory), and an introduction to nonequilibrium quantum statistical mechanics.

plasma dynamics↗
NASA NTRS1971

Linearization of Schwarzschild's line element - Application to the clock paradox.

This article studies the relativistic theory of the motion of a particle in the presence of a uniform acceleration field. The problem is introduced as a linearization of the fundamental line element of general relativity. The linearized line element is a solution of Einstein's field equations. The equations of geodesics corresponding to this line element are solved and applied to the clock paradox problem.-

Broucke, R.↗
NASA NTRS1995

Electric fields and field-aligned currents in polar regions of the solar corona: 3-D MHD consideration

The calculation of the solar rotation electro-dynamical effects in the near-the-Sun solar wind seems more convenient from the non-inertial corotating reference frame. This implies some modification of the 3-D MHD equations generally on the base of the General Theory of Relativity. The paper deals with the search of stationary (in corotating non-inertial reference frame) solutions of the modified 3-D MHD equations for the in near-the-Sun high latitude sub-alfvenic solar wind. The solution is obtained requiring electric fields and field-aligned electric currents in the high latitude near-the-Sun solar wind. Various scenario are explored self-consistently via a number of numerical experiments. The analogy with the high latitude Earth's magnetosphere is used for the interpretation of the results. Possible observational manifestations are discussed.

Pisanko, Yu. V.↗
NASA NTRS1978

An integral solution to a nonlinear diffusion problem

The generalized Galerkin method (or the method of integral relations) is applied to the type of problem described by quasilinear parabolic equations. As an example the problem of nonlinear transient slab diffusion with a general reservoir boundary condition is worked out. The integral relations are given for an arbitrary number of strips, and solutions using up to seven strips have been obtained in order to investigate the convergence of the method.

Hussaini, M. Y.↗
NASA NTRS1988

Ion-ion waves in the auroral region - Wave excitation and ion heating

The properties of the ion-ion mode which is excited in plasmas when two or more cold ion beams are streaming relative to one another are investigated assuming a warm electron distribution at rest (the model consistent with particle distributions observed in auroral plasma cavities). Numerical solutions are derived for the generalized electrostatic dispersion equation for parallel propagation and for oblique propagation of ion waves. It is shown that the relative ion temperature and concentration have significant effects on the stability of the accelerated plasma and the expected ion heating. Finally, a relationship between ion drift and thermal speed is derived using the marginal stability of ion-ion waves for nonzero ion temperature; the relationship was found to compare favorably with DE-1 ion observations at the high-altitude boundary of the auroral cavity.

Dusenbery, P. B.↗