Skip to main content
Article thumbnail
Location of Repository

Solving the transport equation by the use of 6D spectral methods in spherical geometry

By Silvano Bonazzola, Nicolas Vasset and Bruno Peres


We present a numerical method for handling the resolution of a general transport equation for radiative particles, aimed at physical problems with a general spherical geometry. Having in mind the computational time difficulties encountered in problems such as neutrino transport in astrophysical supernovae, we present a scheme based on full spectral methods in 6d spherical coordinates. This approach, known to be suited when the characteristic length of the dynamics is much smaller than the domain size, has the potential advantage of a global speedup with respect to usual finite difference schemes. An analysis of the properties of the Liouville operator expressed in our coordinates is necessary in order to handle correctly the numerical behaviour of the solution. This reflects on a specific (spherical) geometry of the computational domain. The numerical tests, performed under several different regimes for the equation, prove the robustness of the scheme: their performances also point out to the suitability of such an approach to large scale computations involving transport physics for mass less radiative particles.Comment: 18 page

Topics: Physics - Computational Physics, Astrophysics - Solar and Stellar Astrophysics
Year: 2013
OAI identifier:
Download PDF:
Sorry, we are unable to provide the full text but you may find it at the following location(s):
  • (external link)
  • Suggested articles

    To submit an update or takedown request for this paper, please submit an Update/Correction/Removal Request.