Engineering Papers⌕ Search

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.

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

At least 145 records · Page 8

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↗

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

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

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

Simulating magnetized neutron stars with discontinuous Galerkin methods

Discontinuous Galerkin methods are popular because they can achieve high order where the solution is smooth, because they can capture shocks while needing only nearest-neighbor communication, and because they are relatively easy to formulate on complex meshes. We perform a detailed comparison of various limiting strategies presented in the literature applied to the equations of general relativistic magnetohydrodynamics. We compare the standard minmod /ΛΠ N limiter, the hierarchical limiter of Krivodonova, the simple WENO limiter, the HWENO limiter, and a discontinuous Galerkin-finite-difference hybrid method. The ultimate goal is to understand what limiting strategies are able to robustly simulate magnetized Tolman-Oppenheimer-Volkoff stars without any fine-tuning of parameters. Among the limiters explored in the paper, the only limiting strategy we can endorse is a discontinuous Galerkin-finite-difference hybrid method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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

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

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

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

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

The possible role of MHD waves in heating the solar corona

The possible role of waves in the heating of the solar corona has been investigated. A general dispersion relation has been derived for waves propagating in a homogeneous plasma subject to dissipation by viscosity and thermal conduction. The dissipation mechanisms have been incorporated self-consistently into the equations, and no assumptions about the strength of the damping have been made. Solutions of the sixth-order dispersion relation provide information on how the damping of both slow and fast mode waves depends upon the plasma density, temperature, field strength, and angle of propagation relative to the background magnetic field. We provide a detailed comparison to the standard approach, which is to solve for the wave quantities in the absence of dissipation and then to use these quantities in expressions for the heating due to viscosity and thermal conduction.

Porter, Lisa J.↗

Electrical double layers at the oil/water interface

This review presents the historical development and current status of the theory of the electrical double layer at a liquid/liquid interface. It gives rigorous thermodynamic definitions of all basic concepts related to liquid interfaces and to the electrical double layer. The difference between the surface of a solid electrode and the interface of two immiscible electrolyte solutions (ITIES) is analyzed in connection to their electrical properties. The most important classical relationships for the electrical double layer are presented and critically discussed. The generalized adsorption isotherm is derived. After a short review of the classical Gouy-Chapman and Verwey-Niessen models, more recent developments of the double layer theory are presented. These include effects of variable dielectric permittivity, nonlocal electrostatics, hydration forces, the modified Poisson-Boltzmann equation and the ion-dipole plasma. The relative merits of different theories are estimated by comparing them with computer simulation of the ITIES and electrical double layer. Special attention is given to the structure of ITIES and its variation due to adsorption of ions and amphiphilic molecules.

Non-NASA Center↗

Relativistic Navigation: A Theoretical Foundation

