摘 要：本文針對(duì)典型濱海核電廠址開(kāi)展小尺度范圍內(nèi)風(fēng)洞模擬實(shí)驗(yàn)研究，研究了小尺度范圍內(nèi)建筑物群及地形對(duì)局地流場(chǎng)及擴(kuò)散的影響。風(fēng)洞實(shí)驗(yàn)?zāi)Ｐ涂s比為1 ∶ 600，通過(guò)調(diào)整風(fēng)廓線、湍強(qiáng)廓線及排放參數(shù)保證了物理模型與原型的相似。實(shí)驗(yàn)研究了在年均風(fēng)速及小風(fēng)條件下建筑物和地形對(duì)流場(chǎng)的擾動(dòng)，研究結(jié)果表明，在WNW ～ESE 風(fēng)向軸線上，由于建筑物影響引起速度虧損和湍流強(qiáng)度增大，選擇建筑物周圍不同方位的測(cè)量點(diǎn)位置進(jìn)行測(cè)量，可見(jiàn)不同的速度虧損與湍流強(qiáng)度變化。以煙囪和地面不同釋放高度進(jìn)行了年均風(fēng)和小風(fēng)條件的示蹤實(shí)驗(yàn)，年均風(fēng)條件下，煙囪排放時(shí)ESE 風(fēng)向軸向最大彌散因子為4. 05×10-5 s/ m3 ，軸向最大彌散因子的出現(xiàn)距離在240 m;近地面排放時(shí)軸向最大彌散因子為4. 79×10-3 s/ m3 ，軸向最大彌散因子的出現(xiàn)距離在45 m。

關(guān)鍵詞：風(fēng)洞;小尺度;建筑物;小風(fēng)條件

中圖分類號(hào)：X169 文獻(xiàn)標(biāo)志碼：A

0 引言

放射性物質(zhì)的非正常釋放一直是環(huán)境風(fēng)險(xiǎn)和應(yīng)急響應(yīng)研究的重點(diǎn)，大氣擴(kuò)散過(guò)程是核事故中放射性核素和污染物向周圍環(huán)境傳輸轉(zhuǎn)移的主要途徑之一，近距離的擴(kuò)散會(huì)受到周圍建筑物的影響，因此放射性核素的近距離擴(kuò)散模擬一直是人們關(guān)注的重點(diǎn)。目前精確描述建筑物構(gòu)型和局地特征的CFD 模型和簡(jiǎn)單處理建筑物的大氣擴(kuò)散模型[1] 均可用于模擬近距離復(fù)雜流場(chǎng)，但對(duì)于精準(zhǔn)估算污染物在復(fù)雜地形和建筑物周圍的擴(kuò)散還很困難。大氣擴(kuò)散模型的有效性需要實(shí)驗(yàn)數(shù)據(jù)進(jìn)行驗(yàn)證，此類實(shí)驗(yàn)數(shù)據(jù)主要來(lái)自野外大氣示蹤實(shí)驗(yàn)和大氣邊界層風(fēng)洞模擬實(shí)驗(yàn)。大氣邊界層風(fēng)洞模擬實(shí)驗(yàn)是嚴(yán)格按照相似準(zhǔn)則[2] 、在縮小的物理模型上再現(xiàn)實(shí)際物理過(guò)程的物理模擬方法，相對(duì)其它研究手段具有實(shí)驗(yàn)條件可控、可重復(fù)、易于再現(xiàn)、不受實(shí)際天氣條件限制的優(yōu)點(diǎn)[3] ，可進(jìn)行大量基礎(chǔ)性研究和機(jī)制性探索，因此風(fēng)洞模擬實(shí)驗(yàn)在研究復(fù)雜地形和建筑物近場(chǎng)流動(dòng)與擴(kuò)散問(wèn)題方面起著非常重要的作用[4] 。LIU 等[5] 使用PIV（粒子圖像測(cè)速技術(shù)） 對(duì)低層建筑周圍的流動(dòng)進(jìn)行了實(shí)驗(yàn)研究。Finnigan 等[6] 研究了單個(gè)山丘對(duì)大氣流動(dòng)和污染物擴(kuò)散的影響。Gupta 等[7] 使用風(fēng)洞實(shí)驗(yàn)量化了屋頂結(jié)構(gòu)對(duì)羽流擴(kuò)散的下洗效應(yīng)。Zahid等[8] 基于風(fēng)洞實(shí)驗(yàn)及數(shù)值模擬研究了通道對(duì)高層建筑風(fēng)場(chǎng)特性的影響。郭棟鵬等[9] 利用風(fēng)洞實(shí)驗(yàn)并結(jié)合 CFD 數(shù)值模擬技術(shù)模擬了復(fù)雜建筑物群對(duì)周圍流動(dòng)與污染物擴(kuò)散的影響。喬清黨等[10]采用風(fēng)洞實(shí)驗(yàn)和數(shù)值模擬研究了立方體建筑物頂部污染源所排放污染物的流動(dòng)與擴(kuò)散，這些研究都是基于簡(jiǎn)單建筑物模型進(jìn)行的。本文在上述研究的基礎(chǔ)上，以具有復(fù)雜建筑物及地形的實(shí)際核電廠址為研究對(duì)象，建立一定縮比的物理模型，研究污染物在實(shí)際廠址中的擴(kuò)散規(guī)律。

