| 查看: 6049 | 回复: 28 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
fangdian木虫 (正式写手)
|
[求助]
一个老生常谈的话题----亚松弛因子
|
||
|
由于流体力学中要求解非线性的方程,在求解过程中,控制变量的变化是很必要的,这就通过松弛因子来实现的.它控制变量在每次迭代中的变化.也就是说,变量的新值为原值加上变化量乘以松弛因子。 一般来讲,大家都是在收敛不好的时候，采用一个较小的欠松弛因子。 Fluent里面用的是欠松弛，主要防止两次迭代值相差太大引起发散。 松弛因子的值在0～1之间，越小表示两次迭代值之间变化越小，也就越稳定，但收敛也就越慢。 由于FLUENT所解方程组的非线性，我们有必要控制的变化。一般用亚松驰方法来实现控制，该方法在每一部迭代中减少了的变化量。亚松驰最简单的形式为：单元内变量等于原来的值加上亚松驰因子a与变化的积分离解算器使用亚松驰来控制每一步迭代中的计算变量的更新。这就意味着使用分离解算器解的方程，包括耦合解算器所解的非耦合方程（湍流和其他标量）都会有一个相关的亚松驰因子。 注：在FLUENT中，所有变量的默认亚松驰因子都是对大多数问题的最优值。这个值适合于很多问题，但是对于一些特殊的非线性问题（如：某些湍流或者高Rayleigh数自然对流问题），在计算开始时要慎重减小亚松驰因子。 使用默认的亚松驰因子开始计算是很好的习惯。如果经过4到5步的迭代残差仍然增长，你就需要减小亚松驰因子。 有时候，如果发现残差开始增加，你可以改变亚松驰因子重新计算。 但是此时，我有一个问题，就是关于亚松弛因子-----liquid fraction update 我的计算模型主要是3D的VOF和熔化凝固模型。对于亚松弛因子的设定如下：  本来都是默认设定的。liquid fraction update的默认值为0.9. 可是在计算中间，连续性残差不稳，出现突然上升的情况。  此时的Global Courant Number还是小于1的。一个解决办法是减小时间步。可是时间步减小，计算到指定时间所需要的时间步就多了。于是我将liquid fraction update的亚松弛因子由默认的0.9改为0.5接着往下计算，“躲过”了这次的不稳定。  可是，这样直接修改松弛因子接着往下算，有没有问题？对结果有没有影响？ 有的朋友说既然改变了亚松弛因子，就要从头开始算。--------不是我不想从头开始算，主要是我前面的数据已经连续计算5天（120小时）了，重头再来的话，工程浩大。 对于到底取多少松弛因子，有的朋友说多测试，只是问题又来了，我的一个数据计算得上十天，如果一个一个测试，那么时间。。。。。。？ 好了，我总结一下自己的问题吧。 1 亚松弛因子是否可以在计算中途改变，改变前后结果是否有继承性？ 2 对于亚松弛因子的设定，除了按照默认设定，经验设定，还有其他办法吗？ 啰嗦了半天，谢谢大家了 |
回复此楼» 猜你喜欢
2025冷门绝学什么时候出结果
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
fangdian
木虫 (正式写手)
- 应助: 201 (大学生)
- 金币: 2836.4
- 散金: 64
- 红花: 33
- 帖子: 818
- 在线: 206.9小时
- 虫号: 1926821
- 注册: 2012-08-06
- 专业: 多相流热物理学
|
当时不温度的一个文本文件，今天清理的时候发现的，复制上来，留作日后查看。 Updating solution at time level N... Global Courant Number : 0.30 done. physical-dt 2.0000e-05 iter continuity x-velocity y-velocity z-velocity energy time/iter !911520 solution is converged 911520 5.1930e-04 1.4487e-04 1.6887e-04 1.7274e-04 7.8469e-09 0:00:38 70 911530 1.6998e-02 1.7622e-03 1.6115e-03 1.8824e-03 6.7242e-09 0:00:32 60 911540 1.3505e-02 1.5497e-04 1.8145e-04 1.9544e-04 1.5172e-08 0:00:27 50 911550 3.0434e-03 1.9570e-05 2.7743e-05 2.8083e-05 1.3270e-08 0:00:21 40 911560 1.8486e-02 5.8436e-06 7.4465e-06 8.5965e-06 6.6203e-09 0:00:16 30 911570 1.3461e-02 5.0508e-06 6.1176e-06 7.5549e-06 1.5157e-08 0:00:10 20 911580 3.0135e-03 2.3962e-06 3.2288e-06 3.6347e-06 1.3268e-08 0:00:05 10 911590 1.8481e-02 3.7135e-06 3.9322e-06 5.2476e-06 6.6200e-09 0:00:00 0 不分NODE:999999,最高温度300.000000,第一次最高温度时间:1.374859,总最高温度2181.019891,速度0.004167 NODE:999999,最高温度300.000000,总最高温度2181.019891,熔池体积：2.6715e-08,速度0.004167,液面最高点z坐标0.003103,频率53.462099 不分NODE:0,最高温度1924.422149,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:1,最高温度1922.327451,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:2,最高温度2181.019891,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:3,最高温度997.674704,第一次最高温度时间:0.010000,总最高温度2181.019891,速度0.004167 不分NODE:4,最高温度2066.469698,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:5,最高温度1694.314250,第一次最高温度时间:0.010000,总最高温度2181.019891,速度0.004167 不分NODE:6,最高温度1884.067142,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:7,最高温度962.109162,第一次最高温度时间:0.010000,总最高温度2181.019891,速度0.004167 Flow time = 1.863469220648713s, time step = 47222 199878 more time steps Updating solution at time level N... Global Courant Number : 0.30 done. physical-dt 2.0000e-05 iter continuity x-velocity y-velocity z-velocity energy time/iter 911590 1.8481e-02 3.7135e-06 3.9322e-06 5.2476e-06 6.6200e-09 0:00:36 70 911600 2.0701e-02 1.3525e-03 1.2224e-03 1.3139e-03 8.7540e-09 0:00:31 60 911610 6.5794e-03 1.1610e-04 1.3655e-04 1.3304e-04 9.3469e-09 0:00:26 50 911620 2.1620e-02 1.9685e-05 2.5391e-05 2.5167e-05 1.2797e-08 0:00:21 40 911630 1.6690e-02 7.7434e-06 8.9832e-06 1.0102e-05 8.4713e-09 0:00:16 30 911640 6.4028e-03 3.6825e-06 4.4296e-06 5.2573e-06 9.1522e-09 0:00:10 20 911650 2.1597e-02 3.9669e-06 4.1066e-06 5.3098e-06 1.2712e-08 0:00:05 10 911660 1.6677e-02 4.7833e-06 5.0276e-06 6.4548e-06 8.4966e-09 0:00:00 0 不分NODE:999999,最高温度300.000000,第一次最高温度时间:1.374859,总最高温度2181.019891,速度0.004167 NODE:999999,最高温度300.000000,总最高温度2181.019891,熔池体积：2.6725e-08,速度0.004167,液面最高点z坐标0.003103,频率53.462099 不分NODE:0,最高温度1924.422149,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:1,最高温度1922.327451,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:2,最高温度2181.019891,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:3,最高温度997.674704,第一次最高温度时间:0.010000,总最高温度2181.019891,速度0.004167 不分NODE:4,最高温度2066.469698,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:5,最高温度1694.314250,第一次最高温度时间:0.010000,总最高温度2181.019891,速度0.004167 不分NODE:6,最高温度1884.067142,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:7,最高温度962.117140,第一次最高温度时间:0.010000,总最高温度2181.019891,速度0.004167 Flow time = 1.863489220648713s, time step = 47223 199877 more time steps Updating solution at time level N... Global Courant Number : 0.30 done. physical-dt 2.0000e-05 iter continuity x-velocity y-velocity z-velocity energy time/iter 911660 1.6677e-02 4.7833e-06 5.0276e-06 6.4548e-06 8.4966e-09 0:00:36 70 911670 1.4658e-02 1.4186e-03 1.1853e-03 1.3658e-03 5.5972e-09 0:00:30 60 911680 8.1662e-03 1.5166e-04 1.5871e-04 1.6975e-04 1.1039e-08 0:00:25 50 911690 4.2240e-03 2.6927e-05 3.4148e-05 3.5607e-05 5.5006e-09 0:00:21 40 911700 7.5599e-03 7.9561e-06 1.0336e-05 1.0993e-05 1.1048e-08 0:00:16 30 911710 4.1532e-03 5.2922e-06 5.8915e-06 7.1897e-06 5.5026e-09 0:00:10 20 911720 7.5459e-03 3.1028e-06 3.7315e-06 4.4011e-06 1.1048e-08 0:00:05 10 911730 4.1436e-03 4.0764e-06 4.2068e-06 5.3198e-06 5.5027e-09 0:00:00 0 不分NODE:999999,最高温度300.000000,第一次最高温度时间:1.374859,总最高温度2181.019891,速度0.004167 NODE:999999,最高温度300.000000,总最高温度2181.019891,熔池体积：2.6753e-08,速度0.004167,液面最高点z坐标0.003103,频率53.462099 不分NODE:0,最高温度1924.422149,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:1,最高温度1922.327451,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:2,最高温度2181.019891,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:3,最高温度997.674704,第一次最高温度时间:0.010000,总最高温度2181.019891,速度0.004167 不分NODE:4,最高温度2066.469698,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:5,最高温度1694.314250,第一次最高温度时间:0.010000,总最高温度2181.019891,速度0.004167 不分NODE:6,最高温度1884.067142,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:7,最高温度962.125119,第一次最高温度时间:0.010000,总最高温度2181.019891,速度0.004167 Flow time = 1.863509220648713s, time step = 47224 199876 more time steps Updating solution at time level N... Global Courant Number : 0.30 done. physical-dt 2.0000e-05 iter continuity x-velocity y-velocity z-velocity energy time/iter 911730 4.1436e-03 4.0764e-06 4.2068e-06 5.3198e-06 5.5027e-09 0:00:36 70 911740 5.8646e-02 1.5540e-03 1.2590e-03 1.5741e-03 2.9797e-08 0:00:30 60 911750 1.9588e-02 1.7088e-04 1.9684e-04 2.1220e-04 3.3534e-08 0:00:25 50 911760 4.3004e-02 6.6412e-05 8.3443e-05 9.6605e-05 2.2230e-08 0:00:20 40 911770 4.5285e-02 5.5340e-05 7.4715e-05 9.2340e-05 2.9122e-08 0:00:16 30 911780 1.8994e-02 4.7039e-05 6.1886e-05 7.7599e-05 3.3647e-08 0:00:10 20 911790 4.3052e-02 5.1416e-05 6.5150e-05 7.5877e-05 2.2318e-08 0:00:05 10 911800 4.5325e-02 5.2280e-05 7.1187e-05 8.9645e-05 2.9180e-08 0:00:00 0 不分NODE:999999,最高温度300.000000,第一次最高温度时间:1.374859,总最高温度2181.019891,速度0.004167 NODE:999999,最高温度300.000000,总最高温度2181.019891,熔池体积：2.6773e-08,速度0.004167,液面最高点z坐标0.003103,频率53.462099 不分NODE:0,最高温度1924.422149,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:1,最高温度1922.327451,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:2,最高温度2181.019891,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:3,最高温度997.674704,第一次最高温度时间:0.010000,总最高温度2181.019891,速度0.004167 不分NODE:4,最高温度2066.469698,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:5,最高温度1694.314250,第一次最高温度时间:0.010000,总最高温度2181.019891,速度0.004167 不分NODE:6,最高温度1884.067142,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:7,最高温度962.133099,第一次最高温度时间:0.010000,总最高温度2181.019891,速度0.004167 Flow time = 1.863529220648713s, time step = 47225 199875 more time steps Updating solution at time level N... Global Courant Number : 0.30 done. physical-dt 2.0000e-05 iter continuity x-velocity y-velocity z-velocity energy time/iter 911800 4.5325e-02 5.2280e-05 7.1187e-05 8.9645e-05 2.9180e-08 0:00:36 70 911810 3.9618e-02 1.5001e-03 1.2424e-03 1.5815e-03 3.0267e-08 0:00:31 60 911820 7.1477e-02 2.2057e-04 2.7712e-04 3.4080e-04 2.4295e-08 0:00:25 50 911830 5.6085e-02 1.4368e-04 1.5289e-04 2.2288e-04 3.4023e-08 0:00:21 40 911840 3.9570e-02 1.0590e-04 1.3124e-04 1.6023e-04 2.9319e-08 0:00:16 30 911850 7.1773e-02 9.3111e-05 1.5886e-04 1.9296e-04 2.4062e-08 0:00:10 20 911860 5.5813e-02 1.1915e-04 1.2784e-04 1.9573e-04 3.3921e-08 0:00:05 10 911870 3.9409e-02 1.0313e-04 1.2780e-04 1.5650e-04 2.9272e-08 0:00:00 0 不分NODE:999999,最高温度300.000000,第一次最高温度时间:1.374859,总最高温度2181.019891,速度0.004167 NODE:999999,最高温度300.000000,总最高温度2181.019891,熔池体积：2.6848e-08,速度0.004167,液面最高点z坐标0.003103,频率53.462099 不分NODE:0,最高温度1924.422149,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:1,最高温度1922.327451,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:2,最高温度2181.019891,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:3,最高温度997.674704,第一次最高温度时间:0.010000,总最高温度2181.019891,速度0.004167 不分NODE:4,最高温度2066.469698,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:5,最高温度1694.314250,第一次最高温度时间:0.010000,总最高温度2181.019891,速度0.004167 不分NODE:6,最高温度1884.067142,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:7,最高温度962.141079,第一次最高温度时间:0.010000,总最高温度2181.019891,速度0.004167 Flow time = 1.863549220648713s, time step = 47226 199874 more time steps Updating solution at time level N... Global Courant Number : 0.30 done. physical-dt 2.0000e-05 iter continuity x-velocity y-velocity z-velocity energy time/iter 911870 3.9409e-02 1.0313e-04 1.2780e-04 1.5650e-04 2.9272e-08 0:00:36 70 911880 5.0207e-01 1.8983e-03 2.0559e-03 3.2481e-03 2.4034e-08 0:00:31 60 911890 8.2522e+02 1.1072e-01 1.3333e-01 4.6064e-01 1.8953e-08 0:00:26 50 911900 1.6104e+06 1.1591e-01 1.3916e-01 4.9838e-01 1.3101e-08 0:00:21 40 911910 3.1165e+09 1.1573e-01 1.3958e-01 4.9898e-01 2.3098e-08 0:00:16 30 911920 6.1243e+12 1.1577e-01 1.3964e-01 4.9907e-01 1.8810e-08 0:00:10 20 911930 1.1945e+16 1.1590e-01 1.3916e-01 4.9839e-01 1.3080e-08 0:00:05 10 911940 2.3117e+19 1.1573e-01 1.3958e-01 4.9898e-01 2.3089e-08 0:00:00 0 不分NODE:999999,最高温度300.000000,第一次最高温度时间:1.374859,总最高温度2181.019891,速度0.004167 NODE:999999,最高温度300.000000,总最高温度2181.019891,熔池体积：2.6832e-08,速度0.004167,液面最高点z坐标0.003103,频率53.462099 不分NODE:0,最高温度1924.422149,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:1,最高温度1922.327451,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:2,最高温度2181.019891,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:3,最高温度997.674704,第一次最高温度时间:0.010000,总最高温度2181.019891,速度0.004167 不分NODE:4,最高温度2066.469698,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:5,最高温度1694.314250,第一次最高温度时间:0.010000,总最高温度2181.019891,速度0.004167 不分NODE:6,最高温度1884.067142,第一次最高温度时间:1.861049,总最高温度2181.019891,速度0.004167 不分NODE:7,最高温度962.149060,第一次最高温度时间:0.010000,总最高温度2181.019891,速度0.004167 Flow time = 1.863569220648712s, time step = 47227 199873 more time steps Updating solution at time level N... Global Courant Number : 33311284309224931000.00 Primitive Error at Node 0: Global courant number is greater than 250.00 The velocity field is probably diverging. Please check the solution and reduce the time-step if necessary. Primitive Error at Node 1: Global courant number is greater than 250.00 The velocity field is probably diverging. Please check the solution and reduce the time-step if necessary. Primitive Error at Node 2: Global courant number is greater than 250.00 The velocity field is probably diverging. Please check the solution and reduce the time-step if necessary. Primitive Error at Node 3: Global courant number is greater than 250.00 The velocity field is probably diverging. Please check the solution and reduce the time-step if necessary. Primitive Error at Node 4: Global courant number is greater than 250.00 The velocity field is probably diverging. Please check the solution and reduce the time-step if necessary. Primitive Error at Node 5: Global courant number is greater than 250.00 The velocity field is probably diverging. Please check the solution and reduce the time-step if necessary. Primitive Error at Node 6: Global courant number is greater than 250.00 The velocity field is probably diverging. Please check the solution and reduce the time-step if necessary. Primitive Error at Node 7: Global courant number is greater than 250.00 The velocity field is probably diverging. Please check the solution and reduce the time-step if necessary. Interrupting... Interrupting client... |
23楼2012-11-28 17:26:31
fangdian2楼2012-11-06 09:22:40
fangdian
木虫 (正式写手)
- 应助: 201 (大学生)
- 金币: 2836.4
- 散金: 64
- 红花: 33
- 帖子: 818
- 在线: 206.9小时
- 虫号: 1926821
- 注册: 2012-08-06
- 专业: 多相流热物理学
3楼2012-11-06 10:13:25
sh073
木虫 (正式写手)
- 应助: 42 (小学生)
- 金币: 3739.6
- 散金: 51
- 红花: 4
- 帖子: 475
- 在线: 206.8小时
- 虫号: 1695843
- 注册: 2012-03-16
- 专业: 分离过程
4楼2012-11-06 10:54:02