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