張芝永，謝詠哲，陳 紅，潘存鴻

(1.浙江省水利河口研究院(浙江省海洋規(guī)劃設(shè)計(jì)研究院)，浙江 杭州 310020；2.浙江省河口海岸重點(diǎn)實(shí)驗(yàn)室，浙江 杭州 310020；3.河海大學(xué)水利水電學(xué)院，江蘇 南京 210098)

涌潮是外海潮波進(jìn)入逐漸縮窄的寬淺河口后發(fā)生劇烈變形而引起的一種特殊水力學(xué)現(xiàn)象[1]，主要發(fā)生在喇叭形的入海河口。我國(guó)錢塘江涌潮是世界著名的涌潮景觀，其涌潮強(qiáng)度及氣勢(shì)世界罕有。涌潮類似于移動(dòng)的水躍，其水流形態(tài)及流場(chǎng)結(jié)構(gòu)極為復(fù)雜，極強(qiáng)的水流沖擊導(dǎo)致海塘、丁壩等[2]涉水工程破壞嚴(yán)重，同時(shí)也會(huì)導(dǎo)致河床的劇烈沖刷和大量泥沙的懸浮，對(duì)河床沖淤影響極大[3]。因此研究涌潮的局部水動(dòng)力特性及其與結(jié)構(gòu)物、河床的作用機(jī)理有助于為涌潮河段的結(jié)構(gòu)物及河床沖刷防護(hù)方案的實(shí)施提供技術(shù)支撐。

目前對(duì)于涌潮的研究普遍集中于大尺度方面，如對(duì)涌潮在整個(gè)河口區(qū)傳播規(guī)律的研究：潘存鴻等[4]基于KFVS格式開發(fā)了第一個(gè)可以精準(zhǔn)模擬涌潮的大范圍二維數(shù)值計(jì)算模型，謝東風(fēng)等[5-7]先后應(yīng)用FVCOM對(duì)錢塘江涌潮進(jìn)行了三維數(shù)值模擬。以上學(xué)者在數(shù)值模型中均是利用外海邊界的規(guī)則潮波進(jìn)入逐漸束窄河段后發(fā)生潮波變形而自然生成的涌潮，但其模型均因采用基于淺水方程簡(jiǎn)化的二維或三維數(shù)值模型而無(wú)法解決局部建筑物的涌潮作用力及局部沖刷等問(wèn)題。

對(duì)于涌潮與局部建筑物及河床相互作用問(wèn)題的研究，需要借助基于Navier-Stokes方程的CFD數(shù)值模型或涌潮試驗(yàn)水槽。在數(shù)學(xué)模擬或試驗(yàn)?zāi)M過(guò)程中，由于計(jì)算資源或試驗(yàn)?zāi)Ｐ涂臻g的限制，無(wú)法構(gòu)建足夠長(zhǎng)的束窄過(guò)渡段，只能通過(guò)規(guī)則潮波的自然變形來(lái)生成涌潮。為此，如何在邊界處直接生成涌潮成為許多涌潮學(xué)者關(guān)注的問(wèn)題。近年來(lái)，隨著水泵測(cè)控、閘門控制技術(shù)的進(jìn)步，一些學(xué)者相繼提出了在邊界處直接構(gòu)造涌潮的生潮方法，使得涌潮局部尺度問(wèn)題的研究成為可能。如黃靜等[8-9]通過(guò)調(diào)節(jié)水槽末端水泵頻率使流量發(fā)生突變的方法在實(shí)驗(yàn)室內(nèi)實(shí)現(xiàn)了涌潮的準(zhǔn)確模擬，但因流量突變幅度與涌潮高度的相關(guān)性未知，在試驗(yàn)中需要多次調(diào)試才能得到預(yù)期的涌潮水流。而Chanson等[10-12]則基于落閘法得到了反向涌潮波，但其生成的涌潮的底部水流與表層水流運(yùn)動(dòng)方向完全相反，這與河口區(qū)涌潮傳播過(guò)程中的流態(tài)結(jié)構(gòu)完全不同。數(shù)值模擬方面，已有CFD數(shù)值模擬研究中常將涌潮概化為無(wú)漲潮流速的潰壩波[13]，但實(shí)際上涌潮除了涌潮高度外還伴隨有一定的漲潮流速，且漲潮流速不能隨意給定，它是與涌潮高度、潮前水深和落潮流速有關(guān)的參數(shù)。戎貴文等[14-16]直接在邊界處設(shè)置實(shí)測(cè)水位及流速值進(jìn)行涌潮的模擬，取得了較好的效果，但該方法只能針對(duì)特定的涌潮工況。戚藍(lán)等[17-18]借鑒潘存鴻提出的涌潮條件基本方程，推導(dǎo)出與涌潮高度、潮前水深、落潮流速有關(guān)的涌潮漲潮流速計(jì)算方法，并應(yīng)用數(shù)值模型分別進(jìn)行了同時(shí)考慮水深和漲潮流速的涌潮模擬，與理論結(jié)果進(jìn)行對(duì)比驗(yàn)證了邊界生潮法的可靠性。

