作用力求取方法对多相伪势模型的影响.pdfVIP

第 29卷 第6期 计 算 物 理 Vo1．29．No．6 2012年 11月 CHINESEJOURNALOFCOMPUTATIONALPHYSICS NOV．．20l2 文章编号：1001-246X(2012)06-0828-09 作用力求取方法对多相伪势模型的影响 胡安杰 ， 李隆键 ， 曾建邦 (1．重庆大学低 品位能源利用技术及系统教育部重点实验室，重庆 400030； 2．中国科学院可再生能源与天然气水合物重点实验室，广州 510640) 摘 要：基于伪势模型将有效密度和势函数结合计算粒子间相互作用力，提出一种改进描述气液相变的格子 Boltzmann模型．对不同状态方程控制的单组份流体相变过程进行模拟发现，与单纯 由有效密度或势函数计算粒子 间相互作用力的方法相 比，该模型的结果与Maxwell构建更吻合，且最大伪速度更小，扩大了温度的计算范围，提高 了模型的稳定性 ．基于 Laplace定律对界面现象模拟结果进行讨论 ，模拟结果与 Laplace定律吻合 ． 关键词：格子Boltzmann方法；粒子间相互作用力；相变；伪势模型 中图分类号 ：0359 ．1 文献标识码 ：A 0 引言 格子 Bohzmann方法 (LBM)是近二十几年发展起来的一种基于分子动力学的模拟方法 ，其具有清晰的 物理背景 在多相流领域 ，由于该方法在宏观上是离散方法 ，微观上是连续方法 ，因此 ，在模拟复杂流动问 题上具有常规方法所没有的优势，也为研究多相及相变问题提供了一种有效的途径．LBM在多相流领域的 研究主要有：①1991年Gunstensen等提出的着色模型 ；②1993年Shan等提出的伪势模型 ；③Swift等提 出的自由能模型 ；④近几年来提出的几种基于动理学理论的模型 ．其中伪势模型能方便地描述不同相 间的相互作用，自动追踪相界面，且计算程序相对简单，在近十几年来得到充分的发展，在多相流领域已经得 到较为广泛的应用 ．伪势模型中，相间相互作用来源于粒子间相互作用力．传统的粒子间相互作用力的 计算主要有两种 ：①通过有效密度梯度求取 ；②Zhang等人 提出的通过直接利用势函数梯度求取．然而 模拟发现 ，通过 以上两种作用力求取方式所得 的饱和密度 曲线与热力学平衡状态下，描写饱和密度 的 Maxwell方程得出的解存在很大的偏差．因此，为消除这一偏差，提高模型稳定性，本文采用有效密度和势函 数相结合的方式求取粒子间相互作用力，并引入多种具有代表性的实际气体状态方程对相变分离现象进行 模拟，对不同参数下模型的计算结果进行详细的讨论．结果表明与单纯通过有效密度或势函数计算粒子间相 互作用力的方法相 比，该模型结果与Maxwell构建更吻合，且最大伪速度相对更小，扩大了温度的计算范围， 提高了模型的稳定性．同时模型适用于所有状态方程 ，并满足 Laplace定律 ，为模型的应用奠定 了 良好 的 基础 ． 1 基本理论 1．1 格子 Boltzmann模型 伪势模型中单组份多相格子 Boltzmann演化方程可 以表示为 (+eaAt，t+△f)一 (，t)=一( (，t)一／(，t))／r， (1) 式中 为无量纲弛豫时间，由物质的粘性 决定，即 =(2r一1)cAx／6；c=Ax／At为格子速度 ， 为网格步 长，△ 为时间步长； 为各个离散方 向上的速度矢量， =0，1，2，…，Ⅳ代表离散速度矢量方 向． (，t)， (，t)分别是 t时刻位置 处的粒子分布函数和平衡态分布函数；对于D2Q9模型，其Maxwell平衡态分布 收稿 日期：2011—12—20；修 回 日期 ：2012—05—05 基金项 目：国家 自然科学基金资助项 目 作者简介：胡安杰 (1988一)，男，重庆，博士生，主要从事气液两相流的格子Boltzmann方法模拟研究，E-mail：anjie@equ．edu．cn 第 6期 胡安杰等：作用力求取方法对多相伪势模型的影响