[發明專利]彎曲邊界上的低雷諾數不可壓縮流中的壓力差的仿真方法

有效

• 申請號: 201310428248.3
• 申請日: 2013-09-18
• 公開（公告）號: CN103473418B
• 公開（公告）日: 2017-02-08
• 發明（設計）人: 梁朝湘;張超英;劉海燕;黃寄洪;冀肖榆;時偉
• 申請（專利權）人: 梧州學院
• 主分類號: G06F17/50
• 分類號: G06F17/50
• 代理公司: 廣州市越秀區海心聯合專利代理事務所(普通合伙)44295
• 代理人: 黃為,蔡國

1.一種彎曲邊界上的低雷諾數不可壓縮流中的壓力差的仿真方法，包括計算模塊和輸出模塊，其特征在于：在該方法中，計算模塊給定一低雷諾數不可壓縮流，該低雷諾數不可壓縮流在具有曲線邊界的管道內流動，所述的具有曲線邊界的管道為軸對稱的、非平直的并且管道壁為固壁邊界的管道，且該管道的出入口兩端壓差為△p，其中△p?=?p₁－p₂，p₁為入口處的壓力、p₂為出口處的壓力，并采用以x和r分別表示軸向坐標和徑向坐標的柱坐標系，該方法具體包括以下步驟，

步驟S1、計算模塊對于上述給定的低雷諾數不可壓縮流，判斷公式（1）是否成立，

ϵ=δx0<<1,]]>R0x0=O]]>或R₀<<?x₀，L>>?R₀，R_e<1????????????（1）；

其中，δ為壁厚，x₀為任意一點的軸向坐標，R₀為管道平直處的半徑，L為管道的長度，R_e為雷諾數；

步驟S2、若步驟S1中判斷公式（1）成立，計算模塊在模擬流場中生成或繼承粒子，并采用多邊界切線技術處理固壁邊界，根據該技術的需要隨之生成虛擬粒子；

步驟S3、計算模塊根據公式（2）計算每個粒子在任意位置x點處的體積力F_A2，并且入口處的體積力F_A1?由公式（3）確定；

FA2FA1≈∫A1f(r,x)dA∫A2f(r,x)dA---(2);]]>

FA1=Δpρ∫-0.5L0.5L∫A1f(r,x)dA∫A2f(r,x)dAdx---(3);]]>

其中，dA表示面積元，ρ為流體密度，Δp為壓力差；

步驟S4、輸出模塊輸出步驟S3中任意位置x點處的體積力F_A2。