LBM相变传热与流体流动数值分析13.pptVIP

相变传热与流体流动数值分析(第11-14讲)

;主要内容;根据不同的划分规那么，LBM的模型包括：;;;13.1LBGK模型;;;;;;;;;压力入口压力出口;多松弛模型和单松弛模型的主要区别在于它的碰撞过程包含多个松弛时间〔Multiple-Relaxation-Time〕。;上式描述了离散速度分布函数

在速度空间的时空演化过程。;;;;;颜色模型〔Chromodynamicsmodel〕;Gunstensen等人于1991年提出的第一个多组分多相模型。

基于Rothman和keller提出的LGA两相流模型。

Grunau等又将此模型推广到密度和粘性变化的两相流系统中。;以两相流为例：

引入两个分布函数，分别表示红色相和蓝色相流体。混和流体的分布函数的演化方程为;c;确定了后，混和分布函数在每一个时间步内根据方程〔〕进行演化，红蓝两相的分布函数那么通过重新标色〔Recoloring〕过程来获得。这一过程是求解下面的极大化问题;综上，Gunstensen等人的颜色模型可分为三步进行演化：

1.单相碰撞：a)计算混和分布函数；

b)计算和平衡态分布函数；

c)按照单相碰撞规则计算碰撞后分布函数。

2.两相碰撞：a)将表面张力扰动加入碰撞后分布函，；

b)对重新标色，得到。

3.迁移：

;其中分别为每一相的碰撞算子，采用BGK模型;仍然是外表张力引起的扰动，形式为;局限性：;1993年和1994年，Shan和Chen提出了一种能直接地刻画粒子之间相互作用的LBE模型。;SCMP模型;‘a’termduetoattractiveforcesbetweenmolecules[atmL2mol-2]

‘b’termduetofinitevolumeofmolecules[Lmol-1];CO2:P-VSpace;Liquid-VaporCoexistence:CO2,P-VSpace;IncludingDirectionalH-Bonds:Water,298K:P-rSpace;;Gistheinteractionstrength

Yistheinteractionpotential:

Y0andr0arearbitraryconstants

G0forattractionbetweenparticles.

Forceisstrongerwhenthedensityishigher

;;Non-idealComponent;;伪势模型;Martys等提出;伪势模型;伪势模型;;;其中，k表示不同的组分。;参数的绝对值决定了组分之间相互作用的强度，

符号决定了两者之??是相互吸引还是相互排斥，即负值是吸引力而正值为排斥力。;同样的，通过改变平衡态速度来实现的影响，即;在Shan-Doolen模型中，混合流体的速度为;参数的绝对值决定了组分和壁面之间相互作用的强度。

符号决定了作用力是吸附还是排斥，即负值是吸引力，表明壁面对流体时润湿的，而正值为排斥力，表面壁面对流体而言是非润湿壁面。

s为壁面指数，当格子特性为固体时，其值为1，当格子特性为流体时，其值为0。;同样地，在D2Q9中;;伪势模型;;;;LaplaceLaw;;66×132个格子；

半步长反弹格式、周期边界；

恒定的体积力；

密度均为1；

运动黏度均为0.16；;;颜色模型和伪势模型都是基于界面现象的唯象模型，含有一些人为的假设。1995年，Swift等人直接从多相/多组分流体的自由能理论出发，构造了与热动力学理论一致的多相和多组分流体模型。;根据vanderWaals理论，对于包含相界面的单组分非理想流体系统，当处于局部平衡时，其自由能泛函为;根据自由能泛函定义一个非局部压力;为构造符合上述热力学理论的模型，Swift等建议在LBGK模型中使用一个包含密度