Solving the Time-dependent Discrete-ordinates Boltzmann equation on Parallel Architectures

Doug Swesty
State University of New York, Stony Brook

The solution of the time-dependent discrete-ordinates form of the
equation arises in wide variety of transport and radiation-hydrodynamic
problems. In this talk I will describe one approach to the numerical
solution of this equation on massively parallel architectures. Instead
of using a more traditional source-iteration method, we employ Krylov
subspace algorithms to solve the full linear system arising
from the implicitly-discretized Boltzmann equation. This approach
numerous advantages and a few drawbacks. The approach is also easily
extensible, via Newton-Krylov methods, to non-linear cases where
radiation-heating or radiation-hydrodynamic effects are included in the
problem. The success of this full linear system method lies in the
development of preconditioners which aid the convergence of the
Krylov-subspace algorithms. I will discuss some preconditioning
that are useful for broad classes of problems.

Presentation Files (Zip Archive)

Back to Computational Methods in Transport