蘭孝政， 萬　榮， 唐衍力， 黃六一， 趙芬芳， 鮑偉光， 劉長東， 王召磊

(中國海洋大學(xué)水產(chǎn)學(xué)院，山東 青島 266003)

?

圓臺型人工魚礁單體流場效應(yīng)的數(shù)值模擬*

蘭孝政， 萬榮**， 唐衍力， 黃六一， 趙芬芳， 鮑偉光， 劉長東， 王召磊

(中國海洋大學(xué)水產(chǎn)學(xué)院，山東 青島 266003)

應(yīng)用CFD軟件FLUENT開展了圓臺型人工魚礁單體流場效應(yīng)的數(shù)值計算，對比分析了圓臺型人工魚礁與金字塔型、三棱柱型、方型人工魚礁的上升流最大高度、最大速度、平均速度等參數(shù)，討論了圓臺型礁的背渦區(qū)尺寸和流態(tài)。研究顯示，圓臺型礁的上升流最大高度可達(dá)礁體高度的2.14～2.17倍，上升流最大速度可達(dá)來流速度的0.76倍，上升流平均速度可達(dá)來流速度的0.125倍；背渦區(qū)長度約為礁體高度的2.3倍。研究結(jié)果表明，圓臺型人工魚礁流場流態(tài)的大部分指標(biāo)優(yōu)于其它4種常見礁型的流場流態(tài)，且背渦區(qū)漩渦多變、流態(tài)復(fù)雜，具有較好的應(yīng)用前景。

圓臺型人工魚礁； 流場效應(yīng)； 數(shù)值模擬； 水槽實驗

引用格式：蘭孝政， 萬榮， 唐衍力， 等. 圓臺型人工魚礁單體流場效應(yīng)的數(shù)值模擬[J]. 中國海洋大學(xué)學(xué)報(自然科學(xué)版), 2016, 46(8): 47-53.

LAN Xiao-Zheng, WAN Rong, TANG Yan-Li, et al. Numerical simulation of the flow field around the truncated-cone shaped artificial reef [J]. Periodical of Ocean University of China, 2016, 46(8): 47-53.

近十多年來，以人工魚礁建設(shè)、增殖放流為主要內(nèi)容的海洋牧場建設(shè)已成為中國近海漁業(yè)資源修復(fù)與保護的重要手段，并取得了顯著成效。大量實踐表明,通過在適宜的海域投放人工魚礁，可以改善局部水域的流場環(huán)境、生物環(huán)境和魚類棲息環(huán)境，促進(jìn)該海域漁業(yè)生態(tài)環(huán)境的修復(fù)，從而達(dá)到增殖漁業(yè)資源的目的。

在人工魚礁區(qū)設(shè)計與規(guī)劃中要求人工魚礁投放后的流場效應(yīng)要有利于促進(jìn)水域餌料生物的富集和魚類的集聚。由于人工魚礁流場效應(yīng)主要受投放水域底質(zhì)、地形、本底流場和礁體形狀及其配置方式等因素的影響，因此掌握和了解不同類型人工魚礁礁體的水動力學(xué)特性，對于礁型選擇和人工魚礁區(qū)的設(shè)計與規(guī)劃具有重要意義。風(fēng)洞或水槽模型實驗是人工魚礁水動力學(xué)研究的傳統(tǒng)方法，具有實驗條件可控、實驗過程直觀等優(yōu)點，但需要建造昂貴的風(fēng)洞或水槽實驗設(shè)施，實驗費用較高。近年來，隨著計算流體動力學(xué)的迅速發(fā)展及其應(yīng)用領(lǐng)域的拓展，許多研究者開始嘗試將CFD理論與技術(shù)應(yīng)用于人工魚礁流場效應(yīng)的數(shù)值模擬研究。例如，F(xiàn)ujihara采用數(shù)值模擬法分析礁體上升流流態(tài)變化[1]，為人工魚礁的流場研究提供了新思路。潘靈芝對實心方型礁體的二維流場進(jìn)行了數(shù)值模擬分析[2]；劉洪生[3]、李珺[4]借助水槽實驗，驗證了正方體、三棱柱型、金字塔型和米字型人工魚礁流場效應(yīng)的數(shù)值模擬研究的可行性；隨著CFD技術(shù)在人工魚礁流場效應(yīng)研究中應(yīng)用的深入，先后在礁體開孔[5]、迎流面形狀[6]、礁體與水流的擺放角度[7]、礁體布設(shè)間距[8]、投放個數(shù)與排列方式[9]等諸要素對人工礁體流場效應(yīng)的影響研究方面，取得了系列成果，經(jīng)歷了從二維流場到三維流場模擬，從簡單礁型到復(fù)雜礁型的發(fā)展過程。

