Galerkin methods for primary ion transport in inhomogeneous media

Journal article, 2010

This paper concerns the energy deposition of high-energy (e.g., approximate to 50 - 500 MeV) proton and carbon ions and high-energy electrons (of approximate to 50 MeV), in inhomogeneous media. Our goal is to develop a flexible model incorporated with the analytic theory for ions based on bipartition and Fokker-Planck developments. Both procedures are leading to convection dominated convection diffusion equations. We study convergence for semi-discrete and fully discrete approximations of a such obtained equation, for abroad beam model, using the standard Galerkin and streamline diffusion finite element methods. The analytic broad beam model of the light ion absorbed dose were compared with the results of the modified Monte Carlo (MC) code SHIELD-HIT+ and those of Galerkin streamline diffusion approach.