Engineering topics

Webb, G. M.

Publications and source records attributed to Webb, G. M..

At least 19 records

NASA NTRS1995

Symmetries of the TDNLS equations for weakly nonlinear dispersive MHD waves

In this paper we consider the symmetries and conservation laws for the TDNLS equations derived by Hada (1993) and Brio, Hunter and Johnson, to describe the propagation of weakly nonlinear dispersive MHD waves in beta approximately 1 plasmas. The equations describe the interaction of the Alfven and magnetoacoustic modes near the triple umbilic, where the fast magnetosonic, slow magnetosonic and Alfven speeds coincide and a(g)(exp 2) = V(A)(exp 2) where a(g) is the gas sound speed and V(A) is the Alfven speed. We discuss Lagrangian and Hamiltonian formulations, and similarity solutions for the equations.

Webb, G. M.↗

NASA NTRS1995

Multi-dimensional MHD simple waves

In this paper we consider a formalism for multi-dimensional simple MHD waves using ideas developed by Boillat. For simple wave solutions one assumes that all the physical variables (the density rho, gas pressure p, fluid velocity V, gas entropy S, and magnetic induction B in the MHD case) depend on a single phase function phi(r,t). The simple wave solution ansatz and the MHD equations then require that the phase function has the form phi = r x n(phi) - lambda(phi)t, where = n(phi) = Delta phi / (absolute value of Delta phi) is the wave normal and lambda(phi) = omega/k = -phi t / (absolute value of Delta phi) is the normal speed of the wave front. The formalism allows for more general simple waves than that usually dealt with in which n(phi) is a constant unit vector that does not vary along the wave front. The formalism has implications for shock formation for multi-dimensional waves.

Webb, G. M.↗

NASA NTRS1995

The modulational instability for the TDNLS equations for weakly nonlinear dispersive MHD waves

In this paper we study the modulational instability for the TDNLS equations derived by Hada (1993) and Brio, Hunter, and Johnson to describe the propagation of weakly nonlinear dispersive MHD waves in beta approximately 1 plasmas. We employ Whitham's averaged Lagrangian method to study the modulational instability. This complements studies of the modulational instability by Hada (1993) and Hollweg (1994), who did not use the averaged Lagrangian approach.

Webb, G. M.↗

NASA NTRS1995

Contact and pressure balance structures in two-fluid cosmic-ray hydrodynamics

The role of cosmic-ray-modified contact discontinuities and pressure balance structures in two-fluid cosmic-ray hydrodynamics in one Cartesian space dimension are investigated by means of analytic and numerical solution examples, as well as by weakly nonlinear asymptotics. The fundamental wave modes of the two-fluid cosmic-ray hydrodynamic equations in the long-wavelength limit consist of the backward and forward propagating cosmic-ray-modified sound waves, with sound speed dependent on both the cosmic-ray and thermal gas pressures; the contact discontinuity; and a pressure balance mode in which the sum ofthe cosmic ray and thermal gas pressure perturbations is zero. The pressure balance mode, like the contact discontinuity is advected with the background flow. The interaction of the pressure balance mode with the contact discontinuity is investigated by means of the method of multiple scales. The thermal gas and cosmic-ray pressure perturbations satisfy a linear diffusion equation, and entropy perturbations arising from nonisentropic initial conditions for the thermal gas are frozen into the fluid. The contact discontinuity and pressure balance eigenmodes both admit nonzero perturbations in the thermal gas, whereas the cosmic-ray-modified sound waves are isentropic. The total entropy perturbation is shared between the contact discontinuity and pressure balance eigenmodes, and examples are given in which there is a transfer of entropy between the two modes. In particular, N-wave type density disturbances are obtained which arise as a result of the entropy transfer between the two modes. A weakly nonlinear geometric optics perturbation expansion is used to study the long timescale evolution of the short-wavelength entropy wave and the thermal gas sound waves in a slowly varying, large-scale background flow. The weakly nonlinear geometric optics expansion is also used to generalize previous studies of squeezing instability for short-wavelength sound waves in the two fluid model, by including a weakly nonlinear wave steepening term that leads to shock formation, as well as the effect of long time and space dependence of the background flow. Implications of cosmic-ray-modified pressure balance structures and contact discontinuities in models of the interaction of traveling interplanetary shocks and compression and rarefraction waves with the solar wind termination shock are briefly discussed.