本文借助CFD軟件FLUENT，以圓臺型人工魚礁為研究對象，建立人工魚礁流場效應(yīng)的數(shù)值模擬分析平臺，開展不同來流條件下圓臺型人工魚礁流場效應(yīng)的數(shù)值模擬研究，通過計算礁體上升流和背渦流特征，分析其形成機制及主要影響因素，為人工魚礁選型和礁區(qū)設(shè)計與規(guī)劃提供科學(xué)依據(jù)。

1　材料與方法

1.1 實物礁體

如圖1所示，圓臺型人工魚礁礁體高2.4m，坡度320%，上、下表面直徑分別為1.5、3.0m，礁體壁厚0.1m。其中，在礁體壁面上均勻開有8個長方型縫隙(0.1m×1.2m)，縫隙開口寬度0.1m，其上緣距離礁體上表面0.6m。由于本實驗僅限于分析人工魚礁的流態(tài)變化，故不考慮礁體重量、使用年限等因素。

圖1　圓臺型人工魚礁結(jié)構(gòu)圖

1.2 水槽模型實驗

根據(jù)模型實驗水槽的實際條件，確定模型尺度比λL=20，模型實驗按照重力相似準(zhǔn)則進(jìn)行。圓臺型人工魚礁為鋼筋混凝土結(jié)構(gòu)，糙率為0.014。根據(jù)糙率相似原理，魚礁模型的糙率為0.0085，選取有機玻璃作為模型材料[10]。以魚礁底面圓心為原點，建立三維直角坐標(biāo)系，x軸正向為來流方向，y軸正向為鉛直向上，z軸分別垂直于x和y軸。如圖2所示，將模型礁體置于水槽測試區(qū)域，并使其中的一個開縫正對來流，所設(shè)置的49個測速點均分布于礁體的中心縱截面上(z=0截面)。

圖2　人工魚礁模型實驗流速測點分布示意圖

模型實驗在中國海洋大學(xué)“增養(yǎng)殖工程水動力實驗室”的動水槽中進(jìn)行。實驗時將模型礁體固定于光滑薄銅板上，以防礁體出現(xiàn)滑移、滾動。假設(shè)薄銅板為光滑平面，其邊界層對流場的影響可以忽略[9]。設(shè)定水槽實驗的流速為0.18m·s-1(相當(dāng)于實際流速0.80m·s-1)，采用多普勒聲學(xué)測速儀測量各測點的流速。

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

假設(shè)流體為不可壓縮的黏性流體，且不考慮流場溫度變化[11]，則水流流經(jīng)魚礁礁體的控制方程如下：

(1)

(2)

(3)

(4)

1.3.1 計算域圖3為圓臺型人工魚礁流場計算域示意圖。如圖所示，設(shè)定x和y方向計算域長分別為27m(魚礁前9m、后15m)和12m，z軸方向計算域為18m，且呈左右對稱。為計算方便，將計算域再次劃分成A、B兩部分，其中B部分尺寸在x和y方向的長度分別為14和7.2m。

1.3.2 網(wǎng)格劃分和邊界條件采用ICEM分別對A、B兩部分進(jìn)行網(wǎng)格劃分，A部分采用六面體網(wǎng)格劃分，B部分采用四面體網(wǎng)格劃分。六面體網(wǎng)格質(zhì)量高且需要生成的網(wǎng)格數(shù)量相對較少，而四面體網(wǎng)格能很好地貼合復(fù)雜幾何模型。將A、B兩部分網(wǎng)格分別導(dǎo)入FLUENT，通過interface建立聯(lián)接[12]。

其次，進(jìn)行網(wǎng)格獨立性檢驗，最終將A、B兩部分網(wǎng)格最大尺寸分別設(shè)置為0.3和0.2m，靠近魚礁表面網(wǎng)格加密。

