鄂加強，王紅林，劉騰，左威，張慶玲，蘇秀超

?

基于恒風(fēng)速虛擬風(fēng)洞實驗的碟式太陽能聚光器風(fēng)動載荷

鄂加強1, 2，王紅林2，劉騰1, 2，左威2，張慶玲2，蘇秀超2

(1. 湖南大學(xué) 汽車車身先進設(shè)計制造國家重點實驗室，湖南 長沙，410082；2. 湖南大學(xué) 機械與運載工程學(xué)院，湖南 長沙，410082)

針對碟式太陽能聚光器最佳避風(fēng)姿勢問題，采用恒風(fēng)速虛擬風(fēng)洞實驗方法對不同姿態(tài)和不同風(fēng)速下的碟式太陽能聚光器風(fēng)動載荷進行仿真分析。研究結(jié)果表明：在靠近來流方向處碟式太陽能聚光器總壓力較大，而遠離來流方向處其總壓力相對較小，且碟式太陽能聚光器反射面的最大壓強隨恒風(fēng)速的增大而不斷增加；當風(fēng)速由=15 m/s增大到=25 m/s時，碟式太陽能聚光器所受風(fēng)力和風(fēng)力矩的最大幅值至少增加了110%。

聚光器；虛擬風(fēng)洞實驗；風(fēng)動載荷；太陽能

太陽能作為可再生能源與傳統(tǒng)能源相比適用范圍廣且效益顯著[1]，并擁有廣闊的發(fā)展前景[2]。對于碟式太陽能聚光器而言，風(fēng)的擾動作用是最大的不確定載荷，而影響碟式太陽能聚光器風(fēng)動載荷最重要部位是反射面。由于反射面迎風(fēng)面積大，空氣阻力系數(shù)也大，因此導(dǎo)致碟式太陽能聚光器在某些狀態(tài)下的風(fēng)載十分巨大。由于導(dǎo)致碟式太陽能聚光器相對于氣流方向變化的不確定性，故氣流對碟式太陽能聚光器的風(fēng)動載荷不是一個恒定的載荷。而且氣流流過聚光器時，會在反射面之后形成脫體漩渦，相當于給聚光器施加了脈動激振力，容易導(dǎo)致碟式太陽能聚光器振動，甚至引起破壞性的共振。因此，對各種工況下的碟式太陽能聚光器受到的風(fēng)動載荷特性研究顯得十分重要。20世紀90年代以來，國內(nèi)外研究工作者關(guān)于碟式太陽能熱發(fā)電系統(tǒng)的研究主要集中在提高吸熱器能量轉(zhuǎn)化效率[3?6]、聚光器聚光效率[7?9]等方面，而沒有針對碟式太陽能熱發(fā)電系統(tǒng)推廣應(yīng)用的風(fēng)激振動對碟式太陽能聚光器風(fēng)載壓強分布與模態(tài)影響機理進行研究。風(fēng)洞試驗是常用的計算結(jié)構(gòu)風(fēng)載荷的方法，成本高、周期長，而數(shù)值仿真方法可以快速準確地得到模型受風(fēng)作用時的模型風(fēng)效應(yīng)情況，避免風(fēng)洞試驗由于忽略細節(jié)而引起的計算誤差而顯示出蓬勃生機[10?12]。國內(nèi)外學(xué)者也提出多種風(fēng)載荷計算方法[13?15]。因此，本文作者采用簡化的物理模型和高精度的數(shù)值算法，在恒風(fēng)速虛擬風(fēng)洞實驗下的碟式太陽能聚光器進行風(fēng)動載荷分析，獲得碟式太陽能聚光器風(fēng)動載荷分布，為改進碟式太陽能聚光器的結(jié)構(gòu)以及對工程實踐中的最佳避風(fēng)姿勢問題提供理論依據(jù)。

1 風(fēng)動載荷仿真模型

1.1 數(shù)學(xué)模型

碟式太陽能聚光器工作在大氣邊界層的風(fēng)場中，由于氣流在大氣邊界層中的流動屬于湍流，為低速、不可壓縮的黏性牛頓流體。故對碟式太陽能聚光器流體區(qū)域進行數(shù)值模擬時，不考慮熱交換，即不考慮能量方程。

1) 質(zhì)量守恒方程為

2) 動量守恒方程為

3) 湍流方程模型為

式中：v，v和v分別為風(fēng)速度在，和方向的分量，m/s；為時間，s；為空氣動力黏度，Pa·s，且=0+T，0為空氣分子黏性系數(shù)，T=μ2/，μ為與?方程有關(guān)的常數(shù)，為湍動能，為湍動能耗散率；為空氣密度，kg/m3；為絕對壓強，Pa；k為湍流動能的有效Prandtl數(shù)；K=T(?u/?x+?u/?x) (?u/?x)。

式中：ε為湍流動能耗散率的有效Prandtl數(shù)；1和2為與?方程有關(guān)的常數(shù)。

數(shù)學(xué)模型中的有關(guān)常數(shù)如表1所示。

表1 數(shù)學(xué)模型常數(shù)

1.2 力學(xué)模型

風(fēng)力與空氣密度、風(fēng)速、空氣的黏性、碟式太陽能聚光器的直徑等有關(guān)，故風(fēng)力可表示為

用因式分解法可得：

