Modeling fluid instabilities in inertial confinement fusion hydrodynamics codes

The numerical tools typically used to model the evolution of fluid instabilities in inertial confinement fusion hydrodynamics codes are examined, and some are found to have properties which would seem to be incompatible with the accurate modeling of small-amplitude perturbations, i.e., perturbations in the linear stage of evolution. In particular a “differentiability condition” which is satisfied by the physics in such situations is not necessarily satisfied by the numerical algorithms in typical use. It is demonstrated that it is possible to remove much of the nondifferentiability in many cases, and that substantial improvement in one’s ability to accurately model the evolution of small-amplitude perturbations can result. First a simple example involving a nondifferentiable radiation transport algorithm is shown, and then the nondifferentiabilities introduced by the use of upwind and “high resolution” hydrodynamics algorithms are analyzed.

[1]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[2]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[3]  Chi-Wang Shu,et al.  The Runge-Kutta Discontinuous Galerkin Method for Conservation Laws V , 1998 .

[4]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[5]  M.M.R. Williams.  Computing as a language of physics: IAEA, Vienna (1972). 616 pp., £9.00 , 1979 .

[6]  A. Velikovich,et al.  Saturation of perturbation growth in ablatively driven planar laser targets , 1998 .

[7]  R. D. Richtmyer,et al.  Difference methods for initial-value problems , 1959 .

[8]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[9]  P. Lax,et al.  Systems of conservation laws , 1960 .

[10]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

[11]  L. Spitzer.  Physics of fully ionized gases , 1956 .

[12]  A. Velikovich,et al.  Direct observation of mass oscillations due to ablative Richtmyer–Meshkov instability and feedout in planar plastic targets , 2001 .

[13]  I. Bohachevsky,et al.  Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics , 1959 .

[14]  A. Velikovich,et al.  Growth of Pellet Imperfections and Laser Imprint in Direct Drive Inertial Confinement Fusion Targets , 2001 .

[15]  V. Goncharov.  Theory of the Ablative Richtmyer-Meshkov Instability , 1999 .

[16]  John M. Dawson,et al.  Correct values for high-frequency power absorption by inverse bremsstrahlung in plasmas , 1973 .

[17]  Gregory A. Moses,et al.  Inertial confinement fusion , 1982 .

[18]  Phillip Colella,et al.  Efficient Solution Algorithms for the Riemann Problem for Real Gases , 1985 .

[19]  D. Mihalas,et al.  Foundations of Radiation Hydrodynamics , 1985 .

[20]  J. Dahlburg,et al.  Computational modeling of direct-drive fusion pellets and KrF-driven foil experiments , 1998 .

[21]  S. Zalesak.  Introduction to “Flux-Corrected Transport. I. SHASTA, A Fluid Transport Algorithm That Works” , 1997 .