[发明专利]一种数学模型数值计算程序的验证方法

权利要求书

1.一种数学模型数值计算程序的验证方法，其特征在于：所述方法包括以下步骤：

1)确定程序求解的控制方程，假设该控制方程的边界条件；

2)确定步骤1)中控制方程的离散方法的理论收敛阶；

3)确定待测试的数学模型的所有分支程序；并根据边界条件判断控制方程是否适定分支程序，若适定，则进行步骤4)；

4)确定控制方程的虚构解；

5)将虚构解输入控制方程的右端形成源项；

6)将源项加入原控制方程的右端，并完成计算程序的输入，得到控制方程的数值解；

7)判断数值解和虚构解的值是否一致，若一致，则进行步骤8)；

8)进行网格细化来获得不同网格上的误差，并利用数值解计算控制方程的实际收敛阶；

9)将实际收敛阶和理论收敛阶进行对比，若实际收敛阶与理论收敛阶一致，则完成当前测试。

2.根据权利要求1所述的数学模型数值计算程序的验证方法，其特征在于：所述步骤7)中若数值解与虚构解不一致，则进行步骤10)查错并改正错误，查到错误并进行改正后，则返回步骤7)。

3.根据权利要求2所述的数学模型数值计算程序的验证方法，其特征在于：所述步骤9)若实际收敛阶和理论收敛阶不一致，则进行步骤10)查错并改正错误，查到错误并进行改正后，则返回步骤7)。

4.根据权利要求3所述的数学模型数值计算程序的验证方法，其特征在于：所述步骤9)之后还包括步骤11)判断是否需要确定新的虚构解，若是，则返回步骤4)继续。

5.根据权利要求1或2或3或4所述的数学模型数值计算程序的验证方法，其特征在于：所述步骤1)中控制方程的边界条件包括可压或不可压、有粘或无粘。

6.根据权利要求5所述的数学模型数值计算程序的验证方法，其特征在于：所述步骤8)中网格细化的具体方式是：对结构网格整数倍进行加密，对非结构网格对每个三角形单元的各边中点连线进行剖分细化。

7.根据权利要求6所述的数学模型数值计算程序的验证方法，其特征在于：所述步骤10)中实际收敛阶和理论收敛阶不一致的原因包括：测试公式或建立的测试存在的错误、计算精度阶在给定的时间、位置上存在错误、计算误差度量或实际精度存在错误或者收敛不完全。