Relaxation Schemes for the M1 Model with Space-Dependent Flux: Application to Radiotherapy Dose Calculation. (15th January 2016)

• Record Type: Journal Article
• Authors: Pichard, Teddy; Aregba-Driollet, Denise; Brull, Stéphane; Dubroca, Bruno; Frank, Martin
• Is Part Of: Communications in computational physics. Volume 19:Number 1(2016:Jan.)
• Page Start: 168
• Page End: 191
• Publication Date: 2016-01-15
• Subjects: 82C40, -- 35L65, -- 65M25; Radiotherapy, -- Moments models, -- relaxation models, -- method of characteristics; Mathematical physics -- Data processing -- Periodicals; Physics -- Data processing -- Periodicals; 530.150285
• DOI: 10.4208/cicp.121114.210415a
• Languages: English
• ISSNs: 1815-2406

Abstract

Abstract: Because of stability constraints, most numerical schemes applied to hyperbolic systems of equations turn out to be costly when the flux term is multiplied by some very large scalar. This problem emerges with the M 1 system of equations in the field of radiotherapy when considering heterogeneous media with very disparate densities. Additionally, the flux term of the M 1 system is non-linear, and in order for the model to be well-posed the numerical solution needs to fulfill conditions called realizability. In this paper, we propose a numerical method that overcomes the stability constraint and preserves the realizability property. For this purpose, we relax the M 1 system to obtain a linear flux term. Then we extend the stencil of the difference quotient to obtain stability. The scheme is applied to a radiotherapy dose calculation example.