Webb, G. M.↗

NASA NTRS1995

The termination shock: Physical processes

A brief overview of the physics thought to determine the properties of the heliospheric termination shock is given and the central role of hot interstellar pick-up ions and cosmic rays is emphasized.

Zank, G. P.↗

NASA NTRS1994

Green's formula and variational principles for cosmic-ray transport with application to rotating and shearing flows

Green's theorem and Green's formula for the diffusive cosmic-ray transport equation in relativistic flows are derived. Green's formula gives the solution of the transport equation in terms of the Green's function of the adjoint transport equation, and in terms of distributed sources throughout the region R of interest, plus terms involving the particle intensity and streaming on the boundary. The adjoint transport equation describes the time-reversed particle transport. An Euler-Lagrange variational principle is then obtained for both the mean scattering frame distribution function f, and its adjoint f(dagger). Variations of the variational functional with respect to f(dagger) yield the transport equation, whereas variations of f yield the adjoint transport equation. The variational principle, when combined with Noether's theorem, yields the conservation law associated with Green's theorem. An investigation of the transport equation for steady, azimuthal, rotating flows suggests the introduction of a new independent variable H to replace the comoving frame momentum variable p'. For the case of rigid rotating flows, H is conserved and is shown to be analogous to the Hamiltonian for a bead on a rigidly rotating wire. The variable H corresponds to a balance between the centrifugal force and the particle inertia in the rotating frame. The physical interpretation of H includes a discussion of nonrelativistic and special relativistic rotating flows as well as the cases of aziuthal, differentially rotating flows about Schwarzs-child and Kerr black holes. Green's formula is then applied to the problem of the acceleration of ultra-high-energy cosmic rays by galactic rotation. The model for galactic rotation assumes an angular velocity law Omega = Omega(sub 0)(omega(sub 0)/omega), where omega denotes radial distance from the axis of rotation. Green's functions for the galactic rotation problem are used to investigate the spectrum of accelerated particles arising from monoenergetic and truncated power-law sources. We conclude that it is possible to accelerate particles beyond the knee by galactic rotation, but not in sufficient number to adequately explain the observed spectrum.

Webb, G. M.↗

NASA NTRS1994

Time-dependent evolution of cosmic-ray-modified shock structure: Transition to steady state

Steady state solutions to the two-fluid equations of cosmic-ray-modified shock structure were investigated first by Drury and Volk (1981). Their analysis revealed, among other properties, that there exist regions of upstream parameter space where the equations possess three different downstream solutions for a given upstream state. In this paper we investigate whether or not all these solutions can occur as time-asymptotic states in a physically realistic evolution. To do this, we investigate the time-dependent evolution of the two-fluid cosmic-ray equations in going from a specified initial condition to a steady state. Our results indicate that the time-asymptotic solution is strictly single-valued, and it undergoes a transition from weakly to strongly cosmic-ray-modified at a critical value of the upstream cosmic ray energy density. The expansion of supernova remnant shocks is considered as an example, and it is shown that the strong to weak transition is in fact more likely. The third intermediate solution is shown to influence the time-dependent evolution of the shock, but it is not found to be a stable time-asymptotic state. Timescales for convergence to these states and their implications for the efficiency of shock acceleration are considered. We also investigate the effects of a recently introduced model for the injection of seed particles into the shock accelerated cosmic-ray population. The injection is found to result in a more strongly cosmic-ray-dominated shock, which supports our conclusion that for most classes of intermediate and strong cosmic-ray-modified shocks, the downstream cosmic-ray pressure component is at least as large as the thermal gas pressure, independent of the upstream state. As a result, cosmic rays almost always play a significant role in determining the shock structure and dissipation and they cannot be regarded as test particles.

Donohue, D. J.↗

NASA NTRS1993

Particle injection and the structure of energetic-particle-modified shocks

