首页/文章/ 详情

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

8月前浏览6774

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与计算出的截断误差之间的比值。

来源:Fluent学习笔记
ACTFluentSTEPS多相流UDFUM控制
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-02-28
最近编辑:8月前
Fluent学习笔记
博士 签名征集中
获赞 124粉丝 327文章 133课程 3
点赞
收藏
未登录
1条评论
豌@豆
签名征集中
4月前

回复
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