如何验证FiPy求解得到的解满足目标偏微分方程？

一、热方程解的基础验证 #

假设你已通过FiPy求解得到各时间步的解数组Solution和对应的时间序列T，验证步骤如下：

1. 还原热方程与网格配置
   确保方程形式、网格参数与求解阶段完全一致，示例代码：

   from fipy import CellVariable, Grid1D, TransientTerm, DiffusionTerm
   
   # 复用求解时的网格参数
   mesh = Grid1D(nx=你的网格单元数, dx=你的网格步长)
   alpha = 热扩散系数  # 必须和求解时的系数一致
   eqn = TransientTerm() == DiffusionTerm(coeff=alpha)
   

2. 逐时间步计算残差
   遍历每个时间步的解，代入方程验证：

   phi = CellVariable(mesh=mesh)
   # 若时间步长非均匀，需按序列取对应步的间隔
   dt = T[1] - T[0]
   
   for idx, t in enumerate(T):
       # 加载当前时间步的解
       phi[:] = Solution[idx]
       # 计算残差（瞬态问题需传入时间步长dt）
       residual = eqn.residual(phi, dt=dt)
       # 计算残差的L2范数
       residual_norm = residual.norm()
       # 根据精度要求判断解是否满足PDE
       print(f"时间步 {idx}, 时间 {t:.2f}, 残差L2范数: {residual_norm:.6e}")
   

   若残差范数处于1e-6~1e-8量级（可根据问题精度调整），说明该时间步的解满足热方程。

二、复杂PDE的残差验证（用residualVectorAndNorm） #

针对含对流、扩散等项的复杂PDE，且已加载cell-centered解phi_cell，可直接利用eqn.residualVectorAndNorm避免重新求解的开销：

1. 还原复杂PDE配置
   定义与求解时完全一致的方程，示例含对流、扩散的瞬态方程：

   from fipy import ConvectionTerm
   
   # 复用求解时的对流速度、扩散系数等参数
   u = 对流速度变量  # 可以是CellVariable或数组
   D = 扩散系数  # 可以是常数或CellVariable
   eqn = TransientTerm() + ConvectionTerm(coeff=u) == DiffusionTerm(coeff=D)
   

2. 加载已有的cell-centered解
   将外部加载的phi_cell赋值给FiPy的CellVariable：

   phi = CellVariable(mesh=mesh)
   phi[:] = phi_cell  # 确保phi_cell长度与网格单元数匹配
   

3. 计算残差向量与范数
   直接调用方法计算，无需执行求解：

   # 瞬态问题需传入对应时间步的dt
   residual_vec, residual_norm = eqn.residualVectorAndNorm(phi, dt=对应时间步长)
   # 稳态问题无需传入dt
   # residual_vec, residual_norm = eqn.residualVectorAndNorm(phi)
   
   print(f"残差L2范数: {residual_norm:.6e}")
   

关键注意事项 #

• 所有系数（扩散系数、对流速度等）、网格参数必须与求解时完全一致，否则残差计算会失真。
• 瞬态问题中，dt需对应解所在时间步的时间间隔，确保瞬态项计算准确。
• 残差范数的合理阈值需结合问题精度要求判断，通常与求解时设置的收敛公差（如solve()的tolerance参数）处于同一量级即可认为解满足PDE。

内容的提问来源于stack exchange，提问作者user3166083