綜上所述，目前物理模型試驗(yàn)和數(shù)值計(jì)算模型中成熟的涌潮生成方式有所不同，物理模型試驗(yàn)普遍采用流量生潮方法，而數(shù)值計(jì)算模型則是采用邊界生潮方法。因無(wú)法確定流量生潮法和邊界生潮法生成的涌潮水動(dòng)力特性是否存在區(qū)別及數(shù)值模型中湍流模型、網(wǎng)格尺度等參數(shù)對(duì)涌潮精確模擬的影響，有必要對(duì)比現(xiàn)有涌潮生成方法的差別及其內(nèi)在聯(lián)系，探討網(wǎng)格尺度、湍流模型選擇對(duì)涌潮模擬的影響，提出適用于涌潮模擬的生潮方法、網(wǎng)格劃分策略及湍流模型選擇，進(jìn)而為后續(xù)涌潮對(duì)局部結(jié)構(gòu)物及河床沖刷影響的研究提供關(guān)鍵的數(shù)值模擬工具。

1 涌潮特性介紹

在涌潮傳播過(guò)程中，涌潮以一定的水位差H和流速向上游傳播，具體形式見圖1，其行進(jìn)方向的前方稱為潮前，反之稱為潮后。 圖中潮頭行進(jìn)的速度c稱為涌潮傳播速度。潮前水深是涌潮到來(lái)之前河道(水槽)中的初始水深h0，此時(shí)水體可能是靜止的，也可能存在一定的潮前流速(或落潮流速)，即涌潮到來(lái)之前水體的流速v0，一般與涌潮傳播速度相反，代表了徑流在涌潮現(xiàn)象中的作用；h1為涌潮經(jīng)過(guò)之后的水深即潮后水深，或稱為漲潮水深，潮后水深與潮前水深的差值為涌潮高度H，v1為漲潮時(shí)的流速，即潮后流速(或漲潮流速)。

圖1 涌潮潮頭結(jié)構(gòu)Fig.1 Tidal bore structure

根據(jù)涌潮形態(tài)及涌潮波破碎情況的不同，涌潮可分為3種：①潮頭未發(fā)生破碎的波狀涌潮；②潮頭破碎強(qiáng)烈的強(qiáng)旋滾涌潮；③介于兩者之間的稱為弱旋滾涌潮[19]。

涌潮可視為一種移動(dòng)的水躍。通過(guò)空間相對(duì)速度變換并聯(lián)立一維連續(xù)性方程和動(dòng)量方程，可得到涌潮各要素滿足如下方程：

(1)

(2)

由式(1)(2)可知，如已知v0、h0、h1(或H)，則可得c，繼而由c可推求v1。

2 數(shù) 值 模 型

2.1 控制方程

局部尺度涌潮數(shù)值模擬的控制方程為連續(xù)性方程和動(dòng)量方程(Navier-Stokes方程)。

2.2 自由液面追蹤方法

自由面的捕捉采用高精度的VOF方法。VOF法構(gòu)建了一種體積函數(shù)以追蹤單元網(wǎng)格中的流體體積，通過(guò)其函數(shù)值及導(dǎo)數(shù)值確定自由面位置，適用于不可壓縮、黏性、瞬變流體的問(wèn)題。

2.3 湍流模型

