[发明专利]液滴碰壁运动的相界面的获取方法和装置

权利要求书

1.一种液滴碰壁运动的相界面的获取方法，其特征在于，包括以下步骤：

S1、根据计算域中计算单元所在的区域对第一相函数和第二相函数进行初始化，并设定所述计算单元的初始速度；

S2、获取两相流流动的单相和混合相控制方程并求解以获取所述计算单元的速度和压强，并根据所述计算单元的速度和压强对所述计算单元进行更新；

S3、分别对所述第一相函数和所述第二相函数的对流输运方程进行求解，以获取求解后的第一相函数和第二相函数；

S4、判断液滴与壁面是否发生接触；

S5、如果是，则根据壁面润湿模型获取接触角，并根据所述接触角获取相界面的法向向量，以及根据所述法向向量和所述求解后的第二相函数获取所述相界面的位置，其中，所述相界面为气相区域和液相区域的交界面；

S6、如果否，则根据所述求解后的第一相函数获取所述相界面的法向向量，并根据所述法向向量和所述求解后的第二相函数获取所述相界面的位置；

S7、对所述求解后的第一相函数进行重新距离化；以及

重复执行S2-S7，直至达到预设结束条件，并将达到所述预设结束条件时的相界面的位置作为最终结果。

2.如权利要求1所述的液滴碰壁运动的相界面的获取方法，其特征在于，所述第一相函数为Level Set函数，所述第二相函数为VOF函数。

3.如权利要求2所述的液滴碰壁运动的相界面的获取方法，其特征在于，所述两相流流动的单相和混合相控制方程为：

▿u→=0]]>

ρ(φ)(∂u→∂t+u→·▿u→)=-▿p+ρ(φ)g→-σκ(φ)▿H(φ)+▿τ,---(1)]]>

其中，为速度，p为压强，φ为所述Level Set相函数，σ为表面张力系数，κ为相界面曲率，τ为粘性应力张量，ρ为密度，为重力加速度，H(φ)为Heaviside函数。