本文選擇的是典型濱海核電廠址，研究近距離建筑物與地形影響下的流動(dòng)及擴(kuò)散特征。模擬區(qū)域半徑2 km，模擬風(fēng)速包括年均風(fēng)速與小風(fēng)條件，模擬實(shí)驗(yàn)的風(fēng)向結(jié)合反應(yīng)堆廠房、汽輪機(jī)廠房等建筑物組合體的布局確定，模擬釋放高度包括煙囪釋放與地面釋放，模擬實(shí)驗(yàn)內(nèi)容包括垂直及水平流場(chǎng)分布， 空中、地面軸線及地面場(chǎng)濃度分布。

1 實(shí)驗(yàn)方法

1. 1 廠址特征

實(shí)驗(yàn)選取的典型核電廠址為濱海廠址，北、東、南三面環(huán)海，西面與山體相連，廠址周圍的山體與核電廠的反應(yīng)堆安全殼、汽輪機(jī)廠房等建筑物組合體的尾流將影響污染物在廠區(qū)和鄰近地區(qū)的遷移和擴(kuò)散。

該廠址包括6 臺(tái)反應(yīng)堆機(jī)組，屬于典型多堆廠址。根據(jù)廠址工程布局，不同風(fēng)向時(shí)各臺(tái)機(jī)組建筑物之間會(huì)相互影響，引起氣流的空氣動(dòng)力學(xué)形變，從而影響污染物在廠區(qū)的流動(dòng)與擴(kuò)散。各臺(tái)機(jī)組反應(yīng)堆廠房高度70 m，煙囪高度75 m。煙囪的高度略高于反應(yīng)堆廠房的高度。煙囪內(nèi)徑1. 82 m×1. 98 m。廠址地形如圖1 所示。

1. 2 實(shí)驗(yàn)設(shè)備及實(shí)驗(yàn)條件

本實(shí)驗(yàn)在中國(guó)輻射防護(hù)研究院1 號(hào)環(huán)境風(fēng)洞實(shí)驗(yàn)室完成，風(fēng)洞洞體全長(zhǎng)36 m，實(shí)驗(yàn)段長(zhǎng)17 m、寬1. 5 m、高1～ 1. 4 m，實(shí)驗(yàn)段風(fēng)速范圍為0. 2 ～20 m / s。本實(shí)驗(yàn)使用的主要測(cè)量?jī)x器和裝置包括流場(chǎng)探測(cè)儀和一氧化碳示蹤實(shí)驗(yàn)測(cè)量系統(tǒng)，其中流場(chǎng)探測(cè)儀使用丹麥DANTEC 公司生產(chǎn)的Stream-Line 熱線風(fēng)速儀（ 附有自動(dòng)探頭標(biāo)定裝置）、一維和二維熱線探頭。實(shí)驗(yàn)中利用熱線風(fēng)速儀二維（X 型） 探頭進(jìn)行速度測(cè)量，整個(gè)熱線測(cè)量過(guò)程中采樣頻率為200 Hz，平均采樣時(shí)間為5 s;一氧化碳示蹤實(shí)驗(yàn)測(cè)量系統(tǒng)采用一氧化碳（CO）為示蹤劑，整個(gè)測(cè)量系統(tǒng)包括：CO 釋放裝置、樣品采集系統(tǒng)、兩臺(tái)紅外氣體分析儀。采用CO 氣體作為示蹤劑，兩臺(tái)分析儀器配合使用，1 臺(tái)用于常量分析（量程為0 ～ 10%），1 臺(tái)用于微量分析（量程為0～2 000 ppm），可得到從微量濃度到高濃度的脈動(dòng)測(cè)量結(jié)果。

