Couette流动编程中遇到的一个问题
看了那本英文的Anderson 的计算流体力学后按照第九章的思路自己用VC++作了一个Couette流动的程序,迭代的时候发生错误。使用初始条件第一次计算与书上给出的结果完全一样,但是迭代到12次,和书上结果不一样,而且也没有收敛到与解析计算结果相同的解上面去。 请问按照△t的时间步长进行迭代计算是不是就是循环一定的次数啊? 比方说:12△t,就是循环12次。 360△t就是循环360 次? 哪位高人做过给我讲讲,不胜感激。[ 本帖最后由 xinyuxf 于 2007-7-22 16:50 编辑 ] 补充一下,我的程序是这样做的:先做的是△t=1时候,因为这个时候涉及到初始条件的代入。然后从第二次开始才做的按照步长的时间推进,就是一次一次的循环。现在发现第一次△t=1时候的结果是正确的,和书上的结果一样,可最终收敛的结果并不是解析解的结果。差不少。检查很多次了,好像没错误阿。是我方法不对么? 楼主用的是什么算法?
页:
[1]