基于fluent的海洋溢油数值模拟研究.pdfVIP

基于Fluent的海洋溢油数值模拟研究 吕国录刘志萌 (集美大学轮机工程学院福建厦门 361021) 摘要 针对海洋溢油扩散问题．基于Fluent模拟软件，建立二维数值水槽．分析比较不同海流速度和不同含气率条件下油 气混合物的输移扩散特性。计算结果袁明，海流速度越大，油滴尺寸越小，油气混合物的扩散面积相对较大：含气率越大，溢油上 升的速度越大，油气混合物的扩散面积就相对较小。此数值模拟研究可为海上溢油回收工作及环境污染防治提供参考依据。 关键词 海洋溢油 污染防治 海流速度含气率Fluent 中图分类号：X55 文献标识码：A 文章编号：1672—9064(2016)06—030—03 在大力开发海上油气田的过程中．海底管道在输送油气 (2)动量方程。动量守恒方程，也称为N—S方程，表达式 方面发挥了独一无二的作用．但海底管道泄漏的风险不可忽 视。水下油气泄漏对海洋渔业、水产养殖业及海洋环境将造 成致命的伤害．并会严重破坏海洋生态平衡．影响海洋生物 p面du妒h+警+3计pyx+警 生存．因此对海底管道泄漏问题进行一定的研究．将为海洋 溢油的预防及溢油应急处理提供一定的参考价值。从20世 p宏=峨+警+鲁+警 纪70年代起．国外已经开始对海洋溢油问题进行了系统地 p面dw矽。+警+警+0出Pzz 研究．挪威Johansen等建立了DEEPFLOW水下溢油模型…． 并利用实验对模型进行了验证；Yapa和Zheng建立了一个式中：民、R、凡分别为单位质量流体上的质量力在三个方 比较完善的三维数学模型CDOG来模拟水下溢油形成的浮 向的分量；＆为流体内应力张量的分量。 射流【2、。近20年来．国内学者也积极开展了相关的研究．汪 守东在Yapa和Zheng的溢油模型基础上建立海底溢油的数 学模型．并首次考虑了溢油的乳化作用㈨：王晶提出了预测 资料．本文采用标准的K一8模型，其中湍动能K和湍流耗散 海底管道由于焊接裂纹或腐蚀所形成的小孔发生溢油的模 率￡有如下形式[7]： 型．对溢油在水中输移扩散的过程及影响其扩散的主要因素 进行了研究．并分析了各种因素对油滴尺寸的影响[4]。 p警=毒[(p+盯It,。、l肚Okj】+q+Ga一舻一‰ 本文主要利用Fluent软件．采用基于压力求解的PISO 算法，选择VOF模型、标准k-s湍流模型来解决瞬态条件下 压力速度的耦合问题．对较深水域的海洋输油管道泄漏问题 1．2 PISO算法 进行数值模拟。主要研究了不同海流速度、不同含气率以及在 with of PISO(Pressure-ImplicitSplittingOperators)是压 相同含气率不同海流速度时油气混合物在水中的输移特性． 力隐式分割算法．该算法是基于压力速度耦合关系的一种算 为海洋溢油回收工作及海洋环境污染防治提供参考依据。 法。PISO算法的主要思路就是将修正方程中求解阶段中的 1 数学模型 为了便于阐述海洋溢油过程的机理．本文建立二维求解 算法．校正的速度会更加接近流体力学中的连续性和动量方 模型．其基本方程包括非定常不可压黏性流体流动的连续性 程，因此PISO算法在求解瞬态问题中具有明显的优势[8]。 方程、动量方程及湍流方程等基本控制方程，采用PISO算法 2数值模拟计算结果与分析 解决压力与速度耦合的问题；其离散格式均采用默认设置， 油一水一气三相流动采用VOF模型进行数值模拟计算[5]。 底部距离左侧1m处设置一直径为2cm的小孔作为泄漏孔． 1．1基本控制方程 选用四边形结构化网格；并进行边界条件设定：泄漏孔为速 (1)连续性方程。按照质量守恒定律，在控制体内流人的 度入口，速度大小为3m／s，方向垂直向上，左、右边界设置为 质量与流出的质量之差．应该等于控制体内部流体质量的增 速度人口，方向垂直边界