ITSM系统配置详细手册
对于高瑞利数流动(Ra > 10^8),为了得到最佳结果,你应该遵循下面所介绍的某一处理程序:
第一个程序使用定常状态方法:
-
开始解决时使用较低的瑞利数(如:10^7),然后使用一阶格式运行直到收敛。
-
改变有效瑞利数,改变重力加速度的数值(如:从9.8改为0.098来使瑞利数减少两个量级)。
-
使用上面的结果作为高瑞利数流动的初始猜测,然后用一阶格式开始高瑞利数流动的计算。
-
用一阶格式获得解之后,你可以采用高阶格式继续计算。有关更多信息,请参阅一维非定常热传导方程。
第二个程序使用时间相关方法来获取定常解:
-
使用相同或较低瑞利数时得到的定常状态解开始计算。
-
估计时间常数为:TLg L Ra L U L ∆ == − βα τ 21 2 Pr~其中L和U分别是长度和速度。使用时间步长D t:4 τ =∆t如果使用更大的时间步长D t可能会导致发散。
-
当频率f t = 0.05--0.09的振动衰减之后,就达到了定常状态。注意,t是方程1中估计的时间常数,f是单位为Hz的振动频率。一般说来,要达到定常状态一般要进行5000个时间步。关于如何实现,请参考一维非定常热传导方程FTCS格式Fortran程序。
注意:除非使用Boussinesq近似,否则非定常方法不能用于封闭区域。它总是用于具有入口和出口的区域。浮力驱动流动的后处理报告和其它热传导计算的报告一样。详情请参阅热传导的报告一节。
周期性流动和热传导周期流是指流动和热的解具有周期性重复的特点。周期性流动分两种:一种是在周期性平面内没有压降的周期流;第二种是流向周期流。本解讨论流向周期流以及周期性热传导,关于没有压降的周期流请参阅周期性边界条件一节。