Surface element in fenics

[Haibo_Ni]

Hi all,

I would like to solve a diffusion problem on a shell-like surface mesh (3D). Does the legacy fenics or fenicsx support 3D surface elements?

Any suggestions will be deeply appreciated. Thanks

[dokken]

Yes, it does, ref: How are 2D manifolds immersed/embedded in 3D space handled in dolfinx? - #25 by dokken
and Add manifold example with real spaces · Issue #118 · scientificcomputing/scifem · GitHub
and

[SOLVED] Solve PDEs on a Manifold

As far as Im aware these features are in DOLFINx. For instance, you can adapt Figure 3 to work with ffcx: from ufl import Cell, VectorElement, FunctionSpace, Mesh, Coefficient # Define triangle cell embedded in Rˆ3 cell = Cell("triangle", 3) # Define Lagrange vector element coordinate_element = VectorElement("Lagrange", cell, 1) mesh = Mesh(coordinate_element) V = FunctionSpace(mesh, coordinate_element) # Arguments defined over V will have 3 components: u = Coefficient(V) print(u[0], u[1], u…

and maybe particularly interesting for you

FEniCSx-Shells — FEniCSx-Shells 0.10.0.dev0 documentation