樣品采集系統(tǒng)為多點(diǎn)采樣（30 點(diǎn)），采集時(shí)間為5 分鐘，采集結(jié)束后立即用紅外氣體分析儀進(jìn)行測(cè)量，測(cè)量濃度為5 分鐘的平均濃度。

1. 3 相似準(zhǔn)則的實(shí)用處理

（1）幾何相似。以二期工程3# 機(jī)組煙囪為中心，選擇WNW ～ ESE 風(fēng)向進(jìn)行模型制作，山體和反應(yīng)堆建筑物群體的實(shí)驗(yàn)?zāi)Ｐ涂s比為1 ∶ 600，模擬區(qū)域半徑2 km。模型表面具有適當(dāng)粗糙度，維持邊界條件和固壁粗糙條件相似。風(fēng)洞模型示意圖如圖2 所示。

（2）運(yùn)動(dòng)學(xué)相似。根據(jù)運(yùn)動(dòng)學(xué)相似準(zhǔn)則，以該核電廠址實(shí)際氣象觀測(cè)資料為依據(jù)建立風(fēng)洞實(shí)驗(yàn)背景流場(chǎng)結(jié)構(gòu)，取上游平均風(fēng)速廓線、湍流強(qiáng)度及其垂直分布等特征量作為運(yùn)動(dòng)學(xué)相似指標(biāo)。實(shí)驗(yàn)中氣流的速度剖面和湍流強(qiáng)度分布主要由速度車來(lái)控制，并輔以適當(dāng)?shù)牡孛娲植谠?，從而維持運(yùn)動(dòng)學(xué)相似。根據(jù)該廠址氣象觀測(cè)系統(tǒng)數(shù)據(jù)統(tǒng)計(jì)，該廠址D 類穩(wěn)定度的上游平均風(fēng)速廓線指數(shù)為0. 12，圖3 給出了風(fēng)洞模擬試驗(yàn)中來(lái)流的風(fēng)速廓線與湍強(qiáng)廓線。

（3）動(dòng)力學(xué)相似。根據(jù)廠址鐵塔氣象資料，廠址地區(qū)80 m 高度處的年平均風(fēng)速為4. 4 m / s，核電機(jī)組煙囪高度75 m。按照根本茂準(zhǔn)則計(jì)算，小尺度年均風(fēng)速條件下風(fēng)洞中煙囪排放高度的風(fēng)速為0. 49 m / s，小風(fēng)條件下風(fēng)洞中煙囪排放高度的風(fēng)速為0. 32 m / s，年均風(fēng)速條件下模型煙囪的排氣速率為1. 21 m / s，小風(fēng)條件下模型煙囪的排氣速率為1. 20 m / s。

（4）邊界條件相似。模型與原型是幾何相似的，且滿足固壁粗糙條件，調(diào)節(jié)速度車和地面粗糙元使來(lái)流為湍流并滿足μ? ×z0 / υ gt;2. 5，即可滿足雷諾數(shù)自準(zhǔn)的要求。

2 實(shí)驗(yàn)結(jié)果及討論

2. 1 局地流場(chǎng)和湍流特征

核電廠址小尺度模擬范圍內(nèi)主要分布有廠區(qū)建筑物及周邊區(qū)域小部分山體。廠區(qū)建筑物密集，其中反應(yīng)堆廠房高度70 m，其他建筑物高度10 m～50 m 不等。ESE 實(shí)驗(yàn)風(fēng)向下風(fēng)向偏西方位有山體分布，山體最高高度155 m。廠區(qū)復(fù)雜的建筑物布局與山體的存在增加了小尺度范圍內(nèi)流場(chǎng)的復(fù)雜性。氣流流態(tài)受山體和建筑物的幾何形狀和尺寸變化影響，流向與流速發(fā)生改變，湍流強(qiáng)度增大;近距離影響程度大，遠(yuǎn)距離影響程度小，高大建筑物附近的擾動(dòng)影響尤為顯著。

