fluent出现下面的错误(Error at Node 0: floating point exception)发散,主要和连续性方程收敛性不好有关,初始为1e+0。还有Global Courant Number超过1也不稳定,很容易出现发散。该case网格为2.6w。不知道为何,前面的1.5w和后面的3.5w都没有出现发散,这个却出现了发散,其他的都一样,有时候感觉这个CFD计算就是玄学,很难控制,头疼。怪不得做实验的会不相信CFD,就像中医一样,太玄了。有高手可以解决吗?
这个.err输出提示如下。
Error at Node 1: floating point exception
Error at Node 2: floating point exception
Error at Node 3: floating point exception
Error at Node 4: floating point exception
Error at Node 5: floating point exception
Error at Node 6: floating point exception
Error at Node 7: floating point exception
Error at Node 8: floating point exception
Error at Node 9: floating point exception
Error at Node 10: floating point exception
Error at Node 11: floating point exception
Error at Node 12: floating point exception
Error at Node 13: floating point exception
Error at Node 14: floating point exception
Error at Node 15: floating point exception
Error at Node 16: floating point exception
Error at Node 17: floating point exception
Error at Node 18: floating point exception
Error at Node 19: floating point exception
Error at Node 20: floating point exception
Error at Node 21: floating point exception
Error at Node 22: floating point exception
Error at Node 23: floating point exception
Error at Node 24: floating point exception
Error at Node 25: floating point exception
Error at Node 26: floating point exception
Error at Node 27: floating point exception
Error at Node 28: floating point exception
Error at Node 29: floating point exception
Error at Node 30: floating point exception
Error at Node 31: floating point exception
Error at Node 32: floating point exception
Error at Node 33: floating point exception
Error at Node 34: floating point exception
Error at Node 35: floating point exception
Error at Node 36: floating point exception
Error at Node 37: floating point exception
Error at Node 38: floating point exception
Error at Node 39: floating point exception
Error at Node 40: floating point exception
Error at Node 41: floating point exception
Error at Node 42: floating point exception
Error at Node 43: floating point exception
Error at Node 44: floating point exception
Error at Node 45: floating point exception
Error at Node 46: floating point exception
Error at Node 47: floating point exception
Error at Node 48: floating point exception
Error at Node 49: floating point exception
Error at Node 50: floating point exception
Error at Node 51: floating point exception
Error at Node 52: floating point exception
Error at Node 53: floating point exception
Error at Node 54: floating point exception
Error at Node 55: floating point exception
Error at host: floating point exception
Error: floating point exception
Error Object: #f
下面为出现上面错误的fluent输出的out文件,里面有监测的残差值。
Updating solution at time level N...
Global Courant Number [Explicit VOF Criteria] : 1.3063
done.
iter continuity x-velocity y-velocity energy time/iter
366210 2.8559e-04 2.6830e-10 3.2377e-10 3.8046e-10 0:00:01 100
366220 4.6132e+00 4.7100e-06 2.5978e-06 5.1215e-08 0:00:01 90
366230 7.7570e-01 5.6073e-07 7.5951e-07 1.0598e-08 0:00:01 80
366240 2.0724e-01 1.3639e-07 1.6916e-07 4.3341e-09 0:00:01 70
366250 6.7063e-02 3.9896e-08 5.2783e-08 2.4767e-09 0:00:01 60
366260 2.2692e-02 1.0148e-08 1.6518e-08 1.6371e-09 0:00:00 50
366270 9.4862e-03 4.0784e-09 5.8462e-09 1.1559e-09 0:00:00 40
366280 4.0637e-03 1.5603e-09 2.2459e-09 8.4378e-10 0:00:00 30
366290 1.7522e-03 6.2017e-10 8.9443e-10 6.2242e-10 0:00:00 20
366300 7.6186e-04 2.5292e-10 3.6839e-10 4.6061e-10 0:00:00 10
366310 3.3322e-04 1.0607e-10 1.5391e-10 3.4142e-10 0:00:00 0
step flow-time heat-011 heat-010 heat-001 heat-000
7080 7.0800e-02 3.2857e+02 3.2359e+02 3.2227e+02 3.2125e+02
Flow time = 0.07080000000000002s, time step = 7080
92922 more time steps
Updating solution at time level N...
Global Courant Number [Explicit VOF Criteria] : 2.28101
done.
iter continuity x-velocity y-velocity energy time/iter
366310 3.3322e-04 1.0607e-10 1.5391e-10 3.4142e-10 0:00:01 100
366320 4.8766e+00 4.2570e-06 2.9602e-06 4.1813e-08 0:00:01 90
366330 7.6210e-01 9.7182e-07 6.4520e-07 7.6455e-09 0:00:01 80
366340 2.6377e-01 4.3941e-07 3.3542e-07 4.4046e-09 0:00:01 70
366350 8.3969e-02 1.3088e-07 1.2110e-07 2.9210e-09 0:00:01 60
366360 3.3536e-02 4.9237e-08 4.2995e-08 2.0419e-09 0:00:01 50
366370 1.3501e-02 1.6594e-08 9.6834e-09 1.4659e-09 0:00:01 40
366380 5.5080e-03 5.7036e-09 5.3286e-09 1.0664e-09 0:00:00 30
366390 2.3171e-03 2.4907e-09 2.1193e-09 7.8120e-10 0:00:00 20
366400 9.6322e-04 8.2802e-10 7.9221e-10 5.7491e-10 0:00:00 10
366410 4.2601e-04 3.8377e-10 3.3231e-10 4.2444e-10 0:00:00 0
step flow-time heat-011 heat-010 heat-001 heat-000
7081 7.0810e-02 3.2856e+02 3.2358e+02 3.2227e+02 3.2125e+02
Flow time = 0.07081000000000001s, time step = 7081
92921 more time steps
Updating solution at time level N...
Global Courant Number [Explicit VOF Criteria] : 3.48722
done.
iter continuity x-velocity y-velocity energy time/iter
366410 4.2601e-04 3.8377e-10 3.3231e-10 4.2444e-10 0:00:01 100
366420 1.3333e+01 1.3856e-05 2.5209e-05 4.5945e-06 0:00:01 90
366430 2.6263e+01 5.2787e-05 4.6978e-05 3.4613e-06 0:00:01 80
366440 1.6964e+02 2.9984e-04 3.5737e-04 5.3731e-06 0:00:01 70
temperature limited to 1.000000e+00 in 1 cells on zone 28 in domain 1
temperature limited to 1.000000e+00 in 6 cells on zone 28 in domain 1
temperature limited to 1.000000e+00 in 3 cells on zone 28 in domain 1
temperature limited to 1.000000e+00 in 4 cells on zone 28 in domain 1
temperature limited to 1.000000e+00 in 2 cells on zone 28 in domain 1
temperature limited to 1.000000e+00 in 2 cells on zone 28 in domain 1
366450 4.5317e+02 8.1187e-04 6.5988e-04 1.3082e-05 0:00:01 60
temperature limited to 1.000000e+00 in 2 cells on zone 28 in domain 1
temperature limited to 1.000000e+00 in 1 cells on zone 28 in domain 1
temperature limited to 1.000000e+00 in 12 cells on zone 28 in domain 1
366460 3.4402e+03 1.7428e-03 3.5660e-03 4.7011e-06 0:00:01 50
temperature limited to 5.000000e+03 in 1 cells on zone 28 in domain 1
temperature limited to 1.000000e+00 in 3 cells on zone 28 in domain 1
temperature limited to 1.000000e+00 in 1 cells on zone 28 in domain 1
temperature limited to 1.000000e+00 in 12 cells on zone 28 in domain 1
Reversed flow on 97 faces of pressure-outlet 10.
366470 3.9692e+04 7.9946e-04 2.6420e-03 9.8477e-05 0:00:00 40
temperature limited to 1.000000e+00 in 5 cells on zone 28 in domain 1
temperature limited to 1.000000e+00 in 2 cells on zone 28 in domain 1
temperature limited to 1.000000e+00 in 2 cells on zone 28 in domain 1
temperature limited to 1.000000e+00 in 2 cells on zone 28 in domain 1
temperature limited to 1.000000e+00 in 1 cells on zone 28 in domain 1
Reversed flow on 105 faces of pressure-outlet 10.
366480 7.8534e+05 8.4420e-04 2.6200e-03 7.8447e-06 0:00:00 30
temperature limited to 1.000000e+00 in 1 cells on zone 28 in domain 1
temperature limited to 1.000000e+00 in 1 cells on zone 28 in domain 1
Reversed flow on 105 faces of pressure-outlet 10.
366490 1.5482e+07 7.8912e-04 2.5800e-03 2.1943e-05 0:00:00 20
temperature limited to 1.000000e+00 in 2 cells on zone 28 in domain 1
Stabilizing pressure coupled to enhance linear solver robustness.
temperature limited to 1.000000e+00 in 1 cells on zone 28 in domain 1
temperature limited to 1.000000e+00 in 2 cells on zone 28 in domain 1
Reversed flow on 99 faces of pressure-outlet 10.
366500 1.6071e+120 1.0588e-03 6.6939e+55 7.4031e-05 0:00:00 10
Stabilizing pressure coupled to enhance linear solver robustness.
Stabilizing pressure coupled using GMRES to enhance linear solver robustness.
Experiencing convergence difficulties - temporarily relaxing and trying again...
Experiencing convergence difficulties - temporarily relaxing and trying again...
Experiencing convergence difficulties - temporarily relaxing and trying again...
Experiencing convergence difficulties - temporarily relaxing and trying again...
Stabilizing pressure coupled to enhance linear solver robustness.
Stabilizing pressure coupled using GMRES to enhance linear solver robustness.
Experiencing convergence difficulties - temporarily relaxing and trying again...
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Reversed flow on 96 faces of pressure-outlet 10.
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
Divergence detected in AMG solver: pressure coupled
===============Message from the Cortex Process================================
Compute processes interrupted. Processing can be resumed.
==============================================================================
Warning: An error or interrupt occurred while reading the journal file.
Some commands may not have been completed.
> Halting due to end of file on input.
Mon Apr 8 18:01:06 CST 2024
2024/04/08 18:01:08 INFO Demo license file /share/software/apps/ansys_inc/shared_files/licensing/license_files/ansyslmd.lic exists.
下面为一个正常的计算fluent输出的out文件。
Flow time = 0.33899s, time step = 33899
66101 more time steps
Updating solution at time level N...
Global Courant Number [Explicit VOF Criteria] : 0.342811
done.
iter continuity x-velocity y-velocity energy time/iter
!2604010 solution is converged
2604010 5.4795e-06 6.1640e-12 6.0958e-12 6.1643e-16 0:00:01 100
Reversed flow on 42 faces of pressure-outlet 10.
2604020 1.7271e-01 8.4543e-07 7.2380e-07 2.5860e-09 0:00:01 90
Reversed flow on 41 faces of pressure-outlet 10.
2604030 1.7611e-02 8.2713e-08 6.5798e-08 5.5610e-12 0:00:01 80
Reversed flow on 41 faces of pressure-outlet 10.
2604040 3.1860e-03 8.5239e-09 9.6657e-09 2.4607e-13 0:00:01 70
Reversed flow on 41 faces of pressure-outlet 10.
2604050 9.1098e-04 1.8133e-09 1.9492e-09 7.0191e-14 0:00:01 60
Reversed flow on 41 faces of pressure-outlet 10.
2604060 3.2160e-04 5.1264e-10 5.3389e-10 2.5191e-14 0:00:01 50
Reversed flow on 41 faces of pressure-outlet 10.
2604070 1.2701e-04 1.5807e-10 1.6934e-10 1.0070e-14 0:00:01 40
Reversed flow on 41 faces of pressure-outlet 10.
2604080 5.3259e-05 5.1947e-11 5.6763e-11 4.2821e-15 0:00:00 30
Reversed flow on 41 faces of pressure-outlet 10.
2604090 2.3206e-05 1.8130e-11 1.9562e-11 1.9011e-15 0:00:00 20
Reversed flow on 41 faces of pressure-outlet 10.
2604100 1.0389e-05 6.6094e-12 6.9218e-12 8.8240e-16 0:00:00 10
Reversed flow on 41 faces of pressure-outlet 10.
!2604110 solution is converged
2604110 4.7438e-06 2.5105e-12 2.5263e-12 4.3602e-16 0:00:00 0
step flow-time heat-011 heat-010 heat-001 heat-000
33900 3.3900e-01 3.2169e+02 3.2342e+02 3.2207e+02 3.1992e+02
Flow time = 0.339s, time step = 33900
66100 more time steps
Updating solution at time level N...
Global Courant Number [Explicit VOF Criteria] : 0.330329
done.
iter continuity x-velocity y-velocity energy time/iter
!2604110 solution is converged
2604110 4.7438e-06 2.5105e-12 2.5263e-12 4.3602e-16 0:00:01 100
Reversed flow on 41 faces of pressure-outlet 10.
2604120 2.0064e-01 8.5213e-07 7.3077e-07 2.5816e-09 0:00:01 90
Reversed flow on 41 faces of pressure-outlet 10.
2604130 2.0751e-02 7.5992e-08 7.4416e-08 5.2422e-12 0:00:01 80
Reversed flow on 41 faces of pressure-outlet 10.
2604140 3.8967e-03 1.0739e-08 1.2610e-08 1.5227e-13 0:00:01 70
Reversed flow on 41 faces of pressure-outlet 10.
2604150 1.1279e-03 2.2763e-09 2.8113e-09 3.2554e-14 0:00:01 60