[發明專利]基于變量差及相對位移的r型網格自適應移動方法及設備

有效

• 申請號: 202110275460.5
• 申請日: 2021-03-15
• 公開（公告）號: CN112665820B
• 公開（公告）日: 2021-06-04
• 發明（設計）人: 龔小權;李明;馬戎;付云峰;張培紅;崔鵬程;趙輝;何協;唐靜;程鋒;周乃春;賈洪印
• 申請（專利權）人: 中國空氣動力研究與發展中心計算空氣動力研究所
• 主分類號: G01M9/06
• 分類號: G01M9/06
• 代理公司: 成都九鼎天元知識產權代理有限公司 51214
• 代理人: 徐靜

1.一種基于變量差及相對位移的r型網格自適應移動方法，其特征在于，以歸一化網格點的變量差，以及網格點相對位移量這兩類參數作為計算網格點移動的權函數，通過加入網格點相對位移量這一權函數，實現網格點的跟隨移動；歸一化消除所述兩類參數的量值及量綱對網格點移動量影響，以所述兩類參數形成新的組合權值，以此確定激波附近網格點的位移，并將激波附近網格點的位移信息向網格點的相鄰網格點傳遞，實現激波附近的幾排網格點都同步向激波位置移動；

以壓力這一物理量作為特征變量，計算網格中每個網格點的變量值，包括以下步驟：

步驟一：基于距離權計算所有網格點的壓力；

步驟二：統計所有網格點的壓力最大值和壓力最小值；

步驟三：根據所有網格點的壓力最大值和壓力最小值歸一化網格點上的壓力；

步驟四：將網格點與相鄰網格點的壓力之差作為網格點移動位移的第一個權值；

步驟五：將上一迭代步網格點移動位移與網格點之間的初始位移的比值作為網格點移動位移的另一個權值；

步驟六：通過加權當前網格點到周圍所有網格點的矢量獲得網格點位移；

步驟一中，根據格心型二階有限體積方法，網格單元的變量存儲在網格單元體心，將網格點a周圍網格單元體心變量通過距離加權到網格點a：

其中，k為網格點a周圍的網格單元，n為網格點a周圍的網格單元數量，且從0開始計數；p_k為網格單元k體心的壓力，p_a為網格點a的壓力，為網格點a到網格單元k體心的距離。

2.根據權利要求1所述的一種基于變量差及相對位移的r型網格自適應移動方法，其特征在于，步驟二中，所有網格點的壓力最大值p_max和壓力最小值p_min為：

其中，j為所有網格點，即j=a,b,c,d,e；b,c,d,e為網格點a的相鄰網格點。

3.根據權利要求2所述的一種基于變量差及相對位移的r型網格自適應移動方法，其特征在于，步驟三中，歸一化網格點壓力的計算方法為：

其中，為網格點壓力歸一化后的值，scale為網格點壓力影響系數，值為大于等于1的整數。

4.根據權利要求3所述的一種基于變量差及相對位移的r型網格自適應移動方法，其特征在于，步驟四中，網格點a移動位移的第一個權值的計算方法如下：

其中，為網格點a壓力歸一化后的值，為網格點i壓力歸一化后的值，i=b,c,d,e；b,c,d,e為網格點a的相鄰網格點。