FLUENT中实际气体性质模型(UDRGM)在大规模数值模拟中的应用.doc

1前言工程实际中遇到的各类流动问题，当工质热物理性质偏离理想气体时，需要应用实际气体，性质进行求解，才能比较精确地反映流动特性 例如在蒸汽透平低压部分，低压水蒸汽的，性质与理想气体性质偏离较大，如采用理想气体性质进行求解，误差较大IAPWS-IF97标准 中给出了计算水蒸汽气体性质的计算公式，但由于计算量很大，并不适合在大规模三维计算中直接采用因此有必要发展一种高效高精 度的实际气体，性质计算方法本文以水蒸汽实际气体的求解为例，利用FLUENT软件所提供的实际气体模型(UDRGM )实现了大规 模复杂计算中水蒸汽实际气体性质的应用2 FLUENT中的实际气体模型(UDRGM )及实现FLUENT软件中提供了实际气体性质计算的接口 UDRGM (User Defined Real Gas Model)利用该接口，可以将水蒸汽性质表导入 FLUENT中，实现对水蒸汽实际气体性质的计算UDRGM要求对实际气体提供如表1所示的12个函数对水蒸汽定义这12个函数 后，即可实现在FLUENT中调用水蒸汽实际气体性质进行计算表1 UDRGM所定义的函数1巒度p^p（T^p）7动力黏性2燔h = h（T. p）8热传导系3-s=s（T,p）9严)4-定压比熱q二訴Q10空T引&分子量11ch 丽dh恥)&12■ 严釧ioS 切Tr2.1实现方法首先由IAPWS-IF97计算公式生成表1中12个函数的数据表，在启动FLUENT时载入该数据表，在FLUENT中进行水蒸汽实际气体 ，性质计算时，直接利用该数据表进行插值得到所需的函数值。

由于插值输入变量的数目为2 ,因而插值可采用双线性插值或双三次插值 如图1所示，双线性插值根据插值点(xO,y0)周围4个节点上的函数值进行插值，而双三次插值根据插值点(x0,y0)周围16个节点上的 函数值进行插值因此双线性插值的精度低于双三次插值的精度，但节省计算时间按照文献［3］的估计，采用双线性插值比采用理想气 体状态方程进行计算的时间增加13%左右，采用双三次插值比采用理想气体状态方程进行计算的时间增加19%左右这种方法省去了 大量计算水蒸汽实际气体性质的时间，因而适用于大规模的三维CFD计算分析所生成的水蒸汽性质表输入变量等间距分布，如图1所示，可以得到插值公式为：io = int[(xo-xmin)/(xmax-xmin)] + 1 (1)jo = int[(yo-ymin)/(ymax-ymin)] + 1 (2)可以看到，插值所需要的时间与数据表中数据的密度无关因此，只要计算机内存条件允许，可以在数据表中生成足够密度的数据，从 而提高插值的精度2.2误差分析为了检验插值方法的精度，对插值误差进行了分析，以压力p和温度T为输入变量，对表1中函数1的插值精度进行了检验检验 中插值函数用了双线性插值。

结果如图2和图3所示，可以看到与直接由IF-97公式计算值相比，整个区域内插值计算的最大误差小于1%完全满足工程计算的要求图3插值对数误差图2由水蒸汽表插值和IF-97公式计算得到的密度2.3算例分析为了对比采用水蒸汽实际气体性质和理想气体，性质计算的差异，对某Laval喷管中的流动进行了计算计算采用H型网格，网格数为1 90x60计算中分别采用了理想气体状态方程、FLUENT内置的维里型状态方程以及本文发展的基于IAPWS-IF97水蒸汽性质表插值的 方法对该算例计算50步所用的时间对比结果列在表2中可以看到，与采用理想气体状态方程相比较，采用FLUENT内置的维里型状态方 程，计算时间增加了 180%,而采用本文发展的方法，计算时间增加了 27%表2计算50步所用计算时间对比状态方程计算时间相对计算时间理想气休状态方程131.00FLUENT内置维里型状态方程422.80IP97水蒸汽表插值-N k图4给出了三种方法计算得到的压力、温度、密度等值线分布，可以看到计算结果之间仍然存在明显差异采用维里型状态方程与采用 基于IF 97的水蒸汽表插值计算结果基本是吻合的，二者与采用理想气体状态方程得到的结果有明显差异。

考虑到计算所费的时间， 在大规模三维计算中，采用基于IAPWS IF 97的水蒸汽性质表插值比采用FLUENT内置的维里型状态方程可以节省大量计算时间440000 400000 360000 320000 280000 240000 200000 160000 丨2闻0 80000(b)温度(K )二753197531111on^n^oIT L(c)密度(kg/m3 )——理想气体状态方程——维里型状态方程……IF 97水蒸汽表插值图4参数等值线分布3结论基于FLUENT软件提供的UDRGM接口，发展了一种在大规模数值计算中采用实际气体，性质进行计算的方法，对该方法的精度和计算 时间进行了检验，表明具有高效高精度的特点以水蒸汽流动为例，利用所发展的方法，可以采用IAPWS IF 97水蒸汽性质数据高效 地进行三维复杂流动计算。