數(shù)值模擬的邊界條件采用速度入口、壓力出口，礁體底面及側(cè)面采用無滑移壁面邊界條件，上表面為滑移壁面邊界條件[13]。

圖3　魚礁單體計算域及其側(cè)視圖、俯視圖

2　結(jié)果與討論

2.1 實驗驗證

圖4和5分別為實際來流速度0.8m·s-1時，各測點x方向與y方向流速的水槽實測值與數(shù)值計算值的比較。如圖4、5所示，在礁體前部及上方遠(yuǎn)離魚礁處，各點流速幾乎不受魚礁的影響；在礁體近前方附近，呈現(xiàn)x方向流速下降，y方向流速增加的趨勢，其中個別測點處x方向速度接近0，而y方向流速處于峰值。各測點速度的整體變化趨勢實測值(虛線)與數(shù)值模擬結(jié)果(實線)基本一致。方差檢驗表明，在誤差較大的測點處(位于礁體背渦區(qū)的P31、P37)，其不同時刻所測數(shù)據(jù)之間的方差較大，表明實際流態(tài)變化劇烈，紊流強度大，從而也使得其實際測量難度加大。如圖4、5所示，在位于礁體背渦區(qū)的P31、P37處，速度在x、y方向的數(shù)值計算值均為負(fù)值，這是因為在礁體的背渦區(qū)形成了順時針漩渦，此結(jié)果與劉洪生[3]、李珺[4]所得結(jié)果吻合。

圖4　x方向流速水槽實驗與數(shù)值模擬結(jié)果的比較

圖5　y方向流速水槽實驗與數(shù)值模擬結(jié)果的比較

2.2 比較與分析

由于人工魚礁的投放導(dǎo)致魚礁區(qū)流態(tài)的變化(例如產(chǎn)生背渦流、上升流等)，從而引起礁區(qū)生物生境改變，因此，本文結(jié)合數(shù)值模擬結(jié)果，重點分析不同來流速度下圓臺型人工魚礁周圍上升流與背渦流的強度、影響范圍，以探討圓臺型人工魚礁的流場效應(yīng)。

2.2.1 上升流特性為便于分析比較，將y方向速度分量與來流速度之比大于或等于5%的區(qū)域定義為上升流區(qū)[14]。圖6為圓臺型人工魚礁附近z=0m截面上的流速分布隨來流速度變化的數(shù)值模擬結(jié)果，采用不同顏色表現(xiàn)不同區(qū)域的流速差異，比色柱最左端的藍(lán)綠色表示流速接近0，右端的紅色代表最大流速，約為來流速度的1.5倍，從藍(lán)綠色向紅色漸變，即速度不斷增大。如圖6所示，在各流速下礁體迎流面上方均產(chǎn)生了明顯的上升流區(qū)域。

圖6　不同來流速度下在z=0截面上圓臺型人工魚礁的速度分布

表1為圓臺型人工魚礁與金字塔型等4種其它構(gòu)型的人工魚礁[15-18]上升流特征參數(shù)的結(jié)果比較。數(shù)值計算結(jié)果顯示:

(1)隨著來流速度的增大，圓臺型人工魚礁上升流最大高度(Hmax)在5.13～5.22m之間波動，幾乎不隨來流速度的增加而改變，約為魚礁高度(H)的2.14～2.17倍，明顯大于底面邊長和礁高均為3m的金字塔型魚礁(Hmax/H=1.77)，但小于邊長均為3m的三棱柱型魚礁(Hmax/H=2.88)、多孔方型魚礁(Hmax/H=2.68)和無孔方型魚礁(Hmax/H=2.63)。

(2)上升流最大速度Vmax與來流速度V的比值(Vmax/V)為0.76，也幾乎不隨來流速度的增加而變化。Vmax/V值與金字塔型魚礁相同，略大于多孔方型魚礁(Vmax/V=0.74)，明顯大于三棱柱型魚礁(Vmax/V=0.58)和無孔方型魚礁(Vmax/V=0.64)。由此可見，圓臺型人工魚礁可產(chǎn)生較大的最大上升流流速。

