明渠紊流中淹沒柔性植被流固耦合大渦數(shù)值模擬

許?棟，梁心雅，張炳昌，及春寧

明渠紊流中淹沒柔性植被流固耦合大渦數(shù)值模擬

許?棟，梁心雅，張炳昌，及春寧

(天津大學(xué)水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室，天津 300072)

水生植被廣泛存在于天然河道、湖泊及濕地，是河流生態(tài)系統(tǒng)的重要成分．它通過改變明渠水流時(shí)均及紊動(dòng)結(jié)構(gòu)來改變河道水流阻力，并在水流作用下發(fā)生彎曲和擺動(dòng)，形成復(fù)雜的流固耦合作用．研究柔性植被在水流作用下的運(yùn)動(dòng)響應(yīng)，能夠更精確地解析植被的阻水效應(yīng)．本文采用基于嵌入式迭代浸入邊界法的高效能三維紊流計(jì)算程序CgLes，開展流固耦合數(shù)值模擬，研究柔性植被在水流作用下的彎曲及運(yùn)動(dòng)響應(yīng)，分析植被渦激振動(dòng)特性及水流阻力規(guī)律．模擬結(jié)果表明，植被運(yùn)動(dòng)以展向擺動(dòng)為主導(dǎo)，植被剛度對(duì)植被運(yùn)動(dòng)軌跡和擺動(dòng)周期有顯著影響，進(jìn)而影響植被的阻水規(guī)律．在水動(dòng)力條件相同的情況下，植被的擺動(dòng)周期與植被的彈性模量有關(guān)，且不受流場(chǎng)影響．隨著剛度系數(shù)增大，植被運(yùn)動(dòng)模式從準(zhǔn)周期性擺動(dòng)變成周期性擺動(dòng)，最終發(fā)展為剛性情況下的靜止直立狀態(tài)；植被剛度較大且間距較小時(shí)，植被存在逆流向彎折現(xiàn)象；在剛度較大時(shí)，柔性植被隨著水流而擺動(dòng)，抑制流速發(fā)展；但是當(dāng)柔性植被剛度較小時(shí)，植被接近倒伏，明渠的平均流速對(duì)比剛性植被情況增大；本文算例中植被剛度對(duì)平均流速的影響最大約為2.47%．植被間距對(duì)植被的下彎幅度影響較小，但對(duì)展向擺幅作用顯著；平均流速隨植被間距增大而減小，植被間距主要影響植被層內(nèi)部的水流流速，本文算例中植被間距對(duì)平均流速的影響最大約為8.19%．

柔性植被；渦激振動(dòng)；浸入邊界法；流固耦合；數(shù)值模擬

水生植被廣泛存在于天然河道、湖泊及濕地，是河流生態(tài)系統(tǒng)的重要組成部分．水生植被的研究對(duì)實(shí)現(xiàn)河流保護(hù)十分關(guān)鍵．從水動(dòng)力學(xué)的角度而言，水生植被可改變河流水流的阻力，增高河流水位，使其呈現(xiàn)出復(fù)雜的水動(dòng)力特征[1-2]．此外，水生植被能降解水生態(tài)系統(tǒng)中的微生物，影響污染物輸移[3-4]，也影響著營(yíng)養(yǎng)物質(zhì)和重金屬的分布[5-6]．另一方面，水流動(dòng)力特性也對(duì)水生植被的分布和形態(tài)起影響作用[7]．

為了簡(jiǎn)化計(jì)算，水動(dòng)力數(shù)值模擬中往往將植被概化為剛性圓柱[1-2]．然而，水流中的天然植被通常是柔性的，與保持直立的剛性植被相比，植被運(yùn)動(dòng)對(duì)水流結(jié)構(gòu)及水動(dòng)力特征的影響不可忽略[8]．通?？蓪⒀蜎]柔性植被水流分成不同區(qū)域進(jìn)行研究．槐文信等[9]以植被平均倒伏高度為界，將淹沒柔性植被水流分為植被層和植被上層兩個(gè)區(qū)域．Nepf和Vivoni[10]將植被層中的水流分為位于植被頂部附近的垂直交換區(qū)和位于垂直交換區(qū)下方的縱向交換區(qū)．Nikora等[11]提出柔性植被后方垂向流速分布從上到下分別滿足對(duì)數(shù)律、尾流律(受植被影響的自由水層)、雙曲正切分布(植被頂部附近)和均勻分布(植被內(nèi)部)，且區(qū)域的劃分會(huì)根據(jù)植被類型、試驗(yàn)布置、水流流態(tài)等不同而發(fā)生變化．與剛性植被相比，植被的柔性對(duì)水流阻力的影響是增大還是降低，目前尚缺少?gòu)V泛統(tǒng)一的認(rèn)識(shí)．Ghisalberti和Nepf[12]發(fā)現(xiàn)植被和水流之間的耦合作用減小了植被的拖曳力，增加了流速和湍流應(yīng)力.楊婕等[13]通過模型試驗(yàn)指出柔性植被覆蓋下阻力系數(shù)隨著雷諾數(shù)和流量增大均呈減小趨勢(shì)．王曉燕[14]認(rèn)為植被剛度越小，阻力系數(shù)也越?。嘶莸萚15]提出對(duì)不同流量分別存在一個(gè)臨界單株植被剛度，并定義小于該剛度為“植被粗糙區(qū)”，其水頭損失同時(shí)受水流流量和植被剛度的雙重影響，且均為正相關(guān)；大于該剛度為“植被光滑區(qū)”，水頭損失與植被剛度無關(guān)，僅與流量呈正相關(guān)．師冰雪[16]認(rèn)為在柔性植被發(fā)生穩(wěn)定彎曲時(shí)，阻力系數(shù)減小；而植被發(fā)生擺動(dòng)時(shí)，阻力系數(shù)增大．Abdolahpour等[17]研究了兩種柔性淹沒植被對(duì)海岸水流的影響，發(fā)現(xiàn)柔性冠層重分布抑制了冠層中的流速衰減．

