黃 文，何開巖，鐘水庫(kù)，謝 安，曾 威，孫 備

（廣西大學(xué)物理科學(xué)與工程技術(shù)學(xué)院，廣西 南寧 530004）

微藻因其具有極高的光合效率和獨(dú)特的物質(zhì)組成，在食品、飼料、醫(yī)藥、生物固碳和生物燃料等領(lǐng)域有廣泛的應(yīng)用[1-2]。雖然現(xiàn)在已經(jīng)有管式、直筒式和平板式等多種封閉式的光生物反應(yīng)器可供選擇，但是開放式跑道池是最早提出也是目前唯一實(shí)現(xiàn)了商業(yè)化大規(guī)模微藻培養(yǎng)的光生物反應(yīng)器[3-5]。它具有結(jié)構(gòu)簡(jiǎn)單、容易擴(kuò)大規(guī)模和成本較低等優(yōu)點(diǎn)，常被用于培養(yǎng)螺旋藻和小球藻等能忍受高鹽度和極端pH值的藻種。然而開放式跑道池也存在占地面積大、水分蒸發(fā)多、采收成本高、微藻培養(yǎng)密度和單位產(chǎn)量低、CO2和光能利用率低、易受外來(lái)藻種污染等缺點(diǎn)[5]，因此還有待進(jìn)一步的改善。

微藻的生長(zhǎng)受培養(yǎng)液組分、光照、溫度、pH值、溶氧濃度和CO2供給等多種因素的影響，其中光照條件成為了提高光生物反應(yīng)器中微藻生產(chǎn)率的一個(gè)最主要的因素[6]。傳統(tǒng)跑道池由一個(gè)機(jī)械葉輪驅(qū)動(dòng)，使藻液沿著跑道循環(huán)流動(dòng)，但是這樣的循環(huán)方式中死區(qū)的比例較大，能耗較高，湍流也只影響到葉輪附近，其它區(qū)域大都是層流，由于微藻對(duì)光的吸收和散射作用，處于池底層的微藻總是光照不足[7]。針對(duì)以上問(wèn)題，國(guó)內(nèi)外學(xué)者主要從優(yōu)化跑道池的長(zhǎng)寬比例，在跑道池的轉(zhuǎn)角處增加導(dǎo)流板，或采用氣升驅(qū)動(dòng)等方法來(lái)降低能耗和改善流場(chǎng)的均勻性[8-13]。然而這些研究大多是針對(duì)水平方向流場(chǎng)的改善，對(duì)于真正能夠促進(jìn)微藻細(xì)胞實(shí)現(xiàn)光暗循環(huán)的豎直方向上的流場(chǎng)改善則很少報(bào)道。因此，本文作者設(shè)計(jì)了一個(gè)V 形導(dǎo)流槽，在改善光照的同時(shí)也能改善豎直方向的流場(chǎng)分布，混合采用的是比較溫和的氣升方式?？紤]到CFD 方法在設(shè)計(jì)光生物反應(yīng)器中正方興未艾[6]，具有省時(shí)、高效和可視化等優(yōu)點(diǎn)；Fluent 軟件包采用有限體積法，并且具有豐富的湍流模型，在美國(guó)市場(chǎng)的CFD 領(lǐng)域占據(jù)主導(dǎo)地位。因此，采用Fluent6.3 對(duì)該豎直剖面上的速度場(chǎng)分布進(jìn)行了研究，以改善跑道池的混合條件。

1 數(shù)值計(jì)算

1.1 幾何結(jié)構(gòu)與數(shù)學(xué)模型

圖1所示為加導(dǎo)流槽跑道池的三維結(jié)構(gòu)圖，導(dǎo)流槽安裝在跑道池的直線區(qū)，槽內(nèi)充水以平衡液壓。圖2所示為其A-A剖面圖，與圖3 傳統(tǒng)跑道池的剖面圖相比，除了光強(qiáng)隨藻液深度增加而衰減從而分成光區(qū)與暗區(qū)。本研究假設(shè)跑道池的寬度w=2 m，導(dǎo)流槽的寬度w'=1 m，進(jìn)氣邊界寬0.05 m。