(4)圓臺型魚礁的上升流最大高度Hmax與礁體迎流面投影面積S的比值(Hmax/S)為1.60，明顯大于底面邊長和礁高均為3m的金字塔型魚礁(Hmax/S=1.18)、邊長均為3m的三棱柱型魚礁(Hmax/S=0.96)、多孔方型礁(Hmax/S=1.15)和無孔方型礁(Hmax/S=0.88)。

此外，數(shù)值計算結(jié)果還顯示，在各流速條件下，上升流最大速度值均出現(xiàn)于圓臺型魚礁縱截面左側(cè)頂部端點(p25，-0.80,2.45,0)附近，即靠近流動分離點處。

表1　5種礁型的主要流場效應(yīng)指標(biāo)

綜上所述，迎流面投影面積對上升流的影響顯著。此外，礁體迎流面坡度也是主導(dǎo)因素之一。例如，方型魚礁迎流面的垂直結(jié)構(gòu)阻流作用明顯，整體上不利于來流沿豎直方向的加速；三棱柱型礁坡度為173%，對來流的加速起緩解作用，而圓臺型礁與金字塔型礁的迎流面坡度在200%～400%之間，對來流加速的促進(jìn)作用明顯，其Hmax/S優(yōu)于其它3種礁型。其次，礁體開口也是影響上升流的主要因素，多孔方型礁的開孔緩解了礁體的阻流作用，提高了礁體透水性，不利于上升流的產(chǎn)生。但是，適當(dāng)?shù)拈_口尺寸反而有利于上升流最大速度的產(chǎn)生。圓臺型礁與多孔方型礁的開口結(jié)構(gòu)增強了礁體的透水作用，不利于其上升流平均速度的提高，且圓臺型礁帶有弧度的迎流面形狀，使來流更易于向礁體兩側(cè)繞流。三棱柱型與金字塔型礁雖然沒有開口，但其迎流面的坡度結(jié)構(gòu)更易產(chǎn)生較高的上升流平均流速。因此，礁體周圍的上升流高度及其平均速度受到礁體的迎流面坡度、迎流面形狀、迎流面投影面積及礁體開口等的綜合影響。

2.2.2 背渦流特性大量實驗研究表明，人工魚礁不僅產(chǎn)生上升流，而且在其背后部會產(chǎn)生背渦流，形成背渦區(qū)。人工魚礁背渦區(qū)的形成更易促使水層混合，豐富水體中的溶解氧和營養(yǎng)鹽，其影響效應(yīng)與背渦區(qū)范圍及內(nèi)部流態(tài)變化、漩渦大小與強度等因素有關(guān)。

如圖6所示，圓臺型人工魚礁的流動分離點位于上表面的前端。由于邊界層分離，在分離點下游的礁體壁面上不斷有不穩(wěn)定漩渦產(chǎn)生，并從壁面脫落，形成典型的尾跡流，組成速度很小的背渦區(qū)。數(shù)值計算結(jié)果表明，隨著流速的增加，背渦區(qū)形態(tài)變化不大，即背渦流的結(jié)構(gòu)幾乎不隨來流速度增加而改變，其背渦區(qū)長度均約為礁體高度的2.3倍，高度略大于礁體高度，約為礁體高度的1.2倍。

圖7為不同來流速度下在z=0截面上魚礁內(nèi)部及背渦區(qū)的流速矢量圖。由于圓臺型人工魚礁的迎流面傾斜且?guī)в谢《?，結(jié)合礁體表面縫隙的透水作用，其產(chǎn)生的漩渦更為復(fù)雜。如圖7所示，僅魚礁內(nèi)部就形成了3個漩渦，魚礁內(nèi)部的上半部分前端形成一個順時針漩渦，相鄰靠后部分形成一個逆時針漩渦，由于縫隙的透水性，來流進(jìn)入魚礁內(nèi)部，在魚礁下半部形成一個順時針漩渦。魚礁背渦區(qū)也形成了多個漩渦，其中一個順時針的大型漩渦幾乎貫穿整個背渦區(qū)。背渦區(qū)上部靠近礁體頂端處有一順時針小型漩渦，靠近魚礁底部則形成一個逆時針小漩渦。從背渦區(qū)的放大圖可以看出，圓臺型人工魚礁內(nèi)部及背渦區(qū)流態(tài)復(fù)雜，形成多種漩渦，加劇了水層的混合。