柔性植被在水流影響下表現(xiàn)出一定的運(yùn)動(dòng)和彎曲．當(dāng)水流流經(jīng)時(shí)，圓柱后方會(huì)出現(xiàn)交替的渦旋脫落，誘發(fā)渦激振動(dòng)[18]．Vogel[19]的研究表明，瞬時(shí)的、被動(dòng)的和可逆的形狀變化會(huì)導(dǎo)致承受高速水流的柔性植物內(nèi)應(yīng)力顯著減小．對(duì)于這些流動(dòng)誘導(dǎo)變形的所謂適應(yīng)性，Vogel[20]建議使用“重分布”一詞，因?yàn)檫@個(gè)詞有更積極的含義．Leclercq等[21]研究了彈性薄葉片的重構(gòu)以及當(dāng)葉片周圍的流體振蕩時(shí)產(chǎn)生的內(nèi)應(yīng)力，根據(jù)流體顆粒相對(duì)于葉片尺寸的偏移和流動(dòng)振蕩的頻率，確定了4種運(yùn)動(dòng)狀態(tài)．Lei和Nepf[22]研究了波流聯(lián)合作用下的葉片運(yùn)動(dòng)，發(fā)現(xiàn)對(duì)于一個(gè)固定的波浪條件，增加流量會(huì)降低平均偏轉(zhuǎn)高度和葉片運(yùn)動(dòng)范圍．O’Connor和 Revell提出[23]，單株柔性植被在明渠水流下的運(yùn)動(dòng)模式可以分為靜態(tài)模式、擺動(dòng)模式、周期倍增模式和混沌模式4類；多株柔性植被在明渠水流下的運(yùn)動(dòng)模式可以分為靜態(tài)模式、規(guī)則波動(dòng)模式、不規(guī)則波動(dòng)模式和擺動(dòng)模式4類；同時(shí)還發(fā)現(xiàn)多株柔性植被的規(guī)則波動(dòng)模式、不規(guī)則波動(dòng)模式和擺動(dòng)模式下的波動(dòng)頻率隨抗彎剛度比變化的曲線基本重合．Zhang等[24]研究了不同彎曲剛度、質(zhì)量比和雷諾數(shù)對(duì)二維柔性絲的影響，其動(dòng)態(tài)行為模式可以分為3種不同的模式，即倒伏(lodging)、規(guī)則渦激振動(dòng)(VIV)和靜態(tài)重構(gòu)(static reconfiguration)．

在植被對(duì)于水流結(jié)構(gòu)和阻力的影響方面，許多學(xué)者通過使用玻璃棒或塑料棒來模擬剛性植被[25]，用塑料水草[13]模擬柔性植被來進(jìn)行模型試驗(yàn)．同時(shí)，越來越多的學(xué)者開始使用大渦模擬(LES)或格子玻爾茲曼(LBM)等數(shù)值模擬方法，來研究柔性植被在水流作用下的運(yùn)動(dòng)響應(yīng)[23,26]，且在數(shù)值模擬方法中常將柔性植被簡(jiǎn)化為柔性圓柱或柔性板．而在流固耦合方面，部分學(xué)者事先應(yīng)用懸臂梁理論推導(dǎo)出植被彎曲后各點(diǎn)的角度以及彎曲后的高度，然后再依據(jù)不同的植被阻力計(jì)算公式結(jié)合Navier-Stokes(NS)方程建立水流-植被的單向耦合模型[27-28]．在單向耦合模型中，植被不隨水流運(yùn)動(dòng)，與實(shí)際情況偏差較大．Chen和Zou[29]基于浸入邊界法(IBM)和有限體積法，提出了一種能夠同時(shí)在歐拉-拉格朗日混合網(wǎng)格上解析撓度較大的柔性植被運(yùn)動(dòng)和流場(chǎng)的雙向耦合模型．雙向耦合模型考慮了植被運(yùn)動(dòng)對(duì)流場(chǎng)的影響，能更好地反映水流與植被之間的相互作用機(jī)理．