選擇較為常見的標(biāo)準(zhǔn)k-ε和RNGk-ε兩種湍流模型來(lái)比較各模型模擬涌潮的效果。標(biāo)準(zhǔn)k-ε模型由k方程(即湍流動(dòng)能方程)與ε方程(耗散方程)組成，因簡(jiǎn)單高效而應(yīng)用廣泛。RNGk-ε模型全稱為重整規(guī)化群(RNG)k-ε模型[17]，該模型修改了耗散方程并考慮了湍流渦旋，精度有所提高。重整規(guī)化群理論為k-ε方程提供了湍流Prandtl數(shù)，而標(biāo)準(zhǔn)k-ε模型中這一參數(shù)取經(jīng)驗(yàn)值。RNGk-ε模型還增強(qiáng)了模擬低雷諾數(shù)流動(dòng)的能力，使其具有更高的泛用性與可靠性，能夠更好地模擬大曲率流動(dòng)問(wèn)題。

2.4 數(shù)值水槽中生潮邊界構(gòu)造

為了檢驗(yàn)流量生潮法、邊界生潮法的差異及內(nèi)在聯(lián)系，構(gòu)建長(zhǎng)50 m、高0.5 m的數(shù)值水槽，以便觀察涌潮自生潮邊界出現(xiàn)后向上游行進(jìn)、發(fā)展的過(guò)程。邊界生潮法的邊界設(shè)置情況如圖1所示。邊界需同時(shí)給定流速和水深條件，其中水深為潮前水深疊加涌潮高度后的潮后水深h1，流速為v1，在給定v0、h0及H的情況下，可通過(guò)式(1)(2)求解得到。對(duì)于流量生成法，數(shù)值模擬中生潮邊界的構(gòu)造參考了模型試驗(yàn)的布置情況，在長(zhǎng)水槽下游增加一段長(zhǎng)度和深度均為1 m的生潮前室，前室下游側(cè)為固壁，前室底部為流量進(jìn)口邊界，對(duì)應(yīng)特定涌潮工況，生潮流量的大小需在模擬過(guò)程中經(jīng)不斷調(diào)試得到。

2.5 數(shù)值模擬計(jì)算方案

針對(duì)生潮方法、湍流模型及網(wǎng)格尺度對(duì)涌潮模擬的影響，設(shè)計(jì)15組數(shù)值模擬方案，見表1。數(shù)值模擬中潮前水深保持為0.15 m時(shí)，涌潮高度分別取0.05 m、0.10 m和0.15 m，其對(duì)應(yīng)的涌潮形態(tài)分別為波狀涌潮、弱漩滾涌潮和強(qiáng)漩滾涌潮。網(wǎng)格尺度影響方面，在水位波動(dòng)影響區(qū)的網(wǎng)格尺寸d(垂向和縱向網(wǎng)格尺寸一致)分別取0.050H、0.100H、0.125H和0.250H。

表1 數(shù)值模擬研究方案

3 模擬結(jié)果及分析

3.1 生潮方法比較

表2 不同生潮方法對(duì)應(yīng)的邊界條件參數(shù)

在數(shù)值模擬中發(fā)現(xiàn)，對(duì)于邊界生潮法，給定合理的邊界水深和漲潮流速后模擬得到的涌潮高度基本與邊界處一致。而流量生潮法則需要通過(guò)不斷調(diào)試進(jìn)口流量[20]，才能得到與邊界生潮法一致的涌潮高度。表2為應(yīng)用兩種生潮方法時(shí)各涌潮高度對(duì)應(yīng)的邊界條件值，其中v1由式(1)(2)聯(lián)立求解推求出。

針對(duì)涌潮漲潮流速、水面線形態(tài)及涌潮傳播速度對(duì)兩種生潮方法模擬得到的涌潮進(jìn)行對(duì)比，以x=15 m處為監(jiān)測(cè)斷面，選取兩個(gè)典型時(shí)刻：①潮頭到達(dá)監(jiān)測(cè)斷面使縱向流速達(dá)到極值時(shí)，對(duì)于波狀涌潮即首個(gè)波峰，②潮頭經(jīng)過(guò)監(jiān)測(cè)斷面，對(duì)于波狀涌潮即首個(gè)波谷，對(duì)各方法生成的涌潮潮頭形態(tài)、涌潮內(nèi)部縱向流速及傳播速度進(jìn)行對(duì)比。圖2為波狀涌潮(H=0.05 m)及強(qiáng)旋滾涌潮(H=0.15 m)在兩個(gè)典型時(shí)刻的水面線情況。由圖2可見，對(duì)于波狀涌潮，流量生潮法和邊界生潮法生成的波峰、波谷高度基本一致，但波長(zhǎng)略有不同，對(duì)比來(lái)看，流量生潮法得到的波長(zhǎng)要大于邊界生潮法。而對(duì)于強(qiáng)旋滾涌潮來(lái)說(shuō)，兩種方法生成的涌潮潮頭水面線非常接近。從流速最大時(shí)刻來(lái)看，波狀涌潮流速最大值發(fā)生在首峰時(shí)刻，而旋滾涌潮則發(fā)生在涌潮升高的過(guò)程中。