年均風(fēng)速條件下的實(shí)驗(yàn)結(jié)果表明，流場(chǎng)受到廠區(qū)范圍內(nèi)建筑物的影響，風(fēng)速普遍減小至來(lái)流風(fēng)速的30%左右，湍流強(qiáng)度增大到0. 4 左右。如圖4 和圖5 所示，在ESE 風(fēng)向，氣流經(jīng)過(guò)反應(yīng)堆廠房后在下風(fēng)距離130 m 左右速度虧損達(dá)到最大，其速度約為來(lái)流的26%左右，湍流動(dòng)能也達(dá)到最大，湍流強(qiáng)度增大到0. 38。隨后氣流經(jīng)過(guò)一期工程建筑物，擾動(dòng)影響較大，隨著下風(fēng)距離的增加流動(dòng)逐漸接近來(lái)流狀態(tài)。

從歸一化速度與縱向湍流強(qiáng)度可以看出，測(cè)量點(diǎn)位置不同，影響因素與影響程度不同，造成不同的速度虧損與湍流強(qiáng)度。圖6 給出了小尺度風(fēng)洞模擬實(shí)驗(yàn)典型實(shí)驗(yàn)風(fēng)向下，在煙囪下風(fēng)向不同距離處平均速度的觀測(cè)值，和受建筑物和地形高度影響所確定的縱向湍流強(qiáng)度的尾流區(qū)域示意圖，同時(shí)給出了廠址上游受建筑物與地形影響的區(qū)域。

實(shí)驗(yàn)數(shù)據(jù)顯示了該地區(qū)由于地形和建筑物的影響湍流場(chǎng)空間分布不均勻，表現(xiàn)出湍流的非定常效應(yīng)，特別是在回流區(qū)和強(qiáng)湍流區(qū)表現(xiàn)更為明顯。WNW～ ESE 風(fēng)向軸線上，由于流動(dòng)連續(xù)受到一期與二期工程反應(yīng)堆廠房影響，ESE 風(fēng)向還受到部分地形和建筑物擾動(dòng)性影響，相比其他風(fēng)向速度虧損和湍流強(qiáng)度增大明顯。

2. 2 擴(kuò)散特征

為了反映近場(chǎng)建筑物群對(duì)煙囪排放煙羽遷移擴(kuò)散的影響，以煙囪設(shè)計(jì)高度（Hs = 75 m）和近地面（10 m）為示蹤劑釋放點(diǎn)，開(kāi)展了大量的示蹤實(shí)驗(yàn)，包括年均風(fēng)速與小風(fēng)風(fēng)速兩種條件下的實(shí)驗(yàn)。

2. 2. 1 垂向濃度分布

選擇距離釋放點(diǎn)不同距離的代表性點(diǎn)V1 ～V5 進(jìn)行垂向濃度的測(cè)量以反映關(guān)注區(qū)域的濃度分布。代表性點(diǎn)位選擇原則是：山體背風(fēng)坡，以考察山體尾流對(duì)大氣彌散的影響情況;山側(cè)，以考察地形引起的繞流特征;建筑物密集區(qū)域，以考察復(fù)雜建筑物群對(duì)流場(chǎng)和濃度分布的擾動(dòng)。由于有些點(diǎn)位偏離軸線較遠(yuǎn)，故某些點(diǎn)位的垂向濃度近似為零。圖7 給出了ESE 風(fēng)向年均風(fēng)速與小風(fēng)條件下代表性點(diǎn)位的垂向濃度分布。其中，V1、V2、V3、V4、V5 分別代表下風(fēng)向距釋放點(diǎn)軸線300 m、1 000 m、2 000 m、3 600 m 和6 000 m 的5 個(gè)代表性測(cè)量點(diǎn)位。總體上，濃度水平分布近似為高斯型。但是當(dāng)煙羽運(yùn)行到山體和建筑物兩側(cè)時(shí)，氣流紊亂，受其影響濃度水平呈偏態(tài)分布。濃度的垂直分布近似呈高斯分布，但由于空氣動(dòng)力下洗的影響，煙羽軸線明顯下傾。年均風(fēng)速條件下，近地面釋放時(shí)，地面最大濃度高于煙囪釋放將近1個(gè)數(shù)量級(jí)，最大落地點(diǎn)在1 000 m，濃度隨高度的衰減較劇烈。煙囪釋放時(shí)， 最大落地點(diǎn)為3 600 m，濃度隨高度的變化較緩慢。小風(fēng)條件下，兩種釋放高度的濃度廓線變得比較復(fù)雜，近地面釋放的垂直衰減還是強(qiáng)于煙囪釋放，煙囪釋放的垂直擴(kuò)散比年均風(fēng)條件下增強(qiáng)。

