FEniCS:计算点处的函数梯度

Question

设u是函数空间V上用FEniCS求解的给定问题的解

from fenics import *
...
u = Function(V)
solve(a==L, u, bcs)

x是一个Point对象。如何计算u在x上的梯度？

我试过了

g = grad(u)
g(x)

这将引发一个有关维度的错误：

UFLException: Expecting dim to match the geometric dimension, got dim=1 and gdim=2.

在函数空间上投影梯度也不能带来令人鼓舞的结果：

h = project(grad(u), V)

还会引发错误。

UFLException: Shape mismatch.

Answer 1

我猜你在标量函数空间中投影函数时，可能存在维数问题。你可能需要这样做：

 V_vec = VectorFunctionSpace(mesh, "CG", 1)
 h = project(grad(u),V_vec)

你好，莱昂纳多