PISO算法详解.pdfVIP

CFD 模型理论 对于不可压流体的控制方程主要包括连续性方程及动量守恒方程，对于CFD-DEM 耦合模型中的液相，其连续性方程及动量守恒方程如下： 连续性方程：     u 0 (1) t  f   f f  动量守恒方程：   P   F  u   u u  f  g  f f  f  f f   f f f  f (2) t f f f 在 CFD-DEM 模型中我们采用下标f 来表示连续相，下标p 表示离散相。根据 Boussinesq 假设，应力项τ 包括基于分子粘性的平均应力τ 及湍流应力τ 两部分，不可 f m t 压缩流体消掉密度，其计算公式变换为如下所示： 将雷诺应力分解为两个部分：  T  2 偏分量： u  u   u I (46) f ,t  f f  3 f  f  基于分子粘性的平均应力：  T  2   u  u   u I (47) m f  f f  3 f  f  2 各向同性分量： f I (3) 3 将式4-48 归入压力项，令压力： P 2 P   (4) f 3 f 则不可压缩流动的有效应力表示为：   u  T u  2 eff eff    u I (5)  f f  3 efff f 其中有效粘度定义：   (6) eff f f ,t 动量方程中相间动量交换定义如下： Ff  k pf u f up  (7) 相间动量交换液相部分做隐式处理，动量方程根据显隐式划分变换为： 对于液相