在現(xiàn)有的柔性植被數(shù)值模擬研究中，反映植被運(yùn)動(dòng)與三維流場(chǎng)雙向耦合的研究較少．本文采用大渦模擬(LES)和有限元-離散元方法(FDEM)，對(duì)柔性植被覆蓋的明渠水流開展雙向耦合的精細(xì)數(shù)值模擬，研究和分析柔性植被在水流作用下的運(yùn)動(dòng)響應(yīng)，分析渦激振動(dòng)特性及水流阻力規(guī)律．

1?數(shù)學(xué)模型

采用大渦模擬對(duì)復(fù)雜湍流進(jìn)行研究，植被采用浸入邊界法(IBM)描述，不可壓縮黏性流體Navier-Stokes方程為

植被單元的運(yùn)動(dòng)和碰撞由Munjiza開發(fā)的有限元-離散元方法[32]處理，即聯(lián)合有限元-離散元方法(FDEM)．植被單元采用四節(jié)點(diǎn)四面體單元進(jìn)行網(wǎng)格劃分．固體運(yùn)動(dòng)、變形和碰撞用第二牛頓定律對(duì)每個(gè)節(jié)點(diǎn)求解如下：

該方法的詳細(xì)介紹參見Ji等[30]和Uhlmann[31].流體固體求解方法見參考文獻(xiàn)[32]，流固耦合技術(shù)流程見圖1．

2?模擬工況和驗(yàn)證

2.1?模型建立

采用單列植被配合周期邊界來代表植被群，研究其在明渠水流作用下的流固耦合，通過改變植被單元?jiǎng)偠?、植被群密度，共設(shè)置6個(gè)工況，見表1．

表1?工況參數(shù)

Tab.1?Case parameters

圖2?計(jì)算域和邊界條件示意

在所有工況中，采用相同的網(wǎng)格劃分，其中流場(chǎng)部分采用結(jié)構(gòu)化網(wǎng)格，、、向網(wǎng)格數(shù)量分別為192、64和64，在植被直徑方向平均包含約13個(gè)網(wǎng)格，能夠較好地解析植被繞流流場(chǎng)．植被單元采用正四面體網(wǎng)格表述，四面體邊長(zhǎng)為0.08．植被雷諾數(shù)的表達(dá)式如下：

2.2?模型驗(yàn)證

2.2.1?水動(dòng)力模型驗(yàn)證

為了驗(yàn)證植被影響下明渠紊流模擬的可靠性，采用Liu等[2]的淹沒剛性植被試驗(yàn)．試驗(yàn)在矩形明渠中采用交錯(cuò)植被的排布方式，采用LDV測(cè)量了流場(chǎng)中3個(gè)不同位置的流速垂向剖面．測(cè)點(diǎn)的位置如圖3(a)所示．

2.2.2?柔性植被變形模擬驗(yàn)證

柔性植被可以簡(jiǎn)化為懸臂梁模型，可采用歐拉伯努利方程求解懸臂梁彎曲情況．懸臂梁在均布荷載下的最大撓度方程為

圖5?柔性圓柱的變形驗(yàn)證

3?結(jié)果分析

3.1?剛度系數(shù)對(duì)植被單元運(yùn)動(dòng)響應(yīng)的影響

3.1.1?植被單元的運(yùn)動(dòng)軌跡

單根植被在水流作用下呈現(xiàn)出不同的運(yùn)動(dòng)形態(tài)，許多學(xué)者對(duì)其進(jìn)行了模式分類研究．模擬結(jié)果表明，剛度對(duì)植被運(yùn)動(dòng)軌跡和擺動(dòng)周期有顯著影響．經(jīng)過初始的瞬態(tài)運(yùn)動(dòng)后，柔性植被的運(yùn)動(dòng)變?yōu)檎穹皖l率恒定的周期運(yùn)動(dòng)．圖6為工況E8.1-S15中植被單元在一個(gè)周期內(nèi)的運(yùn)動(dòng)軌跡．此時(shí)植被在展向上的運(yùn)動(dòng)軌跡類似圓形封閉曲線，期間植被沿流向高低起伏．圖6觀察到的植被運(yùn)動(dòng)軌跡和周期倍增模式一致．植被頂端在橫流向的最大變形略大于其在順流向上的最大變形，運(yùn)動(dòng)軌跡類似于橢圓形．當(dāng)植被頂端展向位移最大時(shí)，植被的下彎幅度最小．

