自动调节时间步长循环计算某一时间步

[merlinizlf]

各位老师好：

用自动调节时间步的时候经常发现一直在循环计算某一个时间，这是什么原因，是发散了吗？
因为一开始测试为了简单和利于收敛，设置的层流条件，实际上有部分地方应该是湍流，这个会影响吗？

如果降低进口速度，deltaT，有助于改善吗？

[李东岳]

@merlinizlf 在 自动调节时间步长循环计算某一时间步 中说：

一直在循环计算某一个时间

这不是闹鬼了么？时间旅行？一直循环在1-2秒？算完了2秒又回到1秒了？

[merlinizlf]

不好意思老师，是我说的不清楚，发现一直在计算同一个时间，比如一直计算0.1111秒

[merlinizlf]

@李东岳
老师好，是我没有表达清楚，是类似于这样
库朗数和deltaT都非常非常小，计算非常慢，到最后deltaT基本变成0了
我设置的最大库朗数是0.9，但是计算出来才0.0012，不应该自动缩小时间步长呀

请问老师，这个应该怎么解决，如何提高计算速度？

Courant Number mean: 3.58164e-05 max: 0.901511
Interface Courant Number mean: 2.94865e-08 max: 0.00119605
deltaT = 1.7962e-05
Time = 2.01626

MULES: Solving for alpha.Sarin
Phase-1 volume fraction = 6.16252e-05  Min(alpha.Sarin) = 0  Max(alpha.Sarin) = 0.109368
MULES: Solving for alpha.Sarin
Phase-1 volume fraction = 6.16252e-05  Min(alpha.Sarin) = 0  Max(alpha.Sarin) = 0.109368
smoothSolver:  Solving for alpha.Sarin, Initial residual = 5.46198e-10, Final residual = 5.46198e-10, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 0.00020754, Final residual = 9.85112e-08, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.000110383, Final residual = 4.3286e-08, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 0.000191943, Final residual = 1.79758e-06, No Iterations 1
GAMG:  Solving for p_rgh, Initial residual = 0.00248148, Final residual = 1.41595e-05, No Iterations 6
time step continuity errors : sum local = 3.84475e-13, global = 1.89675e-14, cumulative = -3.04519e-09
GAMG:  Solving for p_rgh, Initial residual = 0.00152746, Final residual = 5.47895e-06, No Iterations 7
time step continuity errors : sum local = 1.48826e-13, global = -6.48062e-15, cumulative = -3.04519e-09
ExecutionTime = 680.02 s  ClockTime = 680 s

Courant Number mean: 3.57663e-05 max: 0.90151
Interface Courant Number mean: 2.94446e-08 max: 0.00119381
deltaT = 1.79319e-05
Time = 2.01628

[李东岳]

就是发散了，deltaT太小了，比如1e-15, 1e-16就没有步进

[merlinizlf]

@李东岳

非常感谢老师，虽然库朗数还是很小，但是其实已经发散了是吗？

我用的是层流求解，用的还是upwind格式，下一步是不是只能降低进口流速了？

[OItoCFD]

@merlinizlf 这就是发散了 但是我遇到这种问题时候 经常跟网格有关 我网格加密 然后尽量使他每个方向长度均匀 就好一些

[merlinizlf]

@OItoCFD 非常感谢，我们也准备加密网格，但是具体加密的程度还不知道怎么把握，您这边是使所有网格大小尽量均一吗？

[OItoCFD]

@merlinizlf 我就是加密一点后 然后尽可能让网格均匀 网格别太细长 就好了 当然也有时候是其他原因导致的 但是毕竟改网格是最直观容易的第一个尝试选择

[vien锁]

@李东岳 请问老师，可以实现计算一定步数后，时间步倒退，重新计算一遍嘛？我想尝试计算1000时间步，然后从第500步计算到1500步，再从第1000步计算到2000步，这样可以实现嘛？