We present a theoretical foundation for relativistic astronomical measurements in curved space-time. In particular, we discuss a new iterative approach for describing the dynamics of an isolated astronomical N-body system in metric theories of gravity. To do this, we generalize the Fock-Chandrasekhar method of the weak-field and slow-motion approximation (WFSMA) and develop a theory of relativistic reference frames (RF's) for a gravitationally bounded many-extended-body problem. In any proper RF constructed in the immediate vicinity of an arbitrary body, the N-body solutions of the gravitational field equations are formally presented as a sum of the Riemann-flat inertial space-time, the gravitational field generated by the body itself, the unperturbed solutions for each body in the system transformed to the coordinates of this proper RF, and the gravitational interaction term. We develop the basic concept of a general WFSMA theory of the celestial RF's applicable to a wide class of metric theories of gravity and an arbitrary model of matter distribution. We apply the proposed method to general relativity. Celestial bodies are described using a perfect fluid model; as such, they possess any number of internal mass and current multipole moments that explicitly characterize their internal structures. The obtained relativistic corrections to the geodetic equations of motion arise because of a coupling of the bodies' multiple moments to the surrounding gravitational field. The resulting relativistic transformations between the different RF's extend the Poincare group to the motion of deformable self-gravitating bodies. Within the present accuracy of astronomical measurements we discuss the properties of the Fermi-normal-like proper RF that is defined in the immediate vicinity of the extended compact bodies. We further generalize the proposed approximation method and include two Eddington parameters (gamma, Beta). This generalized approach was used to derive the relativistic equations of satellite motion in the vicinity of the extended bodies. Anticipating improvements in radio and laser tracking technologies over the next few decades, we apply this method to spacecraft orbit determination. We emphasize the number of feasible relativistic gravity tests that may be performed within the context of the parameterized WFSMA. Based on the planeto-centric equations of motion of a spacecraft around the planet, we suggested a new null test of the Strong Equivalence Principle (SEP). The experiment to measure the corresponding SEP violation effect could be performed with the future Mercury Orbiter mission. We discuss other relativistic effects, including the perihelion advance and the redshift and geodetic precession of the orbiter's orbital plane about Mercury, as well as the possible future implementation of the proposed formalism in software codes developed for solar-system orbit determination. All the important calculations are completely documented, and the references contain an extensive list of cited literature.

Turyshev, Slava G.↗

Solution of the N-body problem with recurrent power series.

Adaptation of Steffensen's method of solution of both the restricted and general three-body problems in terms of power series in time to the solution of the equations of motion of N planets around the sun. The key feature of this adaptation lies in the introduction of N new variables and in the use of new relations which contain the time derivatives and can thus be treated as differential equations.

Broucke, R.↗

Low-Dispersion Scheme for Nonlinear Acoustic Waves in Nonuniform Flow

The linear dispersion-relation-preserving scheme and its boundary conditions have been extended to the nonlinear Euler equations. This allowed computing, a nonuniform flowfield and a nonlinear acoustic wave propagation in such a medium, by the same scheme. By casting all the equations, boundary conditions, and the solution scheme in generalized curvilinear coordinates, the solutions were made possible for non-Cartesian domains and, for the better deployment of the grid points, nonuniform grid step sizes could be used. It has been tested for a number of simple initial-value and periodic-source problems. A simple demonstration of the difference between a linear and nonlinear propagation was conducted. The wall boundary condition, derived from the momentum equations and implemented through a pressure at a ghost point, and the radiation boundary condition, derived from the asymptotic solution to the Euler equations, have proven to be effective for the nonlinear equations and nonuniform flows. The nonreflective characteristic boundary conditions also have shown success but limited to the nonlinear waves in no mean flow, and failed for nonlinear waves in nonuniform flow.

Baysal, Oktay↗

Cold-hot coupled waves in a flowing magnetized plasma

Abstract Nonlinear coupling of cold and hot waves in a flowing magnetized plasma is analyzed with the Vlasov equation. An analytical solution is obtained for cold waves of a small amplitude (weak flow) and a long wavelength. The distribution function is obtained by integrating the kinetic equation along a perturbed phase-space trajectory for a time-varying plasma flow. The kinetic description presents a generalized dispersion relation that involves resonances depending on cold and hot wave dispersions. Coherent fluid motion leads to radiation peaks in addition to the cyclotron harmonics, where the wavenumber of the cold wave determines the peak frequencies. The peaks appear narrow when the wave propagates perpendicular to the time-averaged flow while they become broad due to the Doppler effect when the wave propagates parallel to the flow. Fully kinetic particle-in-cell simulations corroborate the theoretical predictions. The dispersion relation and resulting wave spectra provide information about plasma parameters and flow properties.

Physics↗

Implications of chaos, scale-invariance, and fractal statistics in geology

A set of three nonlinear total differential equations (Lorenz equations) exhibiting deterministic chaos is considered, and it is shown that these equations demonstrate that deterministic equations with deterministic initial conditions can yield stocastic solutions with fractal statistics. The logistic map, fractal distributions, and fragmentation are discussed. It is pointed out that well-defined fractal distributions of earthquakes are found both regionally and globally, and that the general applicability of the fractal relation for seismicity can provide the basis for a quantitative seismic hazard assessment. It is suggested that the governing physics of erosional topography is nonlinear and may be related to a fractal distribution of storms and floods that generate and renew erosional feature such as gullies and drainage systems.

Turcotte, D. L.↗

Open-closed string duality, branes, and topological recursion

We consider matrix models exhibiting open-closed string duality in two-dimensional string theories with various amounts of supersymmetry. In particular, a relationship between matrix models in the β = 2 Wigner-Dyson class and models in the (1 + 2Γ, 2) Altland-Zirnbauer class relates the perturbative solutions of the two systems’ string equations. Point-like operator insertions in the closed string theory are mapped to the topological expansion of the free energy in the open string theory. We compute correlation functions of macroscopic loop operators and FZZT branes in a general topological gravity background. The relationship between the topological recursion of moduli space volumes and branes is discussed by analyzing the Virasoro conditions in the matrix models.

2D Gravity↗