[发明专利]面处理的特征无反射边界条件设定方法与系统

权利要求书

1.一种面处理的特征无反射边界条件设定方法，其特征在于，包括步骤：

根据纳维叶－斯托克斯方程，对可压缩流体进行数值模拟，获取可压缩流体控制方程，其中，所述可压缩流体控制方程包括连续方程、动量方程、温度方程、整理方程和状态方程；

根据所述可压缩流体控制方程，获取可压缩流体在三维空间坐标轴X、Y和Z三个方向上的传播特征波方程；

根据所述可压缩流体在三维空间坐标轴X、Y和Z三个方向上的传播特征波方程，忽略横向项、可压缩流体的粘性以及源项，通过数学运算，生成一维无粘特征发射方程，所述一维无粘特征发射方程用于获取所述可压缩流体从计算区域外部进入计算区域的特征波的关系；

获取所述可压缩流体的求解变量数据，根据所述一维无粘特征发射方程，生成所述可压缩流体特征波形式的控制方程，计算所述横向项、所述可压缩流体的粘性以及所述源项，并将计算结果与所述所述可压缩流体特征波形式的控制方程整合，生成可压缩流体进口面求解变量时间推进方程；

根据所述可压缩流体特征波形式的控制方程，获得可压缩流体出口面中求解变量的时间导数方程，根据所述可压缩流体出口面中求解变量的时间导数方程，生成可压缩流体出口面求解变量的时间推进方程；

根据所述压缩流体进口面求解变量时间推进方程和所述出口面求解变量时间推进方程，确定可压缩流体中面边界的条件。

2.根据权利要求1所述的面处理的特征无反射边界条件设定方法，其特征在于，所述可压缩流体特征波形式的控制方程具体如下：

∂ρ∂t+1c2[L2+12(L1+L5)]+V1=0]]>

L₁，L₂，L₅为特征波的振幅，c是声速，V1是所述横向项，ρ为所述可压缩流体的密度，t为时间。

3.根据权利要求1或2所述的面处理的特征无反射边界条件设定方法，其特征在于，所述根据所述压缩流体进口面求解变量时间推进方程和所述出口面求解变量时间推进方程，确定可压缩流体中面边界的条件具体包括步骤：

根据所述压缩流体进口面求解变量时间推进方程和所述出口面求解变量时间推进方程，选取相应的时间推进格式；

采用所述时间推进格式，获得下一个时间上所述压缩流体进口面和出口面的求解变量的更新值；

根据下一个时间上所述压缩流体进口面和出口面的求解变量的更新值，确定可压缩流体中面边界的条件。

4.根据权利要求1或2所述的面处理的特征无反射边界条件设定方法，其特征在于，所述根据所述压缩流体进口面求解变量时间推进方程和所述出口面求解变量时间推进方程，确定可压缩流体中面边界的条件之后还有步骤：

对确定可压缩流体中面边界的条件进行验证。

5.根据权利要求1或2所述的面处理的特征无反射边界条件设定方法，其特征在于，面求解变量包括可压缩流体的速度、密度和温度。