[发明专利]基于物理量梯度修正的紧支基函数多场耦合数据传递方法

权利要求书

1.一种基于物理量梯度修正的紧支基函数多场耦合数据传递方法，其特征在于，包括以下步骤：

S1、对于所有已知物理量网格点进行物理量梯度的计算；

S2、根据梯度值对紧支基函数进行各向异性修正，并根据修正结果对紧支半径计算公式进行更新；

S3、根据更新后的紧支半径计算公式得到修正后的径向基的取值，并按照紧支基函数标准插值流程进行数据传递。

2.根据权利要求1所述的基于物理量梯度修正的紧支基函数多场耦合数据传递方法，其特征在于，所述步骤S1中的网格点包括内部点、网格块角点和边界点。

3.根据权利要求2所述的基于物理量梯度修正的紧支基函数多场耦合数据传递方法，其特征在于，所述步骤S1的具体步骤为：

S11、遍历所有内部点对应的网格块角点和边界点；

S12、根据内部点的相邻网格块角点的物理量和该内部点的物理量计算内部点的梯度；

S13、将边界点的相邻内部点梯度赋值为该边界点的梯度；

S14、将网格块角点最近的边界点梯度赋值为该网格块角点的梯度。

4.根据权利要求3所述的基于物理量梯度修正的紧支基函数多场耦合数据传递方法，其特征在于，所述步骤S12中内部点的梯度的计算公式为：

上式中，U_i为网格点i上物理量的值，U_j为网格点j上物理量的值，Δ(·)_ij＝(·)_j-(·)_i，i、j均为网格点的序号，(x,y,z)为网格点的坐标，N为网格点个数，θ_j为权系数。

5.根据权利要求1所述的基于物理量梯度修正的紧支基函数多场耦合数据传递方法，其特征在于，所述步骤S2的具体步骤为：

S21、计算所有网格点在x,y,z方向上的梯度总和；

S22、根据网格点在x,y,z方向上的梯度总和计算网格点在x,y,z方向上的平均梯度；

S23、对网格点在x,y,z方向上的平均梯度实现几何外形上的尺度归一；

S24、根据尺度归一后的平均梯度对所有网格点的坐标进行缩放；

S25、根据径向基函数公式和缩放后的坐标得到更新后的紧支半径计算公式。

6.根据权利要求5所述的基于物理量梯度修正的紧支基函数多场耦合数据传递方法，其特征在于，所述步骤S22中网格点在x,y,z方向上的平均梯度的计算公式为：

上式中，当i＝1时，Ave_gradx₁＝Ave_gradx，ALL_gradx₁＝ALL_gradx；当i＝2时，Ave_gradx₂＝Ave_grady，ALL_gradx₂＝ALL_grady；当i＝3时，Ave_gradx₃＝Ave_gradz，ALL_gradx₃＝ALL_gradz；Ave_gradx、Ave_grady、Ave_gradz分别为网格点在x、y、z方向上的平均梯度，ALL_gradx、ALL_grady、ALL_gradz分别为网格点在x、y、z方向上的梯度总和，i，j均为网格点的序号，N为网格点个数。

7.根据权利要求6所述的基于物理量梯度修正的紧支基函数多场耦合数据传递方法，其特征在于，所述步骤S23中尺度归一的计算公式为：

上式中，Ave_gradx′_i为尺度归一后的平均梯度，Ave_gradx'、Ave_grady'、Ave_gradz'为尺度归一后x、y、z方向上的平均梯度。