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)sy...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.展开更多
In plasma physics domain,the electron transport is described with the FokkerPlanck-Landau equation.The direct numerical solution of the kinetic equation is usually intractable due to the large number of independent va...In plasma physics domain,the electron transport is described with the FokkerPlanck-Landau equation.The direct numerical solution of the kinetic equation is usually intractable due to the large number of independent variables.That is why we propose in this paper a new model whose derivation is based on an angular closure in the phase space and retains only the energy of particles as kinetic dimension.To find a solution compatible with physics conditions,the closure of the moment system is obtained under a minimum entropy principle.This model is proved to satisfy the fundamental properties like a H theorem.Moreover an entropic discretization in the velocity variable is proposed on the semi-discrete model.Finally,we validate on numerical test cases the fundamental properties of the full discrete model.展开更多
文摘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.
基金supported by EURATOM within the"Keep-in-Touch"activities and was granted access to the HPC resources of CINES under the allocation 2011-056129 made by GENCI(Grand Equipement National de Calcul Intensif).
文摘In plasma physics domain,the electron transport is described with the FokkerPlanck-Landau equation.The direct numerical solution of the kinetic equation is usually intractable due to the large number of independent variables.That is why we propose in this paper a new model whose derivation is based on an angular closure in the phase space and retains only the energy of particles as kinetic dimension.To find a solution compatible with physics conditions,the closure of the moment system is obtained under a minimum entropy principle.This model is proved to satisfy the fundamental properties like a H theorem.Moreover an entropic discretization in the velocity variable is proposed on the semi-discrete model.Finally,we validate on numerical test cases the fundamental properties of the full discrete model.