王志遠(yuǎn) 李煜 張月婷

摘要：核電廠放射性核素在大氣環(huán)境中的遷移與擴(kuò)散備受廣泛的關(guān)注，而了解核電廠正常工況下釋放到環(huán)境中的放射性核素的遷移規(guī)律，對于核電廠的規(guī)劃以及環(huán)境影響評價具有重要意義。復(fù)雜地形與建筑物近場流動及氣載放射性核素擴(kuò)散問題一直是國內(nèi)外非常關(guān)注的課題。本文主要針對風(fēng)洞實驗和數(shù)值模擬兩種研究方法展開介紹，并對兩種方法的結(jié)果進(jìn)行比較分析。結(jié)果表明，數(shù)值模擬結(jié)果與風(fēng)洞物理模擬結(jié)果表現(xiàn)出較一致的變化趨勢，總體上CFD技術(shù)可以作為復(fù)雜地形和建筑物近場流動和擴(kuò)散研究的一種較為有效的工具。

關(guān)鍵詞：核電廠;大氣擴(kuò)散;風(fēng)洞實驗;數(shù)值模擬

引言

對于核電廠放射性核素在大氣中的擴(kuò)散問題在國際上已有多年的研究，放射性核素大氣擴(kuò)散過程的研究方法主要分為數(shù)值模擬、物理模擬和大氣示蹤試驗三種基本手段，實際工作中，通常將三種方法結(jié)合起來，更有效的研究污染物在大氣環(huán)境中的遷移擴(kuò)散規(guī)律。

大氣示蹤實驗作為大氣擴(kuò)散研究最基本的方法之一，能夠通過現(xiàn)場實測獲取最真實的擴(kuò)散情況，在我國《核電廠廠址大氣擴(kuò)散試驗規(guī)范》（NB/T 20202-2013）中有具體要求和規(guī)范[1]，然而僅依靠事故后的觀測結(jié)果會存在覆蓋范圍有限、無法重現(xiàn)、無法提前進(jìn)行安全預(yù)案等問題[2]。風(fēng)洞實驗相較于示蹤實驗可以更加準(zhǔn)確、穩(wěn)定且可重復(fù)地模擬大氣情況，實驗花費也較少，是一種已經(jīng)被證實實驗室環(huán)境下有效的擴(kuò)散分析方法，實驗結(jié)果具有更好的重復(fù)性，從而可以更好地模擬核素擴(kuò)散的情況[3]。數(shù)值模擬技術(shù)，即計算流體力學(xué)技術(shù)（Computational Fluid Dynamics），是研究放射性核素遷移擴(kuò)散的一種有效手段[4]，應(yīng)用數(shù)值模式模擬能快速有效地評估放射性物質(zhì)泄漏的污染范圍和區(qū)域。我國有眾多用于國防、科研以及能源發(fā)電的核設(shè)施，模擬放射性核素的擴(kuò)散。相對于實驗研究，具有其獨特的優(yōu)點：（1）研究成本低、周期短;（2）無實驗儀器干擾;（3）能獲得完整的數(shù)據(jù);（4）能將計算情況在計算機(jī)屏幕上形象地再現(xiàn)[5];（5）數(shù)值模擬技術(shù)通常不受相似準(zhǔn)則等條件的限制，并且CFD可以給出整個模擬區(qū)域內(nèi)流場結(jié)構(gòu)與污染物的擴(kuò)散規(guī)律[6]。同時CFD模型借助于精確的三維建模和網(wǎng)格劃分，可以有效解決有效解決物理模擬和現(xiàn)場試驗存在的不足，并且有國內(nèi)外的研究成果表明，CFD的模擬結(jié)果與風(fēng)洞實驗的結(jié)果吻合較好[7]。

本文主要針對風(fēng)洞實驗和數(shù)值模擬兩種研究方法展開介紹，并對兩種方法的結(jié)果進(jìn)行比較分析。風(fēng)洞實驗以中國輻射防護(hù)研究院2號風(fēng)洞為代表，數(shù)值模擬以流體力學(xué)軟件Fluidyn-PANACHE為代表。

1風(fēng)洞實驗介紹

中國輻射防護(hù)研究院大氣邊界層2號風(fēng)洞實驗室為直流下吹氣風(fēng)洞，洞體全長71.1m，

試驗段截面為5m×3.5m，長28m，如圖1.1所示。