傳統(tǒng)跑道池底面為平面，藻液深為0.15～0.25 m，在水平方向藻液的混合由一個(gè)機(jī)械葉輪驅(qū)動(dòng)，為了實(shí)現(xiàn)藻液上下層之間的有效混合，使微藻細(xì)胞相繼經(jīng)歷光區(qū)與暗區(qū)，暫不考慮機(jī)械葉輪影響，僅研究圖2所示剖面由氣升驅(qū)動(dòng)的豎直方向的混合，以θ、d、vair-in為變量，研究這些參數(shù)對(duì)流場(chǎng)分布的影響。

假設(shè)氣液兩相均為不可壓流體，二相流采用Euler 模型。氣泡對(duì)液體的曳力用Schiller-Naumann模型，因?yàn)樵撃Ｐ陀?jì)算穩(wěn)定，適合于形狀為球形的顆粒，且對(duì)水氣二相流也有較多的實(shí)驗(yàn)數(shù)據(jù)。一般藍(lán)藻和綠藻藻液的黏性系數(shù)為0.9～1.2 mPa?s，密度為1000～1030 kg/m3，與水相差不大[14]，因此將水設(shè)置為液相來(lái)代替實(shí)際的藻液，而空氣為氣相。不考慮傳熱傳質(zhì)，湍流計(jì)算用標(biāo)準(zhǔn)k-ε模型，則連續(xù)方程和動(dòng)量守恒方程分別為式（1）、式（2）[8]。

湍流模型的k方程和ε方程分別為式（3）、式（4）[8]。

式中，C1ε、C2ε、σκ和σε分別為湍流參數(shù)，對(duì)于空氣和水，一般取值為1.44、1.92、1.0和1.3。

1.2 計(jì)算方法與邊界條件

圖2所示剖面為一個(gè)二維流域，從而可以簡(jiǎn)化計(jì)算。圖4所示為其網(wǎng)格劃分和邊界條件示意圖，采用四邊形結(jié)構(gòu)化網(wǎng)格劃分方式，網(wǎng)格大小5 mm，采用速度入口和壓力出口邊界條件，其余邊界全為壁面邊界，入口邊界氣體的相含率為100%，出口邊界相對(duì)壓力為0，回流氣體的相含率為100%。邊界流用標(biāo)準(zhǔn)壁面函數(shù)來(lái)擬合，通過(guò)網(wǎng)格無(wú)關(guān)性檢測(cè)，即網(wǎng)格加密一倍后計(jì)算結(jié)果差別小于5%。氣泡直徑5 mm，考慮重力作用，初始液面分布用Region菜單定義，將流域初始化。

采用時(shí)間推進(jìn)的非穩(wěn)態(tài)求解和Simple 迭代方法，動(dòng)量、湍流強(qiáng)度以及耗散率方程為二階迎風(fēng)離散格式，而相分?jǐn)?shù)為一階迎風(fēng)離散，各松弛因子采用默認(rèn)值。將殘差收斂標(biāo)準(zhǔn)統(tǒng)一設(shè)為1.0×10-5，時(shí)間步長(zhǎng)為1 ms，單步迭代50 次，同時(shí)監(jiān)測(cè)進(jìn)出口的流量和兩相的分布，當(dāng)殘差小于1.0×10-5，或者進(jìn)出口流量基本平衡，相分布不再隨時(shí)間變化時(shí)可以認(rèn)為流場(chǎng)已經(jīng)達(dá)到穩(wěn)定，結(jié)束計(jì)算。

2 結(jié)果分析與討論

2.1 與傳統(tǒng)跑道池的對(duì)比

圖5所示傳統(tǒng)跑道池單點(diǎn)曝氣時(shí)豎直剖面的速度場(chǎng)分布圖，其中進(jìn)氣速度為0.1 m/s，跑道池寬w=2 m，藻液深度d=0.2 m。發(fā)現(xiàn)混合主要是在曝氣點(diǎn)附近的流域。

圖6所示為諸發(fā)超等[9]用CFD 模擬傳統(tǒng)跑道池池底通氣情況下，垂直剖面上速度矢量的分布示意圖，其中藻液深度為0.2 m，而藻液的平均速度為0.026 m/s，在跑道池的直線區(qū)等間距布置了4 排的曝氣點(diǎn)，結(jié)果發(fā)現(xiàn)藻液的混合比較均勻。

