问题标签 [fipy]

在示例文档中，有以下部分：

FiPy 不会绘制或输出任何内容，除非您告诉它：

我了解当前配置会导致使用 matplotlib 或 Mayavi 查看器打开查看器。但是，我希望能够导出 .pvd 或 .xdmf 文件以合并模拟。

谢谢你的帮助！

我正在使用 FiPy 对给定速度场中的化学物质的对流/扩散进行建模，但我在将该场设置为 ConvectionTerm 系数时遇到了问题。

我的速度场由二维数组表示，假设Ugrid水平速度和Vgrid垂直速度，表示每个单元面的值。出于这个原因，我想说将此字段设置coeff为 ConvectionTerm 参数的正确方法是将其分配给 FaceVariable。

但是，我不知道如何传递这两个数组作为valueFaceVariable。显然，我通过使用将字段设置为 CellVariable 没有问题，value = [np.ravel(Ugrid),np.ravel(Vgrid)]并且对流模拟似乎也很有意义，但我认为这是不正确的，正如我上面简要提到的那样。

有什么建议么？

我想在 Spyder IDE 中使用 PDE Solver Fipy 并尝试按照官方安装指南安装它： 现在建议在使用包之前运行 fipy.test() ，这里出现以下错误。

我不明白这在这种情况下意味着什么，这里是这个错误所指的代码：

我想用下面的代码定义的“L 形”几何形状来解决热扩散问题：

现在我想将“顶部”和“右”边界定义为狄利克雷边界，常数值为 1 和 0。我想将“内”和“右”定义为通量等于 0 的诺伊曼边界。我应该怎么走关于设置约束？我用了：

但我得到的只是“AttributeError：'Gmsh2D'对象没有属性'Front'”错误。我怎样才能解决这个问题？非常感谢你。

我刚刚找到了用于描述对流项的数值方案，这些方案用于 FiPy ( https://www.ctcms.nist.gov/fipy/documentation/numerical/scheme.html )。

我目前正在实施类似的方案来解决对流问题。我想研究这些方案与此处介绍的方案相比的性能。

是否有参考资料，我可以在其中阅读有关所提出方案的更多信息？

先感谢您！

我想通过求解一组偏微分方程来模拟化学反应器的吸附床。我想通过空间项的有限体积技术将 eq_1 和 eq_2 简化为一组常微分方程。Ct 、 DL 和 eb 是常量，我很难在 FiPy 代码中找到合适的模型和技术。

基本上，我想用 100 个切片形成一维模拟，其中气体种类 i、j 等从 x=0 开始并移动到列的其余部分。我应该使用什么类型的模型？我已经检查了来自 FiPy 代码网站的示例，但找不到我想要的确切配方。

谢谢，

eq_1

eq_2

我生成了一个 gmsh 网格，我通过 Gmsh2D 将其导入 FiPy。我想处理表面面以设置边界条件，但我不知道如何。

在我在文档中找到的 FiPy 示例中，建议命名某些行以便稍后解决它们。将网格导入 fipy 后如何执行此操作？

-----------------

编辑：

对于简单的表面，这将起作用：我忽略了mesh.exteriorFaces. 对于一个简单的圆圈，这导致了简单的解决方案：

已识别表面细胞

尽管如此，我仍在寻找有关如何从外部访问 gmsh 代码的答案，但这可能对其他人有所帮助:)

我正在尝试解决由 Fipy 中的 gmsh 创建的网格中的 Meinhart 模型。但是，我不确定如何添加零通量边界条件。下面，你可以找到我的代码。我想知道这些行是否：

应激活以确保零通量边界条件。

我读到：https ://github.com/usnistgov/fipy/issues/674 “FiPy 的边界默认情况下是无通量的”。但是，我不知道它是否也适用于由 gmsh 选项创建的网格。

如果不需要激活或在我的代码中添加额外的行来设置零通量边界条件，它也适用于我有更复杂的网格，例如多边形、圆形或其他不规则形状？

谢谢，

我目前正在研究 fipy 文档中的示例，并尝试对其进行调整examples.cahnHilliard.mesh2DCoupled以执行一些模拟。链接是：https ://www.ctcms.nist.gov/fipy/examples/cahnHilliard/generated/examples.cahnHilliard.mesh2DCoupled.html

似乎示例代码声明了 aneq3和 adfdphi_并且没有在其他任何地方使用这些方程，即使在解决方案中也是如此。

从文档中：

在解决方案块中：

我将不胜感激有关eq3和dfdphi_条款的作用的任何见解。作为对此的扩展，求解器如何知道使用 $\phi$ 的旧值来评估泰勒表达式后的线性化形式，因为求解器块似乎没有考虑到这一点。

有没有这样的功能或简单的方法？到目前为止，我发现的唯一功能是mesh.vertexCoords，mesh.faceVertexIDs但我不知道它们是否对我有帮助。