圖7　不同來流速度下在z=0截面上圓臺型人工魚礁背渦區(qū)流速矢量圖

3　結(jié)論

本文基于CFD模擬技術(shù)，探討了圓臺型人工魚礁在不同來流速度下的流場流態(tài)。

(1)通過水槽模型實驗的初步驗證，可初步認(rèn)為數(shù)值模擬結(jié)果能夠基本反應(yīng)流場中礁體周圍流態(tài)變化，滿足圓臺型人工魚礁流場效應(yīng)的研究需要。

(2)圓臺型人工魚礁產(chǎn)生的上升流最大高度、上升流最大高度與礁體迎流面投影面積比、上升流最大速度與來流速度比、上升流平均速度與來流速度比等主要特征參數(shù)幾乎均不隨來流速度的增加而改變。其中，上升流最大高度約為魚礁高度的2.14～2.17倍，上升流最大高度與礁體迎流面投影面積比約為1.60，上升流最大速度與來流速度比約為0.76，上升流平均速度與來流速度比約為0.125。而且不同來流速度下上升流最大速度都出現(xiàn)于(-0.80,2.45,0)附近，靠近流動分離點。魚礁迎流面坡度、投影面積、形狀及開口等都是影響上升流流態(tài)的重要因素。

(3)不同來流速度下，背渦流的結(jié)構(gòu)幾乎不隨來流速度增加而改變。圓臺型魚礁背渦區(qū)長度均約為礁體高度的2.3倍，背渦區(qū)高度約為礁體高度的1.2倍，且流態(tài)復(fù)雜、漩渦多變。

綜上所述，圓臺型人工魚礁在與其它4種主要礁型的相比具有一定優(yōu)勢，可為今后人工魚礁的投放提供一種有效的礁型選擇，但今后有待于對其礁體穩(wěn)定性及組合礁的流場流態(tài)等進(jìn)行更為全面、深入的研究。

[1]Fujihara M, Takeuchi T G, Ohashi. Physical-biological coupled modeling for artificially generated upwelling[J]. Marine Biology, 1997, 189: 69-79.

[2]潘靈芝, 林軍, 張守宇. 鉛直二維定長流中人工魚礁流場效應(yīng)的數(shù)值實驗[J]. 上海水產(chǎn)大學(xué)學(xué)報, 2005, 14(4): 406-412.

PAN Ling-zhi, LIN Jun, ZHANG Shou-yu. A numerical experiment of the efects of artificial reef on vertical 2-dimensional steady flow fiel[J]. Journal of Shanghai Fisheries University, 2005, 14(4): 406-412.

[3]劉洪生, 馬翔, 章守宇, 等. 人工魚礁流場風(fēng)洞實驗與數(shù)值模擬對比驗證[J]. 中國水產(chǎn)科學(xué), 2009, 16(3): 365-371.

LIU Hong-sheng, MA Xiang, ZHANG Shou-yu, et al. Validation and comparison between wind tunnel experiments and numerical simulation of flow field around artificial reefs[J]. Journal of Fishery Sciences of China, 2009, 16(3): 365-371.

[4]李琚, 章守宇. 米字型人工魚礁流場數(shù)值模擬與水槽實驗的比較[J]. 水產(chǎn)學(xué)報, 2010, 34(10): 1587-1594.

LI Jun, ZHANG Shou-yu. The comparison between numerical simulation and water channel experiment on an M i-zi artificial reef[J]. Journal of Fisheries of China, 2010, 34(10): 1587-1594.

[5]李琚, 林軍, 章守宇. 方形人工魚礁通透性及其對礁體周圍流場影響的數(shù)值實驗[J]. 上海海洋大學(xué)學(xué)報, 2010, 19(6): 836-840.

LI Jun, LIN Jun, ZHANG Shou-yu. The numerical experiment on the permeability of a cubic artificial reef and the efect on the flow field around the reef[J]. Journal of Shanghai Ocean University, 2010, 19(6): 836-840.

[6]付東偉, 欒曙光, 張瑞瑾, 等. 人工魚礁開口比和迎流面形狀對流場效應(yīng)影響的雙因素方差分析[J]. 大連海洋大學(xué)學(xué)報, 2012, 27(3): 274-278.

