Code for the rSN method

The code can be found here. Below is the documentation that you also find on the linked website.

We implement the rSN method. This is a modification of the standard SN method that mitigates ray-effects by rotation of the quadrature step after each time step. An interpolation procedure is required which furthermore requires different (non-tensorized) quadrature sets. All these aspects are implemented in this code with Julia.
The two figures below demonstrate the effect that the rotation of the angular quadrature set has on the overall solution quality. In the first picture we see ray-effects for the line-source problem. In the second picture these effects were mitigated due to the rotation.

Continue reading “Code for the rSN method”