基于FLUENT的大型储油罐内原油非稳态传热的数值模拟.doc

基于FLUENT的大型储油罐内原油非稳态传热的数值模拟摘要：本文针对大型浮顶储油罐内原油的温降问题，应用计算流体软件FLUENT对储油罐内原油的温降过程进行了数值模拟，得到了在自然冷却条件下，原油在浮顶储油罐内温降过程中温度场和速度场的分布图。通过数据对比，得到计算流体软件计算的数值与实际测试的数值吻合较好，相对误差在3%以内。关键词：储油罐 温降 FLUENT 数值模拟一、物理模型的建立本文利用FLUENT软件对大型储油罐的温降进行了模拟计算。图1为储油罐的二位几何模型，对此物理模型的相关假设如下：认为罐内横向温度温度分布均匀，所以只对纵向二维温度场进行计算；对计算区域进行简化，取过油罐中心轴的径向切面为研究对象，变三维问题为二维问题； 忽略太阳与油罐间的虑辐射换热；忽略原油内部物理化学因素而产生的内热源，忽略油罐边缘以及底部原油和石蜡的凝固潜热；假设外界环境温度近似为随时间变化的周期性函数；将储油罐浮顶近似处理为单层钢板，同时也将罐壁、罐底假设为均匀的等厚钢板。其中各符号的具体意义为：To为稳罐后储油罐内原油的初始平均温度，为常数；Td为大地恒温层处土壤的温度，为常数；T1为大地温度场的初始温度分布；fw为储油罐壁（保温层外侧）环境的温度函数；ft为储油罐顶部环境温度的分布函数；H1为浮顶储油罐内原油液位的高度（10.1米）；H2为储油罐罐底钢板到大地恒温层的深度（=10米）；R为储油罐径向平面半个平面的宽度；1为保温层厚度，保温层平均厚度（1=80毫米）；2为储油罐罐壁平均厚度（2=10毫米）；3为储油罐顶板厚度（3=5毫米）；4储油罐底板的平均厚度（4=10毫米）；二、数学模型控制方程的直角坐标张量符号形式如下所示质量守恒方程： （1）动量方程： （2）其中p为静压，ij为应力张量，gi和Fi分别为i方向上的重力体积力和外部体积力，式中： （3）能量方程： （4）其中E为内能，keff是有效热传导系数，Jj是组分j的扩散流量。而且hj为同分子物质的焓值，Sh包括了化学反应热以及其它用户定义的体积热源项。式中：，mj是组分j的质量分数。湍动能k方程为： （5）湍动能耗散方程为： （6）； ； ； 在标准k-模型中，C1、C2和C3为经验常数；k和分别是与湍动能k和耗散率对应的Prandtl数；根据Launder等的推荐以及后来的验证，有下式：C1=1.44， C2=1.92，C=0.09，k=1.0，=1.3。三、边界条件1、储油罐顶部与大气直接接触，它们之间的热交换主要形式是对流换热。公式表示为：其中：为钢的导热系数，W/(m)； Tw钢表面温度，；ft大气温度函数，； 边界；h大气与钢板间的对流换热系数，W/(m2)。2、储油罐底部基座为大地，取大地温度恒温层位置=10米，近似认为基座左边界上为绝热。应用公式表示为： 其中：q为热流密度q=0，W/m2；大地的导热系数，W/(m)；边界。3、储油罐基座底边，由罐底向下10米深处为（=10米）恒温。用公式表示为： Td储油罐底向下10米处至无穷远处温度，为常数*为2.2，边界。4、同理可知罐壁为第三类边界条件，对称轴边界上为绝热边界条件。四、结论1、高温区并非集中在油罐中心，而是相对在靠近罐壁附近的原油温度最高。2、在自然对流过程中，原油温降的流动分为两种形式，一是沿着原油上液面对称轴罐底侧壁面这一路进行，这是主要形式；二是在靠近罐壁处原油沿着罐壁向下形相对一形式逆时针形态的涡流。3、对比分析现场实测数据与应用FLUENT模拟计算的结果数据，实验值与数值计算值的相对误差在3%以内，符合工程上的要求。参考文献：1 王明吉、张勇、曹文：立式浮顶原油储罐温降规律J油田地面工程， 2005,24(12)：11122 大庆油