[發(fā)明專利]一種完全基于當(dāng)?shù)亓鲌鰠?shù)對轉(zhuǎn)捩模型可壓縮修正的方法

有效

• 申請?zhí)枺?/td>: 202110669135.7
• 申請日: 2021-06-17
• 公開（公告）號(hào): CN113361173B
• 公開（公告）日: 2022-04-26
• 發(fā)明（設(shè)計(jì)）人: 趙金山;柳森;張志剛;石義雷;陳挺;肖雨;粟斯堯;廖軍好
• 申請（專利權(quán)）人: 中國空氣動(dòng)力研究與發(fā)展中心超高速空氣動(dòng)力研究所
• 主分類號(hào): G06F30/23
• 分類號(hào): G06F30/23;G06F30/28;G06F111/10;G06F119/14
• 代理公司: 北京中濟(jì)緯天專利代理有限公司 11429
• 代理人: 王丹

1.一種完全基于當(dāng)?shù)亓鲌鰠?shù)對轉(zhuǎn)捩模型可壓縮修正的方法，其特征在于，所述的可壓縮修正方法，包括以下步驟：

a. 通過數(shù)值求解NS方程和原始無任何修正的轉(zhuǎn)捩模型，獲得流場基本參數(shù)γ、、；其中，γ表示間歇因子，反映邊界層轉(zhuǎn)捩區(qū)域的特定空間位置在時(shí)間歷程中處于層流或者湍流流態(tài)的概率；表示動(dòng)量厚度雷諾數(shù)；表示原始無任何修正的轉(zhuǎn)捩模型中的壓力梯度因子；

b. 通過引入湍流馬赫數(shù)，采用當(dāng)?shù)赝膭?dòng)能k和當(dāng)?shù)匾羲賏，計(jì)算得到流場當(dāng)?shù)赝牧黢R赫數(shù)；

c. 通過引入可壓縮修正關(guān)系式，計(jì)算出全流場內(nèi)每個(gè)網(wǎng)格點(diǎn)處的修正因子；其中，C₁表示常數(shù)，C₁的范圍為3~4；

d. 對原始無任何修正的轉(zhuǎn)捩模型轉(zhuǎn)捩動(dòng)量厚度雷諾數(shù)輸運(yùn)方程中的壓力梯度因子進(jìn)行修正，，并進(jìn)行全場迭代求解；其中，表示修正后的轉(zhuǎn)捩模型中的壓力梯度因子；如果結(jié)果沒有完全收斂，采用上一迭代歩更新后的流場基本參數(shù)γ、、，返回步驟a繼續(xù)計(jì)算；如果結(jié)果收斂，則結(jié)束計(jì)算。