圖7所示為加V 形導(dǎo)流槽的跑道池豎直剖面上的液體速度場(chǎng)分布圖。從中發(fā)現(xiàn)，氣泡受浮力作用從池底開始上升并夾帶周圍液體往上流動(dòng)。這股上升的氣液混合流體在遇到導(dǎo)流槽的V形尖端阻擋后被對(duì)稱分割成兩股往左右方向流動(dòng)的分流，該分流受導(dǎo)流槽的V 形面所引導(dǎo)，將中心區(qū)上升氣泡流的影響范圍擴(kuò)展到跑道池側(cè)壁區(qū)域，實(shí)現(xiàn)了整個(gè)底層藻液和上層藻液的混合。

與圖5 相比，傳統(tǒng)跑道池池底曝氣只能對(duì)曝氣點(diǎn)附近的流域起到混合作用，遠(yuǎn)離曝氣點(diǎn)的兩側(cè)區(qū)域基本沒(méi)有混合。此外，傳統(tǒng)跑道池中氣泡從池底經(jīng)過(guò)深度為d的藻液后直接進(jìn)入大氣，而加導(dǎo)流槽后，除了液深為d的路程，氣泡還要沿著導(dǎo)流槽的傾斜壁面運(yùn)動(dòng)一段路程才能進(jìn)入大氣，這相當(dāng)于增加了藻液的深度，從而增加了氣泡在藻液中的停留時(shí)間，有利于提高CO2利用效率。

與圖6 相比，在沒(méi)有導(dǎo)流槽的情況下，只有當(dāng)池底曝氣點(diǎn)布置的密度足夠大時(shí)，藻液上下層之間才能夠有效混合，圖6 中跑道池的寬度方向上就布置了4 排的曝氣點(diǎn)，這意味著更大的能量消耗。而增加導(dǎo)流槽后，僅在池中心使用一排的曝氣點(diǎn)就實(shí)現(xiàn)了整個(gè)藻液上下層之間的有效混合，大大減少了能量的消耗。

此外，從圖7 中還可以發(fā)現(xiàn)，曝氣點(diǎn)的兩側(cè)出現(xiàn)了兩個(gè)流線循環(huán)中心，液體圍繞該中心進(jìn)行上下層間的循環(huán)流動(dòng)，使微藻細(xì)胞經(jīng)歷光暗周期。流速較小的區(qū)域出現(xiàn)在循環(huán)中心區(qū)和兩側(cè)池壁附近，而流速較大的區(qū)域出現(xiàn)在池底和導(dǎo)流槽底面附近的流域以及氣泡上升區(qū)。導(dǎo)光槽下表面附近較大的速度分布有利于防止微藻附壁生長(zhǎng)，以保持導(dǎo)流槽良好的透光性。而池底面區(qū)域較大的速度分布則有利于防止微藻沉淀。在實(shí)際養(yǎng)藻過(guò)程中，應(yīng)使光區(qū)的范圍至少達(dá)到流線循環(huán)中心的位置，從而保證所有藻細(xì)胞都有機(jī)會(huì)接收足夠大的光強(qiáng)。

表面積體積比是光生物反應(yīng)器的一個(gè)重要參數(shù)，如式（5）所示，它等于光照面積除以藻液的體積。表面積體積比越大，意味著所能培養(yǎng)的微藻密度越大。

圖8所示為傳統(tǒng)跑道池與增加導(dǎo)流槽后的表面積體積比的比較。下曲面為傳統(tǒng)跑道池的表面積體積比，它與藻液深度呈反比關(guān)系。上曲面為增加導(dǎo)流槽后跑道池的表面積體積比，它隨藻液深度減小而增大，隨傾斜角的增大而增大。不管藻液深度如何，加導(dǎo)流槽后表面積體積比均比傳統(tǒng)跑道池有所增大，但最大值仍小于10。

與傳統(tǒng)跑道池相比，通過(guò)增加導(dǎo)流槽和將平面的池底改造成V形不僅可以改善藻液上下層之間的混合，還增加了反應(yīng)器的表面積體積比。與直接曝氣相比，不僅減少了能量消耗，也提高了CO2在藻液中的停留時(shí)間。

2.2 傾斜率的優(yōu)化

