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)

问候，莱昂纳多