×

zbMATH — the first resource for mathematics

On the performance of SPAI and ADI-like preconditioners for core collapse supernova simulations in one spatial dimension. (English) Zbl 1196.85008
Summary: The simulation of core collapse supernovæ calls for the time accurate solution of the (Euler) equations for inviscid hydrodynamics coupled with the equations for neutrino transport. The time evolution is carried out by evolving the Euler equations explicitly and the neutrino transport equations implicitly. Neutrino transport is modeled by the multi-group Boltzmann transport (MGBT) and the multi-group flux limited diffusion (MGFLD) equations. An implicit time stepping scheme for the MGBT and MGFLD equations yields Jacobian systems that necessitate scaling and preconditioning. Two types of preconditioners, namely, a sparse approximate inverse (SPAI) preconditioner and a preconditioner based on the alternating direction implicit iteration (ADI-like) have been found to be effective for the MGFLD and MGBT formulations. This paper compares these two preconditioners. The ADI-like preconditioner performs well with both MGBT and MGFLD systems. For the MGBT system tested, the SPAI preconditioner did not give competitive results. However, since the MGBT system in our experiments had a high condition number before scaling and since we used a sequential platform, care must be taken in evaluating these results.

MSC:
85-08 Computational methods for problems pertaining to astronomy and astrophysics
85A25 Radiative transfer in astronomy and astrophysics
76M28 Particle methods and lattice-gas methods
PDF BibTeX XML Cite
Full Text: DOI
References:
[1] D’Azevedo, E.F.; Messer, B.; Mezzacappa, A.; Liebendoerfer, M., An adi-like preconditioner for Boltzmann transport, SIAM J. sci. comput., 26, 810-820, (2005) · Zbl 1137.82304
[2] Swesty, F.D.; Smolarski, D.C.; Saylor, P.E., A comparison of algorithms for the efficient solution of the linear systems arising from multi-group flux-limited diffusion problems, Astrophys. J. suppl. ser., 153, 369-387, (2004)
[3] Mezzacappa, A.; Bruenn, S.W., A numerical method for solving the neutrino Boltzmann equation coupled to spherically symmetric stellar core collapse, Astrophys. J., 405, 669-684, (1993)
[4] M. Liebendörfer, Consistent modelling of core-collapse supernovæ in spherically symmetric relativistic space – time, Ph.D. Thesis, The University of Basel, 2000
[5] O.E.B. Messer, Questing for the grail: Spherically symmetric supernova simulations with Boltzmann neutrino transport, Ph.D. Thesis, The University of Tennessee, 2000
[6] Mihalas, D.; Mihalas, B.W., Foundations of radiation hydrodynamics, (1984), Dover Mineola · Zbl 0651.76005
[7] Turner, N.J.; Stone, J.M., A module for radiation hydrodynamic calculations with ZEUS-2D using flux-limited diffusion, Astrophys. J. suppl. ser., 135, 95-107, (2001)
[8] Varga, R.S., Matrix iterative analysis, (1962), Prentice-Hall Englewood Cliffs · Zbl 0133.08602
[9] Alvarado, F.L.; Schreiber, R., Optimal parallel solution of sparse triangular systems, SIAM J. sci. stat. comput., 14, 446-460, (1993) · Zbl 0774.65011
[10] Santos, E.E., On designing optimal parallel triangular solvers, Inform. comput., 161, 172-210, (2000) · Zbl 1045.68981
[11] Raghavan, P., Efficient parallel triangular solution with selective inversion, Parallel process. lett., 8, 29-40, (1998)
[12] Higham, N.J.; Pothen, A., The stability of partitioned inverse approach to triangular solution, SIAM J. sci. stat. comput., 15, 139-148, (1994) · Zbl 0795.65012
[13] H.I.M. Gould, J.A. Scott, On approximate-inverse preconditioners, Technical Report RAL 95-026, Computing and Information Systems Department, Rutherford Appleton Laboratory, Oxfordshire, England, 1995
[14] Benzi, M.; Bertaccini, D., Approximate inverse preconditioning for shifted linear systems, Bit, 43, 231-244, (2003) · Zbl 1037.65043
[15] Bollhöfer, M.; Mehrmann, V., Algebraic multilevel methods and sparse approximate inverses, SIAM J. matrix anal. appl., 24, 191-218, (2002) · Zbl 1035.65030
[16] Bollhöfer, M.; Saad, Y., A factored approximate inverse preconditioner with pivoting, SIAM J. matrix anal. appl., 23, 692-705, (2002) · Zbl 0999.65027
[17] Bollhöfer, M.; Saad, Y., On the relations between ILUs and factored approximate inverses, SIAM J. matrix anal. appl., 24, 219-237, (2002) · Zbl 1017.65019
[18] Chen, K., An analysis of sparse approximate inverse preconditioners for boundary integral equations, SIAM J. matrix anal. appl., 22, 1058-1078, (2001) · Zbl 0985.65038
[19] E.T.-F. Chow, Robust preconditioning for sparse linear systems, Ph.D. Thesis, The University of Minnesota, 1997
[20] Mezzacappa, A.; Messer, O.E.B., Neutrino transport in core collapse supernovae, J. comput. appl. math., 109, 281-319, (1999) · Zbl 0946.76063
[21] Bethe, H.; Wilson, J.R., Revival of a stalled supernova shock by neutrino heating, Astrophys. J., 295, 14-23, (1985)
[22] Grote, M.J.; Huckle, T., Parallel preconditioning with sparse approximate inverses, SIAM J. sci. comput., 18, 838-853, (1997) · Zbl 0872.65031
[23] Golub, G.H.; van Loan, C.F., Matrix computations, (1996), Johns Hopkins University Press Baltimore · Zbl 0865.65009
[24] Saad, Y., Iterative methods for sparse linear systems, (2003), SIAM Philadelphia · Zbl 1002.65042
[25] Young, D.M., Iterative solution of large linear systems, (1971), Academic Press New York · Zbl 0204.48102
[26] Greenbaum, A., Iterative methods for solving linear systems, (1997), SIAM Philadelphia · Zbl 0883.65022
[27] Press, W.H.; Teukolsky, S.A.; Vetterling, W.T.; Flannery, B.P., Numerical recipes in Fortran 77, (1992), Cambridge University Press Cambridge · Zbl 0778.65002
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.