多重网格的加速方法、装置、终端设备及存储介质制造方法及图纸

本申请公开了一种多重网格的加速方法、装置、终端设备和存储介质，通过获取初始流场变量，并在细网格单元上进行迭代计算，得到下一个迭代步的第一流场变量；将第一流场变量和残差从细网格单元传递到粗网格单元上，进行变量限制处理，得到处理后的流场变量；根据处理后的流场变量，在粗网格单元上进行一次迭代计算，得到流场变量的修正量；将流场变量的修正量插值到细网格单元上，并对细网格单元的第一流场变量进行修正，采用径向基函数插值方法进行插值，考虑了周围单元对细网格单元上修正量的影响，提升了修正量的插值精度，能提升多重网格方法的收敛效率，求解精度高，真实物理问题的模拟程度高、从而大幅度地提高了鲁棒性。从而大幅度地提高了鲁棒性。从而大幅度地提高了鲁棒性。

【技术实现步骤摘要】
多重网格的加速方法、装置、终端设备及存储介质


[0001]本申请属于流体力学
，尤其涉及一种多重网格的加速方法、装置、终端设备及存储介质。

技术介绍

[0002]计算流体力学（Computational Fluid Dynamics， CFD），是一门利用计算机，数值求解流动控制NS方程的交叉学科。在工业产品研发过程中，通过借助CFD数值模拟，可大幅降低研发费用，缩短设计周期，提高产品的精细化设计水平，为国防工业和国民经济装备建设发展和数字化转型提供有力支撑。
[0003]随着CFD和计算机技术的发展，CFD模拟的流动问题也越来越复杂，计算网格规模也随之增大。在有限的计算域中，随着网格数量的增加，网格尺度不断减小，在稳定性条件下，局部推进时间步长减小。因此，CFD计算收敛速度将大幅减小，计算时间将大幅增加，常常不能满足工程应用的需要，亟需发展一种高效的加速收敛技术。
[0004]现有技术中可采用多重网格法来进行加速收敛，在CFD中已得到广泛的发展和应用，在计算过程中，通过采用一阶线性插值，将粗网格单元上的值直接对细网格单元进行修正，并未考虑周围其他粗网格的影响，插值精度较低，从而多重网格的加速效果大幅下降，如何能够加速多重网格的加速速率、且提高插值精度，是目前急需解决的问题。

技术实现思路

