Class AAHSweepChunkSpherical

Inheritance Relationships

Base Type

Class Documentation

class AAHSweepChunkSpherical : public opensn::SweepChunk

AAH sweep chunk for 1D spherical geometry.

Solves, for each direction mu_n of a GLProductQuadrature1DSpherical, the PWLD discretization (radial weight r^2) of mu d psi/dr + (1/r) [alpha_{n+1/2}/(w_n tau_n) + 2 mu_n] (psi_n - psi_{n-1/2})

  • sigma_t psi_n = q_n, where psi_{n-1/2} is the angular-edge flux from the previous direction in order of increasing mu and psi_{n+1/2} = (psi_n - (1 - tau_n) psi_{n-1/2}) / tau_n (weighted diamond difference in angle). The recursion starts from the zero-weight direction mu = -1, which has no angular term. The 1/r term uses mass matrices with weight r. Directions must reach each cell in order of increasing mu: the angle aggregation groups them into an inward and an outward set and the inward set precedes the outward set.

A boundary face at r = 0 (the center of a solid sphere) has zero area and receives no incident flux; any other boundary face, including the inner surface of a hollow sphere, uses the boundary condition.

Public Functions

AAHSweepChunkSpherical(DiscreteOrdinatesProblem &problem, LBSGroupset &groupset)
virtual void Sweep(AngleSet &angle_set) override

Sweep chunks should override this.