隨著植被剛度增加，柔性植被單元的運(yùn)動(dòng)軌跡從流向主導(dǎo)轉(zhuǎn)換為展向主導(dǎo)．圖7為工況E32-S15中植被運(yùn)動(dòng)軌跡．受到明渠水流的作用，植被單元不但在順流向發(fā)生彎折，在展向還存在振動(dòng)，其展向運(yùn)動(dòng)軌跡類似于“∞”，且在從中間位置向外搖擺的過程中運(yùn)動(dòng)軌跡弧度更大，在從最大振幅位置向軸線中心搖擺的過程中其運(yùn)動(dòng)軌跡的弧度較?。藭r(shí)植被運(yùn)動(dòng)模式處于周期振動(dòng)模式，但是其在流向方向上的位移并不明顯，若剛度持續(xù)增大，植被將從周期性振動(dòng)變成靜止直立狀態(tài)．

圖6?典型植被運(yùn)動(dòng)軌跡(工況E8.1-S15)

圖7?典型植被運(yùn)動(dòng)軌跡(俯視，工況E32-S15)

圖8為植被頂端中點(diǎn)在展向上的擺動(dòng)位移隨時(shí)間變化圖，植被在展向上周期性擺動(dòng)，不同植被的頂端運(yùn)動(dòng)軌跡相似，植被之間的展向運(yùn)動(dòng)存在固定的相位差．這兩種工況位移圖分別與O’Connor[23]觀察到的準(zhǔn)周期波動(dòng)、周期波動(dòng)圖像一致．隨著剛度增大，植被運(yùn)動(dòng)模式從準(zhǔn)周期波動(dòng)轉(zhuǎn)化為周期波動(dòng)，若剛度持續(xù)增大，植被的運(yùn)動(dòng)模式將最終轉(zhuǎn)化為靜止?fàn)顟B(tài)．

圖8?植被頂端中點(diǎn)位移時(shí)程(流向)

3.1.2?時(shí)均流速和紊動(dòng)強(qiáng)度

在不同剛度植被的作用下，各工況平均流速存在明顯差異．在剛度系數(shù)最大時(shí)，柔性植被隨著水流擺動(dòng)，增大了水體的能量耗散，進(jìn)而起到減低流體動(dòng)能的效果．但是在植被剛度最小的算例中，平均流速反而比剛性工況的要大．這是由于此時(shí)的植被接近于倒伏，植被群有效高度的變化，自由水層高度顯著增大，水體受到的阻力明顯減小，不能起到阻水效果．圖9(b)的流速剖面中可以觀察到明顯的拐點(diǎn)．植被區(qū)域內(nèi)部的流速小，而自由水層的流速大，這是由于植被作用下，水流受到植被群的阻力增大，阻力在植被群與上層水流交界處的不連續(xù)性使流速剖面在該位置形成拐點(diǎn)．

圖9?植被剛度對(duì)積分平均速度分布的影響

先秦兩漢直至魏晉南北朝時(shí)期，人們的“斥巧”態(tài)度，極大地影響到隋唐五代及北宋中期人們的觀念。隋唐時(shí)期，凡是“巧”與政治、道德、人性等相聯(lián)系，則人們的態(tài)度一定是“斥巧”。如李世民云：“朕歷觀前代，讒佞之徒，皆國(guó)蝥賊，巧令朋比?！盵26](P122)姚班亦云：“至于工巧造作，寮史直司，實(shí)為末事，無足勞慮。”[26](P1733)岑文本也說：“去智絕巧，圣人之至德?！盵26](P1526)上述文獻(xiàn)皆說明：政治之“巧”，是敗壞政體根基的禍害，因此務(wù)須“斥巧”。

圖10 植被剛度對(duì)時(shí)均流向速度脈動(dòng)值分布的影響(黑色圓點(diǎn)表示植被橫截面位置)

3.1.3?剛度系數(shù)和間距對(duì)植被運(yùn)動(dòng)和周期的影響

對(duì)比表2中各柔性植被工況發(fā)現(xiàn)，在外力條件相同的情況下植被的擺動(dòng)周期與植被自身材料性質(zhì)有關(guān)，植被剛性越大，周期越?。S著剛度系數(shù)增大，植被的下彎幅度減小，但展向擺幅增大．模擬結(jié)果表明，植被間距不影響植被下彎幅度，但對(duì)展向擺幅作用顯著．

表2?柔性植被運(yùn)動(dòng)周期和幅度

Tab.2?Period and amplitude of flexible vegetation motion

3.2?植被間距對(duì)柔性植被運(yùn)動(dòng)響應(yīng)的影響

