1風(fēng)扇流場分析

1.1案例介紹

風(fēng)扇可以用于發(fā)動機的冷卻等很多場景，合理的風(fēng)扇設(shè)計將極大地提高風(fēng)扇的效率，但由于管道風(fēng)扇內(nèi)部流動非常復(fù)雜，通過理論計算對其流動進行定性分析十分困難，風(fēng)洞試驗雖然可以得到其流動參數(shù)和噪聲特性，但也無法對流場內(nèi)部的流動細(xì)節(jié)進行描述。

本案例演示如何利用Fluent進行風(fēng)扇流動特性和噪聲特性計算。

1.2幾何建模和流場計算域建立

本案例風(fēng)扇外徑為384mm，輪轂直徑為140mm，輪轂比為0.365，8扇葉均勻分布，外流場建模充分考慮到進氣試驗標(biāo)準(zhǔn)，入口區(qū)長度至少為入口處管道直徑的六倍；而出口區(qū)的長度則應(yīng)保證至少為出口位置管道直徑的十倍；至于旋轉(zhuǎn)流體區(qū)，是指包含了風(fēng)扇本體以及周圍流場的圓柱體區(qū)域，應(yīng)當(dāng)保證其尺寸盡量靠近風(fēng)扇葉片的直徑，最終風(fēng)扇模型和外流場模型分別如下圖所示。

1.3模型網(wǎng)格的劃分

網(wǎng)格生成作為仿真計算中的關(guān)鍵環(huán)節(jié)，其結(jié)果直接控制了后續(xù)計算過程的效率與精度。為了保證劃分結(jié)果的質(zhì)量，應(yīng)選擇合適的網(wǎng)格尺寸，防止太疏或太密的網(wǎng)格產(chǎn)生，在流量梯度較大的流動區(qū)域內(nèi)，應(yīng)當(dāng)盡量提高網(wǎng)格質(zhì)量（高細(xì)密度，較小的歪斜度）；至于梯度小的區(qū)域可以在保證精度的基礎(chǔ)上適當(dāng)較少網(wǎng)格數(shù)目。

本案例旋轉(zhuǎn)流體區(qū)由于包含了風(fēng)扇本體且流動情況最為復(fù)雜，為了保證足夠的計算精度，該區(qū)域網(wǎng)格尺寸最小。管道區(qū)網(wǎng)格尺寸較旋轉(zhuǎn)區(qū)略大，最終劃分結(jié)果如下圖。

1.4邊界條件設(shè)定與旋轉(zhuǎn)模型選取

完成網(wǎng)格生成后需進行邊界條件的設(shè)置。在流動的計算過程需要設(shè)定的邊界條件包括：

（1） 流動入口條件：根據(jù)吸氣試驗的要求將流動入口設(shè)置為壓力邊界條件，其中入口處壓值定義為大氣壓力，且氣體沿軸線方向流動；

（2） 流動出口條件：根據(jù)吸氣試驗的要求將流動入口設(shè)置為壓力邊界條件，出口壓力值定義為 0，即出口處沒有外界的作用；

（3） 壁面邊界條件：主要為通流區(qū)的管壁表面。

對風(fēng)扇旋轉(zhuǎn)運動的仿真則是通過 MRF 模型來實現(xiàn)的。Fluent 中常用的多運動坐標(biāo)系模型包括： SMM（滑移網(wǎng)格模型）,MPM(混合面模型)以及 MRF（多重參考系模型）。考慮到風(fēng)扇中氣體運動屬于定常流動，所以選擇計算量相對較少的 MRF 基準(zhǔn)。

作為旋轉(zhuǎn)機械仿真中最常使用的模型，MRF 模型計算思路在于：將算法區(qū)間分成數(shù)個運動相互獨立的子區(qū)間，先在各子部分間對流場方程進行求解，通過各部分間的交界面完成流場信息的傳遞。

