电场作用下两相流的计算方法研究

发布时间：2020-03-21 21:52

【摘要】：以气泡为分散相的气液多介质流体流动广泛存在于能源、石油、动力、化工等实际工程领域。施加外部电场作用提供了一种非接触式主动控制气泡动力学行为的有效途径。然而,在电场作用下气泡的运动行为的模拟研究中存在计算量大,计算效率低、精度不高等问题,所以研究高效率的计算方法来求解电场作用下两相流模型具有非常重要的意义。本研究基于不同计算方法,以气液两相流为研究对象,将不可压缩流体的连续性方程、Navier-Stokes方程、电场方程等控制方程耦合建立了一个研究粘性流体中气泡的动力学模型,并进行数值模拟计算,研究了耦合电场的两相流行为特征,分析电场作用对气泡运动的影响。论文的主要研究内容和研究结论如下:1.采用水平集(Level-set)和相场法(Phase-field)来跟踪记录气泡上升过程中可能出现的形变、破裂等变化,通过质量守恒验证了两种计算方法的准确性和有效性,同时为了计算的精确度对计算区域进行了一定的网格划分,并使用网格独立性检验来计算得到了最优的网格大小。2.相较于水平集方法,相场法在计算效率和计算精度上都有一定提升。其中,在相同网格下相场计算方法的计算效率最大提升了5倍,在垂直电场环境下的计算效率提升了3倍,而且相场法更容易捕获气泡上升过程中气泡发生的细微变化,模拟计算结果的质量也更好,验证了相场计算方法的有效性。3.通过两种计算方法对耦合电场下气泡上升过程的模拟分析,结果表明气泡在静电力、浮力和表面张力的相互作用下,沿电场线方向拉伸成椭球形,其长短轴之比与外加电场强度成正比。气泡上升速度受电场作用影响,垂直电场加速气泡的上升,而水平方向减速。4.利用相场计算方法模拟了计算外部均匀竖直电场环境下一对共轴气泡的运动、融合现象,计算不同电场强度对上升气泡相互作用和相互融合的影响。研究结果表明,均匀强电场加速了流体中双气泡的碰撞融合过程,提高了气泡的上升速度。当电场强度较大时,融合气泡后气泡会发生破裂,破裂时间也随电场强度的增加而减少。将相场方法与同条件下LBM(Lattice Boltzmann method)方法的结果进行了对比,两种方法所得气泡形态变化吻合良好,验证了本文相场计算模型的正确性。
【图文】：

离散模型,格子

图 1.1 D3Q19 LBM 离散模型的计算：0 0 ( 1, 0, 0) , (0 1 0) (0 0 1) ( 1, 1, 0) , ( 1 0 1) (0 1 1) c c c c c c （，），，，，，，，，，，3c1 3 01 18 1 61 36 7 18iiw ii [55]格子 Boltzmann 方法(LBM)究了气。孟辉等[56]利用格子 Boltzmann 法计位置的两个气泡的运动情况，分析两影响。史冬岩等[57]改进了 Swift 的 LB比(1:1000)的气液两相流问题，通过与数值结果精度并验证了其质量守恒。

示意图,网格划分,示意图

同时对界面处经过局部加密网格，得到网格划分图如图 2.3 所示。图 2.3 网格划分示意图2.3.2 数值模拟计算步骤计算流体力学 Computational Fluid Dynamics(CFD)，是为了弥补理论分析和实验研究分析方法的不足而发展起来的，主要是通过计算机和数值方法来求解流体力学流体流动的相关控制方程，从而来描述多介质流体流动等过程中的物理特性。本文基于以上控制方程，将方程进行耦合求解，对自由界面的追踪模拟分别采用水平集方法和相场方法来实现，，针对在外加电场作用下气泡在黏性流体中的上升的模拟，在对流体问题进行求解时，具体计算步骤如下所示：1、分析物理问题，阅读相关参考文献，确定物理模型所使用的控制方程2、创建几何模型，定义物理参数3、设定边界条件和初始值4、生成网格:将求解域分解成单元，生成计算节点，调整网格选取合适的网格大小。5、求解(设定 Solver 参数)离散方程：求解组装得到的大型稀疏矩阵，这是最消耗内存和时间的步骤。6、后处理，结果可视化、动画、输出，后处理的目的是有效或直观的观察和分析计算结果。为更好的处理模拟结果，如气液形状、气泡体积、压力等，可通
【学位授予单位】：兰州理工大学
【学位级别】：硕士
【学位授予年份】：2019
【分类号】：O359

【参考文献】