圖9所示為藻液深度為0.3 m、進(jìn)氣速度為0.1 m/s時(shí)速度場(chǎng)分布隨傾斜率的變化圖。其中黑色圓點(diǎn)代表流線循環(huán)中心的位置，相鄰圖形中的圓點(diǎn)用直線連接起來(lái)，由于兩側(cè)是對(duì)稱的，所以只在右側(cè)標(biāo)示，下面的橫坐標(biāo)表示圓點(diǎn)與中軸的距離L。從中可以發(fā)現(xiàn)，隨著傾斜率增加，L先增加后減小然后又增加，在傾斜率為0.1時(shí)L達(dá)到極大值0.6 m，在傾斜率為0.15時(shí)L達(dá)到極小值0.2 m。此外，當(dāng)傾斜率大于0.1時(shí)，圍繞圓點(diǎn)的橢圓形流線發(fā)生了變形，該橢圓的長(zhǎng)軸隨著傾斜率的增加而往上傾斜，側(cè)壁處的藻液深度也逐漸減小，造成了速度分布的不均勻；而當(dāng)傾斜率小于0.1時(shí)，橢圓形的流線雖然變形最小，但是L小于0.5 m，混合影響的范圍沒(méi)有達(dá)到側(cè)壁邊緣，因而側(cè)壁區(qū)域的流速較小。

圖1 加導(dǎo)流槽跑道池的三維結(jié)構(gòu)圖

圖2 加導(dǎo)流槽跑道池的A-A 剖面圖

圖3 傳統(tǒng)跑道池的剖面圖

圖4 網(wǎng)格劃分和邊界條件示意圖

圖5 傳統(tǒng)跑道池單點(diǎn)曝氣的剖面速度場(chǎng)分布圖

圖6 等間距四排曝氣點(diǎn)通氣時(shí)剖面速度矢量分布圖[9]

圖7 加導(dǎo)流槽的跑道池剖面上的速度場(chǎng)分布圖

圖8 傳統(tǒng)跑道池與加導(dǎo)流槽后的表面積體積比比較

圖9 速度場(chǎng)隨傾斜率變化的分布圖

圖10 速度分布均勻度及氣含率隨傾斜率變化的曲線

速度均勻度反映的是整個(gè)速度場(chǎng)分布均勻性的指標(biāo)，也是定量描述整個(gè)藻液混合效果的一個(gè)重要參數(shù)，該值越大，表明混合的效果越好，其表達(dá)式如式（6）所示，v為藻液速率，為平均速率，Sliquid為藻液所在的截面面積。式（7）為氣含率的表達(dá)式，它等于藻液中氣泡的體積Vair除以藻液的表觀體積Vliquid。氣含率越大，氣液傳質(zhì)系數(shù)也越大。因此，增加氣含率對(duì)提高CO2的利用效率就顯得十分重要。

圖10所示為速度均勻度及氣含率隨傾斜率變化的曲線，從中可以發(fā)現(xiàn)，當(dāng)傾斜率小于0.1時(shí)，速度均勻度和氣含率都隨著傾斜率的增加而增大，且都在傾斜率為0.1時(shí)達(dá)到極大值，分別為0.37和0.03；而當(dāng)傾斜率大于0.1時(shí)，隨著傾斜率的增加，速度均勻度和氣含率都是先減小后增加，但是都沒(méi)有超過(guò)之前的極大值。

因此，從速度均勻度和氣含率最大化的目標(biāo)來(lái)考慮，較優(yōu)的傾斜率為tanθ=0.1。并且從圖9(b)可以發(fā)現(xiàn)，與其它傾斜率下的情況相比，此時(shí)整個(gè)池底以及側(cè)壁區(qū)域的速度都比較大，有利于藻液上下層的混合。

為了進(jìn)一步改善流場(chǎng)分布均勻性，對(duì)優(yōu)化參數(shù)下的反應(yīng)器進(jìn)行了倒角處理，其流場(chǎng)分布如圖11所示。對(duì)比圖11和圖9(b)可以發(fā)現(xiàn)，倒角處理后倒角處的地方流速有了明顯提高。速度均勻度也由原來(lái)的0.376 增加到0.379，而平均速度由原來(lái)的0.205 m/s 減小到0.185 m/s，這兩個(gè)指標(biāo)雖然變化不大，但是卻能消除原來(lái)角落處的流動(dòng)死區(qū)，避免了微藻細(xì)胞在此沉淀。因此，倒角處理可以更進(jìn)一改善反應(yīng)器中藻液的混合效果。

