python fipy的安装及简单使用

这是笔者在安装fipy和使用中做的总结和自己的理解，希望能帮助到读者。

前言

安装fipy自带的库不能解决三维问题的显示和任意边界区域的设定。

• 前者需要额外安装mayavi
• 后者需要额外安装gmsh.

fipy

FiPy是基于标准有限体积（FV）方法的用Python编写的面向对象的偏微分方程（PDE）求解器。

安装

pip install future
pip install fipy

示例代码

examples.diffusion.mesh20x20

from fipy import CellVariable, Grid2D, Viewer, TransientTerm, DiffusionTerm
from fipy.tools import numerix
from fipy import input

nx = 20
ny = nx
dx = 1.
dy = dx
L = dx * nx
mesh = Grid2D(dx=dx, dy=dy, nx=nx, ny=ny)

phi = CellVariable(name = "solution variable",
                   mesh = mesh,
                   value = 0.)

D = 1.
eq = TransientTerm() == DiffusionTerm(coeff=D)                  

valueTopLeft = 0
valueBottomRight = 1

X, Y = mesh.faceCenters
facesTopLeft = ((mesh.facesLeft & (Y > L / 2))
                | (mesh.facesTop & (X < L / 2)))
facesBottomRight = ((mesh.facesRight & (Y < L / 2))
                    | (mesh.facesBottom & (X > L / 2)))



phi.constrain(valueTopLeft, facesTopLeft)
phi.constrain(valueBottomRight, facesBottomRight)

if __name__ == '__main__':
    viewer = Viewer(vars=phi, datamin=0., datamax=1.)
    viewer.plot()


timeStepDuration =