gpt4 book ai didi

python - 使用 Fipy 的非线性 HJB PDE,PDE 转换

转载 作者:行者123 更新时间:2023-12-04 10:06:22 32 4
gpt4 key购买 nike

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

我试过类似的东西:
eqX = TransientTerm() == -DiffusionTerm(coeff=1) + (phi.faceGrad * phi.faceGrad)
它不起作用,因为梯度的平方

我的方程:(du/dt = -\delta u + ||\grad(u)||^2)

FiPy 是否允许求解此类方程?如果没有,是否有使用有限差分的软件包或方法来解决它?
谢谢!

最佳答案

可以将最后一项重铸为扩散项和源项,以便方程可以重写为,

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

这不会出错,但可能不是很稳定

关于python - 使用 Fipy 的非线性 HJB PDE,PDE 转换,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/61569845/

32 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com