山東師大新校區(qū)綜合教學(xué)樓風(fēng)環(huán)境模擬研究VIP

山東師大新校區(qū)綜合教學(xué)樓風(fēng)環(huán)境模擬研究摘要：基于fluent計(jì)算流體力學(xué)平臺(tái)，采用standardk-模型simple算法對(duì)山東師大新校區(qū)綜合教學(xué)樓周?chē)娘L(fēng)環(huán)境進(jìn)行了數(shù)值模擬，結(jié)果表明：當(dāng)計(jì)算域尺度從15h10h6h調(diào)整到15w10d6h（其中w為建筑物長(zhǎng)度，d為建筑物寬度，h為建筑物高度）時(shí)，計(jì)算精度明顯提高。另外當(dāng)入流風(fēng)方向?yàn)槲鞅憋L(fēng)，來(lái)流風(fēng)速為1.6m時(shí)，距建筑物上風(fēng)向200m位置受建筑物影響西北風(fēng)變?yōu)楸憋L(fēng)，在教學(xué)樓區(qū)域形成多湍窩及高風(fēng)區(qū)的復(fù)雜風(fēng)場(chǎng)，在距建筑物下風(fēng)向500600m位置，建筑物對(duì)風(fēng)的影響基本消除。關(guān)鍵詞：小區(qū)風(fēng)場(chǎng)；數(shù)值模擬；計(jì)算流體力學(xué)；建筑繞流中圖分類(lèi)號(hào)：P461.1文獻(xiàn)標(biāo)識(shí)碼：A文章編號(hào)：16749944（2013）020177031引言隨著社會(huì)、經(jīng)濟(jì)的高速發(fā)展，城市的規(guī)模日益擴(kuò)大，城市中高大建筑的數(shù)量和高度也與日俱增，這些建筑顯著改變了城市的風(fēng)環(huán)境。一方面高大密集的建筑群，降低了城市的通風(fēng)、自凈能力，加劇了在低風(fēng)速條件下城市的空氣污染和熱島效應(yīng)；而另一方面在風(fēng)速較大時(shí)，高大建筑周?chē)鷷?huì)產(chǎn)生局地強(qiáng)風(fēng)，影響到行人的舒適與安全，引出行人風(fēng)環(huán)境（PedestrianLevelWindEnvironment）問(wèn)題1，2。對(duì)于建筑物周?chē)L(fēng)場(chǎng)情況的獲取通常來(lái)自三種方法：現(xiàn)場(chǎng)監(jiān)測(cè)、物理風(fēng)洞模擬和利用計(jì)算機(jī)進(jìn)行數(shù)值模擬3?，F(xiàn)場(chǎng)監(jiān)測(cè)費(fèi)時(shí)費(fèi)力，物理風(fēng)洞花費(fèi)較高而且我國(guó)目前物理風(fēng)洞實(shí)驗(yàn)室較少，而隨著計(jì)算機(jī)技術(shù)的不斷發(fā)展，通過(guò)數(shù)值模擬了解建筑物周?chē)L(fēng)場(chǎng)成為趨勢(shì)。山東師大新校區(qū)綜合教學(xué)樓是新校區(qū)的主要教學(xué)樓，教學(xué)樓風(fēng)環(huán)境跟同學(xué)們的日常學(xué)習(xí)和生活息息相關(guān)。建筑主體分為A、B、C、D、E5個(gè)區(qū)，每個(gè)區(qū)域通過(guò)3層以上的走廊連接，相互連通。教學(xué)區(qū)為5層建筑，樓高約20m。教學(xué)樓處于校園西北部，四周寬闊，周?chē)鸁o(wú)明顯遮擋物，因此運(yùn)用數(shù)值模擬的方法研究建筑物局地風(fēng)環(huán)境有良好的適用性（圖1）。圖1教學(xué)區(qū)樓體示意（東北視角）2模擬方案建筑物周?chē)L(fēng)場(chǎng)的數(shù)值模擬是將建筑物致于流動(dòng)的風(fēng)場(chǎng)中，以流動(dòng)風(fēng)的質(zhì)量守恒方程、動(dòng)量守恒方程、湍流輸運(yùn)方程作為控制方程，采用離散化的數(shù)值方法獲得風(fēng)場(chǎng)中離散風(fēng)壓的解，可以確定出風(fēng)場(chǎng)內(nèi)建筑物表面的風(fēng)壓力值，同時(shí)獲得流暢中其他相關(guān)的物理量狀態(tài)值。本文以fluent軟件包為平臺(tái)，對(duì)綜合教學(xué)樓周?chē)L(fēng)流場(chǎng)進(jìn)行數(shù)值模擬。2.1物理建模教學(xué)樓本身體型不規(guī)則，四周較開(kāi)闊，無(wú)明顯遮擋物。在CAD中按照建筑物原尺寸建立三維物理模型，避免尺寸效應(yīng)可能對(duì)結(jié)果產(chǎn)生的影響，并將CAD中的物理模型導(dǎo)入到fluent6.2的前處理模塊gambit中。2.2確定計(jì)算域位于大氣邊界層中的建筑物，當(dāng)風(fēng)對(duì)其繞流時(shí)，處于一個(gè)完全開(kāi)口的流動(dòng)風(fēng)場(chǎng)中，但風(fēng)對(duì)建筑物作用的影響有一定的范圍，在數(shù)值模擬時(shí)可以給定有限的三維計(jì)算區(qū)域，并確定計(jì)算區(qū)域各邊界的邊界條件，以模擬實(shí)際的流動(dòng)風(fēng)場(chǎng)。計(jì)算流域的確定需要考慮到計(jì)算量，計(jì)算流域越大，網(wǎng)格越多，計(jì)算耗時(shí)越長(zhǎng)，對(duì)計(jì)算機(jī)的配置要求越高。然而計(jì)算區(qū)域過(guò)小則會(huì)導(dǎo)致模擬失真。本文方案一入口距建筑物45h（h為建筑物高度）的長(zhǎng)度，出口距建筑物910h，側(cè)面及頂面距建筑物5h的方案確定計(jì)算區(qū)域，總范圍為450m280m120m；方案二入口距建筑物5w（w為建筑物長(zhǎng)度），出口距建筑物10w；側(cè)面距建筑物5d（d為建筑物寬度），頂面距建筑物5h確定為計(jì)算域，總范圍為2400m880m120m。2.3網(wǎng)格劃分為了避免計(jì)算時(shí)計(jì)算機(jī)出現(xiàn)內(nèi)存不足的錯(cuò)誤，將整個(gè)計(jì)算區(qū)域劃分為4個(gè)區(qū)域，分別進(jìn)行網(wǎng)格劃分。網(wǎng)格的類(lèi)型和數(shù)量對(duì)計(jì)算結(jié)果的準(zhǔn)確性有很大的影響，非結(jié)構(gòu)網(wǎng)格具有良好的靈活性和適應(yīng)性，易于網(wǎng)格的自適應(yīng)，因?yàn)榻ㄖ锏捏w型不規(guī)則，因此采用非結(jié)構(gòu)網(wǎng)格對(duì)其周?chē)饔蜻M(jìn)行網(wǎng)格劃分。方案一得到網(wǎng)格總數(shù)1658929，EquiAngleSkew小于0.83，AspectRadio為3.41小于51，網(wǎng)格質(zhì)量良好以交付fluent計(jì)算；方案二得到網(wǎng)格總數(shù)2285771，EquiAngleSkew小于0.82，AspectRadio為3.91小于51，網(wǎng)格質(zhì)量良好可以交付fluent計(jì)算。2.4邊界條件的選取計(jì)算流域的入口邊界條件采用velocity-inlet（速度入口邊界條件），入口速度1.6m/s，入流風(fēng)向?yàn)槲鞅憋L(fēng)；出口采用outflow（完全發(fā)展出流邊界條件），計(jì)算域側(cè)面采用symmetry（對(duì)稱(chēng)邊界層，定義靠近壁面速度為0），建筑物表面及地面故采用無(wú)滑移壁面條件，即壁面處流體速度為0。2.5湍流模型的選擇k-湍流模型是應(yīng)用最為廣泛的湍流模型，standardk-湍流模型采用布西內(nèi)斯克渦湍粘度假定，由Lauder和Spalding在1972年提出的。目前國(guó)內(nèi)外很多學(xué)者成功地運(yùn)用標(biāo)準(zhǔn)k-兩方程模型模擬建筑風(fēng)場(chǎng)，并取得了滿意的結(jié)果4，5。2.6求解參數(shù)的設(shè)置數(shù)值模擬依據(jù)的控制方程是連續(xù)方程（質(zhì)量守恒方程）和N-S方程（動(dòng)量守恒方程），不考慮熱交換，所以屏蔽能量方程。流體介質(zhì)為空氣具有不可壓縮行，密度為常數(shù)，材料參數(shù)使用默認(rèn)缺省值。2013年2月綠色科技第2期石小倩，等：山東師大新校區(qū)綜合教學(xué)樓風(fēng)環(huán)境模擬研究環(huán)境與安全3模擬結(jié)果與分析將不同模擬方案的模擬結(jié)果與2011年春季在教學(xué)樓附近的風(fēng)場(chǎng)監(jiān)測(cè)實(shí)驗(yàn)數(shù)據(jù)對(duì)比分析，發(fā)現(xiàn)模擬的點(diǎn)位風(fēng)速變化趨勢(shì)基本一致，而方案二的模擬精度要高于方案一。用spss分析軟件分別對(duì)方案一及方案二模擬數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)進(jìn)行配對(duì)T檢驗(yàn)，置信區(qū)間為95%，結(jié)果方案一與實(shí)測(cè)數(shù)據(jù)P值為0.013小于0.05，方案二與實(shí)測(cè)數(shù)據(jù)P值為0.092大于0.05，所以說(shuō)明方案一模擬數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)差異性顯著，方案二模擬數(shù)據(jù)與實(shí)測(cè)數(shù)據(jù)差異性不顯著（圖2）。圖2不同點(diǎn)位模擬風(fēng)速與實(shí)測(cè)風(fēng)速對(duì)比圖3教學(xué)樓區(qū)域2m高風(fēng)場(chǎng)模擬教學(xué)區(qū)2m高風(fēng)場(chǎng)模擬結(jié)果見(jiàn)圖3（為了圖片清晰，截取建筑物周?chē)植浚?。從圖3中可以看出，風(fēng)對(duì)教學(xué)樓繞流，形成復(fù)雜風(fēng)場(chǎng)。在上風(fēng)向距建筑物約200m的位置，風(fēng)向由西北風(fēng)轉(zhuǎn)變?yōu)楸憋L(fēng)，在下風(fēng)向500600m位置，教學(xué)樓對(duì)風(fēng)的影響基本消除。在A區(qū)西側(cè)，由于來(lái)流風(fēng)和A區(qū)南側(cè)的通道由東向西的大風(fēng)，形成逆時(shí)針和順時(shí)針兩個(gè)湍窩，風(fēng)向較為復(fù)雜但風(fēng)速都較小。在E區(qū)的東南角（建筑物尾部），因風(fēng)場(chǎng)突然開(kāi)闊形成大風(fēng)，而且風(fēng)向由北向西變化，并在E、D區(qū)南側(cè)形成一個(gè)大的渦流，導(dǎo)致在C區(qū)西部風(fēng)向變?yōu)槟巷L(fēng)或者西南風(fēng)。受E區(qū)南部湍窩的影響，在D區(qū)與C區(qū)之間的通道與C區(qū)和A區(qū)之間的通道風(fēng)向相反，變?yōu)橛晌飨驏|。教學(xué)樓東側(cè)的大風(fēng)由C區(qū)廣場(chǎng)進(jìn)，在C2區(qū)南部形成順時(shí)針?shù)鰷u，容易導(dǎo)致污染物的聚集。在C區(qū)廣場(chǎng)，由于B區(qū)東側(cè)和B、C區(qū)之間的通道的大風(fēng)相互作用，在B區(qū)南側(cè)形成一個(gè)大的順時(shí)針?shù)鰷u和一個(gè)小的逆時(shí)針?shù)鰷u，形成復(fù)雜風(fēng)場(chǎng)。總之在教學(xué)樓C區(qū)廣場(chǎng)及E、D區(qū)南側(cè)分布著范圍不同的空腔區(qū)和小風(fēng)區(qū)，不利于通風(fēng)，會(huì)導(dǎo)致污染物的聚集；在A區(qū)迎風(fēng)尖角部分、C區(qū)B區(qū)之間的通道以及B區(qū)的東北角均會(huì)形成局部的高風(fēng)區(qū)；這與楊德江6等人研究結(jié)果一致?；?980年Visser關(guān)于室外熱舒適的研究結(jié)果，建筑物周?chē)腥藚^(qū)1.5m處風(fēng)速v4結(jié)語(yǔ)（1）對(duì)于大跨度外形建筑物，建筑物長(zhǎng)、寬大于建筑物高度，不能單純的將計(jì)算域尺度按照高度來(lái)確定（15h10h6h），要按照長(zhǎng)度寬度的尺度來(lái)調(diào)整計(jì)算域（15w10d6h，其中w為建筑物長(zhǎng)度，d為建筑物寬度，h為建筑物高度），調(diào)整之后計(jì)算精度明顯提高。（2）當(dāng)入流風(fēng)方向?yàn)槲鞅憋L(fēng)，來(lái)流風(fēng)向?yàn)?.6m時(shí)，風(fēng)向在上風(fēng)向距教學(xué)樓200m位置處又西北風(fēng)變?yōu)楸憋L(fēng)，在教學(xué)樓區(qū)域風(fēng)向、風(fēng)速均發(fā)生變化，形成復(fù)雜風(fēng)場(chǎng)，在C區(qū)廣場(chǎng)前面形成空腔區(qū)，而在迎風(fēng)面的建筑物尖叫部分會(huì)形成局部大風(fēng)區(qū)，在下風(fēng)向距教學(xué)樓約500600m的位置，教學(xué)樓對(duì)風(fēng)的影響基本消除。參考文獻(xiàn)：1IsyumovN.AStudyofthepedestrianlevelWindsfortheCanaryWharfProjectR.London：EngineeringScienceResearchReport，BLWT-SS16，1986.2BottemaM.AmethodforoptimizationofwinddiscomfortcriteriaJ.BuildEnviron，2000（35）：118.3黃瑩.基于fluent軟件的建筑物風(fēng)場(chǎng)數(shù)值模擬D.武漢：華中科技大學(xué)，2005.4JiangXinbo，LiuZehua，LiuYuanquan.AirenvironmentinalaboratoryanimalroomunderdifferentpressuredifferentialsC/TheFifthInternationalWorkshoponEnergyandEnvironmentofResidentialBuildings，Guilin，ChinaEnvironmentSciencePress，2009.收稿日期：20130109作者簡(jiǎn)介：汪洲瓊（1984），女，四川巴中人，主要從事農(nóng)產(chǎn)品質(zhì)量安全工作。5JiangXinbo，LiuZehua，HarunoriYoshida，etal.PressuredifferentialanalysisofalaboratoryanimalroomC/Maximi