【问题标题】:How to represent these non-standard pdes with FiPy?如何用 FiPy 表示这些非标准 pdes?
【发布时间】:2022-01-25 10:54:58
【问题描述】:

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

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

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

感谢您的帮助。提前致谢。

【问题讨论】:

    标签: pde fipy


    【解决方案1】:

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

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

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


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

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

    并明确近似最终的。

    再一次,我们可以走得更远,

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

    【讨论】:

    • 但是f函数依赖于\phi,所以我不能用这个方法...
    • 好的,我已经尝试在编辑中解决这个问题。如果明确评估 f_x,仍然可以使用初始方程。但是,这似乎不如使用其他身份之一那么令人满意。
    猜你喜欢
    • 1970-01-01
    • 2015-01-03
    • 2021-10-03
    • 2021-08-29
    • 2020-03-03
    • 1970-01-01
    • 2015-07-08
    • 2019-11-07
    • 2013-06-01
    相关资源
    最近更新 更多