这可能是一个简单的问题,但是如果我想用梯度之外的空间变化系数求解 PDE,那么正确的 FiPy 语法是什么?到目前为止,我看到的所有示例都只讨论了梯度内部的系数。
例如:
d/dt(Sigma) = (1/r) d/dr (r^0.5 d/dr(nu Sigma r^0.5))
(我忽略了数字因素)
我想求解 Sigma(t,r)。我如何在 d/dr 面前处理 (1/r)?
我知道可以对这个简单的方程进行按摩,这样我就不必担心梯度之外的空间变化系数(或者只是将系数移动到时间导数项内),但我必须为我试图解决的实际问题,这个技巧将不再有效。例如,如果我的方程式如下所示,我该怎么办:
d/dt (var) = f(r) d^2/dr^2 (var) + g(r) d/dr (var)
任何帮助将不胜感激!