Bulk-surface coupling with vector elements

[dokken]

Just to clarify, you would like to implement a H-curl space on \Gamma which is the conforming boundary of a 3D domain \Omega? I don’t see any obvious issues with this.
The only thing that comes to mind is to be careful about global orientations, ref this length post:

Does Dolfinx's solver have an issue or require some preconditioning when handling certain meshes like the Icosahedral mesh?

I am going to try to illustrate the mistake I mean exists in the legacy code with a small example: from dolfin import * class OutwardNormal(UserExpression): def __init__(self, mesh, *arg, **kwargs): super().__init__(*arg, **kwargs) self.mesh = mesh def eval_cell(self, values, x, cell): c = Cell(self.mesh, cell.index) normal = c.cell_normal() scale = 1 if c.orientation == 0 else -1 values[0] = normal[0] * scale values[1] = normal…

with code at: Add manifold example with real spaces · Issue #118 · scientificcomputing/scifem · GitHub