2. 2. 2 地面濃度分布

為了反映不同排放高度及不同風(fēng)速條件下流動(dòng)對(duì)彌散情況的影響，開(kāi)展了大量的地面濃度示蹤實(shí)驗(yàn)，包括了兩個(gè)排放高度煙囪排放（c-1. 0H）與近地面10 m 排放（c-10m）。圖8 給出了ESE風(fēng)向不同風(fēng)速和不同排放高度條件下地面軸線大氣彌散因子的變化情況。圖9、圖10 給出了ESE風(fēng)向不同風(fēng)速和不同排放高度條件下大氣彌散因子的分布情況。

（1）地面最大彌散因子及其出現(xiàn)距離

在年均風(fēng)速、煙囪設(shè)計(jì)排放高度時(shí)，煙囪排放軸向最大彌散因子為4. 05×10-5 s/ m3 ，出現(xiàn)距離在240 m，此處為2# 反應(yīng)堆廠房上方;近地面排放軸向最大彌散因子為4. 79×10-3 s/ m3 ，出現(xiàn)距離在45 m。近地面釋放時(shí)受建筑物影響較大，由于釋放點(diǎn)在建筑物密集區(qū)，各建筑物幾何結(jié)構(gòu)不同，引起流場(chǎng)的空氣動(dòng)力學(xué)畸變，尾流區(qū)及建筑物間的罅隙效應(yīng)，使得濃度場(chǎng)的分布與常規(guī)不同，與局地建筑物特征密切相關(guān)。

在小風(fēng)、煙囪設(shè)計(jì)排放高度時(shí)，煙囪排放在ESE 風(fēng)向下的軸向最大彌散因子為3. 00×10-5 s/m3 ，出現(xiàn)位置同樣在2#反應(yīng)堆廠房上方，下風(fēng)距離240 m 處;近地面排放ESE 風(fēng)向的軸向最大彌散因子為2. 08×10-2 s/ m3 ，出現(xiàn)距離為40 m。小風(fēng)條件下，建筑物對(duì)煙囪釋放的影響與年均風(fēng)幾乎相同，但對(duì)近地面釋放的影響比較大，建筑物的分布及具體釋放位置是污染物分布的重要影響因素。

（ 2）煙囪排放高度對(duì)擴(kuò)散結(jié)果的影響分析

有效排放高度的增加，使得地面建筑物及地形的影響減小，建筑物的抬升作用延長(zhǎng)了出現(xiàn)最大濃度的下風(fēng)距離，近距離的地面濃度大幅度降低，有效地降低了最大落地濃度，從而也改變了軸向彌散因子的分布特征。

煙囪排放與近地面排放時(shí)的地面軸線污染物的分布存在一定的差異，主要表現(xiàn)在：煙囪排放時(shí)軸向彌散因子的分布近似為單峰型，排放污染物在225～ 840 m 左右達(dá)到最大，之后隨著距離增加而減小;近地面排放時(shí)軸向彌散因子的分布近似呈連續(xù)下降的趨勢(shì);在近距離，近地面排放時(shí)軸向彌散因子明顯大于煙囪排放時(shí)的彌散因子;隨著下風(fēng)距離的增加，兩個(gè)排放高度的地面軸向彌散因子變化逐漸趨于一致。

3 結(jié)果討論

本項(xiàng)研究以典型濱海核電廠址煙囪為中心開(kāi)展了小尺度、不同風(fēng)速條件下流場(chǎng)與濃度場(chǎng)的大量風(fēng)洞模擬實(shí)驗(yàn)。根據(jù)實(shí)驗(yàn)結(jié)果可對(duì)廠址小尺度范圍的流場(chǎng)及湍流特征有一個(gè)直觀與定量的認(rèn)識(shí)，為模式驗(yàn)證提供了豐富的數(shù)據(jù)集。本次風(fēng)洞實(shí)驗(yàn)研究可以得出以下結(jié)論：

（1）小尺度范圍內(nèi)流動(dòng)主要受到廠區(qū)范圍內(nèi)建筑物的影響，WNW～ ESE 風(fēng)向軸線上，由于流動(dòng)連續(xù)受到一期與二期工程反應(yīng)堆廠房影響，速度虧損和湍流強(qiáng)度增大較其他風(fēng)向明顯。

