陳發林臺灣大學：應用力學研究所羅安成Ruo, An-ChengAn-ChengRuo2007-11-292018-06-292007-11-292018-06-292005http://ntur.lib.ntu.edu.tw//handle/246246/62418重力流的研究在理論解析方法上，通常是發展一套數學模式(即所謂的淺水方程式shallow-water equation, SWE)，並求其解來模擬分析重力流的相關現象。然而這套方程式並沒有考慮任何來自外圍流體的阻力，使得其解與許多實驗觀察結果不符合，無法正確預測重力流之流動行為。為了使這套模型能被應用合宜，先前學者們便在重力流的前緣處套用一個波前條件來解SWE。雖然這個方法已被廣泛並有效地應用在許多重力流的研究上，但由於波前條件僅適用於二維流場，因此使得淺水方程式一直無法用來預測三維流場的行為，此為其一大缺點。另外一方面，雖然波前條件的使用已經導致理論分析上很大成功，但卻沒有人曾解釋過為何會成功的原因。從數學上看來，使用這個波前條件來解一套雙曲線型的方程式是不嚴謹的，所以「使用波前條件」被視為只是一個邏輯性的妥協，並非正規或嚴謹的方法。 本論文的主要目的，是發展一套由基本統御方程式出發的完整數學模型，並把影響重力流動力行為甚鉅的形狀阻力(form drag)包含在方程式中，使得本模型可以不用波前條件便足以適當的描述重力流之動力行為。因此本模型可被擴充和應用至三維流場的分析，這是重力流理論分析的一大進展。更特別的是，我們利用微擾法中一種很有用的方法－貼合漸近展開法，可將重力流分成靠近頭部的內場區域和遠離頭部的外場區域，分別解得各區域的近似解，最後再貼合起來成為全解。結果顯示，我們所獲得的外場方程式便等於SWE，而我們的內場分析竟然可以回復成波前條件。最終我們便可由此結果清楚地解釋為何波前條件解淺水方程的方法可以奏效的原因。 研究結果亦顯示，重力流前緣的阻力主要來自於將外圍流體排開之形狀阻力。但當密度比很小時，形狀阻力式微，此時黏滯阻力才變得重要。因此對密度比不很小的重力流而言，本理論模型可以適當地描述出重力流主要的行為。To analyze the motion of gravity currents, a common approach is to solve the hyperbolic shallow water equations (SWE) together with the boundary conditions at both the current source far upstream and the current front at the downstream margin. The use of the front condition aims to take the resistance from the ambient fluid into account, however, which is absent in the SWE. In the present study, we re-derive the SWE by taking the ambient resistance into account, and end up with the so-called modified shallow water equations (MSWE) in which the ambient resistance is accounted for by a nonlinear term, so that the use of the front condition becomes unnecessary. These highly nonlinear equations are approximated by perturbation expansions to the leading order, and the resultant singular perturbation equations are solved by an inner-outer expansion. For constant-flux and constant-volume gravity currents, the outer solutions turn out to be exactly the same as the previous ones obtained by solving SWE with the front condition. The inner solution gives both the profile and the velocity of the current head and also leads to the recovery of the front condition while is in a more general form. The combination of inner and outer solutions gives a composite solution for the whole current, which was called by Benjamin (1968) a “formidably complicated” task. To take the turbulent drag on the current into account, we introduce the semi-empirical Chezy drag term into the MSWE and end up with a result comparable with experimental data. The result implies that the ambient resistance is contributed mainly by the inviscid form drag, and the viscous drag dominates only when the density ratio is small. Furthermore, the MSWE can be extended for three-dimensional viscous currents, while will become more complicate that present analytical approach may not be feasible.英文摘要 i 中文摘要 ii 致謝 iii 目錄 iv 表目錄 vii 圖目錄 vii 符號對照表 x 第一章 緒 論 1 1.1 基本現象的描述 2 1.2 流場的特徵與研究上的困難 5 1.2.1 實驗室的觀察方法 6 1.2.2 研究上的困難點 7 1.3 文獻回顧 8 1.3.1 自由表面流的潰壩問題 10 1.3.2 海面下的潰壩問題 11 1.3.3 穩態分析方法 14 1.3.4 其它相關重要文獻 17 1.4 研究動機 20 1.4.1 邏輯性的妥協 21 1.4.2 本文的目標 21 1.5 本文架構 22 第二章 重力流理論分析之初探－穩態分析 23 2.1 穩態的重力流模型 23 2.2 二維勢流理論之應用 27 2.2.1 單層重力流 28 2.2.2 雙層重力流 31 2.2.3 考慮外圍流體之流動 33 第三章 淺水方程之修正式 34 3.1 二維淺水方程式之推導 35 3.2 二維淺水方程之修正 40 3.2.1 未知速度 的分解 42 3.2.2 包含速度勢 的微分項之分解 44 3.2.3 二維淺水方程之修正式 46 第四章 固定入流率的重力流分析 48 4.1 外場分析 49 4.1.1 外場方程式與異性微擾 49 4.1.2 相似變換法 51 4.2 內場分析 54 4.3 使用貼合的技巧獲得完整的解 57 4.3.1 外解間的貼合 58 4.3.2 內解與外解的貼合 59 4.3.3 完整的貼合解 62 4.4 解的合理性與細部討論 63 4.4.1 解的合理性和成立條件 63 4.4.2 解的極限與其物理涵義 64 第五章 固定體積的重力流分析 68 5.1 初始潰散期 (initial slumping phase) 69 5.1.1 外場分析 69 5.1.2 內場的獨立性 72 5.1.3 貼合出完整解 72 5.1.4 應用與討論 74 5.2自我相似期 (Self-similar phase) 75 5.2.1 外場方程式 75 5.2.2 相似解 76 5.2.3 貼合的工作 78 5.2.4 乾區解 (dry-zone solutions) 80 第六章 MSWE之衍生式 85 6.1 黏滯阻力項的修正－柯西阻力 85 6.1.1 外場分析 87 6.1.2 內場分析 87 6.2 三維的MSWE 90 第七章 結果討論與展望 95 7.1 本理論方法與結果之總回顧 95 7.2 非線性本質和相似解的真義 98 7.2.1 非線性本質 99 7.2.2 相似解的真義 100 7.3 理論與實驗誤差之來源剖析 101 7.4 其它特別應用與展望 102 7.5 結語 103 參考文獻 104 附錄 A 相似轉換群之推導 108 附錄B 第一階微擾解為零之佐證 113 附錄 C 軸對稱重力內解獨立性之證明 115 附錄 D 不同釋放條件下所有解之總整理 118 附錄 E 用微擾法解MSWE之流程圖 119 作者簡歷 1204499222 bytesapplication/pdfen-US重力流異重流淺水波方程式波前條件形狀阻力移動性邊界潰壩貼合展開法gravity currentshallow water equationsfront conditionform dragmoving boundaryDam breakmethod of matched asymptotic expansionthesishttp://ntur.lib.ntu.edu.tw/bitstream/246246/62418/1/ntu-94-D89543004-1.pdf