两相流_中文.pdfVIP

w 页码，1/ 16(W) 流体流动模型 : 两相流 两相流 简介 以下例子分析多孔介质中的两相流. 描述流体如何不发生混合同时通过一个孔隙可以回答很多环境和工业问 题. 不幸的是, 多相分析很复杂需要求解一些未知量以及多个相关变量. 多相彼此之间的液体属性与每个液相 压力和饱和级有关. 这个问题验证两相流按照美国环境保护署的实验结果(参考文献1). 这个试验表明实验和数值计算的两相流结 果相符. 通过这些实验, 科研人员估计流体对(空气-水, 空气-石油, 和石油-水)变化的流动, 然后使实验结果 和计算机使用保持和渗透解析表达式数值计算的结果相匹配. 这个讨论根据Lincoln土壤和使用Mualem (参 考文献2)和van Genuchten(参考文献3)给出的流体属性解析公式展开. 图2-42: 在Hopmans和其他(参考文献1)两相流实验几何结构. 这个例子分多步. 第一部分建立水和空气的两相流模型; 方程求解压力. 饱和度随解而变化. 需要强调假设在 整个分析时间里土体中存在一些残余空气和水. 模型追踪空气前锋移动过湿液体通过观察饱和度而不是假设 离散界面. 第二部分修改空气－水属性对空气－石油和空气－水系统进行仿真. 模型定义 在空气和水的实验中, 插入空气到充满水和沙的实验柱体中. 进入的空气 (对该流体对为非湿相)强迫水(湿 相)向柱体基础出口流动. 在入口, 空气压力随时间逐步增加, 并且在柱体顶端没有水存在. 在向出口流动, 水 通过一个空气流动不可渗透的圆盘. 空气和水都不能通过柱体垂向壁面. 出口处的水压随时间变化, 对应着滴 定管中液体高度增加. 柱体总体高度8.34 cm, 半径6cm, 并且圆盘厚度0.74 cm. 实验时间170 hours. mk:@MSITStore:C:\Users\Administrator\Desktop\COMSOL_LibDoc_ES.CHM::/es_... 2010/11/25 w 页码，2/ 16(W) 以下讨论给出两相流仿真计算的步骤. 首先观察控制方程和定义流体保持和渗透的解析关系(参考文献2和参 考文献3). 然后给出一些详细的执行过程和模型数据表. 后面给出结果. 最后在图形用户界面逐步建立模型. 控制方程和边界条件 多孔介质两相流满足湿(w)和非湿(nw)流体的独立方程: (2-4) (2-5) 2 其中θ为总气隙或者饱和体积率; Se为有效饱和函数; t为时间(s); κ 为多孔介质固有渗透率(m ); k 为给 s int r 2 3 定液体的相对渗透函数; η为流体动力粘度(kg/(m·s)); p 为压力(kg/(m·s )); ρ为流体密度(kg/m ); g为 重力加速度; 和D为垂向坐标(例如, x , y , 或者z )(m). 如果流体分布为连续的, 流体在整个分析时间内未完全填充土壤, 给定湿相容积率, θ , 和非湿项, θ . 对湿 w nw 相, θ从0或者很小的残余值θ到总体孔隙率, θ . 有效饱和度, Se, 源自比例系数θ关于θ和θ从0到1变化.