Numerical solution of the radiative transfer equation with a posteriori error bounds

Add to your list(s) Download to your calendar using vCal

• Olga Mula (CEREMADE, Université Paris-Dauphine)
• Thursday 19 April 2018, 15:00-16:00
• MR14, Centre for Mathematical Sciences.

If you have a question about this talk, please contact Clarice Poon.

We propose a new approach to the numerical solution of radiative transfer equations with certified a posteriori error bounds. A key ingredient is the formulation of 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 based on approximately realizing this outer iteration within dynamically updated accuracy tolerances that still ensure convergence to the exact solution. On the one hand, since in the course of this iteration the global scattering operator is only applied, this avoids solving linear systems with densely populated system matrices while only linear transport equations need to be solved. This, in turn, rests on a Discontinous Petrov—Galerkin scheme which comes with rigorous a posteriori error bounds. These bounds are crucial for guaranteeing the convergence of the outer iteration. Moreover, the application of the global (scattering) operator is accelerated through low-rank approximation and matrix compression techniques. The theoretical findings are illustrated and complemented by numerical experiments with a non-trivial scattering kernel.

This talk is part of the Applied and Computational Analysis series.

This talk is included in these lists:

• All CMS events
• All Talks (aka the CURE list)
• Applied and Computational Analysis
• CMS Events
• DAMTP info aggregator
• Featured lists
• Interested Talks
• MR14, Centre for Mathematical Sciences
• My seminars
• Type the title of a new list here
• bld31

Note that ex-directory lists are not shown.