A macroscopic 'self-consistent' nonlinear two-fluid model is developed for energetic-particle-modeled shocks. The model incorporates particle exchange between a thermal gas and an energetic population via a 'thermal leakage' mechanism. In the model, injection is regarded as a process whereby particles in the cooler thermal core of a total particle distribution are energized sufficiently to cross a 'momentum boundary' into the wings of the distribution. Particles in the wings are identified as 'energetic particles' or 'cosmic rays' and are assumed to propagate according to the diffusive transport equation of cosmic ray theory. Thermal particles are energized either as a result of adiabatic compression in a decelerating flow or directly via anomalous heating within a subshock. It is shown that the nature and model of particle injection directly affects the nature and dynamics of the shock, and hence the efficiency of particle acceleration, and that the shock itself regulates particle injection dynamically.

Zank, G. P.↗

NASA NTRS1992

Properties of mass-loading shocks. II - Magnetohydrodynamics

The one-dimensional magnetohydrodynamics of mass-loading shocks is examined. These shocks, which are distinct from MHD shocks of classical nonreacting fluid dynamics and of combustion theory and which are characterized by the addition of mass within the shock transition, are to be found at comets and, depending upon circumstances, at nonmagnetized and weakly magnetized planets such as Venus and Mars. A completely general mass-loading form of the Hugoniot equation is derived, and some of the most important differences between mass-loading and nonreacting classical MHD shocks are identified. Two new types of MHD shocks are described which have no classical MHD analogues.

Zank, G. P.↗

NASA NTRS1992

Wave diffraction in weak cosmic-ray-modified shocks

Weakly multidirectional, long-wavelength cosmic-ray-modified shocks are studied via multiple scale perturbation techniques. The effects of diffraction are discussed in terms of Green's function solutions of the linearized 1 + 3D Burgers and 1 + 3D KdVB equations, and also in terms of solutions with singular Dirac delta initial distributions. The solutions show a monotonic decrease of the wave-front curvature with increasing time owing to the effects of wave diffraction. The shape of the wave surface is discussed in terms of solutions S to the wave eikonal equation corresponding to singular initial conditions. For the fast magnetosonic wave propagating in the positive x-direction, the wave phase surface S = 0 has elliptic cross sections with the planes x = constant and has a convex paraboloidal shape. Plane-wave solutions of the 1 + 3D KdVB equation are discussed.

Webb, G. M.↗

NASA NTRS1992

Energetic particle transport in relativistic flows

A discussion is given of pitch angle dependent and diffusive transport equations for cosmic rays applicable for both special relativistic and general relativistic flows, derived from the relativistic Boltzmann equation. As an example of particle transport in a curved spacetime we give a pitch angle dependent transport equation appropriate for radial accretion onto a Schwarzschild black hole. The roles of fluid shear, acceleration and compression on the energy changes of particles in the diffusive transport equations are emphasized. Also discussed are special flows (e.g. rigidly rotating flows) associated with a Killing vector for which a constant of the motion of the particles can be identified, and for which simplified transport equations can be constructed.

Webb, G. M.↗

NASA NTRS1991

Mass-loading and parallel magnetized shocks

Recent observations at comets Giacobini-Zinner and Halley suggest that simple nonreacting gas dynamics or MHD is an inappropriate description for the bow shock. The thickness of the observed (sub)shock implies that mass-loading is an important dynamical process within the shock itself, thereby requiring that the Rankine-Hugoniot conditions possess source terms. This leads to shocks with properties similar to those of combustion shocks. The paper considers parallel magnetized shocks subjected to mass-loading, describes some properties which distinguish them from classical MHD parallel shocks, and establishes the existence of a new kind of MHD compound shock. These results will be of importance both to observations and numerical simulations of the comet-solar wind interaction.

Zank, G. P.↗

NASA NTRS1990

Diffusive shock acceleration of energetic electrons and the method of generations

