An adaptive nested source term iteration for radiative transfer equations
HTML articles powered by AMS MathViewer
- by Wolfgang Dahmen, Felix Gruber and Olga Mula HTML | PDF
- Math. Comp. 89 (2020), 1605-1646 Request permission
Abstract:
We propose a new approach to the numerical solution of radiative transfer equations with certified a posteriori error bounds for the $L_2$ norm. A key role is played by stable Petrov–Galerkin-type variational formulations of parametric transport equations and corresponding radiative transfer equations. This allows us to formulate an iteration in a suitable, infinite-dimensional function space that is guaranteed to converge with a fixed error reduction per step. The numerical scheme is then based on approximately realizing this iteration within dynamically updated accuracy tolerances that still ensure convergence to the exact solution. To advance this iteration two operations need to be performed within suitably tightened accuracy tolerances. First, the global scattering operator needs to be approximately applied to the current iterate within a tolerance comparable to the current accuracy level. Second, parameter dependent linear transport equations need to be solved, again at the required accuracy of the iteration. To ensure that the stage dependent error tolerances are met, one has to employ rigorous a posteriori error bounds which, in our case, rest on a Discontinuous Petrov–Galerkin (DPG) scheme. These a posteriori bounds are not only crucial for guaranteeing the convergence of the perturbed iteration but are also used to generate adapted parameter dependent spatial meshes. This turns out to significantly reduce overall computational complexity. Since the global operator is only applied, we avoid the need to solve linear systems with densely populated matrices. Moreover, the approximate application of the global scatterer is accelerated through low-rank approximation and matrix compression techniques. The theoretical findings are illustrated and complemented by numerical experiments with non-trivial scattering kernels.References
- Bradley K. Alpert, A class of bases in $L^2$ for the sparse representation of integral operators, SIAM J. Math. Anal. 24 (1993), no. 1, 246–262. MR 1199538, DOI 10.1137/0524016
- Mohammad Asadzadeh, $L_2$-error estimates for the discrete ordinates method for three-dimensional neutron transport, Transport Theory Statist. Phys. 17 (1988), no. 1, 1–24. MR 950623, DOI 10.1080/00411458808230852
- Mohammad Asadzadeh, A finite element method for the neutron transport equation in an infinite cylindrical domain, SIAM J. Numer. Anal. 35 (1998), no. 4, 1299–1314. MR 1620140, DOI 10.1137/S0036142992238119
- Matias Avila, Ramon Codina, and Javier Principe, Spatial approximation of the radiation transport equation using a subgrid-scale finite element method, Comput. Methods Appl. Mech. Engrg. 200 (2011), no. 5-8, 425–438. MR 2749012, DOI 10.1016/j.cma.2010.11.003
- Guillaume Bal, Inverse transport theory and applications, Inverse Problems 25 (2009), no. 5, 053001, 48. MR 2501018, DOI 10.1088/0266-5611/25/5/053001
- M. Blatt, A. Burchardt, A. Dedner, C. Engwer, J. Fahlke, B. Flemisch, C. Gersbacher, C. Gräser, F. Gruber, C. Grüninger, D. Kempf, R. Klöfkorn, T. Malkmus, S. Müthing, M. Nolte, M. Piatkowski, and O. Sander, The Distributed and Unified Numerics Environment, Version 2.4. Archive of Numerical Software, 4(100), pp. 13–29, May 2016. http://dx.doi.org/10.11588/ans.2016.100.26526.
- D. Broersen, W. Dahmen, and R. P. Stevenson, On the stability of DPG formulations of transport equations, Math. Comp. 87 (2018), no. 311, 1051–1082. MR 3766381, DOI 10.1090/mcom/3242
- T. A. Brunner, Forms of Approximate Radiation Transport, Sandia Report SAND2002-1778, Sandia National Laboratories, July 2002. http://dx.doi.org/10.2172/800993.
- Albert Cohen, Wolfgang Dahmen, and Ronald DeVore, Adaptive wavelet methods for elliptic operator equations: convergence rates, Math. Comp. 70 (2001), no. 233, 27–75. MR 1803124, DOI 10.1090/S0025-5718-00-01252-7
- A. Cohen, W. Dahmen, and R. DeVore, Adaptive wavelet methods. II. Beyond the elliptic case, Found. Comput. Math. 2 (2002), no. 3, 203–245. MR 1907380, DOI 10.1007/s102080010027
- Wolfgang Dahmen, Helmut Harbrecht, and Reinhold Schneider, Compression techniques for boundary integral equations—asymptotically optimal complexity estimates, SIAM J. Numer. Anal. 43 (2006), no. 6, 2251–2271. MR 2206435, DOI 10.1137/S0036142903428852
- Wolfgang Dahmen, Chunyan Huang, Christoph Schwab, and Gerrit Welper, Adaptive Petrov-Galerkin methods for first order transport equations, SIAM J. Numer. Anal. 50 (2012), no. 5, 2420–2445. MR 3022225, DOI 10.1137/110823158
- Wolfgang Dahmen and Rob P. Stevenson, Adaptive strategies for transport equations, Comput. Methods Appl. Math. 19 (2019), no. 3, 431–464. MR 3977482, DOI 10.1515/cmam-2018-0230
- Robert Dautray and Jacques-Louis Lions, Mathematical analysis and numerical methods for science and technology. Vol. 6, Springer-Verlag, Berlin, 1993. Evolution problems. II; With the collaboration of Claude Bardos, Michel Cessenat, Alain Kavenoky, Patrick Lascaux, Bertrand Mercier, Olivier Pironneau, Bruno Scheurer and Rémi Sentis; Translated from the French by Alan Craig. MR 1295030, DOI 10.1007/978-3-642-58004-8
- Herbert Egger and Matthias Schlottbom, A mixed variational framework for the radiative transfer equation, Math. Models Methods Appl. Sci. 22 (2012), no. 3, 1150014, 30. MR 2890452, DOI 10.1142/S021820251150014X
- Herbert Egger and Matthias Schlottbom, An $L^p$ theory for stationary radiative transfer, Appl. Anal. 93 (2014), no. 6, 1283–1296. MR 3195888, DOI 10.1080/00036811.2013.826798
- Konstantin Grella and Christoph Schwab, Sparse discrete ordinates method in radiative transfer, Comput. Methods Appl. Math. 11 (2011), no. 3, 305–326. MR 2844781, DOI 10.2478/cmam-2011-0017
- F. Gruber, Adaptive Source Term Iteration: A Stable Formulation for Radiative Transfer, Ph.D. thesis, RWTH Aachen University, 2018. http://dx.doi.org/10.18154/RWTH-2018-230893.
- F. Gruber, A. Klewinghaus, and O. Mula, The DUNE-DPG Library for Solving PDEs with Discontinuous Petrov–Galerkin Finite Elements, Archive of Numerical Software, 5(1), pp. 111–128, 6 March 2017. http://dx.doi.org/10.11588/ans.2017.1.27719.
- Jean-Luc Guermond and Guido Kanschat, Asymptotic analysis of upwind discontinuous Galerkin approximation of the radiative transport equation in the diffusive limit, SIAM J. Numer. Anal. 48 (2010), no. 1, 53–78. MR 2608358, DOI 10.1137/090746938
- L. G. Henyey and J. L. Greenstein, Diffuse Radiation in the Galaxy, The Astrophysical Journal, 93, pp. 70–83, 1941.
- Claes Johnson and Juhani Pitkäranta, Convergence of a fully discrete scheme for two-dimensional neutron transport, SIAM J. Numer. Anal. 20 (1983), no. 5, 951–966. MR 714690, DOI 10.1137/0720065
- Guido Kanschat, Solution of radiative transfer problems with finite elements, Numerical methods in multidimensional radiative transfer, Springer, Berlin, 2009, pp. 49–98. MR 2484899, DOI 10.1007/978-3-540-85369-5_{5}
- M. F. Modest and J. Yang, Elliptic PDE Formulation and Boundary Conditions of the Spherical Harmonics Method of Arbitrary Order for General Three-dimensional Geometries, Journal of Quantitative Spectroscopy and Radiative Transfer, 109(9), pp. 1641–1666, 2008. http://dx.doi.org/10.1016/j.jqsrt.2007.12.018.
- M. Mokhtar-Kharroubi, Mathematical topics in neutron transport theory, Series on Advances in Mathematics for Applied Sciences, vol. 46, World Scientific Publishing Co., Inc., River Edge, NJ, 1997. New aspects; With a chapter by M. Choulli and P. Stefanov. MR 1612403, DOI 10.1142/9789812819833
- J. C. Ragusa, J.-L. Guermond, and G. Kanschat, A robust $S_N$-DG-approximation for radiation transport in optically thick and diffusive regimes, J. Comput. Phys. 231 (2012), no. 4, 1947–1962. MR 2876596, DOI 10.1016/j.jcp.2011.11.017
Additional Information
- Wolfgang Dahmen
- Affiliation: Mathematics Department, University of South Carolina, 1523 Greene Street, Columbia, South Carolina 29208
- MR Author ID: 54100
- Email: dahmen@math.sc.edu
- Felix Gruber
- Affiliation: IGPM, RWTH Aachen, Templergraben 55, 52056 Aachen, Germany
- Email: felgru@gmx.de
- Olga Mula
- Affiliation: CEREMADE, CNRS, UMR 7534, Université Paris-Dauphine, PSL University, 75016 Paris, France
- MR Author ID: 1027794
- Email: mula@ceremade.dauphine.fr
- Received by editor(s): November 20, 2018
- Received by editor(s) in revised form: October 3, 2019, and October 22, 2019
- Published electronically: January 8, 2020
- Additional Notes: The research of the first author was supported by the NSF Grant DMS 1720297, and by the SmartState and Williams-Hedberg Foundation.
The third author is the corresponding author - © Copyright 2020 American Mathematical Society
- Journal: Math. Comp. 89 (2020), 1605-1646
- MSC (2010): Primary 65N12, 65N15, 65N30; Secondary 65N50
- DOI: https://doi.org/10.1090/mcom/3505
- MathSciNet review: 4081913