FU Dong-wei, LUAN Shu-guang, ZHANG Rui-jin , et al. Two-way analysis of variance of efects of cut-opening ratio and surface shape facing flowing in artificial fish-reefs on the flowing field[J]. Journal of Dalian Ocean University, 2012, 27(3): 274-278.

[7]鄭延璇, 關(guān)長濤, 宋協(xié)法, 等. 星體型人工魚礁流場效應(yīng)的數(shù)值模擬[J]. 農(nóng)業(yè)工程學(xué)報, 2012, 28(19): 185-193.

Zheng Yan-xuan, Guan Chang-tao, Song Xie-fa1, et al. Numerical simulation on flow field around star artificial reefs[J]. Transactions of the Chinese Society of Agricultural Engineering, 2012, 28(19): 185-193.

[8]林軍, 章守宇, 葉靈娜. 基于流場數(shù)值仿真的人工魚礁組合優(yōu)化研究[J]. 水產(chǎn)學(xué)報, 2013, 37(7): 1023-1031.

LIN Jun, ZHANG Shou-yu, YE Ling-na. Optimization study of artificial reef assemblage based on the numerical simulation of flow field[J]. Journal of Fisheries of China, 2013, 37(7): 1023-1031.

[9]鄭延璇, 梁振林, 關(guān)長濤, 等. 三種疊放形式的圓管型人工魚礁流場效應(yīng)數(shù)值模擬與PIV試驗研究[J]. 海洋與湖沼, 2014, 45(1): 11-19.

ZHENG Yan-Xuan, LIANG Zhen-Lin, GUAN Chang-Tao, et al. Numerical simulation and experimental study on flow field of artificial reefs in three tube-stacking layouts[J]. Oceanologia et Limnologia Sinica, 2014, 45(1): 11-19.

[10]唐衍力, 王 磊, 梁振林, 等. 方型人工魚礁水動力性能試驗研究[J]. 中國海洋大學(xué)學(xué)報(自然科學(xué)版), 2007, 37(5): 713-716.

TANG Yan-Li, WANG Lei, LIANG Zhen-Lin, et al. Test of the hydrodynamic performance of square artificial reefs[J]. Periodical of Ocean University of China, 2007, 37(5): 713-716.

[11]王福軍. 計算流體動力學(xué)分析[M]. 北京: 清華大學(xué)出版社, 2004.

WANG Fu-jun. Computational Fluid Dynamics Analysis[M]. Beijing: Tsinghua University Press, 2004.

[12]丁源, 王清. Ansys Icem CFD 從入門到精通[M]. 北京: 清華大學(xué)出版社, 2012.

DING Yuan, WANG Qing. Ansys Icem CFD from Entry to Master[M]. Beijing: Tsinghua University Press, 2012.

[13]Huggins D L, Piedrahita R H, Rumsey T. Analysis of sediment transport modeling using computational fluid dynamics (CFD) for aquaculture raceways[J]. Aquacultural Engineering, 2004, 31(3-4): 277-293.

[14]黃遠(yuǎn)東, 付登楓, 何文榮. 人工魚礁開口比對流場效應(yīng)影響的三維數(shù)值模擬研究[J]. 水產(chǎn)資源與水工程學(xué)報, 2014, 25(4): 39-43.

HUANG Yuan-dong, FU Deng-feng, HE Wen-rong. Three dimensional numerical simulation on influence of cut-opening ratio of artificial reef on flow effect[J]. Journal of Water Resources & Water Engineering, 2014, 25(4): 39-43.

[15]何文榮, 黃遠(yuǎn)東, 黃黎明, 等. 金字塔型人工魚礁繞流的三維CFD模擬研究[J]. 水資源與水工程學(xué)報, 2013, 24(5): 71-76.

HE Wen-rong, HUANG Yuan-dong, HUANG Li-ming, et al. Simulation of three-dimensional CFD of water flow at pyramid artificial reef[J]. Journal of Water Resources & Water Engineering, 2013, 24(5): 71-76.

[16]黃遠(yuǎn)東, 龍 催, 鄧濟通. 三棱柱型人工魚礁繞流流場的CFD 分析[J]. 水資源與水工程學(xué)報, 2013, 24(1): 1-3.

