格子Boltzmann方法在模拟液滴冲击液面中的应用 Application of Lattice Boltzmann method to simulation of impact of droplet on liquid surface.pdfVIP

水 利 学 报 2008年12月 SHUILIXUEBAO 第39卷第12期 文章编号：0559．9350(2008)12．1316．05 格子Boltzmann方法在模拟液滴冲击液面中的应用 张 华1，郝智文2，冯志军3 (1．华北电力大学可再生能源学院暨电站设备状态监测与控制教育部重点实验室，北京102206； 2．山西省气象科学研究所，山西太原030002；2．新华水利水电投资公司，北京100053) 摘要：应用格子Boltzmann两相流方程，结合水平集界面重构技术，提出两相流的数值模型。采用五阶WENO格式和 三阶TVD Runge-Kutta格式来求解水平集函数的输运方程。在不计重力的条件下，对液滴平动进行了数值模拟，在 计算了800个时间步后，平动的液滴仍然保持原来的形状。在考虑重力的条件下，对液滴在空气中自由下落、液滴 冲击液面的全过程进行了数值模拟，精细地刻画了液滴和液面界面的变化过程。数值模拟结果与物理现象一致。 关键词：液滴；液面；格子Boltzmann；水平集；数值模拟 中图分类号：IVl35．2 文献标识码：A 液滴冲击液面是一种常见的自然现象，这种现象在工程中也常常出现，如水利工程中的泄洪雾化现 象。液滴冲击液面可形成空腔，引起水面波动、空气夹带等界面迁移现象，它是一个复杂的液气两相流 边界元方法对液滴冲击同种液体的无限大液面现象做了数值模拟研究。本文试采用格子Boltzmann方 法b】，结合水平集界面重构技术，对液滴平动，以及液滴在空气中自由下落、液滴冲击液面的全过程进行 仿真研究。 1格子Bohzmann模型 述，即LBGK(Lattice Bhatnagar-Gross．Krook)方程： (1) 式中，Z为流体A粒子分布函数与流体B粒子分布函数之和，Ci为粒子速度矢量，F为平衡分布函数， C；II2，t：为系数，可表示为： r为松弛时间，F为外力，P=II t：=4／3，P=o；f；=1／3，P=1；tp。=1／12，P=2(2) 流场离散采用二维9点正方形格子(D2Q9)，其9个方向的粒子速度矢量Ci可表示为b】： c。=(o，1) 江o f i=1—4 {C。=(cosol，sinoi)，0j-(i一1)hi2 (3) L i=5—8 cj=√2cosOj，sinOi)，见=(i一5)，r，2+丌，4 密度P和动量J『可表示为： 收稿日期：2008-03-08 基金项目：国家自然科学基金委员会、二滩水电开发有限责任公司雅砻江水电开发联合研究基金项目 作者简介：张华(1962一)，男，山西平遥人，博士，教授，主要从事水工水力学等研究。E-mail：zhanghua@necpu．edu．en —1316一 万方数据 p(r，t)=∑Z(，，￡) (4) ‘；0 S ， ． J(川)=J+告F，J=∑z(，，￡)C； (5) 利用渐近理论，通过对式(1)进行Chapman．Enskog多