3.2.1?植被間距對(duì)植被運(yùn)動(dòng)軌跡的影響

在植被剛度系數(shù)保持不變的情況下(工況E32-S15、E32-S7.5和E32-S5)，逐漸減小植被間距，發(fā)現(xiàn)當(dāng)植被間距減小到5時(shí)，出現(xiàn)了植被逆流向彎折現(xiàn)象，即植被向上游方向倒伏(見圖11(a))．

隨著植被間距的減小，順流向植被密度增大，植被之間的相互影響更強(qiáng)烈，流場(chǎng)紊動(dòng)加?。谥脖婚g距較小的情況下植被發(fā)生周期性逆流向彎折現(xiàn)象，其統(tǒng)計(jì)周期與植被振動(dòng)周期相同．但在植被間距較大時(shí)并不會(huì)發(fā)生植被逆流向彎折現(xiàn)象，這是由于當(dāng)植被間距較小時(shí)，上游植被單元尾流強(qiáng)度較大，進(jìn)而對(duì)下游植被單元運(yùn)動(dòng)產(chǎn)生影響，當(dāng)植被間距增大后，上游植被尾流強(qiáng)度減弱，不足以驅(qū)動(dòng)下游植被產(chǎn)生逆流向彎折．

圖11?植被姿態(tài)和對(duì)應(yīng)流速云圖(工況E32-S5)

圖11(b)是同一時(shí)刻工況E32-S5的流速場(chǎng)示意圖，黑色圓點(diǎn)表示植被底端的位置．可以看出：植被由于受到慣性和前方植被尾渦的影響，不僅逆流向彎折的植被后方存在回流區(qū)，其前方也存在回流區(qū)．植被尾渦相互耦合，水流結(jié)構(gòu)復(fù)雜．尾流區(qū)域壓強(qiáng)較低，由于植株尾流壓強(qiáng)高于上游側(cè)壓強(qiáng)，出現(xiàn)了逆流彎折現(xiàn)象．

3.2.2?相鄰植被運(yùn)動(dòng)響應(yīng)

以E32-S7.5中的4號(hào)植被一個(gè)運(yùn)動(dòng)周期為基準(zhǔn)，相鄰植被在該時(shí)間范圍內(nèi)的相對(duì)位置變化見圖12．由于模型在流向上是周期邊界的明渠水流，在順流向上，不同位置的植被運(yùn)動(dòng)存在相位差，植被單元的運(yùn)動(dòng)在下一段時(shí)間后總會(huì)按它上游植被的運(yùn)動(dòng)軌跡進(jìn)行．其展向運(yùn)動(dòng)軌跡呈現(xiàn)出波動(dòng)圖像．植被所在的這一條線附近區(qū)域流速普遍較小，而近岸流速較大，這樣在交界處存在周期性大渦結(jié)構(gòu)．植被受到周期性水流力的作用，其受迫擺動(dòng)具有周期性．在植被間距不同的工況中均觀察到此現(xiàn)象．

圖12?一個(gè)周期內(nèi)植被的展向的運(yùn)動(dòng)軌跡

3.2.3?時(shí)均流場(chǎng)

垂向平均流速呈現(xiàn)出隨植被密度增大而減小的規(guī)律(圖13)．相對(duì)于剛度而言，植被間距對(duì)縱向流速剖面的彎曲度和拐點(diǎn)影響較為明顯，在水深/＝9處水流平均流速受植被密度影響最大，正好從植被順流向的最大彎折高度附近．而植被上方自由水層中植被間距對(duì)流速影響不大，說明植被密度對(duì)水流流速的影響主要在植被層內(nèi)部，這與Poggi等[35]的結(jié)論一致．

圖13?植被間距對(duì)垂向平均流速分布的影響

4?結(jié)?論

本文采用大渦模擬和有限元-離散元方法(FDEM)，對(duì)柔性植被覆蓋的明渠水流開展雙向耦合的精細(xì)數(shù)值模擬，研究和分析了柔性植被在水流作用下的運(yùn)動(dòng)響應(yīng)，渦激振動(dòng)特性及水流阻力規(guī)律，試驗(yàn)結(jié)果表明以下幾點(diǎn)．

(1) 剛度系數(shù)不同導(dǎo)致柔性植被存在不同的運(yùn)動(dòng)模式和運(yùn)動(dòng)軌跡，進(jìn)而影響植被的阻水規(guī)律．隨著剛度系數(shù)增大，植被運(yùn)動(dòng)模式從準(zhǔn)周期性擺動(dòng)變成周期性擺動(dòng)，最終發(fā)展為剛性情況下的靜止直立狀態(tài)．植被的運(yùn)動(dòng)軌跡以展向擺動(dòng)為主導(dǎo)，隨著剛度系數(shù)增大，展向擺動(dòng)的主導(dǎo)作用更為明顯．

