发现增加网格,时间步长如果不降低的话,很容易导致这个CFL数超过1,很容易出现发散。这里出现了一个计算到0.9s还出现发散的案例,small 2d-3.5w网格,这个大概率是时间步长没有降低引起的。
还会出现下面的报错。
Compute processes interrupted. Processing can be resumed.
Updating solution at time level N...
Global Courant Number [Explicit VOF Criteria] : 1.60472
done.
iter continuity x-velocity y-velocity energy time/iter
7580790 1.4505e-05 6.4254e-12 2.8694e-12 2.3550e-12 0:00:01 100
Reversed flow on 62 faces of pressure-outlet 10.
7580800 4.6609e-01 1.6361e-06 1.8882e-06 6.0528e-09 0:00:01 90
Reversed flow on 62 faces of pressure-outlet 10.
7580810 5.5459e-02 2.0586e-07 1.5314e-07 9.6171e-11 0:00:01 80
Reversed flow on 62 faces of pressure-outlet 10.
7580820 7.2041e-03 2.6118e-08 2.0835e-08 5.1357e-11 0:00:01 70
Reversed flow on 62 faces of pressure-outlet 10.
7580830 1.2327e-03 3.9607e-09 2.4394e-09 3.6777e-11 0:00:01 60
Reversed flow on 62 faces of pressure-outlet 10.
7580840 3.2143e-04 7.1081e-10 4.7496e-10 2.6282e-11 0:00:01 50
Reversed flow on 62 faces of pressure-outlet 10.
7580850 1.0830e-04 1.9333e-10 1.1204e-10 1.7741e-11 0:00:01 40
Reversed flow on 62 faces of pressure-outlet 10.
7580860 4.5116e-05 5.4563e-11 2.8222e-11 1.1468e-11 0:00:00 30
Reversed flow on 62 faces of pressure-outlet 10.
7580870 2.0378e-05 1.4963e-11 8.3435e-12 7.4702e-12 0:00:00 20
Reversed flow on 62 faces of pressure-outlet 10.
!7580880 solution is converged
7580880 9.7816e-06 4.6136e-12 2.9584e-12 5.0214e-12 0:00:00 10
step flow-time heat-011 heat-010 heat-001 heat-000
93361 9.3361e-01 3.2229e+02 3.2283e+02 3.2179e+02 3.2106e+02
Flow time = 0.9336100000000001s, time step = 93361
6639 more time steps
Updating solution at time level N...
Global Courant Number [Explicit VOF Criteria] : 1.52448
done.
iter continuity x-velocity y-velocity energy time/iter
!7580880 solution is converged
7580880 9.7816e-06 4.6136e-12 2.9584e-12 5.0214e-12 0:00:01 100
Reversed flow on 62 faces of pressure-outlet 10.
7580890 1.0640e+00 2.2608e-06 3.9072e-06 5.8772e-09 0:00:01 90
Reversed flow on 63 faces of pressure-outlet 10.
7580900 9.2296e-02 3.1741e-07 5.0238e-07 3.6264e-11 0:00:01 80
Reversed flow on 63 faces of pressure-outlet 10.
7580910 1.5719e-02 3.7620e-08 5.6747e-08 2.2916e-11 0:00:01 70
Reversed flow on 63 faces of pressure-outlet 10.
7580920 3.7151e-03 6.9992e-09 1.0067e-08 1.8502e-11 0:00:01 60
Reversed flow on 63 faces of pressure-outlet 10.
7580930 1.1597e-03 1.8276e-09 2.7730e-09 1.3806e-11 0:00:01 50
Reversed flow on 63 faces of pressure-outlet 10.
7580940 4.1138e-04 5.0089e-10 8.9288e-10 1.0004e-11 0:00:01 40
Reversed flow on 63 faces of pressure-outlet 10.
7580950 1.5426e-04 1.5255e-10 3.0166e-10 7.1869e-12 0:00:00 30
Reversed flow on 63 faces of pressure-outlet 10.
7580960 6.1797e-05 5.1330e-11 9.2432e-11 5.1652e-12 0:00:00 20
Reversed flow on 63 faces of pressure-outlet 10.
7580970 2.6470e-05 1.7624e-11 3.1436e-11 3.7279e-12 0:00:00 10
Reversed flow on 63 faces of pressure-outlet 10.
7580980 1.2003e-05 7.1216e-12 1.1102e-11 2.7058e-12 0:00:00 0
step flow-time heat-011 heat-010 heat-001 heat-000
93362 9.3362e-01 3.2229e+02 3.2283e+02 3.2179e+02 3.2106e+02
Flow time = 0.9336200000000001s, time step = 93362
6638 more time steps
Updating solution at time level N...
Global Courant Number [Explicit VOF Criteria] : 2.07392
done.
iter continuity x-velocity y-velocity energy time/iter
7580980 1.2003e-05 7.1216e-12 1.1102e-11 2.7058e-12 0:00:01 100
Reversed flow on 63 faces of pressure-outlet 10.
7580990 8.4343e-01 2.9223e-06 2.2536e-06 5.9416e-09 0:00:01 90
Reversed flow on 63 faces of pressure-outlet 10.
7581000 6.9794e-02 2.7402e-07 3.5720e-07 3.0543e-11 0:00:01 80
Reversed flow on 63 faces of pressure-outlet 10.
7581010 1.1744e-02 4.1272e-08 1.0056e-07 1.4124e-11 0:00:01 70
Reversed flow on 63 faces of pressure-outlet 10.
7581020 4.7136e-03 1.5134e-08 4.4993e-08 9.5775e-12 0:00:01 60
Reversed flow on 63 faces of pressure-outlet 10.
7581030 2.1203e-03 6.9628e-09 2.2293e-08 6.4517e-12 0:00:01 50
Reversed flow on 63 faces of pressure-outlet 10.
7581040 1.1231e-03 3.5059e-09 1.1520e-08 4.3930e-12 0:00:00 40
Reversed flow on 63 faces of pressure-outlet 10.
7581050 6.0654e-04 1.8433e-09 6.1044e-09 3.0336e-12 0:00:00 30
Reversed flow on 63 faces of pressure-outlet 10.
7581060 3.2631e-04 9.7742e-10 3.2652e-09 2.1246e-12 0:00:00 20
Reversed flow on 63 faces of pressure-outlet 10.
7581070 1.7593e-04 5.2401e-10 1.7559e-09 1.5079e-12 0:00:00 10
Reversed flow on 63 faces of pressure-outlet 10.
7581080 9.4849e-05 2.8213e-10 9.4623e-10 1.0831e-12 0:00:00 0
step flow-time heat-011 heat-010 heat-001 heat-000
93363 9.3363e-01 3.2229e+02 3.2283e+02 3.2179e+02 3.2106e+02
Flow time = 0.9336300000000001s, time step = 93363
6637 more time steps
Updating solution at time level N...
Global Courant Number [Explicit VOF Criteria] : 2.19464
done.
iter continuity x-velocity y-velocity energy time/iter
7581080 9.4849e-05 2.8213e-10 9.4623e-10 1.0831e-12 0:00:01 100
Reversed flow on 63 faces of pressure-outlet 10.
7581090 2.1846e+00 3.4272e-06 3.3522e-06 7.2802e-09 0:00:01 90
Reversed flow on 63 faces of pressure-outlet 10.
7581100 2.7836e+00 1.9524e-06 6.2758e-06 4.7048e-09 0:00:01 80
Reversed flow on 59 faces of pressure-outlet 10.
7581110 1.5298e+01 2.9112e-06 3.0553e-05 1.5473e-09 0:00:01 70
Reversed flow on 62 faces of pressure-outlet 10.
7581120 2.8817e+01 6.0194e-06 5.5918e-05 4.0947e-09 0:00:01 60
Reversed flow on 67 faces of pressure-outlet 10.
7581130 1.1000e+02 3.0529e-05 2.0806e-04 5.6595e-09 0:00:01 50
Reversed flow on 91 faces of pressure-outlet 10.
7581140 3.5056e+01 3.2985e-05 6.5153e-05 1.8169e-09 0:00:00 40
Reversed flow on 91 faces of pressure-outlet 10.
7581150 4.7813e+02 9.9444e-05 7.8565e-04 4.7113e-09 0:00:00 30
Reversed flow on 92 faces of pressure-outlet 10.
7581160 4.7954e+02 1.5305e-04 7.3062e-04 5.9567e-09 0:00:00 20
Reversed flow on 92 faces of pressure-outlet 10.
7581170 2.6745e+03 3.1434e-04 2.4911e-03 1.9006e-09 0:00:00 10
Reversed flow on 93 faces of pressure-outlet 10.
7581180 2.4283e+03 3.7562e-04 7.0808e-04 4.8893e-09 0:00:00 0
step flow-time heat-011 heat-010 heat-001 heat-000
93364 9.3364e-01 3.2229e+02 3.2283e+02 3.2179e+02 3.2106e+02
Flow time = 0.93364s, time step = 93364
6636 more time steps
===============Message from the Cortex Process================================
Compute processes interrupted. Processing can be resumed.
==============================================================================
Updating solution at time level N...
Global Courant Number [Explicit VOF Criteria] : 2567.17
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.
Wed Apr 10 16:34:31 CST 2024
2024/04/10 16:34:36 INFO Demo license file /share/software/apps/ansys_inc/shared_files/licensing/license_files/ansyslmd.lic exists.