2.3 藻液深度的優(yōu)化

圖12所示為傾斜率tanθ=0.1、進(jìn)氣速度為0.1 m/s時(shí)速度場(chǎng)分布隨藻液深度的變化圖。圖中黑色圓點(diǎn)表示流線循環(huán)中心的位置。從中可以發(fā)現(xiàn)，當(dāng)d≤0.25 m時(shí)，黑色圓點(diǎn)與池中心區(qū)的距離L小于0.125 m，此時(shí)圍繞黑色圓點(diǎn)的流線所形成的橢圓變形較大，因而側(cè)壁附近速度較小，混合不均勻。當(dāng)d≥0.3 m時(shí)，黑色圓點(diǎn)與池中心區(qū)的距離L大于0.5 m，橢圓形流線的變形也較小，因而側(cè)壁附近的流速也較大，速度場(chǎng)分布也更加趨向于均勻。此外，隨著藻液深度的增加，整個(gè)流域內(nèi)的流速均有所增大，在池中心區(qū)氣泡垂直上升區(qū)域，最大速度也有所增大，這可能是因?yàn)殡S著藻液深度的增加，氣泡上升的路徑增加，所帶動(dòng)的藻液的速度也因此而有所增加，將流線循環(huán)中心推向離池中心區(qū)更遠(yuǎn)的距離。

圖13所示為液體速度分布均勻度及氣含率隨藻液深度變化的曲線，從中可以看到，速度均勻度隨著d的增加而增大，從液深0.2 m時(shí)的0.11 增加到液深0.35 m時(shí)的0.44，整整提高了3倍。而氣含率則隨著d的增加先減小后增加，最后又減小，在液深0.3 m時(shí)達(dá)到極大值0.03。而在液深0.35 m時(shí)卻只有0.008，這可能是因?yàn)楫?dāng)藻液深度大于0.3 m后，不僅藻液的體積隨著d而直線增加，而且藻液流速的增大也加快了氣泡逃逸，使藻液中的氣泡數(shù)量減小，從而造成了氣含率的減小。

藻液的深度除了影響速度均勻度和氣含率外，對(duì)光強(qiáng)的分布也有重要影響。根據(jù)光衰減的Lambert-Beer 方程，研究人員對(duì)微藻提出了線性模型[15-16]來(lái)擬合，見式（8）。

式中，a、b為實(shí)驗(yàn)確定的常數(shù)；OD為用消光度表示的微藻密度；x為藻液厚度；I0為初始入射光強(qiáng)。當(dāng)微藻密度小于1.7 gDW/L時(shí)，其相關(guān)系數(shù)R2＞0.96，因?yàn)榕艿莱氐奈⒃迕芏炔淮螅蚨刹捎镁€性模型。假設(shè)太陽(yáng)光的散射部分是全天空均勻的，而直射部分隨時(shí)間和天氣而變化，當(dāng)光線不是垂直液面入射時(shí)在側(cè)壁處會(huì)有遮光效應(yīng)，但對(duì)于一天中早上10 點(diǎn)～下午4 點(diǎn)這個(gè)主要生產(chǎn)時(shí)間段，遮光的比例較小，與垂直入射區(qū)別不大。此外，還采用了月平均光強(qiáng)來(lái)代表入射光強(qiáng)I0。

對(duì)每一種微藻都有一個(gè)補(bǔ)償光強(qiáng)I1，當(dāng)外界光強(qiáng)等于I1時(shí)細(xì)胞光合作用生產(chǎn)的物質(zhì)剛好等于自身代謝的消耗。將藻液中I≥I1的區(qū)域定義為光區(qū)，而將I≤I1的區(qū)域定義為暗區(qū)。傳統(tǒng)跑道池中藻液出現(xiàn)上下分層，為滿足所有藻細(xì)胞的光照條件，全部藻液都要處于光區(qū)，則要求：(a+b×OD)d≤ln(I0/I1)，從中可以發(fā)現(xiàn)，微藻密度與藻液深度呈反比關(guān)系，因此，在以上假設(shè)前提下，反應(yīng)器所能承受的最大產(chǎn)量為常數(shù)。而增加導(dǎo)流槽后，由于藻液上下層之間能夠有效混合，允許底層存在一定比例的暗區(qū)，暗區(qū)中的藻細(xì)胞通過(guò)混合進(jìn)入光區(qū)進(jìn)行光合作用。如圖12(c)所示，為保證所有微藻的光照條件，光區(qū)的平均厚度dlight至少要使光區(qū)覆蓋到流線循環(huán)中心的位置，則要求：(a+b×OD)dlight≤ln(I0/I1)。其中dlight＜d，說(shuō)明增加導(dǎo)流槽可以增加微藻的總產(chǎn)量。因此，可以從混合均勻和氣含率來(lái)選擇d=0.3 m 作為設(shè)計(jì)參數(shù)，而微藻的密度則要求滿足上述關(guān)系的限制。

