这里指出 FiPy 不擅长处理双曲偏微分方程。然而这个论坛是 2016 年写的,所以也许事情已经发生了一些变化。
我想实现 2D Burgers 方程:$$
u_t + (u \partial_x + v \partial_y) u = - 你* abla^{2} u \
v_t + (u \partial_x + v \partial_y) v = - 你* abla^{2} v
$$
其中 $u$ 和 $v$ 是流体流动的 x 和 y 分量 $U = (u, v) $。
为了让我的生活更轻松一些,我使用链式法则将方程重写为
$$
u_t + 阿布拉