1.多重网格的加速方法，其特征在于，所述方法包括：获取初始流场变量，并在细网格单元上进行迭代计算，得到下一个迭代步的第一流场变量；其中，所述初始流场变量至少包括速度、密度或压力中的一种或多种；将所述第一流场变量和残差从细网格单元传递到粗网格单元上，进行变量限制处理，得到处理后的流场变量；根据处理后的流场变量，在粗网格单元上进行一次迭代计算，得到流场变量的修正量；将所述流场变量的修正量插值到细网格单元上，并对细网格单元的第一流场变量进行修正；所述将所述流场变量的修正量插值到细网格单元上，并对细网格单元的第一流场变量进行修正，包括：将所述流场变量的修正量采用径向基函数差值方法插值到细网格单元上，得到细网格单元上的误差修正量；所述细网格单元的误差修正量为： 式中： 为细网格上的误差修正量， 为粗网格上的误差修正量， 为多变量插值的径向基函数的插值算子；通过细网格单元上的修正量对所述第一流场变量进行修正，修正公式为： 式中： 为在细网格上求解NS方程获得的流场变量， 为细网格上的误差修正量， 为细网格上修正之后的流场变量；所述将所述流场变量的修正量采用径向基函数差值方法插值到细网格单元上，得到细网格单元上的误差修正量，包括：对每个细网格单元，查找其所在的粗网格单元及邻居粗网格单元，构造得到N个粗网格单元形成的集合；构造插值系数方程，在N个粗网格单元形成的集合上，对每个粗网格单元进行处理，得到如下方程组： 其中：M是粗网格集形成的基函数矩阵，X、Y、Z是控制结果，Δ表示修正量， 是待求解的插值系数向量, 为多变量插值的径向基函数, N为参与插值的总的粗网格个数；求解插值系数方程，得到 三个插值系数向量；根据求解出来得到的插值系数 ，得到细网格单元上的误差修正量。
2.根据权利要求1所述的多重网格的加速方法，其特征在于，所述获取初始流场变量，并在细网格单元上进行迭代计算，得到下一个迭代步的第一流场变量，包括：根据无粘通量和粘性通量组成的强迫函数，计算细网格上的第一流场变量；其中，所述强迫函数为： ；其中：h为细网格，Q为流场变量， 为由无粘通量和粘性通量组成的强迫函数，L为粗网格向细网格上的插值算子。
3.根据权利要求2所述的多重网格的加速方法，其特征在于,所述将所述第一流场变量和残差从细网格单元传递到粗网格单元上，进行变量限制处理，得到处理后的流场变量，包括：采用体积加权平均的方式获得粗网格单元上的第二流场变量，将细网格上的第一流场变量限制到粗网格上；具体公式如下： ；式中：H表示粗网格， 表示将流场变量Q由细网格传递到粗网格上的限制算子，N表示粗网格单元包含N个细网格单元， 表示粗网格单元包含的第i个细网格单元的体积值，i为自然数； 为求解粗网格上方程获得的流场变量； 为细网格上修正之后的流场变量；采用所有单元残差直接求和的方式，将细网格单元上的残差限制到粗网格单元上； 式中： Res代表残差； 为粗网格上的残差， 为细网格上的残差； 是将残差由细网格传递到粗网格上的限制算子。
4.根据权利要求3所述的多重网格的加速方法，其特征在于，所述根据处理后的流场变量，在粗网格单元上进行一次迭代计算，得到流场变量的修正量，包括：在粗网格单元上迭代计算一次NS方程，得到粗网格单元上的离散方程；所述粗网格单元上的离散方程为： ；根据所述粗网格单元上的离散方程中的从细网格限制得到的残差部分 和利用从细网格上限制得到的流场变量 ，计算相应的残差 ；根据所述相应的残差，计算粗网格单元上的误差修正量；所述误差修正量的计算公式为： 其中： 为粗网格上的误差修正量， 为求解粗网格上方程获得的流场变量， 为由细网格限制获得的粗网格上的流场变量。
5.多重网格的加速装置，其特征在于，所述装置包括：获取模块，用于获取初始流场变量，并在细网格单元上进行迭代计算，得到下一个迭代步的第一流场变量；其中，所述初始流场变量至少包括速度、密度或压力中的一种或多种；处理模块，用于将所述第一流场变量和残差从细网格单元传递到粗网格单元上，进行变量限制处理，得到处理后的流场变量；计算模块，用于根据处理后的流场变量，在粗网格单元上进行一次迭代计算，得到流场变量的修正量；修正模块，用于将所述流场变量的修正量插值到细网格单元上，并对细网格单元的第一流场变量进行修正;所述修正模块具体用于：将所述流场变量的修正量采用径向基函数差值方法插值到细网格单元上，得到细网格单元上的误差修正量；所述细网格单元的误差修正量为： 式中： 为细网格上的误差修正量， 为粗网格上的误差修正量， 为多变量插值的径向基函数的插值算子；通过细网格单元上的修正量对所述第一流场变量进行修正，修正公式为： 式中： 为在细网格上求解NS方程获得的流场变量， 为细网格上的误差修正量， 为细网格上修正之后的流场变量；所述将所述流场变量的修正量采用径向基函数差值方法插值到细网格单元上，得到细网格单元上的误差修正量，包括：对每个细网格单元，查找其所在的粗网格单元及邻居粗网格单元，构造得到N个粗网格单元形成的集合；构造插值系数方程，在N个粗网格单元形成的集合上，对每个粗网格单元进行处理，得到如下方程组： 其中：M是粗网格集形成的基函数矩阵，X、Y、Z是控制结果，Δ表示修正量， 是待求解的插值系数向量, 为多变量插值的径向基函数, N为参与插值的总的粗网格个数；求解插值系数方程，得到 三个插值系数向量；根据求解出来得到的插值系数 ，得到细网格单元上的误差修正量。
6.根据权利要求5所述的多重网格的加速装置，其特征在于，所述获取模块用于：根据无粘通量和粘性通量组成的强迫函数，计算细网格上的第一流场变量；其中，所述强迫函数为： ；其中：h为细网格，Q为流场变量， 为由无粘通量和粘性通量组成的强迫函数，L为粗网格向细网格上的插值算子。
7.一种终端设备，其特征在于，包括：至少一个处理器和存储器；所述存储器存储计算机程序；所述至少一个处理器执行所述存储器存储的计算机程序，以实现权利要求1-4中任一项所述的多重网格的加速方法。
8.一种计算机可读存储介质，其特征在于，该计算机可读存储介质中存储有计算机程序，所述计算机程序被执行时实现权利要求1-4中任一项所述的多重网格的加速方法。
