如何在Fipy中指定Neumann (固定通量法向面)边界条件？

from fipy import *

nx = 50  # 网格的x方向节点数
ny = 50  # 网格的y方向节点数
dx = 1.0  # 网格的x方向步长
dy = 1.0  # 网格的y方向步长
mesh = Grid2D(nx=nx, ny=ny, dx=dx, dy=dy)

phi = CellVariable(name="phi", mesh=mesh)  # 定义变量phi
D = 1.0  # 扩散系数
eq = TransientTerm() == DiffusionTerm(coeff=D)  # 定义扩散方程

valueLeft = 0.0  # 左边界的值
fluxRight = 1.0  # 右边界的通量值

# 定义左边界条件
phi.constrain(valueLeft, mesh.facesLeft)

# 定义右边界条件
phi.faceGrad.constrain([fluxRight], mesh.facesRight)

phi.setValue(0.0)  # 设置初始条件为0

dt = 0.1  # 时间步长
steps = 100  # 总时间步数

for step in range(steps):
    eq.solve(var=phi, dt=dt)