作者wowow2005 (耶稣是主)
看板Cad_Cae
标题[问题]想请问关於叠代次数的问题
时间Sat Oct 11 10:00:46 2008
我目前使用icepak模拟
若是稳态模拟,随着网格数的增加,收敛的残值会越来越小
但收敛速度很慢,要几万或十万次以上才可能达到收敛
而网格数较少的话,残值到某ㄧ定值就不会再收敛了
至於暂态模拟,每一个时间点的叠代次数几乎两三次或四五次就收敛了
这是我觉得蛮奇怪的地方
我是进行温度场的模拟的,收敛标准就照预设的1e-07
而稳态模拟顶多到(1.xxxx)e-06
我看theory部分有介绍AMG法,是因为这个关系才导致收敛这麽快吗?
还是我其它设定的问题需要改进?
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 61.228.154.12
1F：推 largesperm:最後的 flux 有合理吗 10/11 17:49
2F：→ wowow2005:阁下是说解出来的温度值还是HEAT FLUX? 10/11 19:55
3F：推 largesperm:heat flux~ 最後的净通量要比流入的还要小很多才对 10/11 20:57
4F：→ wowow2005:我在计算域内有热源,单位是W/m^3,计算域边壁有四面设 10/11 21:14
5F：→ wowow2005:opening,另外两面为石绵绝热层,所以热会由opening那四面 10/11 21:15
6F：→ wowow2005:出去,阁下所说的净通量是指计算完後出现的小空白视窗里 10/11 21:16
7F：→ wowow2005:面的数据吗? 10/11 21:16
8F：推 largesperm:哈~ 我没用过 icepak 不清楚 10/11 21:17
9F：→ largesperm:觉得结果有问题 都会先去看质量 能量有没有守恒 10/11 21:18
10F：→ largesperm:暂态的残值也有到 1e-7 吗 10/11 21:19
11F：→ wowow2005:有 10/11 21:25
12F：→ wowow2005:稳态则是几乎没到1E-07 10/11 21:29
13F：→ wowow2005:净通量理论上应该要是0吗? 10/11 21:30
14F：推 largesperm:不是 不会完全是 0, 跟流入流出量比起来小很多就行了 10/11 21:35
15F：→ wowow2005:可是我的热能是由系统内固体区域发出的ㄝ... 10/11 21:44
16F：→ wowow2005:流体区域则是AIR 10/11 21:44
17F：推 largesperm:系统内固体发出 热应该会存在空气中 或是流出系统吧 10/11 21:45
18F：→ largesperm:发出的热+在空气的热-流出系统的热 应该快接近 0 吧 @@ 10/11 21:46
19F：→ wowow2005:如果是由系统内发出,怎麽还会有流入的情形? 10/11 21:46
20F：→ largesperm:我觉得那个就算是流入的量吧~ 10/11 21:47
21F：→ wowow2005:是啊,既然这样,根据能量守恒,理论上不应该是0吗? 10/11 21:47
22F：→ largesperm:数值会有截断误差 所以净通量不可能为 0 10/11 21:48
23F：→ largesperm:当然理论上是要为 0 没错 10/11 21:48
24F：→ wowow2005:哦 原来是这样子 了解 10/11 21:49
25F：→ wowow2005:不过 什麽是截断误差啊 10/11 21:49
26F：推 largesperm:truncation error~ 方程式是非线性的 所以才用数值 10/11 21:50
27F：→ largesperm:把方程式离散一小块一小块求解 但是这也是近似而已 10/11 21:51
28F：→ wowow2005:哦 这个我有听过 了解 谢谢! 10/11 21:52
29F：→ largesperm:round error 也是误差的原因之一 电脑的数字有限 10/11 21:52
30F：推 meaning0728:网格愈多..不一定会愈准..要注意..有个上限 10/12 10:58
31F：→ meaning0728:AMG应该指Algebraic Multi-Grid.是一种加速收敛的方法 10/12 11:00
32F：→ wowow2005:对,因为在CFDACE+中,可以选择要不要使用AMG,而在ICEPAK 10/12 15:16
33F：→ wowow2005:里,就没有这个选项,所以我在想是不是因为AMG才导致收敛 10/12 15:16
34F：→ wowow2005:这麽快的(或是叠代次数这麽少的?) 10/12 15:17
35F：推 meaning0728:时间和空间都要做网格独立测试 10/14 01:35
36F：→ wowow2005:楼上不懂你的意思? 10/17 20:50
37F：→ johan1026:因为AMG会视收敛情形自动调整Courant number 10/25 21:02
38F：推 largesperm:楼上帅哥专业 10/26 00:32