圖3對(duì)應(yīng)圖2，為不同生潮方法下x=15 m斷面不同典型時(shí)刻的縱向流速u垂向分布情況。由圖3可見，當(dāng)波狀涌潮首鋒及旋滾涌潮縱向流速最大時(shí)，縱向流速垂向分布可分為兩個(gè)區(qū)域：潮前水深(0.15 m)以上區(qū)域，流速迅速增大，并在自由表面處達(dá)到極值；潮前水深以下區(qū)域縱向流速則略有減小，其流速變化率明顯小于潮前水深以上部分。整體來(lái)看，此時(shí)涌潮縱向流速垂向分布呈Γ型，各類涌潮均呈現(xiàn)相同規(guī)律。從流速大小看，涌潮高度越大，其自由表面處流速及潮前水深以下區(qū)域的平均流速相應(yīng)也會(huì)增大。而在涌潮經(jīng)過(guò)x=15 m后，潮前水深以上區(qū)域縱向流速有所減小。對(duì)于波狀涌潮工況，縱向流速僅是其最大縱向流速的1/4。由于兩種方法生成的潮頭波長(zhǎng)有所區(qū)別，波谷時(shí)刻流速也有所差異，差值在0.1 m/s左右。而對(duì)于旋滾涌潮，潮頭經(jīng)過(guò)后的斷面最大縱向流速減小幅度約40%。對(duì)比圖3中兩種生潮方法的流速垂向分布情況，除了波狀涌潮波谷時(shí)刻流速分布略有差異外，其他時(shí)刻兩種方法生成的涌潮內(nèi)部流速結(jié)構(gòu)基本接近。

圖3 不同生潮方法典型時(shí)刻縱向流速垂向分布Fig.3 Vertical distribution of longitudinal velocity using different bore generation methods

利用涌潮自x=10 m傳播至x=15 m位置的時(shí)間，計(jì)算得到各組次涌潮傳播速度[21]，并與式(1)或式(2)計(jì)算得到的理論值進(jìn)行對(duì)比。由圖4可見，涌潮高度0.05 m、0.10 m和0.15 m對(duì)應(yīng)的傳播速度理論值分別是1.513 m/s、1.808 m/s和2.101 m/s，涌潮傳播速度隨著涌潮高度的增大而增大。邊界生潮法和流量生潮法模擬得到的涌潮傳播速度與理論值之間的平均相差幅度分別為0.8%和2.3%，基本接近。針對(duì)涌潮高度分別為0.05 m、0.10 m和0.15 m的3種涌潮工況在水槽中開展模型試驗(yàn)，測(cè)量各工況下的涌潮傳播速度，其與數(shù)值模擬結(jié)果的對(duì)比結(jié)果見圖4?？傮w來(lái)說(shuō)，水槽試驗(yàn)結(jié)果略小于理論值和數(shù)值模擬結(jié)果，其原因主要是水槽試驗(yàn)中流量突變是依靠水泵頻率的急速提升實(shí)現(xiàn)的，而水泵造成的流量抬升在實(shí)際情況中必將經(jīng)歷一個(gè)短時(shí)的過(guò)渡過(guò)程，并非數(shù)學(xué)模擬中生潮流量的瞬間變化，這造成了涌潮動(dòng)力的減弱，因此水槽試驗(yàn)得到的涌潮傳播速度相對(duì)略小。盡管如此，水槽試驗(yàn)結(jié)果與理論值的偏差基本在10%以內(nèi)，模擬效果較好。

圖4 涌潮傳播速度對(duì)比Fig.4 Comparisons of bore propagation velocity

