Adaptive solution of the domain decomposition+L2-jumps method applied to the neutron diffusion equation on structured meshes

Abstract

At the core scale, neutron deterministic calculations are usually based on the neutron diffusion equation. Classically, this equation can be recast in a mixed variational form, and then discretized by using the Raviart-Thomas-Nédélec Finite Element. The goal is to extend the Adaptive Mesh Refinement (AMR) strategy previously proposed in [1] to the Domain Decomposition+L2 jumps which allows non conformity at the interface between subdomains. We are able to refine each subdomain independently, which eventually leads to a more optimal refinement. We numerically investigate the improvements made to the AMR strategy.