[发明专利]一种改进的三维二物质欧拉界面处理算法

权利要求书

1.一种改进的三维二物质欧拉界面处理算法，整个算法由以下步骤组成：由混合网格周围26个网格的体积份额来确定物质平面的法线方向，根据该混合网格中的某物质体积份额确定平面的位置，计算出通过该混合网格边界的输运量，其特征在于增加物质输运次序的计算内容。

2.根据权利要求1所述的三维二物质欧拉界面处理算法，其特征在于增加物质输运次序的计算内容，其详细内容如下：

当沿x轴正向输运时，设V_l1＝V_ratio(id-1，j，k，1)，V_r1＝V_ratio(ic，j，k，1)，V_l2＝V_ratio(id-1，j，k，2)，V_r2＝V_ratio(ic，j，k，2)，V_l1、V_r1、V_l2和V_r2分别为左边相临网格1物质、右边相临网格1物质、左边相临网格2物质和右边相临网格2物质体积比，某物质的体积比是指网格中某物质的体积除以网格总体积所得到的百分比；如果V_l1大于Epsv则I_l1＝1，否则，I_l1＝0；如果V_r1大于Epsv则I_r1＝1，否则，I_r1＝0；如果V_l2大于Epsv则I_l2＝1，否则，I_l2＝0；如果V_r2大于Epsv则I_r2＝1，否则，I_r2＝0，Epsv为体积比误差限；定义如下两个变量(sgn为符号函数)：L₁＝I_l1·sgn(V_r1-V_l1)；R₁＝I_r1·sgn(V_r1-V_l1)；L₂＝I_l2·sgn(V_r2-V_l2)；R₂＝I_r2·sgn(V_r2-V_l2)；根据L₁、R₁、L₂和R₂的值，作出如表1所示的五种输运构型，这五种构型的输运优先级顺序为从高到低，即：①＞②＞③＞④＞⑤，例如，若1物质属于构型③，2物质输运构型⑤，则1物质的输运优先级大于2物质，因此在输运时，优先输运1物质的输运量，不足部分以2物质进行补充；当沿x轴负向输运时，仿照上面可得到相应的输运构型，沿y、z轴正负方向的计算方法与x轴类似。

表1物质输运构型表