2.4 進(jìn)氣速度的優(yōu)化

圖14所示為傾斜率tanθ=0.1、藻液深度為0.3 m時(shí)的速度場(chǎng)分布隨進(jìn)氣速度的變化。從中可以發(fā)現(xiàn)，隨著進(jìn)氣速度增加，整個(gè)流域的速度都有所提升，當(dāng)進(jìn)氣速度為0.01 m/s時(shí)，最大速度不超過(guò)0.1 m/s；當(dāng)進(jìn)氣速度增加到0.4 m/s時(shí)，最大流速已經(jīng)超過(guò)0.65 m/s，中心區(qū)最大速度隨進(jìn)氣速度的增加而明顯增大。流線循環(huán)中心與跑道池中心區(qū)的距離L隨著進(jìn)氣速度增加先增大后減小，然后又增大，在進(jìn)氣速度為0.1 m/s時(shí)L達(dá)到極大值0.6 m，此時(shí)速度場(chǎng)分布比較均勻。

圖11 跑道池倒角處理后的速度場(chǎng)分布

圖12 速度場(chǎng)隨藻液深度變化的分布圖

圖13 速度分布均勻度及氣含率隨藻液深度變化的曲線

圖14 速度場(chǎng)隨進(jìn)氣速度變化的分布圖

圖15所示為速度均勻度及平均速率隨進(jìn)氣速度變化的曲線。從中可以發(fā)現(xiàn)，當(dāng)vair?in＜0.1 m/s時(shí)，隨著vair?in增大，速度均勻度緩慢減小，而平均速率直線快速增加；當(dāng)vair?in＞0.1 m/s時(shí)，隨著vair?in增大，速度均勻度減小較快，而平均速率則比較緩慢的增加。而且，研究發(fā)現(xiàn)流速為0.1～0.6 m/s時(shí)混合效果較好[13]，因?yàn)楫?dāng)流速小于0.1 m/s時(shí)，微藻容易發(fā)生沉淀；當(dāng)流速大于0.6 m/s時(shí)，流體剪切力會(huì)損傷微藻細(xì)胞，并且能耗也很大。綜合考慮混合均勻度以及能量消耗后，認(rèn)為較優(yōu)的進(jìn)氣速度為0.1 m/s，此時(shí)速度均勻度較高，能量消耗少，平均速度也滿足混合的要求。

圖15 速度均勻度及平均速率隨進(jìn)氣速度變化的曲線

通過(guò)以上理論分析和各項(xiàng)模擬后發(fā)現(xiàn)，藻液的速度均勻度隨傾斜率增加是先增大后減小，最后又增大，在tanθ=0.1時(shí)達(dá)到極大值；隨藻液深度增加而增大；隨進(jìn)氣速度增加而減小。氣含率隨傾斜率增加是先增大后減小，最后又增大，在tanθ=0.1時(shí)達(dá)到極大值；隨藻液深度增加是先減小后增大，最后減小，在d=0.3 m時(shí)達(dá)到極大值。平均速度隨進(jìn)氣速度增加先快速增加后緩慢增加，轉(zhuǎn)折點(diǎn)發(fā)生在vair?in=0.1 m/s。增加導(dǎo)流槽可以提高反應(yīng)器的表面積體積比，但其最大值仍小于10。此外，倒角處理可以增加角落區(qū)的混合，消除流動(dòng)死區(qū)。

3 結(jié)論

