Engineering topics
Publications and source records attributed to Jones, Frank C..
We consider the possibility that distinctive features of the local cosmic ray spectra and composition are influenced by the Solar system being embedded within the cavity of an ancient superbubble. Shifts in the measured cosmic ray composition between 10(exp 11) and 10(exp 20) eV as well as the "knee" and "second knee" may be understood in this picture.
Measuring the chemical composition of the cosmic rays in the energy region of greater than or equal to 10(exp 12)eV would be highly useful in settling several nagging questions concerning the propagation of cosmic rays in the galaxy. In particular an accurate measurement of secondary to primary ratios such as Boron to Carbon would gibe clear evidence as to whether the propagation of cosmic rays is determined by a diffusion coefficient that varies with the particle's energy as E(sup 0.5) or E(sup 0.3). This would go a long ways in helping us to understand the anistropy (or lack thereof) of the highest energy cosmic rays and the power requirements for producing those particles at approximately equal to 10(exp 18) eV which are believed to be highest energy particles produced in the Galaxy. This would be only one of the benefits of a mission such as ACCESS to perform direct particle measurements on very high energy cosmic rays.
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.
In reviewing this book I am going to say some harsh things about it. From this one might think that I did not like the book. In fact, one would be wrong to think so. For this reason I wish to state up front that I did like the book. It was an ambitious undertaking and I think that the author has, for the most part, pulled it off. This extremely broad and rich book will allow those non-specialists with a sufficiently long attention span to obtain a good understanding of what cosmic rays are and what cosmic ray research is all about. I have pointed out things that bothered me because I believe that the book is good enough for a second edition and I hope that these matters can be addressed at that time.
There has recently been interest in the role of inverse bremsstrahlung, the emission of photons by fast suprathermal ions in collisions with ambient electrons possessing relatively low velocities, in tenuous plasmas in various astrophysical contexts. This follows a long hiatus in the application of suprathermal ion bremsstrahlung to astrophysical models since the early 1970s. The potential importance of inverse bremsstrahlung relative to normal bremsstrahlung, i.e. where ions are at rest, hinges upon the underlying velocity distributions of the interacting species. In this paper, we identify the conditions under which the inverse bremsstrahlung emissivity is significant relative to that for normal bremsstrahlung in shocked astrophysical plasmas. We determine that, since both observational and theoretical evidence favors electron temperatures almost comparable to, and certainly not very deficient relative to proton temperatures in shocked plasmas, these environments generally render inverse bremsstrahlung at best a minor contributor to the overall emission. Hence inverse bremsstrahlung can be safely neglected in most models invoking shock acceleration in discrete sources such as supernova remnants. However, on scales approximately > 100 pc distant from these sources, Coulomb collisional losses can deplete the cosmic ray electrons, rendering inverse bremsstrahlung, and perhaps bremsstrahlung from knock-on electrons, possibly detectable.
We have modeled the injection and acceleration of pickup ions at the solar wind termination shock and investigated the parameters needed to produce the observed Anomalous Cosmic Ray (ACR) fluxes. A non-linear Monte Carlo technique was employed, which in effect solves the Boltzmann equation and is not restricted to near-isotropic particle distribution functions. This technique models the injection of thermal and pickup ions, the acceleration of these ions, and the determination of the shock structure under the influence of the accelerated ions. The essential effects of injection are treated in a mostly self-consistent manner, including effects from shock obliquity, cross- field diffusion, and pitch-angle scattering. Using recent determinations of pickup ion densities, we are able to match the absolute flux of hydrogen in the ACRs by assuming that pickup ion scattering mean free paths, at the termination shock, are much less than an AU and that modestly strong cross-field diffusion occurs. Simultaneously, we match the flux ratios He(+)/H(+) or O(+)/H(+) to within a factor approx. 5. If the conditions of strong scattering apply, no pre-termination-shock injection phase is required and the injection and acceleration of pickup ions at the termination shock is totally analogous to the injection and acceleration of ions at highly oblique interplanetary shocks recently observed by the Ulysses spacecraft. The fact that ACR fluxes can be modeled with standard shock assumptions suggests that the much-discussed "injection problem" for highly oblique shocks stems from incomplete (either mathematical or computer) modeling of these shocks rather than from any actual difficulty shocks may have in injecting and accelerating thermal or quasi-thermal particles.
The Fermi shock acceleration mechanism may be responsible for the production of high-energy cosmic rays in a wide variety of environments. Modeling of this phenomenon has largely focused on plane-parallel shocks, and one of the most promising techniques for its study is the Monte Carlo simulation of particle transport in shocked fluid flows. One of the principal problems in shock acceleration theory is the mechanism and efficiency of injection of particles from the thermal gas into the accelerated population. The Monte Carlo technique is ideally suited to addressing the injection problem directly, and previous applications of it to the quasi-parallel Earth bow shock led to very successful modeling of proton and heavy ion spectra, as well as other observed quantities. Recently this technique has been extended to oblique shock geometries, in which the upstream magnetic field makes a significant angle Theta(sub B1) to the shock normal. Spectral resutls from test particle Monte Carlo simulations of cosmic-ray acceleration at oblique, nonrelativistic shocks are presented. The results show that low Mach number shocks have injection efficiencies that are relatively insensitive to (though not independent of) the shock obliquity, but that there is a dramatic drop in efficiency for shocks of Mach number 30 or more as the obliquity increases above 15 deg. Cosmic-ray distributions just upstream of the shock reveal prominent bumps at energies below the thermal peak; these disappear far upstream but might be observable features close to astrophysical shocks.
We discuss the fundamental ideas of particle acceleration in plasma shocks with emphasis on those features that are required to produce the 'universal' power-law spectrum. We compare shock acceleration with the more familiar second-order or stochastic acceleration and see that they are not too different in many respects. We discuss the features of shock acceleration that make it appealing and some of its problems as well.
Important problems are discussed concerning the effects of shock geometry on the injection of thermal particles and the process of standard Fermi particle acceleration in general. In pursuit of a self-consistent model that simultaneously describes shock structure and particle acceleration in shocks of arbitrary obliquity, a Monte Carlo simulation developed for parallel shocks is generalized to address oblique geometries. Attention is given to initial results concerning the ways in which injection and acceleration efficiency varies with obliquity and Mach number.
Ways in which energy change terms arise in the transport equation and how the various terms relate to the modes of energy exchange between the particles and plasma are shown. It is argued that the transport equation cannot be used to describe the initial acceleration of thermal particles by plasma shocks or relativistic shocks where the energetic particle speeds are never much greater than the flow speeds. In most other situations, it describes almost any acceleration process that can be caused by a moving plasma. It describes shock acceleration for both parallel shocks and oblique ones, and stochastic acceleration by the turbulent motion of the scatterers as well as by their motion across the magnetic field.
Previous applications of the Monte Carlo technique at the quasiparallel earth bow shock has motivated the extension of this technique to oblique shock geometries typical of those found in most astrophysical shock environments. In addition, such a generalization will permit the thorough examination of theoretical predictions of rapid acceleration times at quasi-perpendicular shocks. Therefore, we have embarked on the modification of our existing Monte Carlo code and in this paper outline the major technical aspects involved in developing a simulation of cosmic-ray acceleration at modified oblique nonrelativistic shocks.
The history and theory of shock acceleration is reviewed, paying particular attention to theories of parallel shocks which include the backreaction of accelerated particles on the shock structure. The work that computer simulations, both plasma and Monte Carlo, are playing in revealing how thermal ions interact with shocks and how particle acceleration appears to be an inevitable and necessary part of the basic plasma physics that governs collisionless shocks is discussed. Some of the outstanding problems that still confront theorists and observers in this field are described.
The weighted-slab method is modified so that, although it is still not exact, it gives the 'best' approximation in a minimization sense when energy loss cannot be neglected. In this approximation the species-dependent energy-change term that operates on the path-length distribution is 'averaged' over the slab-model solution for that particular energy and path length.
The dynamical halo of the Galaxy offers a natural explanation for the form of the variation of cosmic-ray path length with energy. The variation above 1 GeV per nucleon can be understood as due to the variation of the diffusion coefficient, and hence the resident time in the galaxy, with energy. The flattening of the curve below 1 GeV per nucleon is seen to mark a transition to a convection dominated regime where the variation of the diffusion coefficient is no longer a determining factor. It is possible that the random motion of the cosmic rays about the galaxy that prevents us from seeing their sources in a clear manner may enable us to extract information about the galaxy at large and learn something about its large scale motions.
Starting from the Boltzmann equation, a transport equation is derived for energetic particles in a moving magnetized plasma in which the scattering centers that keep the particles quasi-isotropic are moving with a velocity that is not necessarily the same as that of the plasma. The scattering is characterized by three very loose constraints: (1) there is a rest frame for each scatterer in which the particles scatter elastically; (2) in this frame the scattering will not disturb an isotropic distribution; and (3) the momentum transfer in an average collision may be described by a tensor operating on the particles original momentum. Since the strength of the scattering is not specified, the derivation should be as valid for plasma microturbulence as for hard-sphere scattering. The results show clearly which phenomena are responsible for tying the particles to the plasma in the transport equation.
Monte Carlo calculations of test particle spectra and acceleration times are presented from first-order Fermi particle acceleration for parallel shocks with arbitrary flow velocities and compression ratios r up to seven, shock velocities u1 up to 0.98c, and injection energies ranging from thermal to highly superthermal. Far above the injection energy, the spectra are well-approximated by a power law and the spectra are always harder than for nonrelativistic shocks. Approximate analytic expression are given for the spectral slope as a function of u1 and r. The acceleration time as a function of particle energy is less than for nonrelativistic shocks by a factor that increases with u1 and is about three for u1 = 0.98c. It is confirmed that the spectrum for pitch-angle diffusion is considerably steeper than for large-angle scattering for the same shock parameters.
The electron heating required if protons scatter elastically in a parallel, collisionless shock is calculated. Near-elastic proton scattering off large amplitude background magnetic field fluctuations might be expected if the waves responsible for the shock dissipation are generated by the fire hose instability. The effects of an electrostatic potential jump in the shock layer are included by assuming that the energy lost by protons in traversing the potential jump is converted into electron thermal pressure. It is found that the electron temperature increase is a strong function of the potential jump. Comparison is made to the parallel shock plasma simulation of Quest (1987).
A puzzling aspect of the structure of collisionless shock waves concerns the existence of a component of magnetic field that lies out of the coplanarity plane in the transition layer. It is shown that the existence of this component may be understood and its value calculated from a simple expression derived from first principles. The value of this field component may be derived without reference to the transformation properties of the fields or the changes in energy that the ions and electrons undergo on traversing the electrostatic potential of the shock interface.