Bulk-surface virtual element method for systems of PDEs in two-space dimensions

Frittelli, Massimo, Madzvamuse, Anotida and Sgura, Ivonne (2021) Bulk-surface virtual element method for systems of PDEs in two-space dimensions. Numerische Mathematik, 147. pp. 305-348. ISSN 0029-599X

[img] PDF - Accepted Version
Restricted to SRO admin only

Download (881kB)
[img] PDF - Published Version
Available under License Creative Commons Attribution.

Download (2MB)


In this paper we consider a coupled bulk-surface PDE in two space dimensions. The model consists of a PDE in the bulk that is coupled to another PDE on the surface through general nonlinear boundary conditions. For such a system we propose a novel method, based on coupling a virtual element method (Beirão Da Veiga et al. in Math Models Methods Appl Sci 23(01):199–214, 2013. https://doi.org/10.1051/m2an/2013138) in the bulk domain to a surface finite element method (Dziuk and Elliott in Acta Numer 22:289–396, 2013. https://doi.org/10.1017/s0962492913000056) on the surface. The proposed method, which we coin the bulk-surface virtual element method includes, as a special case, the bulk-surface finite element method (BSFEM) on triangular meshes (Madzvamuse and Chung in Finite Elem Anal Des 108:9–21, 2016. https://doi.org/10.1016/j.finel.2015.09.002). The method exhibits second-order convergence in space, provided the exact solution is H2+1/4 in the bulk and H2 on the surface, where the additional 14 is required only in the simultaneous presence of surface curvature and non-triangular elements. Two novel techniques introduced in our analysis are (i) an L2-preserving inverse trace operator for the analysis of boundary conditions and (ii) the Sobolev extension as a replacement of the lifting operator (Elliott and Ranner in IMA J Num Anal 33(2):377–402, 2013. https://doi.org/10.1093/imanum/drs022) for sufficiently smooth exact solutions. The generality of the polygonal mesh can be exploited to optimize the computational time of matrix assembly. The method takes an optimised matrix-vector form that also simplifies the known special case of BSFEM on triangular meshes (Madzvamuse and Chung 2016). Three numerical examples illustrate our findings.

Item Type: Article
Schools and Departments: School of Mathematical and Physical Sciences > Mathematics
SWORD Depositor: Mx Elements Account
Depositing User: Mx Elements Account
Date Deposited: 24 Nov 2020 08:52
Last Modified: 25 Feb 2022 14:15
URI: http://sro.sussex.ac.uk/id/eprint/95271

View download statistics for this item

📧 Request an update
Project NameSussex Project NumberFunderFunder Ref
InCeM: Research Training Network on Integrated Component Cycling in Epithelial Cell MotilityG1546EUROPEAN UNIONUnset
New predictive mathematical and computational models in experimental sciencesG1949ROYAL SOCIETYWM160017
Unravelling new mathematics for 3D cell migrationG1438LEVERHULME TRUSTRPG-2014-149