A Newton-Krylov solver for implicit solution of hydrodynamics in core collapse supernovae

This paper describes an implicit approach and nonlinear solver for solution of radiation-hydrodynamic problems in the context of supernovae and proto-neutron star cooling. The robust approach applies Newton-Krylov methods and overcomes the difficulties of discontinuous limiters in the discretized equations and scaling of the equations over wide ranges of physical behavior. We discuss these difficulties, our approach for overcoming them, and numerical results demonstrating accuracy and efficiency of the method.

[1]  G. Sod.  A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws , 1978 .

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

[3]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[4]  Hirata,et al.  Observation of a neutrino burst from the supernova SN1987A. , 1988, Physical review letters.

[5]  Park,et al.  Observation of a neutrino burst in coincidence with supernova 1987A in the Large Magellanic Cloud. , 1987, Physical review letters.

[6]  C. Kelley.  Iterative Methods for Linear and Nonlinear Equations , 1987 .

[7]  P. Brown.  A local convergence theory for combined inexact-Newton/finite-difference projection methods , 1987 .

[8]  S. Woosley.  The Birth of Neutron Stars , 1987 .

[9]  Yousef Saad,et al.  Hybrid Krylov Methods for Nonlinear Systems of Equations , 1990, SIAM J. Sci. Comput..

[10]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[11]  M. Norman,et al.  ZEUS-2D: A radiation magnetohydrodynamics code for astrophysical flows in two space dimensions. I - The hydrodynamic algorithms and tests. II - The magnetohydrodynamic algorithms and tests , 1992 .

[12]  F. Swesty,et al.  THE ROLE OF THE EQUATION OF STATE IN THE PROMPT PHASE OF TYPE II SUPERNOVAE , 1994 .

[13]  F. Swesty.  Implicit general relativistic hydrodynamic methods for modeling the late-time explosion mechanism in core collapse supernovae , 1995 .

[14]  H. Janka,et al.  Ledoux Convection in Protoneutron Stars—A Clue to Supernova Nucleosynthesis? , 1996, astro-ph/9610203.

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

[16]  A. S. Umar,et al.  The Interplay between Proto-Neutron Star Convection and Neutrino Transport in Core-Collapse Supernovae , 1997, astro-ph/9709184.

[17]  J. Lattimer,et al.  Evolution of Proto-Neutron Stars , 1998, astro-ph/9807040.

[18]  D. Keyes,et al.  Jacobian-free Newton-Krylov methods: a survey of approaches and applications , 2004 .

[19]  Carol S. Woodward,et al.  Enabling New Flexibility in the SUNDIALS Suite of Nonlinear and Differential/Algebraic Equation Solvers , 2020, ACM Trans. Math. Softw..

[20]  F. Douglas Swesty,et al.  A NUMERICAL ALGORITHM FOR MODELING MULTIGROUP NEUTRINO-RADIATION HYDRODYNAMICS IN TWO SPATIAL DIMENSIONS , 2006 .