Python中pde有限差分求解器的交替方向隐式方法

Question

我正在努力实施Alternating direction implicit method来解决FitzHugh–Nagumo反应扩散模型。我在blog中找到了一个Python实现示例，但我认为方法中存在错误 - 在此处提供的模板中： 不应该是半个时间步长乘以反应项 f ？

Answer 1

用差分商替换差异商，得到

U_t = D/2 * U_xx + D/2 * U_yy + Δt*f

在两个实例中，这不是等式

U_t = D * (U_xx + U_yy) + f

这是最初提出的任务。

所以系数1/(Δt/2)应该是U_t，D/(Δp^2)位于U_pp，p=x,y和1位于{{1 }}。

似乎公式是具有不同商数的公式和下一个级别乘以f的混合。

在下一个公式中，我不需要新的常量Δt/2，α_p=σ_p，然后你是正确的p=x,y因子应该是f。< / p>