热油管道FLUENT模拟

石油工程应用实例7.1 热油管路温降分析7.1.1 案例简介与思路分析我国所产原油80％以上为凝点较高的含蜡原油。

目前含蜡原油多采用加热管输的方式。

在热油管道运行过程中，不可避免地会因计划性或事故性的原因而发生停输。

原油管道的停输再启动，一直是管道安全运行的核心技术问题。

由于热油管道大部分都采用埋地敷设，因此在实际生产、运行中对于停输后热油管道以及土壤温度场分布的研究，具有重要的意义和作用。

【案例简介】本例的输油管为Φ500×6mm钢管，管中心埋深为1.5m，管外保温层厚度10mm，沥青防腐层厚度6mm。

由于涉及到与土壤环境的传热问题，需要考虑大地恒温层及管道热力影响区的水平边界来划定土壤区域范围，本例中选择土壤垂深10m，水平方向的热力半径为10m，如图7-126所示。

图7-1 管路温降几何示意图【思路分析】在模拟的过程中分别划分出流体区域与固体区域，然后在选择求解器时采用非耦合式隐式算法，打开能量方程，由于存在温度梯度，因此需要编写UDF来定义边界温度。

整个模拟计算的流程如图7-127所示。

图7-2 井下节流阀流动模拟流程图【光盘文件】——参见附带光盘中的“START\Ch7\7-6.msh”文件。

——参见附带光盘中的“END\Ch7\7-6.cas，6-6.dat”文件。

——参见附带光盘中的“AVI\Ch7\7-6.avi”文件。

[1]启动GAMBIT，单击→→，依次创建坐标（0，0，0）、(0.25，0，0)、（0.256，0，0）、（0.266，0，0）、（0.272，0，0）的5个点，如图7-128所示。

图7-3 创建的点[2]单击→→，打开Move/Copy Vertices对话框，选中之前创建的5个点，选择Copy选项，在Operation选项卡中选择Rotate，在旋转角度Angle中输入90，旋转轴Axis 选项保持默认，即围绕z轴旋转，单击按钮即可。

在旋转角度Angle中输入-90，其他选项保持不变，单击按钮，得到的点如图7-129所示。

图7-4 创建的点图7-5 创建的圆弧[3]单击→→，打开Create Straight Edge面板，创建点1（Vertix.1）与点2（Vertix.2）、点1与点6、点6与点7、点7与点8、点8与点9、点1与点10、点10与点11、点11与点12、点12与点13，共9条线段；右键单击按钮，在下拉选项中选择按钮，打开Create Circular Arc面板，在Center选项右侧选中点1，End-Points选项右侧选中点2、点6，保持Arc选项为默认，单击按钮即可得到所需弧线。

同理，可得以点1为Center，点3与点7、点4与点8、点5与点9、点2与点10、点3与点11、点4与点12、点5与点13所形成的圆弧，如图7-130所示。

[4]将线段2、线段9、线段10、线段11、线段12组成一个面，线段1、线段7、线段8、线段9、线段10和线段3、线段4、线段5、线段11、线段12分别组成两个面，单击→→，在Create Face from Wireframe面板的Edge黄色输入栏中选取需要围成面的线段，单击按钮生成几何面。

本例生成的三个面如图7-131所示。

图7-6 创建的点图7-7 移动后的面[5]单击→→，打开Move/Copy Faces对话框，选中之前创建的5个面，选择Move选项，在Operation选项卡中选择Translate，在坐标Global 中y方向右侧输入-1.5，其他选项保持默认，单击按钮即可，如图7-132所示。

[6]单击→→，依次创建坐标（0，0，0)、(10，0，0)、（10，-10，0）、（0，-10，0）的4个点；单击→→，打开Create Straight Edge面板，创建点5（Vertix.5）与点14（Vertix.14）、点14与点15、点15与点16、点16与点17、点17与点13，共5条线段；单击→→，在Create Face from Wireframe面板的Edge黄色输入栏中选取需要围成面的线段，依次分别为线段14、线段19、线段20、线段21、线段22、线段23、线段18，单击按钮生成几何面，如图7-133所示。

图7-8 绘制好的土壤区域7.1.3 网格划分[1]单击→→，在Mesh Edges黄色输入栏中选中线段1、线段2、线段6、线段11，线段12，线段13，线段14，线段15，线段16，线段17，线段18，设置Ratio的节点距离比例为1，选择Interval Count分段方式，在左侧输入60，单击按钮生成该线段的线网格；选中线段3、线段5、线段9、线段7，设置Ratio的节点距离比例为1，选择Interval Count分段方式，在左侧输入3，单击按钮生成该线段的线网格；选中线段4、线段8，设置Ratio的节点距离比例为1，选择Interval Count分段方式，在左侧输入4，单击按钮生成该线段的线网格；，划分好的线网格如图7-134所示。

