如何用 FiPy 表示这些非标准的 pdes?

How to represent these non-standard pdes with FiPy?

我正在尝试使用 FiPy 求解一组 PDE,但不确定如何表示方程:see the PDEs here。

给我带来麻烦的术语是红色矩形和蓝色矩形中的那些。红色矩形中的那个是对流项乘以取决于 x、T1 和 T2 的函数;蓝色矩形中的那个是瞬态项乘以另一个函数,具体取决于 x、T1 和 T2。

我找不到包含此类术语的示例。你能给我一些建议,或者给我一个例子吗?

感谢任何帮助。提前致谢。

这个问题类似于cellvariable*Diffusion in fipy,但是是对流而不是扩散。解决方案是相同的。使用,

将红色项转换为系数为 f 和源项的对流项。

对于蓝色项,执行完全相同的操作以获得具有系数和源项的瞬态项。


编辑:如果我们假设 那么我们仍然可以使用

并明确近似(在 FiPy 中将 f 作为变量并使用 grad)。但是,我们可以进一步使用,

并明确地近似最终的

同样,我们可以更进一步,

最后一项可以明确解决。根据 f 的形式,最后一个 应该 变得更加微不足道,因此,明确性不再是一个问题。