1.一种计算流体力学软件采用的隐式方法，所述方法包括：步骤102，读取模拟对象的网格信息，更新所有网格单元的数据，进行网格初始化，生成网格的壁面距离；步骤104，判断所述壁面距离是否大于计算域阈值；步骤106，若所述壁面距离大于所述计算域阈值，则判定为内部区域，并采用块LU-SGS方法对NS方程进行时间离散；所述步骤106包括：步骤202，将网格单元中心的通量采用插值方法进行分裂，计算谱半径；步骤204，块LU-SGS初始化，对后续扫描过程中需使用的中间变量赋零值；步骤206，将NS方程离散化后的离散方程的左端矩阵分解成三个矩阵之和，即前扫矩阵L、矩阵D及后扫矩阵U；步骤208，利用上述前扫矩阵L、矩阵D以及后扫矩阵U，依次执行向前扫描、通讯中间量、向后扫描的操作，得到中间预测值以全场通量值；所述步骤206包括：步骤302，所述前扫矩阵L取值为：“小于当前网格单元序号的相邻网格单元”的无粘通量系数雅可比矩阵与该网格单元无粘通量谱半径之差，再加上该网格单元粘性通量系数雅可比矩阵与该网格单元粘性通量谱半径之差；步骤304：所述矩阵D取值为：当前网格单元的无粘通量系数雅可比矩阵与该网格单元粘性通量系数雅可比矩阵的之差，再加上该网格单元粘性通量的谱半径与粘性通量谱半径的差值；步骤306：所述后扫矩阵U取值为：“大于当前网格单元序号的相邻网格单元”的无粘通量系数雅可比矩阵与该网格单元无粘通量谱半径之差，再加上该网格单元粘性通量系数雅可比矩阵与该网格单元粘性通量谱半径之差；步骤108，若所述壁面距离小于或等于所述计算域阈值，则判定为物面边界附近区域，并采用点LU-SGS方法对NS方程进行时间离散；所述步骤108包括：步骤502，将网格单元中心的通量采用插值方法进行分裂，计算谱半径；步骤504，点LU-SGS初始化，对后续扫描过程中需使用的中间变量赋初值；步骤506，将NS方程离散化后的离散方程左端矩阵分解成三个矩阵之和，即前扫矩阵L’、对角矩阵D’及后扫矩阵U’；步骤508，利用上述前扫矩阵L’、对角矩阵D’以及后扫矩阵U’，依次执行向前扫描、通讯中间量、向后扫描的操作，得到中间预测值以全场通量值；所述步骤506包括：步骤602：对角矩阵D’取值为：当前网格单元无粘通量与粘性通量的谱半径之和，即当前网格单元无粘通量与粘性通量所组成合通量的系数雅可比矩阵的谱半径；步骤604：前扫矩阵L’取值为：“小于当前网格单元序号的相邻网格单元”的无粘通量系数雅可比矩阵与该网格单元无粘通量的谱半径之和，再加上该网格单元粘性通量的谱半径；步骤606：后扫矩阵U’取值为：“大于当前网格单元序号的相邻网格单元”的无粘通量系数雅可比矩阵与该网格单元无粘通量谱半径之差，再减去该网格单元粘性通量的谱半径；步骤110，更新流场变量的数值解，得到下一迭代时间步上的通量值。
2.根据权利要求1所述的方法，其特征在于，所述步骤208包括：步骤402，L步向前扫描运算，即计算某时间迭代步的残差值并作为所述离散方程的右端项，执行L块矩阵运算，进行全场逐点向前扫描得到中间预测值；步骤404：U步向后扫描运算，即将所述中间预测值作为所述离散方程的右端项，执行U块矩阵运算，进行全场向后扫描得到通量增量值。
3.根据权利要求1所述的方法，其特征在于，所述步骤508包括：步骤702：L’块向前扫描运算，即计算某时间层上的残差值，执行前扫矩阵L’块算子运算，进行全场逐点向前扫描得到中间预测值；步骤704：U’块向后扫描运算，即将所述中间预测值作为所述离散方程的右端项，执行后扫矩阵U’块算子运算，进行全场向后扫描得到通量增量值。