（2）風(fēng)速與湍流強(qiáng)度受影響情況相似，差異主要表現(xiàn)在：小風(fēng)條件下風(fēng)速相對(duì)于來(lái)流廓線速度虧損較小，湍流強(qiáng)度在兩種風(fēng)速條件下無(wú)顯著差異;整體來(lái)看，年均風(fēng)速條件下尾流影響高度略高一些。

（ 3）在年均風(fēng)速、煙囪設(shè)計(jì)排放高度排放時(shí)，煙囪排放時(shí)軸向最大彌散因子為4. 05×10-5 s/ m3 ，軸向最大彌散因子的出現(xiàn)距離在240 m 之間，此處為2#反應(yīng)堆廠房上方;近地面排放時(shí)軸向最大彌散因子在4. 79×10-3 s/ m3 ，軸向最大彌散因子的出現(xiàn)距離在45 m。

（4）煙囪排放與近地面排放時(shí)的地面軸線污染物的分布存在一定的差異，主要表現(xiàn)在：煙囪排放時(shí)軸向彌散因子近似為單峰型分布，排放污染物在225～ 840 m 左右達(dá)到最大，之后隨著距離增加而減小;近地面排放時(shí)軸向彌散因子的分布近似呈連續(xù)下降的變化趨勢(shì);在近距離近地面排放時(shí)的軸向彌散因子明顯大于煙囪排放時(shí)的軸向彌散因子;隨著下風(fēng)距離的增加，兩個(gè)排放高度的地面軸向彌散因子的變化逐漸趨于一致。

參考文獻(xiàn)：

[ 1 ] 歐陽(yáng)琰， 蔣維楣， 胡非， 等. 城市小區(qū)環(huán)境流場(chǎng)及污染物擴(kuò)散的風(fēng)洞實(shí)驗(yàn)研究[J]. 南京大學(xué)學(xué)報(bào)（自然科學(xué)版），2003（6）： 770-780.

[ 2 ] 陳謨. 論相似準(zhǔn)則、風(fēng)洞尺寸與數(shù)據(jù)精、準(zhǔn)度與關(guān)系[J]. 宇航學(xué)報(bào)， 1997（1）： 41-47.

CHEN Mo. On the relationship of the similarity criterion wind-tunnel scales and data accuracy [ J ]. Journal ofAstronautics， 1997（1）： 41-47.

[ 3 ] 凌坤. 空氣動(dòng)力學(xué)實(shí)驗(yàn)基礎(chǔ)[M]. 北京： 北京科學(xué)教育出版社， 1961： 7-90.

[ 4 ] Graingfr C， Meroney R N. Inversted flow Wind-Tunnel simulation of stably stratified atmospheric boundary layer flow[J].Atmospheric Environment， 1994， 28（11）， 1887-1893.

[ 5 ] LIU Zhixiang， YU Zhixiang， CHEN Xiaoxiao， et al. An investigation on external airflow around low-rise building withvarious roof types： PIV measurements and LES simulations[J]. Building and Environment， 2020， 106583： 106583.

[ 6 ] Finnigan J， Ayotte K， Harman I， et al. Boundary-layer flow over complex topography[J]. Boundary-layer Meteorology，2020， 177（2）： 247-313.

[ 7 ] Gupta A， Stathopoulos T， Saathoff P. Wind tunnel investigation of the downwash effect of a rooftop structure on plumedispersion[J]. Atmospheric Environment， 2012， 46： 496-507.

[ 8 ] Zahid Iqbal Q M， CHAN Al S. Pedestrian level wind environment assessment around group of high-rise cross-shapedbuildings： Effect of building shape， separation and orientation[J]. Building and Environment， 2016， 101： 45-63.

[ 9 ] 郭棟鵬， 姚仁太， 喬清黨. 建筑物對(duì)污染物擴(kuò)散影響的數(shù)值與風(fēng)洞模擬研究[J]. 實(shí)驗(yàn)流體力學(xué)， 2010， 24（6）： 16-21， 26.

[10] 喬清黨， 郭棟鵬， 姚仁太， 等. 建筑物周圍流場(chǎng)與污染物擴(kuò)散的數(shù)值與風(fēng)洞模擬研究[J]. 輻射防護(hù)， 2011， 31（4）： 202-208， 237.

QIAO Qingdang， GUO Dongpeng， YAO Rentai， et al. Numerical and wind-tunnel simulations of flow and dispersionaround a building[J]. Radiation Protection， 2011， 31（4）： 202-208， 237.