因為碟式太陽能聚光器的2正比于碟式太陽能聚光器的特征面積(即2∝)，且雷諾數(shù)=/，故式(8)可表示為

式中：F為風(fēng)力系數(shù)，F(xiàn)=()。

對于碟式太陽能聚光器，當風(fēng)向不同時，風(fēng)力也不同。為計算方便，無論風(fēng)向如何，式(9)中的取相同值，例如取高度角=0o、方位角=0o時的最大面積，而僅僅是風(fēng)力系數(shù)隨風(fēng)向而變。

圖1所示為碟式太陽能聚光器受風(fēng)載荷作用后的受力圖。一般情況下，碟式太陽能聚光器所受風(fēng)力可分解為側(cè)力F、阻力F和升力F：

F=0.52AC(10)

F=0.52AC(11)

F=0.52AC(12)

式中：C為側(cè)力系數(shù)；C為阻力系數(shù)；C為升力系數(shù)。

圖1 碟式太陽能聚光器受力圖

此外，碟式太陽能聚光器受風(fēng)力作用時還會產(chǎn)生風(fēng)力矩。類似地，碟式太陽能聚光器所受風(fēng)力矩可表示為

式中：C為風(fēng)力矩系數(shù)；為碟式太陽能聚光器半徑，m。

當?shù)教柲芫酃馄魇苎剌S風(fēng)向作用時，風(fēng)力矩在一般情況下可分解為俯仰力矩、翻滾力矩和方位力矩：

=0.5Cρv2(14)

=0.5Cρv2(15)

=0.5Cρv2(16)

1.3 物理模型

碟式太陽能聚光器物理模型如圖2所示。當風(fēng)對碟式太陽能聚光器進行繞流時，碟式太陽能聚光器相當于處于完全開口的流動風(fēng)場中。而風(fēng)對碟式太陽能聚光器的作用存在于一定的范圍內(nèi)，因此在建立碟式太陽能聚光器計算流體域時給定與實際相對應(yīng)的邊界條件。碟式太陽能發(fā)電系統(tǒng)聚光器直徑約為18 m，厚度為0.28 m，為保證氣流在流體域內(nèi)的流動達到充分發(fā)展的狀態(tài)，流體域的長度約為聚光器的10倍，流體域的寬度和高度約為碟式太陽能發(fā)電系統(tǒng)聚光器的5倍，經(jīng)過多次建模仿真，確定流體域長×寬×高為 170 m×80 m×80 m，模型中心距地面高度為10 m，風(fēng)的入口距模型中心距離為55 m，如圖3所示。

圖2 碟式太陽能聚光器物理模型

圖3 碟式太陽能聚光器計算域模型

對模型進行設(shè)定，將空氣看作連續(xù)介質(zhì)、無壓縮黏性流體。網(wǎng)格的類型和網(wǎng)格的數(shù)量對計算結(jié)果的準確性有很大的影響。對該模型的流體區(qū)域采用的是多面體網(wǎng)格，非結(jié)構(gòu)網(wǎng)格具有很好的靈活性和適應(yīng)性，易于進行網(wǎng)格自適應(yīng)。為了更準確地模擬聚光器周圍流場特性，對聚光器周圍區(qū)域有遠及近逐步加密。

1.4 仿真條件設(shè)定

設(shè)置邊界條件前需對模型的流體域進行設(shè)置，流體域中所選流體為空氣，氣溫為常溫且絕熱，參考壓強為100 kPa，且采用3種以下方式確定模擬風(fēng)洞實驗初始條件和邊界條件：

1) 入口初始條件和邊界條件：該區(qū)域的流體為不可壓縮流體流動，入口初始風(fēng)速分為恒定風(fēng)速，該區(qū)域的流體流動為不可壓縮流體流動，入口風(fēng)速分別為15，20和25 m/s。

2) 出口邊界條件：出口處的邊界條件采用壓強出口邊界條件，壓強為100 kPa。

3) 壁面條件：流體域表面和聚光器表面，wall是一種限定流體域和固體域的邊界條件。壁面條件的粗糙度設(shè)置為光滑。對于黏性流體，采用黏附條件，即認為壁面處流體速度與壁面該處的速度相同，無滑移壁面的速度為0 m/s，壁面處流體速度為0 m/s。聚光器表面和地面固定不動，不發(fā)生移動，所以采用無滑移的壁面條件；而流體域的頂面和前后面采用滑移邊界條件。

1.5 控制方程離散化與求解

采用有限差分法對碟式太陽能聚光器風(fēng)動載荷模型的守恒方程組進行離散化。根據(jù)以上確定的初始條件和邊界條件，采用SIMPLEST思想對動量方程進行離散化，對壓力方程采用全域求解，其他守恒方程采用線性迭代法求解。

2 仿真結(jié)果與分析

2.1 流速分布

當=0°和=45°時，流體域速度分布如圖4所示。由圖4(a)可見：當風(fēng)速為15 m/s時，在碟式太陽能聚光器后表面會形成一個大的回流區(qū)，將使碟式太陽能聚光器后表面所受風(fēng)力與前表面所受部分風(fēng)力抵消。另外，在風(fēng)通過碟式太陽能聚光器后表面縫隙時，流出速度有所上升，且風(fēng)在碟式太陽能聚光器上下端面的速度方向沿中心向外，且下端速度大于上端速度。