问使用Fipy，PDE转换的非线性HJB PDE
EN

Stack Overflow用户

提问于 2020-05-03 13:02:42

回答 1查看 83关注 0票数 2

我正在尝试使用FiPy解决一个非线性PDE HJB方程，但是我在将PDE转换为正确的FiPy语法时遇到了一些困难：

我尝试了这样的方法: eqX = TransientTerm() == -DiffusionTerm(coeff=1) + (phi.faceGrad * phi.faceGrad)，但由于梯度的平方而不起作用

我的等式：(du/dt =- \delta + ||\grad(u)||^2)

FiPy允许求解这种方程吗？如果没有，有没有一个软件包或者使用有限差分来解决它的方法呢？谢谢！

发布于 2020-05-04 05:47:45

可以将最终项重写为扩散项和源项，从而可以将方程重写为，

eqn = TransientTerm() = DiffusionTerm(u - 1) - u * u.faceGrad.divergence

这不会产生错误，但可能不是很稳定

票数 2

页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持

原文链接：

https://stackoverflow.com/questions/61569845

复制

相似问题

问使用Fipy，PDE转换的非线性HJB PDEEN