根據(jù)相似理論，兩個流動系統(tǒng)相似必須滿足一套確定的相似準(zhǔn)則，除要求模型與原型之間實現(xiàn)幾何相似外，還要求運動相似和動力相似以及邊界條件相似。實驗中，同樣滿足幾何相似、運動相似、邊界條件相似，其中運動相似的平均速度廓線遵從指數(shù)律，動力相似符合根本茂準(zhǔn)則。

流場測量裝置是丹麥DANTEC公司生產(chǎn)的Multichannel CTA多通道同步熱線風(fēng)速儀，二維熱線探頭，測量速度范圍：0.2m/s～300m/s，精確度：±1.5%或±0.02m/s。濃度測量儀器包括自動采樣器和氣象色譜儀。

2數(shù)值模擬方法介紹

Fluidyn-PANACHE軟件是由法國Transoft公司與法國環(huán)境與能源署（ADEME）協(xié)作開發(fā)，并結(jié)合現(xiàn)代軟件工程技術(shù)，通過三維有限體積法（FVM）求解模擬空氣運動的Navier-Stokes方程。PANACHE包括一個內(nèi)置的自動三維網(wǎng)格生成器，可以創(chuàng)建結(jié)構(gòu)化（矩形）或非結(jié)構(gòu)化（三角形）網(wǎng)格，網(wǎng)格圍繞障礙物或起伏不定地形進(jìn)行貼體計算，如圖2.1所示，同時考慮所有地形地貌的影響，諸如地形起伏、植被冠層和城市冠層等。

Fluidyn-PANACHE的基本求解方程為N-S方程，同時求解物種的濃度、質(zhì)量和能量守恒方程。對于顆粒及氣溶膠的擴(kuò)散，還要求解拉格朗日描述下的顆粒運動軌跡方程。對于湍流模式，PANACHE求解上述方程的雷諾平均形式。雷諾應(yīng)力模型使用線性渦黏性模型（LEVM）方程：

（4）能量守恒

式中，τ為黏性應(yīng)力張量;為變形率張量;ρ為密度;U為速度矢量;ym為物種m的質(zhì)量分?jǐn)?shù);T為溫度;p為壓力;Dm為物種m的有效擴(kuò)散系數(shù);μ為有效黏性系數(shù);k為膨脹黏性系數(shù)（k=0為Stokes流體）;δ為單位張量;Cp為定壓比熱;q為熱通量向量，q=-k▽T;k為有效熱傳導(dǎo)系數(shù);Sm為物種m濃度守恒方程的源項;Sρ為連續(xù)性方程的源項;SU為動量守恒方程的源項;ST為能量守恒方程的源項。

PANACHE使用修正浮力和可壓縮的標(biāo)準(zhǔn)高雷諾數(shù)格式[8]，其求解湍流動能k的輸運方程及其耗散率ε方程為：

式中，SεV為植被冠層產(chǎn)生的源項。

3風(fēng)洞實驗與數(shù)值模擬結(jié)果比較

本次工作主要指針對某核電廠址，同步開展風(fēng)洞實驗和數(shù)值模擬研究，并對兩種方法的流場及地面濃度分布情況進(jìn)行比較。該廠址5km范圍內(nèi)有起伏地形及高大建筑物分布。

風(fēng)洞實驗中，進(jìn)口的氣流條件為：（1）近地層平均風(fēng)廓線冪指數(shù)P=0.16，模型頂部1200mm處的風(fēng)速U=1.0m/s;（2）實驗?zāi)Ｐ涂s比為1：1000。數(shù)值模擬計算采用非結(jié)構(gòu)化網(wǎng)格和有限體積法來求解N-S方程組，模擬湍流流動使用k-ε模型，開展基于實際地形的數(shù)值模擬，并將模擬結(jié)果與風(fēng)洞實驗結(jié)果進(jìn)行比較。計算時邊界層入口邊界條件風(fēng)廓線、湍流廓線、模型頂部風(fēng)速以及污染源排放條件均與風(fēng)洞實驗相同。

3.1流場特征

根據(jù)該廠址氣象資料分析結(jié)果，并結(jié)合實際地形特征，選取該廠址年均主導(dǎo)風(fēng)向下，取煙囪所在軸線上不同距離處的風(fēng)廓線結(jié)果，比較風(fēng)洞實驗和數(shù)值模擬結(jié)果各關(guān)注點歸一化速度隨高度的變化情況，如圖3.1所示。

