《安得森计算流体力学基础及其应用》第九章使用simple算法求解不可压库艾特流

paopaoai11

sdfh根据前人写的fortran源代码和安得森计算流体力学基础及其应用第六章和第九章相关内容用matlab编写了使用simple算法求解不可压库艾特流的程序。程序虽然可以跑通，但是有一个很奇怪的现象，就是只有当v(15,5)<=0.4时，程序可以得到相对正确的解，其他情况都无法收敛。希望感兴趣的同志可以帮忙看一下问题所在。

[ 本帖最后由 paopaoai11 于 2014-4-23 21:25 编辑 ]

hust_illusion

你好，v(15,5)初始化为0.5也是可以计算的。程序默认就是0.5