simple编的非稳态传热程序跟时间步长有关会是什么错误呢？

apfel

用西交大simple程序编的非稳态传热程序，改变时间步长会有较大的影响，会是什么错误呢？
用simple时间算是比较长，但是总感觉理解不透彻，总会出现这样那样的问题，那位用的熟的，指点一下，谈谈体会心得？

anyone

你是算真正非稳态问题, 解和时间的关系不对; 还是用非稳态方法算稳态问题最终解不对?
如果是前者, 很正常. SIMPLE 时间大概用的隐式一阶精度, 时间精度并不是很好, 原则上要把时间步长放小, 直到解和时间步长无关.
如果是后者, 不太正常, 当然, 也许解还没有完全收敛.

apfel

下面引用由anyone在 2005/06/09 11:33pm 发表的内容：
你是算真正非稳态问题, 解和时间的关系不对; 还是用非稳态方法算稳态问题最终解不对?
如果是前者, 很正常. SIMPLE 时间大概用的隐式一阶精度, 时间精度并不是很好, 原则上要把时间步长放小, 直到解和时　...

多谢，anyone：）
是算非稳态问题，时间步长要放多小啊，0.5秒我觉得够小了，没有再往下试，那我再试试看。simple的时间精度可以改善吗？
顺便问一下，“非稳态方法算稳态问题”指什么啊？

anyone

"非稳态方法算稳态问题"是这样的, 如果边界(和其他)条件是定常的, 不管初始条件如何, 大部分情况会最终收敛到一个稳定的解. 求解时你可以一开始就把它看成稳态问题,用稳态方法来求. 但也可以把它看成非稳态问题, 从某个初始条件出发, 经过时间推移, 最后得到同样的稳定解. 这种算法相当于对稳态方法进行某种亚松弛, 但这种亚松弛有物理意义. 有一段时间很流行.
忘了说了, 收敛不好也会造成问题. 每个时间步你可以多算一些ITERATION看看会不会好一点.

fallye

你说的0.5秒太不确切了。对于地基的热传导问题，可能10秒都嫌少；对于微流体，也许总时间才一秒，你直接取步长就0.5秒，那是不可能的。取多少完全与你的问题有关，SIMPLE是隐式的，按理来说多长的时间步都没问题。但是SIMPLE采用的时间离散一般全是后差欧拉型的，这个的精度只有一阶精度，如果时间步长太大的话，那计算的结果肯定不准。对于你说的时间精度改善，我想是可以的，但是比较困难。因为传统的SIMPLE方法对于时间来说是两层格式，如果要提高时间的精度，那肯定就出现多层格式，这样的难度应该相当大。