ABSTRACT In this article, least‐squares weak Galerkin finite element methods are presented to solve the and ‐elliptic interface problems with non‐homogeneous jump conditions in both primal and flux variables. The proposed methods enable the incorporation of discontinuous functions within mesh element partitions composed of general polytopal elements. Optimal‐order error estimates have been derived for both flux and primal unknowns. A series of numerical tests for the proposed problems has been conducted to corroborate our theoretical findings, encompassing interfaces of varying smoothness and complexity.
Raman Kumar (2025) studied this question.