发现增加网格,时间步长如果不降低的话,很容易导致这个CFL数超过1,很容易出现发散。这里出现了一个计算到0.9s还出现发散的案例,small 2d-3.5w网格,这个大概率是时间步长没有降低引起的。

文章描述了在使用ANSYS等软件进行数值模拟时遇到的问题,发散现象源于时间步长过大导致CFL数超过1。文中提到通过更新解决方案和调整时间步长来控制计算过程,同时提到了计算过程中出现的中断及其处理方法。
摘要由CSDN通过智能技术生成

发现增加网格,时间步长如果不降低的话,很容易导致这个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.

  • 22
    点赞
  • 15
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值