(2) 植被剛度較大且間距較小時(shí)，植被存在逆流向彎折現(xiàn)象．這是因?yàn)橹脖皇艿綉T性和前方植被尾渦的影響，尾流區(qū)域壓強(qiáng)較低，植株尾流壓強(qiáng)高于上游側(cè)壓強(qiáng)．流向上不同位置的植被運(yùn)動(dòng)存在相位差，呈現(xiàn)出波動(dòng)現(xiàn)象．

(3) 隨著剛度系數(shù)增大，植被的下彎幅度減小，但展向擺幅增大．模擬結(jié)果表明，植被間距不影響植被下彎幅度，但對(duì)展向擺幅作用顯著．在本文算例中，植被的展向擺幅最大為2.3倍植被直徑．

(4) 在外力條件相同的情況下植被的擺動(dòng)周期與植被自身材料性質(zhì)有關(guān)，且不受流場(chǎng)影響，植被剛性越大，周期越?。?/p>

(5) 植被剛度較小時(shí)，植被彎折減小迎流面積，導(dǎo)致平均流速大于完全剛性算例的平均流速；剛度系數(shù)較大時(shí)，倒伏不明顯，柔性植被隨著水流擺動(dòng)，渦激振動(dòng)作用抑制流速發(fā)展．平均流速呈現(xiàn)出隨植被密度增大而減小的規(guī)律，植被密度對(duì)水流流速的影響主要在植被層內(nèi)部．本文算例中植被剛度對(duì)平均流速的影響最大約為2.47%，植被間距對(duì)平均流速的影響最大約為8.19%．

［1］ Le Bouteiller C，Venditti J G. Sediment transport and shear stress partitioning in a vegetated flow[J]. Water Resources Research，2015，51(4)：2901-2922.

［2］ Liu D，Diplas P，F(xiàn)airbanks J D，et al. An experimental study of flow through rigid vegetation[J]. Journal of Geophysical Research，2008，113：F04015.

［3］ 韓?璐. 柔性植被河道水流特試驗(yàn)研究[D]. 南京：河海大學(xué)，2006.

Han Lu. Research on Hydraulic Characteristics in Flexible Vegetation Channel[D]. Nanjing：Hohai University，2006(in Chinese).

［4］ 王?忖，王?超. 含挺水植物和沉水植物水流紊動(dòng)特性[J]. 水科學(xué)進(jìn)展，2010，21(6)：816-822.

Wang Cun，Wang Chao. Turbulent characteristics in open-cannel flow with emergent and submerged macrophytes[J]. Advances in Water Science，2010，21(6)：816-822(in Chinese).

［5］ 郭和蓉. 水生植物凈化環(huán)境與被的修復(fù)[J]. 長(zhǎng)江大學(xué)報(bào)，2005，2(5)：65-68.

Guo Herong. Purification and restoration of aquatic plants[J]. Journal of Yangtze University，2005，2(5)：65-68(in Chinese).

［6］ Bornette G，Puijalon S. Response of aquatic plants to abiotic factors：A review[J]. Aquatic Sciences，2011，73(1)：1-14.

［7］ Nepf H. Hydrodynamics of vegetated channels[J]. Journal of Hydraulic Research，2012，50(3)：262-279.

［8］ 張?娜，張金鳳，張慶河，等. 規(guī)則波作用下柔性植物蜂巢護(hù)岸水動(dòng)力特性試驗(yàn)研究[J]. 天津大學(xué)學(xué)報(bào)(自然科學(xué)與工程技術(shù)版)，2021，54(5)：441-449.

Zhang Na，Zhang Jinfeng，Zhang Qinghe，et al. Experimental study of the hydrodynamics of revetment protection using honeycomb-type，flexible plants under the action of regular waves[J]. Journal of Tianjin University(Science and Technology)，2021，54(5)：441-449 (in Chinese).

［9］ 槐文信，韓?杰，曾玉紅，等. 淹沒柔性植被明渠恒定水流水力特性的試驗(yàn)研究[J]. 水利學(xué)報(bào)，2009(7)：25-31.

Huai Wenxin，Han Jie，Zeng Yuhong，et al. Experimental stuly on hydraulic behaviors of steady uniform flow in open channel with submerged flexible vegetation[J]. Journal of Water Conservancy，2009(7)：25-31(in Chinese).

［10］ Nepf H M，Vivoni E R. Flow structure in depth-limited，vegetated flow[J]. Journal of Geophysical Research，2000，105(C12)：28547-28557.

［11］ Nikora N，Nikora V，O’donoghue T. Velocity profiles in vegetated open-channel flows：Combined effects of multiple mechanisms[J]. Journal of Hydraulic Engineering，2013，139(10)：1021-1032.

