本发明公开了面向复杂飞行器全速域数值模拟的非线性数值迭代求解方法，包括：基于度量张量生成计算域内各向异性网格，构建可压缩流动的数值离散格式与RANS湍流模型，并采用Gauss–Seidel方法对RANS偏微分方程组的非线性预算子进行迭代计算，减小了各向异性网格在迭代过程中可能出现的高频误差。为了快速访问网格节点并保证节点插值解的准确性，实施节点模板增强策略和高效的着色遍历算法；并采用广义共轭残差法进行前后迭代解静态残差评估，使得计算效率和鲁棒性大为提高。本发明擅长处理各向异性混合网格，既能保证大宽高比边界层网格上数值解的迭代鲁棒性，又能保证在大尺寸粗网格上数值解的收敛效率，进而保证整个非线性方程求解迭代过程的高效和稳定。