Function opensn::ApplyDSACorrectionToDelayedBoundaryFlux

Function Documentation

void opensn::ApplyDSACorrectionToDelayedBoundaryFlux(DiscreteOrdinatesProblem &do_problem, const LBSGroupset &groupset, const std::vector<double> &delta_phi_local)

Applies a diffusion-synthetic scalar-flux correction to the lagged (delayed) incoming angular fluxes of opposing reflecting boundaries, as the isotropic angular correction delta psi = delta phi (quadrature weights sum to one).

Source iteration carries the lagged boundary angular fluxes as unknowns alongside phi. A DSA update that corrects phi but not these fluxes leaves the boundary error to re-enter on the next sweep, which can make source iteration with DSA diverge.

delta_phi_local has the layout of the local flux-moment vector.