七十七、Fluent自适应时间步长操作详解1

1. 自适应时间步长概念

对于瞬态计算来说，时间步长的选择一直是个难点。时间步长较大可以加快计算速度，但可能会降低收敛性，甚至发散。时间步长较小，虽然收敛性会变好，但计算速度会变化。(有时，时间步长较小反而会发散)。

关于时间步长的选取，可参考文章：二十九、Fluent瞬态时间步长与迭代步数的讨论

能不能自动调整时间步长呢？使求解比较精细时选用较小的时间步长，而计算接近稳态时选用较大的时间步长。这样既能保证计算的收敛性，还能提高计算速度。因此出现了自适应时间步长。

2. 自适应时间步长的方法

Fluent提供了三种自适应时间步长的方法：

a. Adaptive界面自适应时间步长

b. 表达式调节

c. UDF调节

其中第一种方法是Fluent内置的，也是最简单的，只需要在Run calculation界面设置即可。

第二种和第三种方法需要写代码来调节，有一定的入门要求。这个我们会在下篇文章介绍

3. Adaptive界面自适应时间步长

Fluent提供了自适应时间步长的设置，对于瞬态计算，在Type处选择Adaptive即可设置自适应时间步长。

3.1 Method

Method:选择设置时间步长的依据。对于Adaptive自适应时间步长有三种选择。

a. CFL-Based：根据库朗数来调节时间步长，CFL-Based在基于密度显示格式时不适用，也不适用于多相流模型。

选择CFL-Based需要设置库朗数，Fluent将根据设置的库朗数自动调整时间步长，使计算出的库朗数小于设置值。默认Courant Number为1。

b. Error-Based：根据计算截断误差来调节时间步长，不适用基于密度的显式格式，也不适用显式格式的欧拉模型

选择Error-Based需要设置误差，如果截断误差小于设置的误差，则时间步长的增大；如果截断误差大于设置的误差，则时间步长的减小。对于一般问题，Error-Based默认值为0.01是比较适合的。

c. Multiphase-Specific：对于多相流模型，自适应时间步长可以选择Multiphase-Specific，此时需要设置全局库朗数。这里不多介绍

注：对于CFL-Based和Error-Based两种方法，下面的设置都相同

3.2 Duration Specification Method

Duration Specification Method：设置计算总时间。对于Adaptive和User-Defined Function方式设置时间步长，此时无法直接确定总的计算时间，因此需要通过Duration Specification Method来设置。

共有四种设置方法：

a. Total Time：设置计算的总时间，包括之前的计算时间。比如目前保存的工况在之前已经计算了3s，设置的Total Time为4s，那么只会再计算1s就会停止计算。

b. Total Time Steps：设置计算的总时间步，这里和Total Time类似，包括了之前计算的时间步。

由于时间步长不确定，因此Total Time Steps无法确定计算的时间，只能通过Fluent控制台查看，或者输入文本控制命令查看。

也可通过Report Simulation Status查看。

单击Report Simulation Status，会弹出下框，包含当前的计算状态。

c. Incremental Time：增量时间，以此时作为计时起点计算的时间，无论当前的dat已经计算了多少时间，都会在此基础上再计算Incremental Time时间。

d. Incremental Time Steps：增量时间步，和Incremental Time类似，以当前的dat作为时间步起点，在此基础上再进行计算更多的时间步。

3.3 Parameters参数设置

a. Number of Fixed Time Steps：固定的时间步，在初始时间步长改变之前必须要执行的时间步数。比如，Number of Fixed Time Steps=100，初始时间步长为0.001s，则必须先使用初始时间步长0.001s计算100步后，时间步长才会改变。

b. Initial Time Step Size：初始时间步长，即第一个时间步长，必须以这个时间步长执行Number of Fixed Time Steps步后时间步长才会改变。

Initial Time Step Size必须介于最小和最大时间步长大小之间。一般应选择使库朗数=1的时间步长作为初始时间步长。

c. Max Iterations/Time Step：每个时间步最大的迭代次数。如果在此值之前就已经收敛，则进入下个时间步的计算，如果没有收敛，则达到最大迭代次数后也会进入下个时间步。

d. Reporting Interval：报告时间间隔，如设置为2，则表示每隔2个时间步在控制台输出一次文本，绘制一次图形。当计算速度比较慢时，可将此值设置为较大值，来减少文本输出频率，增加计算速度。

e. Profile Update Interval：DEFIEN_PROFILE函数的更新频率，设置为1，表示每个时间步，DEFIEN_PROFILE函数都会更新。

f. Time Step Size Update Interval：时间步长的更新间隔，如设置为2，表示每两个时间步才会自动计算并更新时间步长。

3.4 详细设置

点击Setting，会弹出Adaptive Time Stepping

a. Minimum/Maximum Time Step Size：限制的最小和最大时间步长，设置和计算过程中的时间步长必须在此范围内。

b. Minimum/Maximum Step Change Factor：限制时间步长的变化程度。这种限制可以使时间步长不会产生较大的突变，从而造成计算收敛性问题。

比如前一个时间步长为0.01s，但是计算出的当前时间步长为1s，两者差了100倍，这种时间步长的突然增大可能会造成计算问题，因此对其进行限制。

对于CFL-Based和 Multiphase-Specific方法，此处的值表示当前时间步长与前一个时间步长的比值。如上图的值，表示当前的时间步长被限制在前一个时间步长的0.5-2倍之间。

对于Error-Based的方法，此值为设置的error tolerance与计算出的截断误差之间的比值。

只发干货内容，如果您觉得文章内容有用，点赞、在看、收藏、转发是对我最大的鼓励！！！！！