精确二维固相体积分率方案(2).pdfVIP

精确二维固相体积分率方案 为了精确计算二维固相体积分数，将颗粒与网格的部分重叠进一步细分为 四种情形。图 3.1 中给出了各种重叠情形下颗粒与网格的相对位置关系示意。 计算网格的二维固相体积分数时，部分重叠的四种情形还需要进一步区分 网格顶点与颗粒的位置关系。 本章模拟采用正方形网格进行空间离散， 这里令颗 粒的坐标为 ( X ,Y ) ，网格宽度和厚度分别为 dx, dy ，网格的几何中心为 ( x, y) ，网 格的四个顶点分别为 ( x dx / 2, y dy / 2) ，(x dx / 2, y dy / 2) ，(x dx / 2, y dy / 2 ) 和 (x dx / 2, y dy / 2) ，由此可计算颗粒与四个顶点的距离 d , d , d 和 d 。这些距 1 2 3 4 离中，当 d ,d , d 和 d 分别为最小时，单颗粒所占网格的二维固相体积分数的表 1 2 3 4 达形式存在差别， 不失一般性，这里只给出当 d d ，d d 和 d d 同时成立时， 1 2 1 3 1 4 各种重叠情形下单颗粒所占网格的二维固相体积分数的表达式。 图 3.1 颗粒与控制体的五种重叠情形 图 3.2 颗粒与网格重叠时颗粒质心的区域划分 图 3.2 给出了当 d1 d2 ， d1 d3 和 d1 d4 同时成立时，各种重叠情形下颗粒 质心与网格左下顶点的相对位置关系， 与部分重叠的四种情形相对应， 这里部分 重叠时颗粒与顶点的位置关系也划分为四种情形。 图 3.2 还显示了颗粒与网格发 生重叠时， 颗粒质心所处位置的精确范围。 以下将分类讨论各种重叠情形下单个 颗粒对网格的二维固相体积分数的贡献。 r (X ,Y) r 1 1 (x ,y ) (a) Cas 1 (b) Cas 2 图 3.3 颗粒与网格左下顶点最近的重叠情形（ Cas 1 与 Cas 2） 图 3.3 给出了颗粒与网格左下顶点距离最近的重叠情形 1 和 2。这两种情形 下，颗粒所占网格的二维固相体积分数的表达形式是唯一确定的，不存在变例。 对于情形 1，颗粒质心与网格的位置关系为： x1 r X x 且 y1 r Y y 。则单 个满足重叠情形 1 的颗粒 i 所占目标网格的二维固相体积分数为 1 2 f i r (3.3)