A HDG method for elliptic problems with integral boundary condition: theory and applications

In this paper, we address the study of elliptic boundary value problems in presence of a boundary condition of integral type (IBC) where the potential is an unknown constant and the flux (the integral of the flux density) over a portion of the boundary is given by a value or a coupling condition. We first motivate our work with realistic examples from nano-electronics, high field magnets and ophthalmology. We then define a general framework stemming from the Hybridizable Discontinuous Galerkin method that accounts naturally for the IBC and we provide a complete analysis at continuous and discrete levels. The implementation in the Feel++framework is then detailed and the convergence and scalability properties are verified. Finally, numerical experiments performed on the real-life motivating applications are used to illustrate our methodology.