Hi everyone!

I would like to implement a bulk-surface equation. The model problem I am considering to begin with is the one presented in 1.1 in the following paper: Finite element analysis for a coupled bulk–surface partial differential equation | IMA Journal of Numerical Analysis | Oxford Academic. It reads:

I tried to navigate the forum and find how to deal with having two global fem spaces, one defined on the bulk and one on the boundary of the same bulk:

- Boundary Finite Element Spaces - Kunena
- How to define a codimention-1 finite element space on the interface - Kunena

I wonder if it is, currently, how I should treat this kind of problems or if you had suggestions on how to proceed on setting up the FEM spaces?

Thanks a lot for the help and all the functionality your software provides