meiyongyuandeze 发表于 2011-3-15 09:36

请教编程求积二维不可压缩N-S方程时的问题！

本文自己采用微分求积方法（DQM）在Fortan90中编程求解二维不可压缩N-S方程的程序，现在程序已经全部调试通过，没有BUG存在，残差收敛设为10-6。但问题是：给定一个初值，程序计算便会收敛到一个解，初值不同，最后结果页不同。检查了很长时间，很困惑。请教曾经自己编写过程序的高手给点建议，该类问题应该从哪方面入手进行检查，可能是什么原因导致的，是算法还是程序。

meiyongyuandeze 发表于 2011-3-18 10:59

难道没有建议吗！哎，编程序真是太恼火了，疯了{:{19}:}

ChaChing 发表于 2011-3-26 15:27

是所谓CFD程序吧!
外行人问下, 採用的编程算法是否严格收敛!? 是否有唯一解!?

meiyongyuandeze 发表于 2011-3-26 16:03

这个算法其实也是一个比较新的算法，国外有些学者在利用这个算法计算偏微分问题，但收敛性的问题一直没有得到很好的证明。我现在做的就是根据一位新加坡学者C.Shu的理论来应用求解我的问题，但是即使用那位学者的理论，相同的边界条件，什么设置都也一样，同一个模型，算出来的结果也相差很大。我不想对其它的做过多的评论，但只是想弄清楚原因。我估计有可能是该种离散数值方法的收敛性可能存在某些问题。现在我们实验室有个去年新进的数学博士正在做相关的研究。
一个程序，花了我1年得时间，一定要弄个明白！{:{23}:}

ChaChing 发表于 2011-3-26 23:56

回复 4 # meiyongyuandeze 的帖子

为何会那样提醒, 是因为以前曾请教过CFD方面的博士朋友这个问题, 若记的没错, 许多编程算法都尚未有严格收敛唯一解证明的

meiyongyuandeze 发表于 2011-3-27 00:17

恩，谢谢关注，如有新的进展会第一时间在论坛和各位分享
页: [1]
查看完整版本: 请教编程求积二维不可压缩N-S方程时的问题！