Multi-dimensional radiation transport for modeling axisymmetric Z pinches: Ray tracing compared to Monte Carlo solutions for a two-level atom

Radiation plays a critical role in the dynamical evolution and energetics of Z pinches and a dominant role for those loads designed as radiation sources in the keV photon energy region. Therefore, in modeling such plasmas it is essential to employ an accurate and realistic calculation of radiation generation and transport. This study presents a straightforward method to carry out three-dimensional radiation transport calculations in an azimuthally symmetric cylindrical geometry for use in magnetohydrodynamic simulations of high energy Z-pinch plasma radiation sources. The present approach is based on the integral solution of the radiative transfer equation using long characteristics with equal weight angular quadratures, and uses either a multifrequency grid or probability-of-escape-based zone-to-zone coupling. Using previously published Monte Carlo calculations of a two-level atom in cylindrical geometry as a benchmark, the accuracy of the approach is assessed as a function of the frequency grid, number of rays employed, and spatial grid. An analytic expression for the frequency averaged escape probability of an emission line is found to be as accurate for the population kinetics as the multifrequency method. While the escape probability description sacrifices information on the line profile, it is far more computationally efficient yet sufficiently accurate for line dominated plasmas such as Z-pinch radiation sources.

[1]  J. Weisheit.  Photon escape probabilities for stark-broadened lyman series lines , 1979 .

[2]  L. Rudakov,et al.  Model of enhanced energy deposition in a Z-pinch plasma , 2000 .

[3]  P. Gouttebroze.  Radiative transfer in cylindrical threads with incident radiation V. 2D transfer with 3D velocity fields , 2008 .

[4]  A. Skumanich,et al.  Radiative transport in finite homogenous cylinders by the Monte Carlo technique , 1969 .

[5]  Gordon L. Olson,et al.  Short characteristic solution of the non-LTE line transfer problem by operator perturbation—II. The two-dimensional planar slab , 1987 .

[6]  J. Apruzese,et al.  Direct solution of the equation of transfer using frequency- and angle-averaged photon escape probabilities. Interim report , 1979 .

[7]  C. Coverdale,et al.  Enhanced energy coupling and x-ray emission in Z-pinch plasma implosions , 2004 .

[8]  G. Moses,et al.  Non-LTE radiation transport in moderate-density plasmas , 1990 .

[9]  J. Apruzese,et al.  Investigation of K‐shell emission from moderate‐Z, low‐η (‐velocity), Z‐pinch implosions , 1996 .

[10]  J. Kilkenny,et al.  X‐ray spectroscopy of high‐energy density inertial confinement fusion plasmas , 1993 .

[11]  L. Auer,et al.  Short characteristic integration of radiative transfer problems: Formal solution in two-dimensional slabs , 1988 .

[12]  D. Hillier,et al.  2D non-LTE modeling for axisymmetric winds : Method and test cases , 2006 .

[13]  M. Carlsson.  3D radiative transfer in stellar atmospheres , 2008 .

[14]  E. Lewis,et al.  Computational Methods of Neutron Transport , 1993 .

[15]  J. Giuliani,et al.  Non-local radiation transport via coupling constants for the radially inhomogeneous Hg–Ar positive column , 2005 .

[16]  I. Hubeny,et al.  Multidimensional Non-LTE Radiative Transfer. I. A Universal Two-dimensional Short-Characteristics Scheme for Cartesian, Spherical, and Cylindrical Coordinate Systems , 2002 .

[17]  J. Kub'at,et al.  Radiative transfer in moving media : II. Solution of the radiative transfer equation in axial symmetry , 2005, astro-ph/0506378.

[18]  G. Mellema,et al.  Hybrid Characteristics: 3D radiative transfer for parallel adaptive mesh refinement hydrodynamics , 2005, astro-ph/0505213.

[19]  D. Ryutov,et al.  The physics of fast Z pinches , 1998 .

[20]  A. Velikovich,et al.  One- and two-dimensional modeling of argon K-shell emission from gas-puff Z-pinch plasmas , 2007 .

[21]  B. Wandelt,et al.  Adaptive ray tracing for radiative transfer around point sources , 2001, astro-ph/0111033.

[22]  A new exact method for line radiative transfer , 2005, astro-ph/0510616.

[23]  D. Osterbrock,et al.  Astrophysics of Gaseous Nebulae , 1976 .

[24]  J. Apruzese.  An analytic Voigt profile escape probability approximation , 1985 .

[25]  Laser target model , 1975 .

[26]  R. Mancini,et al.  Escape factors for Stark-broadened line profiles , 1987 .

[27]  David G. Hummer,et al.  An accelerated lambda iteration method for multilevel radiative transfer. II : Overlapping transitions with full continuum , 1992 .

[28]  T. Holstein.  Imprisonment of Resonance Radiation in Gases , 1947 .

[29]  C. Coverdale,et al.  Ion viscous heating in a magnetohydrodynamically unstable Z pinch at over 2 x 10(9) Kelvin. , 2006, Physical review letters.

[30]  J. Levine,et al.  Energetics of a long-implosion-time, 12-cm-diameter argon-gas-puff Z pinch at 6.5 MA , 2007, 2007 16th IEEE International Pulsed Power Conference.

[31]  G. Rybicki,et al.  An accelerated lambda iteration method for multilevel radiative transfer. I - Non-overlapping lines with background continuum , 1991 .

[32]  A. Velikovich,et al.  An efficient tabulated collisional radiative equilibrium radiation transport model suitable for multidimensional hydrodynamics calculations , 2001 .

[33]  A. Fisher,et al.  Ion-kinetic-energy measurements and energy balance in a Z-pinch plasma at stagnation. , 2007, Physical Review Letters.