作為 CFD 模型中唯一運動的旋轉(zhuǎn)流體域，將其邊界條件設(shè)定為 Fluid（流動域）。在 Fluent 軟件中選擇 MRF，并且將風(fēng)機轉(zhuǎn)速定義為坐標(biāo)系的轉(zhuǎn)速。

1.5計算方程選擇與仿真參數(shù)設(shè)置

對于風(fēng)扇內(nèi)部的穩(wěn)態(tài)流動，采用定常計算模型進行仿真，且計算過程中不考慮重力的影響。利用 SIMPLE 方程完成速度與壓力的解耦，將湍流模型定義為 RNG k-epsilon;的雙方程模型。

在 Fluent 軟件對參數(shù)進行設(shè)定時，根據(jù)實際工況將流體材料定義為空氣且認(rèn)為風(fēng)扇內(nèi)部流體不可壓縮；由于流動過程中沒有熱能的交換所以不對能量守恒方程進行求解，只考慮流體連續(xù)方程以及動量方程。對于控制方程中的湍動能耗散項以及動量項等使用二階迎風(fēng)的離散格式，在迭代過程中使用欠松馳因子以加速收斂。

1.6風(fēng)扇流場計算結(jié)果分析

用Fluent軟件對轉(zhuǎn)速為2000rpm的風(fēng)扇進行計算，得到包括速度矢量圖、壓力云圖結(jié)果如下所示。

2風(fēng)扇氣動噪聲分析

2.1噪聲分析步驟

在 Fluent 中對于風(fēng)機噪聲的仿真是分為兩個部分先后完成的：

（1） 首先使用大渦模擬模型（LES）對風(fēng)扇流場中的瞬態(tài)控制方程求解獲得流場的動態(tài)穩(wěn)定值，通過計算結(jié)果得到風(fēng)扇的噪聲源（即風(fēng)扇葉片上的動態(tài)載荷）；

（2） 接下來則是通過求解 FW-H 模型的方法對風(fēng)機載荷進行分析并得到噪聲值。

2.2瞬態(tài)流場仿真邊界條件設(shè)定

聲場仿真過程中由于其 CFD 模型與流場極為相似因此不再另行建立模型，而是對原有流場模型的邊界條件進行修改。由于噪聲特性的仿真屬于非定常計算，雖然同樣將旋轉(zhuǎn)流體域設(shè)為唯一的運動區(qū)域，但是改用滑移網(wǎng)格模型對風(fēng)扇的動葉片與靜止區(qū)域進行耦合以保證瞬態(tài)計算的精確度。

在控制方程的離散過程中使用PISO 算法代替原來的 SIMPLE 方法，相比較而言 PISO 算法在原有“預(yù)測-修正”方法的基礎(chǔ)上添加了一個再修正過程，對原有計算結(jié)果進行了二次改進，有效的提高了計算精度與方程的收斂速度。至于迭代過程中參數(shù)的設(shè)置，將時間步長設(shè)定為0.0001s，而計算的截止頻率取 6000Hz，在每個時間步長內(nèi)計算 40 次，迭代次數(shù)為 1000。

2.3氣動噪聲邊界條件設(shè)定與后處理

將2.2節(jié)中的計算結(jié)果與 FW-H 方程相結(jié)合在葉片表面使用二重積分獲得隨需要的聲壓值信號，在進行噪聲參數(shù)設(shè)定時，以風(fēng)扇本體為噪聲源，而監(jiān)測點則按照 GB/T2888-2008《風(fēng)機和羅茨鼓風(fēng)機噪聲測量方法》中的規(guī)定設(shè)置，取風(fēng)機前 1m 處噪聲結(jié)果作為分析。

最終可以得到1m處噪聲值計算結(jié)果為78dB(A)，完成聲場計算過程后得到的數(shù)據(jù)是時域信號，還需使用Fluent 軟件后處理功能中的快速傅立葉變換模塊（FFT）完成時頻轉(zhuǎn)換獲得聲壓級頻譜圖，如下所示。