精确二维固相体积分率方案

精确二维固相体积分率方案为了精确计算二维固相体积分数，将颗粒与网格的部分重叠进一步细分为四种情形。图3.1中给出了各种重叠情形下颗粒与网格的相对位置关系示意。计算网格的二维固相体积分数时，部分重叠的四种情形还需要进一步区分网格顶点与颗粒的位置关系。本章模拟采用正方形网格进行空间离散，这里令颗粒的坐标为),(YX，网格宽度和厚度分别为dydx,，网格的几何中心为),(yx，网格的四个顶点分别为)2/,2/(dyydxx，)2/,2/(dyydxx，)2/,2/(dyydxx和)2/,2/(dyydxx，由此可计算颗粒与四个顶点的距离321,,ddd和4d。这些距离中，当321,,ddd和4d分别为最小时，单颗粒所占网格的二维固相体积分数的表达形式存在差别，不失一般性，这里只给出当1d2d，1d3d和1d4d同时成立时，各种重叠情形下单颗粒所占网格的二维固相体积分数的表达式。图3.1颗粒与控制体的五种重叠情形图3.2颗粒与网格重叠时颗粒质心的区域划分图3.2给出了当1d2d，1d3d和1d4d同时成立时，各种重叠情形下颗粒质心与网格左下顶点的相对位置关系，与部分重叠的四种情形相对应，这里部分重叠时颗粒与顶点的位置关系也划分为四种情形。图3.2还显示了颗粒与网格发生重叠时，颗粒质心所处位置的精确范围。以下将分类讨论各种重叠情形下单个颗粒对网格的二维固相体积分数的贡献。(a)Cas1(b)Cas2图3.3颗粒与网格左下顶点最近的重叠情形（Cas1与Cas2）图3.3给出了颗粒与网格左下顶点距离最近的重叠情形1和2。这两种情形下，颗粒所占网格的二维固相体积分数的表达形式是唯一确定的，不存在变例。对于情形1，颗粒质心与网格的位置关系为：xXrx1且yYry1。则单个满足重叠情形1的颗粒i所占目标网格的二维固相体积分数为rr(x1,y1)(X,Y)21ryxfi(3.3)对于情形2，颗粒质心与网格的位置关系为：rxXx11，ryYy11，且ryYxX2121)()(。此情形下if计算为)arccos2sin(21)arccos2sin(21arccosarccos1121212122ryYrrxXrryYrrxXrryxfi(3.4)图3.4给出了颗粒与网格左下顶点距离最近的重叠情形3。对于情形3，颗粒质心与网格的位置关系存在四种变例，但均满足ryYxX2121)()(。尽管情形3存在四种变例，但经过推导，颗粒所占网格的二维固相体积分数的表达形式是统一的。根据图3.4容易得到)()(1212xXyYra(3.5))()(1212yYxXrb(3.6)则颗粒i所占目标网格的二维固相体积分数为)(21)(2122arccos2111122222xXbyYarbarryxfi(3.7)图3.4颗粒与网格左下顶点最近的重叠情形（Cas3）图3.5给出了颗粒与网格左下顶点距离最近的重叠情形4。情形4存在两种变例，并且颗粒所占网格的二维固相体积分数的表达形式也存在差异。对于情形4的变例1，颗粒质心与网格的位置关系所满足的条件为：rxXx11，1yY，且ryYxX2121)()(。此变例下if计算为)arccos2sin(21arccos112122ryYrryYrryxfi(3.8)对于情形4的变例2，颗粒质心与网格的位置关系为：1xX，ryYy11，且ryYxX2121)()(。此变例下if计算为abaaabb(x1,y1)(X,Y)br(a)Variation1(b)Variation2(c)Variation3(d)Variation4rrr)arccos2sin(21arccos112122rxXrrxXrryxfi(3.9)图3.5颗粒与网格左下顶点最近的重叠情形（Cas4）图3.6颗粒与网格左下顶点最近的重叠情形（Cas5）图3.6给出了颗粒与网格左下顶点距离最近的重叠情形5。情形5也存在两种变例，并且颗粒所占网格的二维固相体积分数的表达形式也存在差异。对于情形5的变例1，颗粒质心与网格的位置关系满足条件：11xXrx，1yY，且ryYxX2121)()(。此变例下有)arccos2sin(21arccos11212rYyrrYyryxfi(3.10)(x1,y1)(X,Y)r(x1,y1)r(X,Y)(a)Variation1(b)Variation2(x1,y1)(X,Y)r(x1,y1)r(X,Y)(a)Variation1(b)Variation2(x1,y1)对于情形5的变例2，颗粒质心与网格的位置关系为：1xX，11yYry，且ryYxX2121)()(。此变例下有)arccos2sin(21arccos11212rXxrrXxryxfi(3.11)在上述各种情形及其变例中，if的计算以1d2d，1d3d和1d4d同时成立为前提。类似地，不难导出当2d、3d和3d分别为最小时的各种if表达式。将if带入式(3.1)计算精确的s2和D2，由此给出了精确二维空隙率方案。