python - 使用 FiPY 在 Python 中求解复杂的偏微分方程

Question

我正在尝试使用 FiPy 解决复杂的 PDE。FiPY 似乎不支持复杂的参数，一般建议似乎是拆分 PDE，分别求解，然后重新组合。现在还是这样吗？

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

谢谢你。

Answer 1

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

请参阅此耦合 Cahn-Hilliard 模拟的后半部分。

当然可以设置非均匀的初始条件。请参阅：examples/cahnHilliard/mesh2DCoupled.py，了解如何将初始条件设置为随机场或示例/diffusion/electrostatics.py，以及示例/phase/中的任何内容，了解如何根据坐标设置初始条件。正如您在这些示例中所看到的，有几种不同的语法可以使用，但 setValue() 是一种完全合法的方法。