图7-9 创建的线网格图7-10 管道区域的面网格[2]单击→→，打开Mesh Faces面板，选中已经完成线网格划分的面1、面2，运用Quad单元，Tri Primitive方法对这2个面进行划分，其他设置保持默认即可；选中已经完成线网格划分的面3、面4、面5，运用Quad单元，Map方法对这3个面进行划分，其他设置保持默认，其生成的面网格如图7-135所示。

[3]单击→→，在Mesh Edges黄色输入栏中选中线段20、线段22，设置Ratio的节点距离比例为1，选择Interval Count分段方式，在左侧输入75，单击按钮生成土壤区域上下两条边的线网格；选中线段21，设置Ratio的节点距离比例为1，选择Interval Count分段方式，在左侧输入80，单击按钮生成土壤区域右侧线段的线网格；选中线段19，线段19的红色箭头方向由上指向下，设置Ratio的节点距离比例为0.95，选择Interval Count分段方式，在左侧输入40，单击按钮生成土壤区域轴线处管道上方线段的线网格；选中线段23，线段23的红色箭头方向由下指向上，设置Ratio的节点距离比例为0.95，选择Interval Count分段方式，在左侧输入75，单击按钮生成土壤区域轴线处管道下方线段的线网格，如图7-136所示。

图7-11 土壤区域的线网格[4]单击→→，打开Mesh Faces面板，选中已经完成线网格划分的面6，运用Tri单元，Pave方法对这1个面进行划分，其他设置保持默认即可，其生成的面网格如图7-137所示。

图7-12 热油管路传热问题的网格[5]单击→，在Specify Boundary Types面板中将地面（edge.20）定义为壁面（Wall），名称为di-biao；将底边（edge.22）定义为壁面（Wall），名称为di-bian；将右侧边（edge.21）定义为壁面（Wall），名称为right；将防腐层（edge.14、edge.18）定义为壁面（Wall），名称为fangfuceng；将保温层（edge.13、edge.17）定义为壁面（Wall），名称为baowenceng；将钢管外壁（edge.12、edge.16）定义为壁面（Wall），名称为gangguan-waibi；将钢管内壁（edge.11、edge.15）定义为壁面（Wall），名称为gangguan-neibi；将对称轴（edge.2、edge.3、edge.4、edge.5、edge.6、edge.7、edge.8、edge.9、edge.19、edge.23）定义为壁面（Wall），名称为duichengzhou；，定义好的边界如图7-138所示。

图7-13 定义好的边界类型名称列表[6]单击→，在Specify Continuum Types面板中Name中输入yuanyou，在Type 选项中选择FLUENT，在Entity中选择面1、面2，单击按钮，完成对这一区域的定义；在Name中输入gangguan，在Type选项中选择SOLID，在Entity中选择面3，单击按钮；和定义gangguan方法一样，依次定义new_baowenceng（面4）、new_fangfuceng（面5）、turang（面6），定义好的介质类型如图7-139所示。

图7-14 定义好的介质类型[7]执行File→Save As…命令，将工程文件保存为“7-6”。

同时，执行Export→Mesh命令，输出网格文件7-6.msh。

7.1.4 Fluent中的求解设置[1]由于沿土壤方向存在温度梯度，故需要定义y方向上的温度变化规律。

本例的程序如下：#include"udf.h"DEFINE_PROFILE(WALL_youbian,thread,index){face_t f;real x[ND_ND];real y;begin_f_loop(f,thread){F_CENTROID(x,f,thread);y=x[1];F_PROFILE(f,thread,index)=273.15+2.2-0.03*y;}end_f_loop(f,thread)}程序编写好后，保存文件名为d.c，并放至工作目录下。

[2]双击桌面Fluent 14.5图标，在弹出的Fluent Launcher窗口中选择2D计算器，单击OK启动Fluent。

[3]执行File→Read→Mesh…命令，读入划分好的网格文件“7-6.msh”。

网格读入之后，在左侧控制树中点击Solution Setup下的General，于中间General控制面板中点击按钮检查网格的信息，待文本窗口中最后一行出现“Done”语句。

[4]点击General控制面板中的按钮，弹出Scale Mesh面板，本例不需要重新定义网格的尺寸。

[5]在General操作面板的Solver下方勾选Gravity，并在Y(m/s2)中输入-9.8，即考虑在Y轴负向存在重力加速度，见图7-140。

图7-15 设置重力加速度[6]点击操作面板Problem Setup的Models，双击Energy，弹出图7-141所示的Energy 对话框，选中Energy Equation，单击OK。

回到Models面板，双击Viscous-Laminar，弹出Viscous Model对话框（图7-142），所有选项保持默认，单击OK。

图7-16 选择能量方程图7-17 Model控制面板[7]由于材料库中没有原油这一材料，需要用户自己定义。

单击Create/Edit Materials操作面板右侧的按钮，弹出的图7-143所示的Open Database操作面板。