我需要在 FEniCS 中求解以下方程:
据我了解,弱形式将如下所示:
我在 FEniCS 中表达这个方程时遇到问题。我想我需要这样的东西:
u = TrialFunction(V)
v = TestFunction(V)
g = as_vector((u.dx(0), u.dx(1))) # <-- modify here?
a = dot(g, grad(v))*dx
u = Function(V)
solve(a == Constant(0) * v * dx, u, bc)
如何添加取决于坐标的东西,到上面的变量g?
我还必须写Constant(0) * v * dx正确的部分吗?可以只有0吗?