【发布时间】:2022-01-25 10:54:58
【问题描述】:
我正在尝试使用 FiPy 求解一组 PDE,但不确定如何表示方程:see the PDEs here。
给我带来麻烦的术语是红色矩形和蓝色矩形中的术语。红色矩形中的一个是对流项乘以取决于 x、T1 和 T2 的函数;蓝色矩形中的一个是瞬态项乘以另一个取决于 x、T1 和 T2 的函数。
我找不到任何包含此类术语的示例。你能给我一些建议,或者分享一个例子吗?
感谢您的帮助。提前致谢。
【问题讨论】:
我正在尝试使用 FiPy 求解一组 PDE,但不确定如何表示方程:see the PDEs here。
给我带来麻烦的术语是红色矩形和蓝色矩形中的术语。红色矩形中的一个是对流项乘以取决于 x、T1 和 T2 的函数;蓝色矩形中的一个是瞬态项乘以另一个取决于 x、T1 和 T2 的函数。
我找不到任何包含此类术语的示例。你能给我一些建议,或者分享一个例子吗?
感谢您的帮助。提前致谢。
【问题讨论】:
这个问题类似于cellvariable*Diffusion in fipy,但针对的是对流而不是扩散。解决方案是相同的。使用,
将红色项转换为系数为 f 和源项的对流项。
对于蓝色项执行完全相同的操作以获得具有系数和源项的瞬态项。
编辑:如果我们假设 那么我们仍然可以使用
并显式近似(在 FiPy 中将 f 作为变量并使用 grad)。但是,我们可以更进一步并使用它,
并明确近似最终的。
再一次,我们可以走得更远,
再次,最后一项可以明确解决。根据f 的形式,最后他们应该变得更加微不足道,因此,明确性不再是问题。
【讨论】: