高手指教编程计算流场停止条件？

zhaocc · 发表于 2006-5-22 15:51:11

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？注册

x

1、请问大家在编程计算流场时，都是通过什么方法控制停止迭代的？
2、我想在前后两次迭代过程中速度场变化小于0.1%时停止迭代，请问是所有离散点的的速度变化都小于0.1%，还是对于某一点参考或者某一截面速度场变化小于0.1%就行？

请高手指教！！！谢谢先！

zhaocc · 发表于 2006-5-22 16:02:25

谁有计算流体力学的MD（分子动力学）程序能发给我，我将不胜感激！

anyone · 发表于 2006-5-23 23:45:15

一般用残差来判断.
比如方程是
AX+BY=C
你的解是X=X'; Y=Y';
残差 R=C-(AX';+BY)

下面引用由zhaocc在 2006/05/22 03:51pm 发表的内容：
1、请问大家在编程计算流场时，都是通过什么方法控制停止迭代的？
2、我想在前后两次迭代过程中速度场变化小于0.1%时停止迭代，请问是所有离散点的的速度变化都小于0.1%，还是对于某一点参考或者某一截面速度场变化小于0.1%就行？

请高手指教！！！谢谢先！

zhaocc · 发表于 2006-5-24 11:09:28

谢谢楼上先。
我现在就使用残差来控制迭代停止的。可是大部分时候残差很难满足小于某个极小值的条件。都是用迭代次数（10000次左右）来控制迭代停止的，这样的话对于不同的参数计算选择同样的迭代停止次数，感觉不是很合理。
听说过可以用速度场在先后两次迭代过程中变化不小于0.1%的方法来控制迭代次数。望大家指教怎么实施。

CFDNHT · 发表于 2006-5-24 17:21:19

控制容积的残差与参考流量的比值,小于一个很小的数,如
SMAX/Flowin < 1.0E-5
或者每隔n次迭代后,压差或摩擦因子两者相差的误差小于一个很小的数,如
abs( dp(t)-dp(t+20)/dp(t+20) ) < 0.001

zhaocc · 发表于 2006-5-25 17:20:04

谢谢斑竹，请问这里的压差是某一个截面的压差，还是某一个参考点的压差？

anyone · 发表于 2006-5-25 23:36:16

使用残差是正招. 一般可以用相对残差. 就是开始计算的残差和结束计算的残差之比.(一般说残差掉了几个数量级)这样就可以大部分解决不同物理模型残差标准应该相应变化的问题.
有时问题本身用现在的CFD技术很难收敛, 残差掉不下来(比如只掉1-2个数量级甚至更少)也是没办法的事情. 但很多情况下,特别是初学者, 是程序写的不好,或计算条件设置不合理. 如果用改变收敛判断来掩盖是没有意义的. 而且也错失了发现自己错误的机会.

zhaocc · 发表于 2006-5-29 09:55:25

谢谢，anyone！

firelet · 发表于 2006-6-11 10:19:51

下面引用由anyone在 2006/05/25 11:36pm 发表的内容：
使用残差是正招. 一般可以用相对残差. 就是开始计算的残差和结束计算的残差之比.(一般说残差掉了几个数量级)这样就可以大部分解决不同物理模型残差标准应该相应变化的问题.
有时问题本身用现在的CFD技术很难收敛 ...

经验之谈，学习了：）

		自动登录	找回密码
密码			注册

高手指教编程计算流场停止条件？

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。