用Python解决Python中复杂的PDE

Question

我试图用FiPy解决一个复杂的PDE。似乎FiPY不支持复杂的论点，一般的建议似乎是将PDE分开，单独解决然后重新组合。这仍然是这样吗？

其次，我看不出任何方式来设置非恒定的初始条件。假设我希望我的网格在每个单元格中都没有相同的值，那么除了使用setValue之外，还有更简单的方法吗？

谢谢。

Answer 1

FiPy不支持复杂的参数，但它确实支持向量的求解，这可能会简化您的任务。

查看this coupled Cahn-Hilliard simulation的后半部分。

当然可以设置非均匀的初始条件。请参阅：examples/cahnHilliard/mesh2DCoupled.py了解如何将初始条件设置为随机字段或examples/diffusion/electrostatics.py和 以及examples/phase/中关于如何根据坐标设置初始条件的任何内容。正如你在这些例子中看到的那样，有几种不同的语法可以工作，但setValue（）是完全合法的方法。