非靜壓水波模型通常被視為對淺水長波方程所滿足的靜壓模型的一種修正,在經(jīng)典水波理論及其數(shù)學模型的論述中鮮有提及。本書內(nèi)容可概括為三部分,首先從水波運動的理論角度,梳理非靜壓水波模型相關(guān)的理論基礎(chǔ),辨明其與淺水長波的靜壓模型的差異性。其次,針對非靜壓水波模型通常采用的預(yù)估校正數(shù)值方法,闡述非靜壓模型與淺水長波的靜壓模型在數(shù)值求解過程中的內(nèi)在關(guān)聯(lián)。后從應(yīng)用角度闡述及驗證了非靜壓模型的適用性,其不僅適用于海岸帶的波浪運動的數(shù)值模擬,也適用于復雜邊界約束條件下帶自由表面的水流運動、流固耦合運動等的數(shù)值模擬。
本書力圖闡明非靜壓水波模型的理論基礎(chǔ)、驗證數(shù)值方法的可靠性、分析模型的適用性,希望為非靜壓水波模型的拓廣應(yīng)用添磚加瓦。
水波理論以其發(fā)展歷史為脈絡(luò),經(jīng)歷了線性水波、非線性水波研究的進階之路。水波問題的研究成果豐碩,水波理論和模型多樣。各種水波模型的差異性來自于不同的建模理論基礎(chǔ),而不同的理論框架限制著不同水波模型的適用性。非靜壓模型是靜壓模型的推廣,即舍棄了靜壓假定這一限制條件。非靜壓模型從提出伊始,即作為靜壓模型應(yīng)用范圍的拓廣,主要應(yīng)用在水波運動的數(shù)值模擬,如海岸帶波浪演化等問題的研究。本書結(jié)合作者多年來關(guān)于非靜壓模型開發(fā)的研究工作,系統(tǒng)地總結(jié)了非靜壓模型20余年的發(fā)展,梳理了理論基礎(chǔ),概括了模型的適用條件,主要內(nèi)容包括:水波模型簡介,非靜壓水波模型的理論基礎(chǔ)、數(shù)值求解、校驗、數(shù)值求解問題分析和在高精度數(shù)值模擬中的應(yīng)用等方面,以期為相關(guān)研究者提供參考。
綠樹掩映下,微風吹拂湖面,生成片片漣漪,熱帶氣旋掠過海面,引起了驚濤駭浪,雖然感觀各異,但這些水流運動均可歸結(jié)為帶自由表面的水波運動。水波運動的研究歷史悠長,吸引了眾多數(shù)學家、物理學家、力學家為之奮斗。水波運動的研究也助推了某些相關(guān)科學問題的研究,其中著名的當屬孤立波的研究。孤立波的研究首先來源于水波運動的觀察及探究,繼而推動了非線性科學的發(fā)展,即孤立子理論。孤立子起源于孤立波,并已在非線性光學、生物學、等離子體及光纖孤立子通信等一系列高科技領(lǐng)域有了令人矚目的應(yīng)用。
水波理論以其發(fā)展歷史為脈絡(luò),經(jīng)歷了線性水波、非線性水波研究的進階之路。水波問題的研究成果豐碩,其中不乏經(jīng)典著作。水波理論多種多樣,理論模型各異,往往使初學者眼花繚亂。分析各水波模型的差異性,大體上表現(xiàn)在兩個方面,其一為對波動非線性特征的描述精度,其二為對波動色散性特征的描述精度。水波各模型的差異性也受到當時科技水平的限制,包括數(shù)學分析手段、計算技術(shù)的先進性等。水波模型的發(fā)展及應(yīng)用除受基礎(chǔ)理論的發(fā)展水平限制以外,為滿足實際工程應(yīng)用的需求,水波模型的提出及發(fā)展還具有針對性。工程需求往往帶有針對性、時效性,在認知水平有限、技術(shù)條件掣肘的情況下,為了適應(yīng)急迫的工程需求,研究者發(fā)展出的水波模型也具有針對性。因此,各類水波模型往往需要注重計算效率與計算精度的有機結(jié)合。
淺水長波模型因其特定的理論基礎(chǔ)及廣闊的適用條件,在眾多水波模型中受到了廣泛的關(guān)注及應(yīng)用。淺水長波模型的缺陷明顯,為提高模型的模擬能力,眾多研究者提出了若干修正模型,其中著名的當屬布西內(nèi)斯克(Boussinesq)方程。同樣以提高淺水長波模型對波動色散性的模擬能力為目標,非靜壓模型的建模思路不同于布西內(nèi)斯克方程,它理論上屬于納維斯托克斯(Navier-Stokes)方程,數(shù)值求解采用壓強分裂技術(shù),即將流動壓強分為靜壓和非靜壓兩部分。非靜壓模型的提出是繼靜壓模型之后,其較早的數(shù)值求解方法借鑒了靜壓模型的相關(guān)技術(shù),故稱其為非靜壓模型,凸顯了其模型特點及其發(fā)展歷程。
非靜壓模型主要是數(shù)值求解過程,并不能稱其為一種嚴格的理論模型,稱為水波的數(shù)值求解模式似乎更為恰當。非靜壓模型的理論基礎(chǔ)是納維斯托克斯方程,較之傳統(tǒng)的納維斯托克斯方程的數(shù)值求解,差別僅在于對壓強的分裂?倝簭姺至殉龅撵o壓成分具有特定的物理含義,而不僅僅是計算過程中的預(yù)估流場值。同時,在自由表面捕捉等方面,借助垂向坐標變換等技術(shù),計算效率有所提高。非靜壓模型雖然起源于靜壓水波模型,但其應(yīng)用范圍可以拓展至靜壓假定不再成立的任何自由表面水流運動的數(shù)值模擬,如地形陡變條件下的明渠流動、水流繞結(jié)構(gòu)物流動、水下射流、變密度流動等的數(shù)值模擬。總體而言,非靜壓模型屬于一種小眾的水波模型,有其特殊的適用性。合理發(fā)揮模型的長處,可為海洋、海岸、河口、河流、湖泊等帶自由表面水流運動的研究提供一種數(shù)值模擬的技術(shù)手段。
本書旨在盡可能系統(tǒng)地總結(jié)非靜壓模型二十余年的發(fā)展,梳理其理論基礎(chǔ),概括模型的適用條件,明確其優(yōu)勢所在,并結(jié)合筆者多年來關(guān)于非靜壓模型開發(fā)的研究經(jīng)驗,將其形諸筆端。相關(guān)研究工作受到了如下項目的資助:國家自然科學基金海上風力機支撐結(jié)構(gòu)局部沖刷的動力機制及基于流動控制的防沖刷技術(shù)研究;國家自然科學基金重點項目海嘯力學及其在南中國海的應(yīng)用;國家重點研發(fā)計劃超大型海上風電機組整機、部件、基礎(chǔ)一體化設(shè)計技術(shù)。本書的出版得到了上海交通大學劉樺教授的大力支持,筆者由衷地表示感謝。多年來筆者所指導的研究生們也在相關(guān)工作中付出了辛勤的汗水,在此表示真摯的謝意。
本書并非水波理論基礎(chǔ)的論著,僅限于一類帶自由表面波動的數(shù)值模擬方法的研究。若能以綿薄之力為相關(guān)研究工作添磚加瓦,將甚感欣慰。鑒于筆者學識有限,認識及實踐工作中難免存在偏頗之處,敬請讀者指正。
第1章 緒論
1.1 水波模型簡介
1.2 靜壓水波模型
1.3 非靜壓水波模型
1.4 小結(jié)
第2章 非靜壓水波模型的理論基礎(chǔ)
2.1 流體中的壓力
2.2 水波理論簡述
2.2.1 波動方程
2.2.2 微幅波簡述
2.3 淺水長波
2.3.1 淺水長波特性
2.3.2 淺水長波控制方程
2.4 Boussinesq方程
2.5 非靜壓模型的描述
2.5.1 納維斯托克斯方程及湍流模擬概述
2.5.2 非靜壓模型的建立
2.5.3 動壓項引入的色散效應(yīng)
2.6 非靜壓模型的數(shù)值求解
2.6.1 控制方程組的封閉
2.6.2 壓力求解的迭代修正法
2.6.3 非靜壓模型數(shù)值求解的預(yù)估校正法
第3章 非靜壓模型的數(shù)值求解
3.1 自由表面捕捉方法
3.1.1 自由表面捕捉的MAC數(shù)值方法
3.1.2 自由表面捕捉的VOF法
3.1.3 自由表面捕捉的間接模擬方法
3.2 坐標變換
3.2.1 變動計算域的模擬方法簡述
3.2.2 坐標變換
3.2.3 坐標變換的討論
3.2.4 坐標系中水平擴散項的計算
3.3 非靜壓模型的數(shù)值求解
3.3.1 計算網(wǎng)格及變量布置
3.3.2 控制體界面的局部坐標系及其度量
3.3.3 非結(jié)構(gòu)網(wǎng)格基礎(chǔ)上的散度計算
3.3.4 非結(jié)構(gòu)網(wǎng)格的TVD數(shù)值離散格式
3.3.5 梯度計算的小二乘法
3.3.6 格心(CC)格式的單元界面的速度修正
3.3.7 坐標系下非靜壓模型的控制方程
3.3.8 非靜壓模型數(shù)值求解的預(yù)估校正法
第4章 非靜壓模型數(shù)值求解方法校驗
4.1 計算網(wǎng)格生成
4.2 數(shù)值算法驗證
4.2.1 半隱式數(shù)值離散格式校驗
4.2.2 TVD數(shù)值離散格式驗證
4.2.3 行波模擬能力驗證
4.3 非靜壓模型有效性校驗
4.3.1 深水駐波驗證
4.3.2 孤立波驗證
4.3.3 非線性波越潛壩的數(shù)值模擬
4.3.4 地形陡變條件下帶自由表面流動的模擬
4.3.5 孤立波繞方柱的波面演化
4.3.6 彎曲河道環(huán)流的模擬驗證
第5章 非靜壓模型的不同形式及數(shù)值求解問題分析
5.1 垂向不同分層數(shù)的非靜壓模型
5.1.1 垂向的多層模型
5.1.2 水深積分的單層模型
5.1.3 水深積分的垂向兩層模型
5.2 動壓場的數(shù)值求解
5.2.1 垂向網(wǎng)格粗化的動壓求解方法
5.2.2 動壓場離散化控制方程組的數(shù)值求解
5.2.3 區(qū)域劃分的非靜壓模型數(shù)值求解
5.3 破碎波的數(shù)值模擬
5.3.1 簡述
5.3.2 滾波模型
5.3.3 人工黏性耗散模型(artificial eddy viscosity)
5.3.4 非靜壓模型的模式分裂法模擬破碎波
5.3.5 流動間斷捕捉法模擬破碎波
5.3.6 非靜壓模型模擬破碎波的算例驗證
第6章 非靜壓模型在高精度數(shù)值模擬中的應(yīng)用
6.1 帶自由表面水流繞結(jié)構(gòu)物的數(shù)值模擬
6.2 非靜壓模型結(jié)合分離渦湍流模式的高精度模擬
6.2.1 分離渦湍流模型
6.2.2 分離渦湍流模型的灰區(qū)問題
6.2.3 高精度數(shù)值格式的影響
6.3 非靜壓模型在高精度湍流運動模擬中的應(yīng)用
6.3.1 平板邊界層數(shù)值模擬
6.3.2 非平底明渠水流的高精度模擬
6.3.3 非靜壓模型用于結(jié)構(gòu)物局部繞流的精細化模擬
6.4 小結(jié)
參考文獻
索引
后記