Skip to main content

Implicit Asymptotic Preserving Method for Linear Transport Equations

  • Qin Li (a1) (a2) and Li Wang (a3)

The computation of the radiative transfer equation is expensive mainly due to two stiff terms: the transport term and the collision operator. The stiffness in the former comes from the fact that particles (such as photons) travel at the speed of light, while that in the latter is due to the strong scattering in the optically thick region. We study the fully implicit scheme for this equation to account for the stiffness. The main challenge in the implicit treatment is the coupling between the spacial and angular coordinates that requires the large size of the to-be-inverted matrix, which is also ill-conditioned and not necessarily symmetric. Our main idea is to utilize the spectral structure of the ill-conditioned matrix to construct a pre-conditioner, which, along with an exquisite split of the spatial and angular dependence, significantly improve the condition number and allows a matrix-free treatment. We also design a fast solver to compute this pre-conditioner explicitly in advance. Our method is shown to be efficient in both diffusive and free streaming limit, and the computational cost is comparable to the state-of-the-art method. Various examples including anisotropic scattering and two-dimensional problems are provided to validate the effectiveness of our method.

Corresponding author
*Corresponding author. Email addresses: (Q. Li), (L. Wang)
Hide All
[1] Ashby, S., Brown, P., Dorr, M. and Hindmarsh, A., A linear algebraic analysis of diffusion synthetic acceleration for the Boltzmann transport equation, SIAM J. Numer. Anal., 32 (1995), 128178.
[2] Adams, M. and Larsen, E., Fast iterative methods for discrete-ordinate particle transport calculations, Prog. Nucl. Eng., 40 (2002), 3159.
[3] Azmy, Y., Unconditionally stable and robust adjacent-cell diffusive preconditioning of weighted-difference particle transport methods is impossible, J. Comput. Phys., 182 (2002), 213.
[4] Brantley, P. and Larsen, E., The simplified P3 approximation, Nucl. Sci. Eng., 134 (2001), 121.
[5] Boscarino, S., Pareschi, L. and Russo, G., Implicit-explicit Runge-Kutta scheme for hyperbolic systems and kinetic equations in the diffusion limit, SIAM J. Sci. Comput., 35 (2013), 2251.
[6] Brown, P., A linear algebraic development of diffusion synthetic acceleration for three-dimensional transport equations, SIAM J. Numer. Anal., 32 (1995), 179214.
[7] Fleck, J. A. and Cummings, J. D., An implicit Monte Carlo scheme for calculating time and frequency dependent nonlinear radiation transpor, J. Comput. Phys., 8 (1971), 313342.
[8] Frank, M., Klar, A., Larsen, E. and Yasuda, S., Time-dependent simplified pn approximation to the equations of radiative transfer, J. Comput. Phys., 226 (2007), 22892305.
[9] Frank, M. and Seibold, B., Optimal prediction for radiative transfer: a new perspective on moment closure, Kinet. Relat. Models, 31 (2011), 717733.
[10] Guthrie, B., Holloway, J. and Patton, B., GMRES as a multi-step transport sweep accelerator, Trans. Theory Stat. Phys., 289(1) (1999), 83102.
[11] Golse, F., Jin, S. and Levermore, D., The convergence of numerical transfer schemes in diffusive regimes i: The discrete-ordinate method, SIAM J. Numer. Anal., 36 (1999), 13331369.
[12] Hauck, C., High-order entropy-based closures for linear transport in slab geometries, Commun. Math. Sci., 9 (2011), 187205.
[13] Hauck, C. and Lowrie, R., Temporal regularization of the pn equations, Multiscale Model. Simul., 7 (2009), 14971524.
[14] Jin, S. and Levermore, D., Fully discrete numerical transfer in diffusive regimes, Transp. Theory Stat. Phys., 22 (1993), 739791.
[15] Jin, S., Pareschi, L. and Toscani, G., Uniformly accurate diffusive relaxation schemes for multiscale transport equations, SIAM J. Numer. Anal., 38 (2000), 913936.
[16] Kupper, K., Frank, M. and Jin, S., An asymptotic preserving 2D staggered grid method for multiscale transport equations, submitted, 2015.
[17] Klar, A., An asymptotic-induced scheme for nonstationary transport equations in the diffusive limit, SIAM J. Numer. Anal., 35(6) (1998), 1097–1094.
[18] Larsen, E., Diffusion theory as an asymptotic limit of transport theory for nearly critical systems with small mean free paths, Ann. Nucl. Energy, 7 (1980), 249255.
[19] Lewis, E. Jr., and Miller, W., Computational Methods of Neutron Transport John Wiley and Sons, 1983.
[20] Lemou, M. and Mieussens, L., New asymptotic preserving scheme based on micro-macro formulation for linear kinetic equations in the diffusion limit, SIAM J. Sci. Comput., 31 (2008), 334368.
[21] Larsen, E. and Morel, J., Asymptotic solutions of numerical transport problems in optically thick, diffusive regimes ii, J. Comput. Phys., 83 (1989), 212236.
[22] Mieussens, L., On the asymptotic preserving property o fate unified gas kinetic scheme for the diffusion limit of linear kinetic modle, J. Comput. Phys., 253 (2013), 138156.
[23] Morel, J., Wareing, T., Lowrie, R. and Parsons, D., Analysis of ray-effect mitigation techniques, Nucl. Sci. Eng., 144(1) (2003), 122.
[24] Olbrant, E., Hauck, C. and Frank, M., A realizability-preserving discontinuous Galerkin method for the m1 model of radiative transfer, J. Comput. Phys., 231 (2012), 56125639.
[25] Olbrant, E., Larsen, E., Frank, M. and Seibold, B., Asymptotic derivation and numerical investigation of time-dependent simplified PN equations, J. Comput. Phys., 238(1) (2012), 315336.
[26] Pomraning, G., Asymptotic and variational derivations of the simplified PN equations, Ann. Nucl. Energy, 20 (1993), 623637.
[27] Sun, W., Jiang, S. and Xu, K., An asymptotic preserving unified gas kinetic scheme for gray radiative transfer equations, J. Comput. Phys., 285(5) (2015), 265279.
[28] Tomasevic, D. and Larsen, E., The simplified P2 approximation, Nucl. Sci. Eng., 122 (1996), 309325.
[29] Wareing, T. A., McGhee, J., Morel, J. and Pautz, S., Discontinuous finite element Sn methods on three-dimensional unstructured grids, Nucl. Sci. Eng., 138(3) (2001), 256268.
[30] Warsa, J., Wareing, T. and Morel, J., Krylov iterative methods and the degraded effectiveness of diffusion synthetic acceleration for multidimensional Sn calculations in problems with material discontinuities, Nuclear Math. Comput. Sci., 147 (2004), 218248.
Recommend this journal

Email your librarian or administrator to recommend adding this journal to your organisation's collection.

Communications in Computational Physics
  • ISSN: 1815-2406
  • EISSN: 1991-7120
  • URL: /core/journals/communications-in-computational-physics
Please enter your name
Please enter a valid email address
Who would you like to send this to? *


MSC classification


Full text views

Total number of HTML views: 0
Total number of PDF views: 0 *
Loading metrics...

Abstract views

Total abstract views: 0 *
Loading metrics...

* Views captured on Cambridge Core between <date>. This data will be updated every 24 hours.

Usage data cannot currently be displayed