The diffusive shock acceleration of energetic electrons at a plane shock is studied by means of the method of generations developed by Petukhov et al. (1986), where the particles are binned according to the number of particle acceleration cycles. The method emphasizes the cyclic character of the shock acceleration process, and takes into account the particle energy increment per cycle, the energy loss per cycle, and the escape probability of particles downstream of the shock. In the limit of a small fractional momentum gain per cycle, the method reproduces the known steady-state semianalytic solution for constant diffusion coefficients on both sides of the shock. It is demonstrated that the method may also be applied to time-dependent shock-acceleration problems.

Webb, G. M.↗

NASA NTRS1990

Application of the sine-Poisson equation in solar magnetostatics

Solutions of the sine-Poisson equations are used to construct a class of isothermal magnetostatic atmospheres, with one ignorable coordinate corresponding to a uniform gravitational field in a plane geometry. The distributed current in the model (j) is directed along the x-axis, where x is the horizontal ignorable coordinate; (j) varies as the sine of the magnetostatic potential and falls off exponentially with distance vertical to the base with an e-folding distance equal to the gravitational scale height. Solutions for the magnetostatic potential A corresponding to the one-soliton, two-soliton, and breather solutions of the sine-Gordon equation are studied. Depending on the values of the free parameters in the soliton solutions, horizontally periodic magnetostatic structures are obtained possessing either a single X-type neutral point, multiple neural X-points, or solutions without X-points.

Webb, G. M.↗

NASA NTRS1989

On the stability of a class of isothermal, magnetostatic atmospheres

The stability of a class of isothermal magnetostatic atmospheres with one ignorable coordinate corresponding to a uniform gravitational field in a plane geometry is explored. Special attention is given to the stability of nonpotential equilibria that have the same magnetic field geometry as the potential equilibrium of the solution family. It is found that equilibria in which the source currents for the potential field are buried at large distances below the photospheric base are stable. Also considered is the stability of configurations in which the source currents for the potential field are located at a finite depth below the base.

Webb, G. M.↗

NASA NTRS1989

On wave stability in relativistic cosmic-ray hydrodynamics

Wave stability of a two-fluid hydrodynamical model describing the acceleration of cosmic rays by the first-order Fermi mechanism in relativistic, cosmic-ray-modified shocks is investigated. For a uniform background state, the short- and long-wavelength wave speeds are shown to interlace, thus assuring wave stability in this case. A JWKB analysis is performed to investigate the stability of short-wavelength thermal gas sound waves in the smooth, decelerating supersonic flow upstream of a relativistic, cosmic-ray-modified shock. The stability of the waves is assessed both in terms of the fluid velocity and density perturbations, as well as in terms of the wave action. The stability and interaction of the short-wavelength cosmic-ray coherent mode with the background flow is also studied.

Webb, G. M.↗

NASA NTRS1989

The diffusion approximation and transport theory for cosmic rays in relativistic flows

Equations describing the transport of cosmic rays in relativistic flows in the diffusion approximation are obtained. The analysis is based on the zeroth, first, and second differential moment equations of the relativistic Boltzmann equation with a BGK collision term. A perturbation solution of the moment equations in the diffusion approximation yields both the co-moving frame particle current and viscous stresses. The resultant cosmic-ray continuity equation contains three readily recognized energy change terms: the adiabatic energy change term; the viscous shear energy change term; and a term proportional to the scalar product of the acceleration vector of the scattering frame and the heat flux.

Webb, G. M.↗

NASA NTRS1988

Cosmic-ray hydrodynamics at relativistic shocks

A test particle description of first-order Fermi acceleration of cosmic rays at plane-parallel relativistic shocks is presented in which the pitch angle distribution is obtained explicitly and used to determine the pressure anisotropy or variable Eddington factor as a function of position. The Eddington factor is in turn used to find the hydrodynamical cosmic-ray sound speed, taking the anisotropy into account. For upstream flows speeds above c/sq rt 3, it is shown that the cosmic-ray sound speed exceeds the shock speed in a region just upstream of the shock and decreases montonically with distance further upstream, reaching the value c/sq rt 3 far from the shock. Thus, the shock succeeds in maintaining an anisotropy in the cosmic-ray distribution which is sufficient to keep the flow subsonic with respect to the cosmic-ray sound speed both downstream and in a finite region upstream. The astrophysical implications of these findings are discussed.

Kirk, J. G.↗