Dear NGSolvers,
A student and I are working on a (surface) shape optimization problem. As part of the problem we need to solve a linear elasticity problem on a surface, and as we might approach the incompressible case in our application, we want to employ a mixed formulation based on a the HDivSurface space plus weakly enforcement of the tangential continuity, very similar to what Wang/Ye did in their work (“New finite element methods in computational fluid dynamics by H (div) elements”, SINUM, 2007).
Therefore we are wondering how to correctly compute the surface divergence, and in particular the surface gradient and its correct surface transpose (or equivalently, the symmetric derivative ) for the test and trial functions from the HDivSurface space? Are these diff operators already implemented? Unfortunately, we haven’t managed to find any documentation regarding this space…
Best,
André