[0005]本申请意在提供一种多重网格的加速方法、装置、终端设备和存储介质，以解决现有技术中存在的不足，本申请要解决的技术问题通过以下技术方案来实现。
[0006]第一个方面，本申请实施例提供一种多重网格的加速方法，所述方法包括：获取初始流场变量，并在细网格单元上进行迭代计算，得到下一个迭代步的第一流场变量；其中，所述初始流场变量至少包括速度、密度或压力中的一种或多种；将所述第一流场变量和残差从细网格单元传递到粗网格单元上，进行变量限制处理，得到处理后的流场变量；根据处理后的流场变量，在粗网格单元上进行一次迭代计算，得到流场变量的修正量；将所述流场变量的修正量插值到细网格单元上，并对细网格单元的第一流场变量进行修正。
[0007]可选地，所述获取初始流场变量，并在细网格单元上进行迭代计算，得到下一个迭代步的第一流场变量，包括：根据无粘通量和粘性通量组成的强迫函数，计算细网格上的第一流场变量；其中，所述强迫函数为：L
h
(Q)=f
h
；其中： h为细网格，Q为流场变量，f
h
由无粘通量和粘性通量组成的强迫函数，L为粗网格向细网格上的插值算子。
[0008]可选地，所述将所述第一流场变量和残差从细网格单元传递到粗网格单元上，进行变量限制处理，得到处理后的流场变量，包括：采用体积加权平均的方式获得粗网格单元上的第二流场变量，将细网格上的第一流场变量限制到粗网格上；具体公式如下： ；式中：H表示粗网格，表示将流场变量 Q由细网格传递到粗网格上的限制算子，N表示粗网格单元包含N个细网格单元， V
i
表示粗网格单元包含的第i个细网格单元的体积值，i为自然数；采用所有单元残差直接求和的方式，将细网格单元上的残差限制到粗网格单元上；式中： Res代表残差；Res
H
为粗网格上的残差，Res
h
为细网格上的残差；是将残差由细网格传递到粗网格上的限制算子。
[0009]可选地，所述根据处理后的流场变量，在粗网格单元上进行一次迭代计算，得到流场变量的修正量，包括：在粗网格单元上迭代计算一次NS方程，得到粗网格单元上的离散方程；所述粗网格单元上的离散方程为：根据所述粗网格单元上的离散方程中的从细网格限制得到的残差部分和利用从细网格上限制得到的流场变量L
h
(Q
H
)，计算相应的残差；根据所述相应的残差，计算粗网格单元上的误差修正量；所述误差修正量的计算公式为： 其中 v
H
为粗网格上的误差修正量， Q
H
为求解粗网格上方程获得的流场变量，为由细网格限制获得的粗网格上的流场变量。
[0010]可选地，所述将所述流场变量的修正量插值到细网格单元上，并对细网格单元的第一流场变量进行修正，包括：将所述流场变量的修正量采用径向基函数差值方法插值到细网格单元上，得到细网格单元上的误差修正量；所述细网格单元的误差修正量为： 式中：v
h
为细网格上的误差修正量，v
H
为粗网格上的误差修正量，为从RBF插值算子；通过细网格单元上的修正量对所述第一流场变量进行修正，修正公式为： 式中： 为S1中在细网格上求解NS方程获得的流场变量，v
h
为细
网格上的误差修正量，为细网格上修正之后的流场变量。
[0011]可选地，所述将所述流场变量的修正量采用径向基函数差值方法插值到细网格单元上，得到细网格单元上的误差修正量，包括：对每个细网格单元，查找其所在的粗网格单元及邻居粗网格单元，构造得到N个粗网格单元形成的集合；构造插值系数方程，在N个粗网格单元形成的集合上，对每个粗网格单元进行处理，得到如下方程组：得到如下方程组：其中：M是粗网格集形成的基函数矩阵，X、Y、Z是控制结果，Δ表示修正量，A
x
、A
y
、A
z
是待求解的插值系数向量，
ϕ
为RBF插值基函数，N为参与插值的总的粗网格个数；求解插值系数方程，到A
x
、A
y
、A
z
三个插值系数向量；根据求解出来得到的插值系数A
x
、A
y
、A
z
，到细网格单元上的误差修正量。
[0012]第二个方面，本申请实施例提供一种多重网格的加速装置，所述装置包括：获取模块，用于获取初始流场变量，并在细网格单元上进行迭代计算，得到下一个迭代步的第一流场变量；其中，所述初始流场变量至少包括速度、密度或压力中的一种或多种；处理模块，用于将所述第一流场变量和残差从细网格单元传递到粗网格单元上，进行变量限制处理，得到处理后的流场变量；计算模块，用于根据处理后的流场变量，在粗网格单元上进行一次迭代计算，得到流场变量的修正量；修正模块，用于将所述流场变量的修正量插值到细网格单元上，并对细网格单元的第一流场变量进行修正。
[0013]可选地，所述获取模块用于：根据无粘通量和粘性通量组成的强迫函数，计算细网格上的第一流场变量；其中，所述强迫函数为：L
h
(Q)=f
h
；其中： h为细网格，Q为流场变量，f
h
由无粘通量和粘性通量组成的强迫函数，L为粗网格向细网格上的插值算子。
[0014]可选地，所述处理模块用于：采用体积加权平均的方式获得粗网格单元上的第二流场变量，将细网格上的第一流场变量限制到粗网格上；具体公式如下： ；式中：H表示粗网格，表示将流场变量 Q由细网格传递到粗网格上的限制算子，N表示粗网格单元包含N个本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.多重网格的加速方法，其特征在于，所述方法包括：获取初始流场变量，并在细网格单元上进行迭代计算，得到下一个迭代步的第一流场变量；其中，所述初始流场变量至少包括速度、密度或压力中的一种或多种；将所述第一流场变量和残差从细网格单元传递到粗网格单元上，进行变量限制处理，得到处理后的流场变量；根据处理后的流场变量，在粗网格单元上进行一次迭代计算，得到流场变量的修正量；将所述流场变量的修正量插值到细网格单元上，并对细网格单元的第一流场变量进行修正。2.根据权利要求1所述的多重网格的加速方法，其特征在于，所述获取初始流场变量，并在细网格单元上进行迭代计算，得到下一个迭代步的第一流场变量，包括：根据无粘通量和粘性通量组成的强迫函数，计算细网格上的第一流场变量；其中，所述强迫函数为：L
h (Q)=f
h
；其中： h为细网格，Q为流场变量，f
h
由无粘通量和粘性通量组成的强迫函数， L为粗网格向细网格上的插值算子。3.根据权利要求2所述的多重网格的加速方法，其特征在于,所述将所述第一流场变量和残差从细网格单元传递到粗网格单元上，进行变量限制处理，得到处理后的流场变量，包括：采用体积加权平均的方式获得粗网格单元上的第二流场变量，将细网格上的第一流场变量限制到粗网格上；具体公式如下：；式中： H表示粗网格， 表示将流场变量 Q由细网格传递到粗网格上的限制算子，N表示粗网格单元包含N个细网格单元， V
i
表示粗网格单元包含的第i个细网格单元的体积值，i为自然数；Q
H
为求解粗网格上方程获得的流场变量；Q
h
为细网格上修正之后的流场变量；采用所有单元残差直接求和的方式，将细网格单元上的残差限制到粗网格单元上； 式中：Res代表残差；Res
H
为粗网格上的残差，Res
h
为细网格上的残差；是将残差由细网格传递到粗网格上的限制算子。4.根据权利要求3所述的多重网格的加速方法，其特征在于，所述根据处理后的流场变量，在粗网格单元上进行一次迭代计算，得到流场变量的修正量，包括：在粗网格单元上迭代计算一次NS方程，得到粗网格单元上的离散方程；所述粗网格单元上的离散方程为：；根据所述粗网格单元上的离散方程中的从细网格限制得到的残差部分和利用
从细网格上限制得到的流场变量L
h
(Q
H
)，计算相应的残差；根据所述相应的残差，计算粗网格单元上的误差修正量；所述误差修正量的计算公式为： 其中： v
H
为粗网格上的误差修正量，Q
H
为求解粗网格上方程获得的流场变量，为由细网格限制获得的粗网格上的流场变量。5.根据权利要求4所述的多重网格的加速方法，其特征在于，所述将所述流场变量的修正量插值到细网格单元上，并对细网格单元的第一流场变量进行修正，包括：将所述流场变量的修正量采用径向基函数差值方法插值到细网格单元上，得到细网格单元上的误差修正量；所述...

【专利技术属性】
技术研发人员：赵钟，万云博，何磊，何先耀，徐刚，陈坚强，
申请(专利权)人：中国空气动力研究与发展中心计算空气动力研究所，
类型：发明
国别省市：