［12］ Ghisalberti M，Nepf H M. Shallow flows over a permeable medium：The hydrodynamics of submerged aquatic canopies[J]. Transport in Porous Media，2009，78(2)：309-326.

［13］ 楊?婕，張寬地，楊?帆. 柔性植被和剛性植被水流水動(dòng)力學(xué)特性研究[J]. 人民黃河，2017，39(8)：85-89.

Yang Jie，Zhang kuandi，Yang Fan. Experimental study on hydraulic characteristics of overland flow under rigid and flexible vegetation coverage[J]. People’s Yellow River，2017，39(8)：85-89(in Chinese).

［14］ 王曉燕. 植被剛度對(duì)水流阻力特性影響的研究[D]. 南京：河海大學(xué)，2007.

Wang Xiaoyan. The Effect of Vegetation Rigidity on Flow Resistance[D]. Nanjing：Hohai University，2007 (in Chinese).

［15］ 潘?惠，丁?磊，吳龍華，等. 單株淹沒植被剛度及其阻力特性研究[J]. 湖北農(nóng)業(yè)科學(xué)，2013，52(18)：4515-4518.

Pan Hui，Ding Lei，Wu Longhua，et al. Study on effect of submerged plant vegetation rigidity on flow resistance [J]. Hubei Agricultural Sciences，2013，52(18)：4515-4518(in Chinese).

［16］ 師冰雪. 淹沒植被作用下流速及含沙量垂線分布試驗(yàn)研究[D]. 鄭州：鄭州大學(xué)，2018.

Shi Bingxue. Study on the Vertical Distribution of Water Flow Velocity and Sediment Concentration in Water Permeable Barrier[D]. Zhengzhou：Zhengzhou University，2018(in Chinese).

［17］ Abdolahpour M，Ghisalberti M，Mcmahon K，et al. The impact of flexibility on flow，turbulence，and vertical mixing in coastal canopies[J]. Limnology and Oceanography，2018，63(6)：2777-2792.

［18］ 徐萬海，欒英森，李志彪，等. 傾斜柔性抑制圓柱的渦激振動(dòng)特性[J]. 天津大學(xué)學(xué)報(bào)(自然科學(xué)與工程技術(shù)版)，2017，50(6)：637-642.

Xu Wanhai，Luan Yingsen，Li Zhibiao，et al. Vortex-Induced vibrations of an inclined flexible straked cylinder[J]. Journal of Tianjin University(Science and Technology)，2017，50(6)：637-642(in Chinese).

［19］ Vogel S. Drag and flexibility in sessile organisms[J]. American Zoologist，1984，24(1)：37-44.

［20］ Vogel S. Drag and reconfiguration of broad leaves in high winds[J]. Journal of Experimental Botany，1989，40(8)：941-948.

［21］ Leclercq T，de Langre E. Reconfiguration of elastic blades in oscillatory flow[J]. Journal of Fluid Mechanics，2018，838：606-630.

［22］ Lei J，Nepf H. Blade dynamics in combined waves and current[J]. Journal of Fluids and Structures，2019，87：137-149.

［23］ O’Connor J，Revell A. Dynamic interactions of multiple wall-mounted flexible flaps[J]. Journal of Fluid Mechanics. 2019，870：189-216.

［24］ Zhang X，He G，Zhang X. Fluid-structure interactions of single and dual wall-mounted 2D flexible filaments in a laminar boundary layer[J]. Journal of Fluids and Structures，2020，92：102787.

［25］ 宋瀅汀，景何仿，張?凱，等. 含淹沒剛性植被明渠水流特性的試驗(yàn)研究[J]. 水力發(fā)電學(xué)報(bào)，2019，38(4)：21-32.

Song Yingting，Jing Hefang，Zhang Kai，et al. Experimental study on characteristics of open channel flows with submerged rigid vegetation[J]. Journal of Hydropower，2019，38(4)：21-32(in Chinese).

［26］ Sundin J，Bagheri S. Interaction between hairy surfaces and turbulence for different surface time scales[J]. Journal of Fluid Mechanics，2019，861：556-584.

［27］ 王偉杰. 明渠植被水流流速分布解析解與阻力特性研究[D]. 武漢：武漢大學(xué)，2016.

Wang Weijie. Analytical Model of Velocity Distribution and Resistance Characteristics of Open Channel Flow within Vegetation[D]. Wuhan：Wuhan University，2016(in Chinese).

［28］ 張?凱. 含柔性植被明渠水流實(shí)測(cè)及數(shù)值模擬研究[D]. 銀川：北方民族大學(xué)，2018.

Zhang Kai. Research on Numerical Simulation Measurement of Open Channel Flow with Flexible Vegetation[D]. Yinchuan：North Minzu University，2018(in Chinese).