綜上所述，流量生潮法和邊界生潮法生成的涌潮在涌潮形態(tài)、水流流速及傳播速度上均較為接近，說(shuō)明兩者在適當(dāng)邊界條件下可生成特征相同的涌潮水流。而在相同涌潮條件下，流速分布及水面線形態(tài)的一致性表明，兩種方法在單位時(shí)間內(nèi)的過(guò)水流量應(yīng)是一致的。對(duì)比表2發(fā)現(xiàn)，應(yīng)用流量生潮法和邊界生潮法模擬3種涌潮工況時(shí)生潮邊界處的單寬流量基本相同，因此確定流量生潮法的流量邊界條件時(shí)可利用邊界生潮法的計(jì)算公式。即：

q=v1h1

(3)

式(3)中v1由式(1)(2)確定。故采用流量生潮法時(shí)可采用式(3)估算生潮流量，從而避免了試驗(yàn)中需要不斷調(diào)試流量的問(wèn)題。值得注意的是，在水槽試驗(yàn)中水泵頻率突變后流量提升所需要的時(shí)間也是影響涌潮生成質(zhì)量的一個(gè)主要因素，過(guò)長(zhǎng)的流量提升時(shí)間將會(huì)減弱生成涌潮的強(qiáng)度。因此物理模型試驗(yàn)中必須配置性能優(yōu)、反應(yīng)速度快的變頻水泵才能更好地復(fù)演涌潮過(guò)程。

3.2 網(wǎng)格尺度影響

對(duì)于數(shù)值模擬計(jì)算來(lái)說(shuō)，模型網(wǎng)格尺寸對(duì)數(shù)值模擬結(jié)果同樣會(huì)有較大影響。適當(dāng)、精細(xì)的網(wǎng)格能得到精確的結(jié)果，而另一方面網(wǎng)格尺寸將直接影響計(jì)算時(shí)間，其一個(gè)數(shù)量級(jí)的變化會(huì)導(dǎo)致計(jì)算時(shí)間大幅增加。 對(duì)于涌潮這一強(qiáng)間斷水流來(lái)說(shuō)，如何選擇合適的網(wǎng)格尺寸以準(zhǔn)確刻畫其潮頭特征是數(shù)值模擬中較為重要的問(wèn)題。為此選用不同的相對(duì)網(wǎng)格尺寸進(jìn)行研究，取各級(jí)網(wǎng)格尺寸d分別為H的0.050、0.100、0.125和0.250倍。圖5為不同網(wǎng)格尺度下旋滾涌潮潮頭水面線特征。由圖5可以看出，各網(wǎng)格尺寸下涌潮高度是基本接近的，但在潮頭水面坡度線(陡度)方面則略有差別，具體表現(xiàn)為網(wǎng)格尺度越小，涌潮潮頭刻畫越為精細(xì)，其潮頭平均坡度會(huì)越陡，在網(wǎng)格尺寸小于0.100H時(shí)，模擬結(jié)果變化基本很小。

圖5 各網(wǎng)格尺寸下的漩滾涌潮潮頭水面線Fig.5 Water surface profiles of breaking tidal bore under different mesh sizes

3.3 湍流模型的影響

對(duì)于涌潮這一強(qiáng)間斷流來(lái)說(shuō)，湍流模型的選擇至關(guān)重要。為此針對(duì)不同涌潮，分別應(yīng)用標(biāo)準(zhǔn)k-ε模型和RNGk-ε模型進(jìn)行模擬分析。圖6對(duì)比了不同湍流模型下兩種涌潮水面線的模擬情況。總體來(lái)看兩種湍流模型得到的涌潮水面線基本接近，相比而言，采用標(biāo)準(zhǔn)k-ε模型的模擬結(jié)果其潮頭水位略低于RNGk-ε模型模擬結(jié)果。而從旋滾涌潮(H=0.15 m)工況中還可以看出，標(biāo)準(zhǔn)k-ε模型的水面線坡度相對(duì)RNGk-ε模型的模擬結(jié)果更緩，說(shuō)明RNGk-ε模型更適用于模擬高強(qiáng)度、高陡度的涌潮。

圖6 各湍流模型下的涌潮潮頭形態(tài)Fig.6 Tidal bore profile under different turbulent models

3.4 原型尺度涌潮的數(shù)值模擬驗(yàn)證

為了驗(yàn)證本文數(shù)值模擬方法是否可以用于原型尺度的涌潮模擬，利用2010年海寧鹽官涌潮觀測(cè)資料進(jìn)行模型驗(yàn)證。2010年實(shí)測(cè)涌潮特征如表3所示。