HUANG Yuan-dong, LONG Cui, DENG Ji-tong. Analysis of water flow field based on CFD at a three-prism artificial reef[J]. Journal of Water Resources & Water Engineering, 2013, 24(1): 1-3.

[17]黃遠(yuǎn)東, 趙樹夫, 姜劍偉, 等. 多孔方型人工魚礁繞流的數(shù)值模擬研究[J]. 水資源與水工程學(xué)報, 2012, 23(5): 15-18.

HUANG Yuan-dong, ZHAO Shu-fu, JIANG Jian-wei, et al. Study on numerical simulation of water flow past a porous square artificial reef[J]. Journal of Water Resources & Water Engineering, 2012, 23(5): 15-18.

[18]黃遠(yuǎn)東, 姜劍偉, 趙樹夫. 方型人工魚礁周圍水流運動的數(shù)值模擬研究[J]. 水資源與水工程學(xué)報, 2012, 23(3): 1-3.

HUANG Yuan-dong, JIANG Jian-wei, ZHAO Shu-fu. Study on numerical model of water flows past a square artificial reef[J]. Journal of Water Resources & Water Engineering, 2012, 23(3): 1-3.

責(zé)任編輯朱寶象

Numerical Simulation of the Flow Field Around the Truncated-cone Shaped Artificial Reef

LAN Xiao-Zheng, WAN Rong, TANG Yan-Li, HUANG Liu-Yi, ZHAO Fen-Fang,BAO Wei-Guang, LIU Chang-Dong, WANG Zhao-Lei

(College of Fisheries, Ocean University of Chain, Qingdao 266003, China)

Artificial reefs, built in the proper sea area, change the flow pattern of local waters thus improving the biological environment and the environment of fish habitat, which promote the restoration of fishery ecological environment and achieve the proliferation of fishery resources. Therefore, the effect on the flow field made by artificial fishery reefs is the basis of their effect on the ecosystem. The flow field is mainly affected by the shape of the artificial reefs. It is important to have a good understand on the hydrodynamic characteristics of various reefs for the shape-selection in the designing plan of artificial reefs. In the present work, the flow field around the truncated-cone shaped artificial reef was numerically simulated to investigate its hydrodynamic characteristics by means of CFD software FLUENT (a special software package for the computational fluid mechanics) under five incoming velocities (0.1, 0.2, 0.4, 0.6 and 0.8 meter per second respectively). And then, flume experiments were performed with the models of selected shape. Comparison between numerical and experimental results showed that the relative error was within twenty percent for more than eighty percent of measuring points. It may be concluded that the results of numerical simulation were reliable. Analyses of the maximum velocity of the upwelling current, the maximum upwelling height, average velocity of the upwelling current and other parameters of the truncated-cone shaped artificial reef, single-bodied, immersed in the flow field, were comparable with those of pyramid , three-prism , porous square and square artificial reefs. We discussed and analyzed also the wake vortices of the truncated-cone shaped artificial reef. The generation and pattern of flow separation, upwelling and vortex development were observed as the water flows past the artificial reefs. It was found that this was a typical flow around a bluff body including complicate nonlinear procedure. The numerical results obtained showed that the maximum upwelling height was 2.14～2.17 times of the reef height; the maximum and average velocities of the upwelling current were 0.76 and 0.125 times of the incoming velocity, respectively; and the length of the back vortex region was about 2.3 times of the reef height. In comparison with the other four types of artificial reefs including pyramid-shape, triangular prism and a square column with or without openings, it was found that the truncated-cone shaped artificial reef had a fairly good performance in most aspects of flow parameters. The wake vortices vary swiftly, which may lead to a complicated flow pattern in the wake of truncated-cone shaped artificial reef. It was expected that this type of artificial reef has a good prospect in the future application.

truncated-cone shaped artificial reef; hydrodynamic characteristics; numerical simulation; flume experiment

山東省科技發(fā)展計劃項目(2011GHY11524)資助

2015-04-20；

2015-11-01

蘭孝政(1988-)，男，碩士生。E-mail: ls-lxz@163.com

**通訊作者：E-mail：rongwan@ouc.edu.cn

S911

A

1672-5174(2016)08-047-07

10.16441/j.cnki.hdxb.20150149

Supported by the Science and Technology Development Project of Shandong Province(2011GHY11524)