[发明专利]基于残差摄动法的边界变化对流场影响差量的计算方法

权利要求书

1.基于残差摄动法的机身局部鼓包对流场影响程度差量的计算方法，其特征在于，包括：

步骤A1：分别为有鼓包飞机的机身与无鼓包飞机的机身建立计算网格，将两种构型的全机计算网格分为C1区域与C2区域，其中两种构型的机身的C1区域的计算网格完全相同，C2区域的计算网格不同；

步骤A2：针对无鼓包飞机建立流场矩阵方程Ax+b₁＝0；其中A为控制方程的矩阵；b₁为向量，表示根据无鼓包飞机建立的边界条件；x为向量，表示待求的流场参数；

步骤A3：计算使流场矩阵方程Ax+b₁＝0收敛于允许残差R₁的向量x的解x'₁；即x'₁满足Ax'₁+b₁＝R₁；

步骤A4：将边界条件更新为b₂，并改写流场矩阵方程，得到(Ax-R₁)+b₂＝0；b₂为根据有鼓包飞机建立的边界条件；

步骤A5：计算使新的流场矩阵方程(Ax-R₁)+b₂＝0收敛于摄动残差R′的向量x的解x'₂；即x'₂满足R′＝(Ax'₂-R₁)+b₂；

步骤A6：计算解x'₁与解x'₂之间的差值便得到飞机机身有鼓包时对流场的影响差量。

2.根据权利要求1所述的基于残差摄动法的机身局部鼓包对流场影响程度 差量的计算方法，其特征在于，步骤A5中，使用迭代方法求解流场矩阵方程(Ax-R₁)+b₂＝0时，将解x'₁中的无鼓包构型的C1区域的解作为有鼓包构型的C1区域的迭代初始值，将解x'₁中无鼓包构型的C2区域的解进行插值计算得到有鼓包构型的C2区域的迭代初始值。

3.基于残差摄动法的飞行器喷流对流场影响差量的计算方法，其特征在于，包括：

步骤B1：根据飞行器建立流场矩阵方程Ax+b₁＝0；其中A为控制方程的矩阵；b₁为向量，表示根据飞行器无喷流时建立的边界条件；x为向量，表示待求的流场参数；

步骤B2：计算使流场矩阵方程Ax+b₁＝0收敛于允许残差R₁的向量x的解x'₁；即x'₁满足Ax'₁+b₁＝R₁；

步骤B3：将边界条件更新为b₂，并改写流场矩阵方程，得到(Ax-R₁)+b₂＝0；b₂为根据飞行器有喷流时建立的边界条件；

步骤B4：计算使新的流场矩阵方程(Ax-R₁)+b₂＝0收敛于摄动残差R′的向量x的解x'₂；即x'₂满足R′＝(Ax'₂-R₁)+b₂；

步骤B5：计算解x'₁与解x'₂之间的差值便得到飞行器有喷流时对流场的影响差量。

4.根据权利要求3所述的基于残差摄动法的飞行器喷流对流场影响差量的计算方法，其特征在于，步骤B4中，使用迭代方法求解流场矩阵方程(Ax-R₁)+b₂＝0时，将无喷流时的流场参数解x'₁作为迭代的初始值。