［29］ Chen H，Zou Q. Eulerian-Lagrangian flow-vegetation interaction model using immersed boundary method and OpenFOAM[J]. Advances in Water Resources，2019，126：176-192.

［30］ Ji C，Munjiza A，Williams J J R. A novel iterative direct forcing immersed boundary method and its finite volume applications[J]. Journal of Computational Physics，2012，231(4)：1797-1821.

［31］ Uhlmann M. An immersed boundary method with direct forcing for the simulation of particulate flows[J]. Journal of Computational Physics，2005，209(2)：448-476.

［32］ Munjiza A，Owen D R J，Bicanic N. A combined finite-discrete element method in transient dynamics of fracturing solids[J]. Engineering Computation，1995，12(2)：145-174.

［33］ Zhang B，Xu D，Ji C，et al. Numerical investigation on the incipient motion of non-spherical sediment particles in viscous sublayer of open channel flows[C]//Proceeding of 14th International Symposium on River Sedimentation. Chengdu，China，2019：16-19.

［34］ 及春寧，劉?爽，楊立紅，等. 基于嵌入式迭代的高精度浸入邊界法[J]. 天津大學(xué)學(xué)報(bào)(自然科學(xué)與工程技術(shù)版)，2014，47(5)：377-382.

Ji Chunning，Liu Shuang，Yang Lihong，et al. An accurate immersed boundary method based on built-in iterations[J]. Journal of Tianjin University(Science and Technology)，2014，47(5)：377-382(in Chinese).

［35］ Poggi D，Porporato A，Ridolfi L，et al. The effect of vegetation density on canopy sub-layer turbulence[J]. Boundary-Layer Meteorology，2004，111(3)：565-587.

Large Eddy Simulation of Submerged Flexible Vegetation Fluid-Solid Coupling in Open Channel Turbulent Flow

Xu Dong，Liang Xinya，Zhang Bingchang，Ji Chunning

(State Key Laboratory of Hydraulic Engineering Simulation and Safety，Tianjin University，Tianjin 300072，China)

Aquatic vegetation is widely found in natural rivers，lakes，and wetlands，and is an important component of river ecosystems. It increases the influence of channel flow resistance by changing the time-averaged and turbulent structure of open channel flow，and bending and oscillating under the action of water flow，forming a complex flow-solid coupling effect. Studying the motion response of flexible vegetation under the action of water flow is helpful to analyze the water blocking effect of vegetation more accurately. This paper used the high-performance three-dimensional turbulence calculation program CgLes based on the built-in iterative immersed boundary method to simulate the flow-solid coupling large eddy，studied the bending and motion response of flexible vegetation under the action of water flow，analyzed the characteristics of vegetation vortex-induced vibration and the law of water flow resistance. The simulation results show that the vegetation movement is dominated by spanwise oscillation. The stiffness has a significant impact on the vegetation movement track and swing period，thus affects the water blocking laws of vegetation. Under constant hydrodynamic conditions，the swing period of the vegetation is related to elastic modulus of the vegetation and is not affected by the flow. As the stiffness increases，the vegetation motion changes from quasi-period to period-doubling and finally reaches into the static motion. With the large vegetation stiffness and the small spacing，the reverse-flow bending phenomenon of vegetation may occur. When the stiffness is large，the flexible vegetation swings with the flow，which inhibits the development of flow velocity. However，when the stiffness is small enough，the vegetation is close to lodging，and the average velocity of the open channel increases compared with it in the rigid vegetation condition. The maximum effect of vegetation stiffness on the average velocity is about 2.47% in this paper. Vegetation spacing has little effect on the bending amplitude of vegetation，but significantly affects the spanwise swing. The average velocity decreases with the increase of vegetation spacing，and the vegetation spacing mainly affects the flow velocity in the vegetation layer. In this paper，the maximum effect of vegetation spacing on the average velocity is about 8.19%.

flexible vegetation；vortex-induced vibration；immersed boundary method；fluid-solid coupling；numerical simulation

TV133.2

A

0493-2137(2022)03-0313-09

10.11784/tdxbz202104027

2021-04-15；

2021-07-01.

許?棟（1980—??），男，博士，副教授，xudong@tju.edu.cn.

及春寧，cnji@tju.edu.cn.

國(guó)家重點(diǎn)研發(fā)計(jì)劃重點(diǎn)專項(xiàng)資助項(xiàng)目(2018YFC0407505)；國(guó)家自然科學(xué)基金資助項(xiàng)目(5217090155，51979186，51779175).

Supported by the National Key Research and Development Program(No.2018YFC0407505)，the National Natural Science Foundation of China(No.5217090155，No.51979186，No.51779175).

(責(zé)任編輯：樊素英)