Wayne G. Roberge

Three-Dimensional X-ray Microtomography

The new technique of x-ray microtomography nondestructively generates three-dimensional maps of the x-ray attenuation coefficient inside small samples with approximately 1 percent accuracy and with resolution approaching 1 micrometer. Spatially resolved elemental maps can be produced with synchrotron x-ray sources by scanning samples at energies just above and below characteristic atomic absorption edges. The system consists of a high-resolution imaging x-ray detector and high-speed algorithms for tomographic image reconstruction. The design and operation of the microtomography device are described, and tomographic images that illustrate its performance with both synchrotron and laboratory x-ray sources are presented.

Time-dependent, Multifluid, Magnetohydrodynamic Shock Waves with Grain Dynamics. I. Formulation and Numerical Tests

This is the first in a series of papers on the effects of dust on the formation, propagation, and structure of nonlinear MHD waves and MHD shocks in weakly ionized plasmas. We model the plasma as a system of nine interacting fluids, consisting of the neutral gas, ions, electrons, and six types of grain fluids comprised of very small grains or polycyclic aromatic hydrocarbons and classical grains in different charge states. We formulate the governing equations for perpendicular shocks under approximations appropriate for dense molecular clouds, protostellar cores, and protoplanetary disks. We describe a code that obtains numerical solutions using a finite difference method, and we establish its accuracy by comparing numerical and exact solutions for special cases.

MOLECULAR LINE EMISSION FROM MULTIFLUID SHOCK WAVES. I. NUMERICAL METHODS AND BENCHMARK TESTS

We describe a numerical scheme for studying time-dependent, multifluid, magnetohydrodynamic shock waves in weakly ionized interstellar clouds and cores. Shocks are modeled as propagating perpendicular to the magnetic field and consist of a neutral molecular fluid plus a fluid of ions and electrons. The scheme is based on operator splitting, wherein time integration of the governing equations is split into separate parts. In one part, independent homogeneous Riemann problems for the two fluids are solved using Godunovs method. In the other, equations containing the source terms for transfer of mass, momentum, and energy between the fluids are integrated using standard numerical techniques. We show that, for the frequent case where the thermal pressures of the ions and electrons are magnetic pressure, the Riemann problems for the neutral and ion-electron fluids have a similar mathematical structure which facilitates numerical coding. Implementation of the scheme is discussed and several benchmark tests confirming its accuracy are presented, including (1) MHD wave packets ranging over orders of magnitude in length- and timescales, (2) early evolution of multifluid shocks caused by two colliding clouds, and (3) a multifluid shock with mass transfer between the fluids by cosmic-ray ionization and ion-electron recombination, demonstrating the effect of ion mass loading on magnetic precursors of MHD shocks. An exact solution to an MHD Riemann problem forming the basis for an approximate numerical solver used in the homogeneous part of our scheme is presented, along with derivations of the analytic benchmark solutions and tests showing the convergence of the numerical algorithm.