A Mixed Extended Virtual Element Method for Elliptic Interface Problems on Polygonal Meshes

arXiv:2606.08526v1 Announce Type: new Abstract: We propose a lowest-order \(H(\operatorname{div})\)-conforming mixed extended virtual element method for elliptic interface problems on interface-unfitted polygonal meshes. The flux and pressure are approximated by subdomain-wise extended \(H(\operatorname{div})\)-VEM spaces and by piecewise constants, respectively. On cut elements, the computable polynomial projection is defined on the whole background element and then restricted to the two subdomains. Compared with BDM-type polynomial spaces, the mixed VEM space contains a non-polynomial component, which gives rise to additional consistency terms on cut elements. To control these terms, we use an enhanced kernel stabilization on cut elements and an interface normal-flux average in the mixed coupling. A corrected interface-flux penalty and a local divergence ghost penalty are added to obtain cut-position-independent stability without using a volume div-div augmentation. We prove continuity, a discrete inf-sup condition, and an optimal first-order error estimate in a mesh-dependent norm. The constants are independent of the mesh size and of the position of the interface relative to the background mesh, but may depend on the coefficient contrast.