Purgatorio—a new implementation of the Inferno algorithm

Abstract

An overview of PURGATORIO, a new implementation of the INFERNO 1 equation of state model, is presented. The new algorithm emphasizes a novel decimation scheme for automatically resolving the structure of the continuum density of states, circumventing limitations of the pseudo-R matrix algorithm previously utilized.

Equation of state, occupation probabilities and conductivities in the average atom Purgatorio code

“…Capturing the resonances is essential for a self-consistent calculation and is ensured either by the refinement described in Ref. [1] or by tracking changes in the phase shifts [6].…”

“…Previous calculations [1] used an adaptive energy grid method based on Gaussian quadrature to resolve detailed resonance structure in the density of states.…”

“…Prototypical results from Purgatorio, including the principal Hugoniots of Be and Al and pressure vs. energy isochores for low-density Al, have been published recently [1].…”

“…The Purgatorio code [1] is a recent implementation of Liberman's Inferno model [2], which belongs to a family of historically successful ion-in-cell or neutral-pseudo-atom (NPA) approaches to the electron-thermal contribution to the EOS [see, for example, Refs. [3][4][5][6]].…”

Capsule modeling of high foot implosion experiments on the National Ignition Facility

“…The bulk of the simulations described above are 2D simulations using ALE hydrodynamics [68], multi-group radiation transport in the diffusion approximation [69], tabular equation of state [70][71][72] and opacity data [73], and Monte Carlo transport of fusion products. These characteristics raise the necessity of demonstrating convergence of the simulations with respect to the computational mesh used to resolve the hydrodynamic motion, the number of radiation groups used to transport x-ray radiation, and the number of Monte Carlo particles used to model the slowing down of fusion products.…”