A surface integral based momentum advection scheme for neutron transport in moving materials Conference Paper uri icon

abstract

  • A new momentum advection scheme for neutron transport in moving materials is developed and tested for 1-D spatial geometries. In previous work, the momentum advection term was discretized by taking integrals over velocity cell volumes. The new scheme replaces these integrals over cell volumes with surface integrals by applying the divergence theorem. The surface integral approach is shown to be accurate and numerically stable and is additionally easier to implement and extend to multi-dimensional geometries than the old scheme. Preconditioning strategies that exploit the structure of the momentum advection matrix equation and can be applied to either advection scheme are also developed and shown to accelerate the convergence of the GMRES solver. The new advection scheme and preconditioning strategies are tested on a problem of interest, namely the r-process in collapsing proto-neutron stars.

published proceedings

  • Mathematics and Computations, Supercomputing in Nuclear Applications and Monte Carlo International Conference, M and C+SNA+MC 2015

author list (cited authors)

  • Fichtl, E. D., Baker, R. S., Dahl, J. A., & Morel, J. E.

complete list of authors

  • Fichtl, ED||Baker, RS||Dahl, JA||Morel, JE

publication date

  • January 2015