表3 海寧鹽官2010年實(shí)測(cè)涌潮資料

利用垂向二維涌潮數(shù)值水槽，根據(jù)涌潮模擬策略對(duì)表3中的3種實(shí)測(cè)大、中及小涌潮進(jìn)行了模擬。圖7為小潮及大潮情況下涌潮傳播過(guò)程中縱向流速分布情況。由圖7可見，對(duì)于涌潮高度為0.8 m的小潮，由于涌潮高度較小，水動(dòng)力強(qiáng)度稍弱，其涌潮潮頭呈現(xiàn)波狀涌潮特征；而涌潮高度為2.74 m的大潮潮頭旋滾非常強(qiáng)烈，數(shù)值模擬得到的涌潮形態(tài)與實(shí)際的涌潮形態(tài)基本一致。圖8對(duì)比了數(shù)值模擬涌潮與實(shí)測(cè)涌潮的傳播速度，可見原型尺度與實(shí)驗(yàn)室尺度下的涌潮傳播速度規(guī)律相同，涌潮高度越大，傳播速度越快，而且數(shù)值模擬得到的涌潮傳播速度與實(shí)測(cè)值基本接近，平均偏差在9%以內(nèi)。

圖7 涌潮縱向流速及形態(tài)數(shù)值模擬結(jié)果Fig.7 Contours of longitudinal velocity and bore shapes from numerical simulation

圖8 涌潮傳播速度實(shí)測(cè)值與數(shù)值模擬結(jié)果對(duì)比Fig.8 Comparisons of numerical results and field measured results for propagation velocity

綜上所述，通過(guò)對(duì)比涌潮形態(tài)及傳播速度，證明了本文數(shù)值模擬方法在模擬原型尺度局部范圍涌潮時(shí)的可靠性。但值得注意的是，本文邊界生潮公式(1)和(2)是在斷面規(guī)則的平直河道(水槽)情況下，基于一維連續(xù)性和動(dòng)量守恒方程推導(dǎo)得出的，公式并未考慮地形、岸線變化對(duì)涌潮水流的影響，而實(shí)際情況中一些河段的地形、岸線變化極大，針對(duì)復(fù)雜河段內(nèi)的涌潮水流運(yùn)動(dòng)進(jìn)行CFD數(shù)值模擬時(shí)，建議以實(shí)測(cè)的潮前水深、漲潮流速、涌潮高度賦予模型邊界條件。

4 結(jié) 論

a.適當(dāng)邊界條件下流量生潮法和邊界生潮法可模擬得到各項(xiàng)特質(zhì)均較為一致的涌潮水流，其中旋滾涌潮的模擬結(jié)果最為接近。兩種生潮方法在邊界構(gòu)造上的不同導(dǎo)致模擬得到的波狀涌潮首峰之后水面線形態(tài)及流速分布略有差異。在傳播速度上，邊界生潮法模擬結(jié)果更接近于理論值。通過(guò)對(duì)比發(fā)現(xiàn)，流量生潮法單寬流量與邊界生潮法在生潮邊界處的計(jì)算流量一致，因此可通過(guò)邊界生潮法公式估算流量生潮法的邊界條件，這為物理模型試驗(yàn)中涌潮生潮流量的確定提供了參考依據(jù)。

b.網(wǎng)格尺度越小，模擬得到的涌潮平均陡度越陡，當(dāng)網(wǎng)格尺寸不大于0.100H時(shí)，網(wǎng)格尺度對(duì)計(jì)算結(jié)果的影響較小，因此進(jìn)行涌潮數(shù)值模擬時(shí)建議網(wǎng)格尺寸不應(yīng)大于0.100H。

c.相比RNGk-ε模型，標(biāo)準(zhǔn)k-ε模型模擬得到的潮頭高度略低，潮頭陡度稍緩。對(duì)于高陡度及高卷曲的涌潮水流來(lái)說(shuō)，RNGk-ε模型更適合于涌潮的CFD數(shù)值模擬。

d.本文涌潮局部數(shù)值模擬方法適用于岸線及地形基本穩(wěn)定的河段內(nèi)的涌潮水流。針對(duì)地形、岸線變化較大的河段，建議以實(shí)測(cè)的潮前水深、漲潮流速、涌潮高度等涌潮參數(shù)作為邊界條件進(jìn)行涌潮CFD數(shù)值模擬。