由圖可以看出，基本上數(shù)值模擬結(jié)果與風(fēng)洞實驗的歸一化速度隨高度的變化趨勢吻合較好?？傮w上，風(fēng)洞實驗的結(jié)果略大于CFD結(jié)果。在近地面處，數(shù)值模擬的結(jié)果略高于風(fēng)洞實驗結(jié)果，考慮受到該廠址廠區(qū)范圍內(nèi)近地面高大建筑物及起伏地形的影響，同時本次數(shù)值模擬所采用的k-ε湍流模型低估了湍流動能。

3.2濃度場特征

分別取下風(fēng)向500m、1000m、2000m和3000m處地面橫向濃度，比較了風(fēng)洞實驗結(jié)果與數(shù)值模擬計算結(jié)果的地面濃度分布情況，如圖3.2所示。

由上圖可知，總體上，數(shù)值模擬結(jié)果與風(fēng)洞實驗結(jié)果能表現(xiàn)出一致的變化趨勢，且風(fēng)洞實驗與數(shù)值模擬表現(xiàn)出的地面濃度分布情況吻合較好。

4結(jié)論

總體上，風(fēng)洞實驗結(jié)果高于數(shù)值模擬結(jié)果，說明本次數(shù)值模擬所使用的k-ε模型低估了湍流動能，使得計算結(jié)果的濃度散步較小，垂向風(fēng)廓變化趨勢較陡。流場驗證結(jié)果表明，本次案例所使用的CFD技術(shù)的模擬結(jié)果與風(fēng)洞實驗結(jié)果所表現(xiàn)出的趨勢基本一致，但對于近地面由于復(fù)雜地形引起的空氣動力學(xué)畸變特點未能表現(xiàn)出來。濃度場驗證結(jié)果表明，本次案例所使用的CFD技術(shù)能夠較好地模擬出污染物的輸送與擴(kuò)散規(guī)律，其結(jié)果與風(fēng)洞實驗結(jié)果的變化趨勢一致，但近地面兩者結(jié)果還存在一定的差異;但對于地面濃度的分布情況與風(fēng)洞實驗?zāi)茌^好的吻合。

通過比較數(shù)值模擬及風(fēng)洞實驗的流場特征及濃度分布特征，發(fā)現(xiàn)數(shù)值模擬結(jié)果與風(fēng)洞實驗結(jié)果表現(xiàn)出較一致的變化趨勢，總體上CFD技術(shù)可以作為復(fù)雜地形和建筑物近場流動和擴(kuò)散研究的一種較為有效的工具，但對于近地面，兩者存在一定的差異性。

參考文獻(xiàn)：

[1]《核電廠廠址大氣擴(kuò)散試驗規(guī)范》（NB/T 20202-2013）.

[2]邰揚，申世飛等，基于Lagrangian-Eulerian耦合的放射性物質(zhì)大氣擴(kuò)散模型研究，2021，43（1）：1-8.

[3]喬清黨，姚仁太，郭占杰，等.內(nèi)陸核電廠冷卻塔對周圍大氣流動和污染物擴(kuò)散影響的風(fēng)洞實驗研究[J].輻射防護(hù)，2011，13（3）：141-149.

[4]葛寶珠，陸芊芊，陳學(xué)舜等.放射性核素大氣擴(kuò)散數(shù)值模擬研究綜述[J].環(huán)境科學(xué)學(xué)報，2021，41（5）：1599-1609.

[5]王博，郭瑞萍，張瓊等.內(nèi)陸核電廠冷卻塔對大氣擴(kuò)散影響的CFD模擬[J].科技導(dǎo)報，2013，31（32）：34-41.

[6]郭棟鵬，閆函，姚仁太，等.復(fù)雜建筑物對近場擴(kuò)散影響的數(shù)值與風(fēng)洞模擬的比較分析[J].實驗流體力學(xué)，2013，27（2）：56-62.

[7]Gromke C，Buccolieri R，Sabatino S D，et al.Dispersion study in a streetcanyon with tree planting by means of wind tunnel and numericalinvestigations-Evaluation of CFD data with experimental data [J].Atmospheric Environment，2008，42（37）：8640-8650.

[8]Ferziger J H，Peric M.Computational methods for fluid dynamics[M].Berlin：Springer，2002.

作者簡介：王志遠(yuǎn)（1970-），男，本科，畢業(yè)于西北建筑工程學(xué)院（現(xiàn)為長安大學(xué)）供熱通風(fēng)空調(diào)專業(yè)。兼職設(shè)總及商務(wù)經(jīng)理，研究方向為全面設(shè)計管理、技經(jīng)管理和商務(wù)管理。