Search NASA⌕ Search

SEARCH · Search NASA

Results for “relativistic fluid dynamics”

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.

14 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Fluid dynamical description of relativistic nuclear collisions

On the basis of both a conventional relativistic nuclear fluid dynamic model and a two fluid generalization that takes into account the interpenetration of the target and projectile upon contact, collisions between heavy nuclei moving at relativistic speeds are calculated. This is done by solving the relevant equations of motion numerically in three spatial dimensions by use of particle in cell finite difference computing techniques. The effect of incorporating a density isomer, or quasistable state, in the nuclear equation of state at three times normal nuclear density, and the effect of doubling the nuclear compressibility coefficient are studied. For the reaction 20Ne + 238U at a laboratory bombarding energy per nucleon of 393 MeV, the calculated distributions in energy and angle of outgoing charged particles are compared with recent experimental data both integrated over all impact parameters and for nearly central collisions.

Nix, J. R.↗
NASA NTRS1994

Cosmic bullets as particle accelerators and radio sources

We have simulated in two dimensions the dynamical evolution of dense gas clouds(`cosmic bullets') moving supersonically through a uniform low-density medium. The diffusive shock acceleration of relativistic protons (cosmic rays) and their dynamical feedback on the background flow are included by the two-fluid model for this process. The acceleration of relativistic electrons is approximated by a test-particle model, and a passive magnetic field is followed by a simple advection scheme. Strong bow shocks, with Mach numbers similar to that of a bullet's motion, are the most important particle accelerators in the flow, while tail shocks and shocks inside the bullets do not play generally significant roles in this regard. For our simulation parameters, approximately greater than 10% of the initial bullet kinetic energy is converted to a combination of internal energy of gas and cosmic-ray protons by the time the bullets begin to be disrupted. Characteristically, the cosmic rays gain several percent of the available kinetic energy. Bullet destruction on timescales only a little larger than the ram pressure bullet crushing time begins in response to Kelvin-Helmholtz and especially to Rayleigh-Taylor instabilities along the forward bullet surface. For dense bullets this happens before the bullet is stopped by ram pressure. According to our simple model for synchrotron emission from relativistic electrons accelerated and transported within the flows, that emission increases rapidly as the bullet begins to fragment, when it is strongly dominated by field enhancement in sheared flows. Synchrotron emission from the acceleration region within the bow shock is, by contrast, much weaker.

Jones, T. W.↗
NASA NTRS2003

Finite Element Method for Capturing Ultra-relativistic Shocks

While finite element methods are used extensively by researchers solving computational fluid dynamics in fields other than astrophysics, their use in astrophysical fluid simulations has been predominantly overlooked. Current simulations using other methods such as finite difference and finite volume (based on finite difference) have shown remarkable results, but these methods are limited by their fundamental properties in aspects that are important for simulations with complex geometries and widely varying spatial and temporal scale differences. We have explored the use of finite element methods for astrophysical fluids in order to establish the validity of using such methods in astrophysical environments. We present our numerical technique applied to solving ultra-relativistic (Lorentz Factor Gamma >> 1) shocks which are prevalent in astrophysical studies including relativistic jets and gamma-ray burst studies. We show our finite element formulation applied to simulations where the Lorentz factor ranges up to 2236 and demonstrate its stability in solving ultra-relativistic flows. Our numerical method is based on the Flowfield Dependent Variation (FDV) Method, unique in that numerical diffusion is derived from physical parameters rather than traditional artificial viscosity methods. Numerical instabilities account for most of the difficulties when capturing shocks in this regime. Our method results in stable solutions and accurate results as compared with other methods.

Richardson, G. A.↗
NASA NTRS1976

Collected Works of Robert T. Jones

Fluid dynamic principles are applied to airfoil stability and control problems, blood pump studies, and in relativistic kinematics.

Aeronautics↗
NASA NTRS1980

Spontaneous formation of knots in relativistic flows - A model for variability in compact synchrotron sources

A model is proposed in which regions of enhanced synchrotron emission are formed via radiative thermal instabilities in relativistic flows. A perturbed volume in which the magnetic field is enhanced by a few tens of a percent relative to the steady value collapses due to the increased cooling rate and consequent loss in pressure. The collapse further increases the magnetic field, leading to even shorter radiative lifetimes. The instability progresses in this way until either the external electron pressure is balanced by the perturbed magnetic pressure or the overall expansion of the flow becomes important. A fluid-dynamical treatment of the instability is performed, and the results are shown to apply to flux variations and formation of knots in the relativistic flows which are thought to occur in quasars and active galactic nuclei. The knots possess the same kinematic properties as the collimated flow and can thus be responsible for the apparent superluminal motions observed in compact radio sources.

Marscher, A. P.↗
NASA NTRS1994

Structure of relativistic shocks in pulsar winds: A model of the wisps in the Crab Nebula

We propose a model of a optical 'wisps' of the Crab Nebula, features observed in the nebular synchrotron surface brightness near the central pulsar, as manifestations of the internal structure of the shock terminating the pulsar wind. We assume that this wind is composed of ions and a much denser plasma of electrons and positrons, frozen together to a toroidal magnetic field and flowing relativistically. We construct a form of solitary wave model of the shock structure in which we self-consistently solve for the ion orbits and the dynamics of the relativistically hot, magnetized e(+/-) background flow. We ignore dispersion in the ion energies, and we treat the pairs as an adiabatic fluid. The synchrotron emission enhancements, observed as the wisps, are then explained as the regions where reflection of the ions in the self-consistent magnetic field causes compressions of the e(+/-).

Gallant, Yves A.↗
NASA NTRS1990

Dynamical effects of annihilation in pair-dominated winds

The steady, spherically symmetric flow of an ideal fluid dominated by photons and ultrarelativistic electron-positron pairs is analyzed. A new wind equation and a set of critical point conditions are obtained which describe the relativistic flow of an annihilation gas in which the flow velocity exceeds the diffusion velocity of the photons. Numerical results are reported which suggest the possible existence of trapped, pure-pair winds driven by the combined pressure of the pairs and the photons. Most of the annihilation occurs below the critical radius in trapped flows, and a substantial fraction of the total energy of the injected pairs is converted into kinetic energy and radiation. Accurate numerical solutions for the flow velocity and the positron loss rate in optically thin, Newtonian winds are obtained, and a useful approximate analytic relation between the positron loss rate and the flow velocity is derived which suggests that a large number of pairs may survive the annihilation region, ultimately escaping the potential well.

Becker, Peter A.↗
NASA NTRS2003

Magnetohydrodynamic Jump Conditions for Oblique Relativistic Shocks with Gyrotropic Pressure

Shock jump conditions, i.e., the specification of the downstream parameters of the gas in terms of the upstream parameters, are obtained for steady-state, plane shocks with oblique magnetic fields and arbitrary flow speeds. This is done by combining the continuity of particle number flux and the electromagnetic boundary conditions at the shock with the magnetohydrodynamic conservation laws derived from the stress-energy tensor. For ultrarelativistic and nonrelativistic shocks, the jump conditions may be solved analytically. For mildly relativistic shocks, analytic solutions are obtained for isotropic pressure using an approximation for the adiabatic index that is valid in high sonic Mach number cases. Examples assuming isotropic pressure illustrate how the shock compression ratio depends on the shock speed and obliquity. In the more general case of gyrotropic pressure, the jump conditions cannot be solved analytically with- out additional assumptions, and the effects of gyrotropic pressure are investigated by parameterizing the distribution of pressure parallel and perpendicular to the magnetic field. Our numerical solutions reveal that relatively small departures from isotropy (e.g., approximately 20%) produce significant changes in the shock compression ratio, r , at all shock Lorentz factors, including ultrarelativistic ones, where an analytic solution with gyrotropic pressure is obtained. In particular, either dynamically important fields or significant pressure anisotropies can incur marked departures from the canonical gas dynamic value of r = 3 for a shocked ultrarelativistic flow and this may impact models of particle acceleration in gamma-ray bursts and other environments where relativistic shocks are inferred. The jump conditions presented apply directly to test-particle acceleration, and will facilitate future self-consistent numerical modeling of particle acceleration at oblique, relativistic shocks; such models include the modification of the fluid velocity profile due to the contribution of energetic particles to the momentum and energy fluxes.

Double, Glen P.↗
NASA NTRS1985

Nonlinear shock acceleration. III - Finite wave velocity, wave pressure, and entropy generation via wave damping

The nonlinear theory of shock acceleration developed in earlier papers, which treated the waves as being completely frozen into the fluid, is generalized to include wave dynamics. In the limit where damping keeps the wave amplitude small, it is found that a finite phase velocity (V sub ph) of the scattering waves through the background fluid, tempers the acceleration generated by high Mach number shocks. Asymptotic spectra proportional to 1/E sq are possible only when the ratio of wave velocity to shock velocity is less than 0.13. For a given asymptotic spectrum, the efficiency of relativistic particle production is found to be practically independent of the value of V sub ph, so that earlier results concerning its value remain valid for finite V sub ph. In the limit where there is no wave damping, it is shown that for modest Alfven Mach numbers, approximately greater than 4 and less than 6, the magnetic field is amplified by the energetic particles to the point of being in rough equipartition with them, as models of synchrotron emission frequently take the field to be. In this case, the disordering and amplification of field energy may play a major role in the shock transition.

Eichler, D.↗
NASA NTRS1994

Low-frequency modes and nonbarotropic effects in pseudo-Newtonian accretion disks

A recently developed formalism is used to reexamine the question of the existence of hydrodynamical modes that pulsate with very low frequencies in the inner regions of accretion disks. The formalism is valid in an exact sense for the adiabatic pulsations of rotating Newtonian fluids that are generally nonbarotropic (such as those with 'nonadiabatic temperature gradients,' for example), and hence its application in the present context represents an improvement over previous analyses that are more approximate. The formalism is applied to thin non-self-gravitating disks, with the gravitational potential of the central source modified in the usual way in order to simulate relativistic effects. In the barotropic limit, the analyses indicate that in many cases nearly Keplerian disks exhibit nonaxisymmetric modes of pulsation that are trapped in the inner disk regions, with pulsation periods much longer than the dynamical timescale. These results are similar to those of earlier calculations that assume disks pulsate without changing the temperature distribution. A method is developed for including lowest order nonbarotropic effects. Previous analyses have been incapable of accurately treating the nonbarotropic regime. The application of the present method to the low-frequency modes reveals that, due to unexpected cancellations among terms, the nonbarotropic correction to the pusation frequency omega is only of order tilde-omega(sub BV exp 2) omega, where tilde-omega(sub BV) is the appropriate dimensionless Brunt-Vaisala frequency. This correction is much smaller than the expected correction of order tilde-omega(sub BV) Omega, where Omega is the rotation angular velocity. The important conclusion drawn from this is that nonbarotropic corrections are generally small and hence that low-frequency modes persist into the nonbarotropic regime. For disk temperatures appropriate to X-ray emission, the adiabatic frequencies of trapped modes are of the same order as the frequencies associated with quasi-periodic oscillations (QPOs) of certain low-mass X-ray binaries and other X-ray sources. It is suggested that further consideration be given to the possiblity that modulation of emission processes by these modes can explain at least some of the QPO phenomena.

Ipser, James R.↗
NASA NTRS1996

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