通過(guò)CFD 模擬了跑道池增加導(dǎo)流槽后的二維簡(jiǎn)化流場(chǎng)分布以及理論分析表明，增加導(dǎo)流槽可以提高反應(yīng)器的表面積體積比，但增幅不大；主要是改善了藻液垂直方向上的混合，從而能夠增加微藻的總產(chǎn)量；此外，還能延長(zhǎng)氣泡在藻液中的停留時(shí)間，有利于提高CO2利用效率。在模型池寬2 m、導(dǎo)流槽寬1 m的假設(shè)下，優(yōu)化的設(shè)計(jì)參數(shù)為tanθ=0.1，d=0.3 m，vair?in=0.1 m/s。 此時(shí)速度場(chǎng)分布比較均勻，能量消耗也較小，平均速率適中，能夠有效抑制微藻附壁現(xiàn)象并改善藻液的混合。

[1]劉明升，魏群，蔡元妃，等.六種微藻固定CO2實(shí)驗(yàn)研究[J].廣西大學(xué)學(xué)報(bào)：自然科學(xué)版，2012，37（3）：544-548.

[2]鄒樹平，吳玉龍，楊明德，等.微藻的綜合開發(fā)利用[J].水產(chǎn)科學(xué)，2007，26（3）：179-181.

[3]Wang Bei，Lan Christopher Q，Horsman Mark，et al.Closed photobioreactors for production of microalgae biomasses[J].Biotechnology Advances，2012，30（4）：904-912.

[4]李志勇，郭祀遠(yuǎn)，等.不同類型光生物反應(yīng)器的螺旋藻培養(yǎng)特性研究[J].廣西大學(xué)學(xué)報(bào)：自然科學(xué)版，1998，23（3）：290-293.

[5]劉娟妮，胡萍，姚領(lǐng)，等.微藻培養(yǎng)中光生物反應(yīng)器的研究進(jìn)展[J].食品科學(xué)，2006，27（12）：772-777.

[6]Bitog J P，Lee I B，Lee C G，et al.Application of computational fluid dynamics for modeling and designing photobioreactors for microalgae production：A review[J].Computers and Electronics in Agriculture，2011，76（2）：131-147.

[7]Huggins Dania L，Piedrahita Raul H，Rumsey Tom，et al.Analysis of sediment transport modeling using computational fluid dynamics(CFD)for aquaculture raceway[J].Aquaculture Engineering，2004，31（3-4）：277-293.

[8]李昱喆，齊祥明，劉天中，等.帶波紋隔板的平板式光生物反應(yīng)器流動(dòng)特性[J].過(guò)程工程學(xué)報(bào)，2010，10（5）：849-855.

[9]諸發(fā)超，黃建科，陳劍佩，等.敞開式跑道池光生物反應(yīng)器的CFD模擬與優(yōu)化[J].化工進(jìn)展，2012，31（6）：1184-1192.

[10]Liffman Kurt，Paterson David A，Liovic Petar，et al.Comparing the energy efficiency of different high rate algal raceway pond designs using computational fluid dynamics[J].Chemical Engineering Research and Design，2013，91（2）：221-226.

[11]Hekmat Atefeh，Amooghin Abtin Ebadi，Moraveji Mostafa Keshavarz，et al.CFD simulation of gas-liquid flow behaviour in an air-lift reactor：determination of optimum distance of the draft tube[J].Simulation Modeling Practice and Theory，2010，18（7）：927-945.

[12]Chiaramonti David，Prussi Matteo，Casini David，et al.Review of energy balance in raceway ponds for microalgae cultivation ：Rethinking a traditional system is possible[J].Applied Energy，2013，102：101-111.

[13]Ketheesan B，Nirmalakhandan N.Development of a new airlift-driven raceway reactor for algal cultivation[J].Applied Energy，2011，88（10）：3370-3376.

[14]Petkov G D.Adsorber tower as a photobioreactor for microalgae[J].Russian Journal of Plant Physiology，2000，47（6）：786：788.

[15]康瑞娟，蔡昭鈴，施定基.光強(qiáng)在魚腥藻7120 培養(yǎng)液中的衰減及其對(duì)藻細(xì)胞生長(zhǎng)的影響[J].化工冶金，2000，21（4）：384-388.

[16]陳燁，駱其君，馬斌，等.角毛藻光衰減特性及其在光生物反應(yīng)器培養(yǎng)中的應(yīng)用[J].寧波大學(xué)學(xué)報(bào)：理工版，2004，17（2）：143-146.