河勢貼體河道平面二維正交網(wǎng)格生成方法的研究及應用_第1頁
河勢貼體河道平面二維正交網(wǎng)格生成方法的研究及應用_第2頁
河勢貼體河道平面二維正交網(wǎng)格生成方法的研究及應用_第3頁
河勢貼體河道平面二維正交網(wǎng)格生成方法的研究及應用_第4頁
河勢貼體河道平面二維正交網(wǎng)格生成方法的研究及應用_第5頁
已閱讀5頁,還剩2頁未讀 繼續(xù)免費閱讀

下載本文檔

版權說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權,請進行舉報或認領

文檔簡介

1、河勢貼體河道平面二維正交網(wǎng)格生成方法的研究及河勢貼體河道平面二維正交網(wǎng)格生成方法的研究及應用應用摘要:基于河勢概念和 Hermite 三次插值函數(shù),提出了河勢貼體河道平面二維正交四邊形網(wǎng)格的生成 方法 ;與邊界層坐標系下水深平均流體力學控制方程和 SIMPLER 算法聯(lián)合使用,建立了河道平面二維數(shù)學模型;進行了葛洲壩樞紐至磨盤溪河段二維網(wǎng)格生成及水流泥沙數(shù)學模型的實例 研究 ;討論了該網(wǎng)格生成方法的主要特征。關鍵詞:正交 網(wǎng)絡 數(shù)學模型 河勢 Hermite 函數(shù) 邊界層 SIMPLER 算法1河勢貼體網(wǎng)格河道平面二維數(shù)學模型網(wǎng)格生成方法研究中兩個關鍵 問題 是:(1)網(wǎng)格與河道擬合的貼體問題

2、;(2)二維網(wǎng)格、控制方程和數(shù)值方法三者之間的匹配問題。 目前 ,常用的河道二維正交網(wǎng)格生成方法是邊界擬合坐標系方法,即河道 Thompson 法1,它主要是通過物理平面(天然河道平面)與變換平面(數(shù)模 計算 平面)之間 Poisson 方程邊值問題數(shù)值解實現(xiàn)二維正交網(wǎng)格的生成。式中:x,y 為物理平面網(wǎng)格坐標;,為變換平面網(wǎng)格坐標;D 為物理平面區(qū)域。河道 Thompson 法存在的主要問題有:(1)二維網(wǎng)絡主要考慮與河道岸線的擬合,多數(shù)情況下與河道的河勢或主流線之間的擬合不理想;(2)復雜洲、灘及岸線河道邊界、非恒定流動岸和數(shù)值求解引起的動邊界等情況下,網(wǎng)格與河道岸線之間的擬合同樣會出現(xiàn)實

3、際偏離;(3)控制方程變換引起的數(shù)值求解困難和數(shù)值誘發(fā)耗散問題;(4)對于常見的寬、窄相間的河道平面形態(tài),二維網(wǎng)格不均勻間距可能導致的數(shù)值計算精度問題;(5)與河道一維斷面數(shù)學模型的嵌套和聯(lián)合存在接口困難等等。為了盡可能地避免河道 Thompson 法的上述問題和困難,本文提出了河勢貼體河道平面二維正交網(wǎng)格生成方法。采用這一網(wǎng)格生成方法及相應的河道二維數(shù)學模型,不僅可以進行天然河道大多數(shù)邊界條件下的水深平均平面二維水流、泥沙及溫排和污排等計算 分析 及 應用 研究,而且還為長河段河道二維及一、二維嵌套數(shù)學模型的研究和應用提供了新思路。2網(wǎng)格生成方法張瑞瑾教授2認為:河道水流的流態(tài)(或河勢)具有

4、很廣的含義,一切標志河道水流總體傾向的現(xiàn)象,都被納入這一概念之中。將河流動力學中基本和核心的河勢概念,引入到河道二維網(wǎng)格生成方法及其河道二維數(shù)學模型的研究之中,具有更加堅實的 理論 依據(jù)和物理基礎。河勢貼體網(wǎng)格生成方法的總體思想為:放棄傳統(tǒng)的網(wǎng)格生成方法中嚴格要求網(wǎng)格與河道岸線相擬合的思路,采用 Hermite 三次插值函數(shù)3,生成河道沿程縱向與河勢或主流線相擬合的河勢擬合線(曲線),并使得河寬方向橫向網(wǎng)格線(直線)與網(wǎng)格控制斷面相吻合,從而構造出平面二維正交四邊形網(wǎng)格。由此生成的二維網(wǎng)格,一方面避開了河道復雜和變動的洲、灘及河岸岸線,另一方面體現(xiàn)了控制河道水流運動的河勢概念。2.1Hermi

5、te 插值函數(shù)Hermite 插值函數(shù)不僅要求插值節(jié)點上的函數(shù)值相等,而且還要求節(jié)點處一階甚至高階導數(shù)相等。兩個插值節(jié)點情況下的 Hermite 三次插值函數(shù)可表述如下3:需要滿足的節(jié)點條件為函數(shù)表達式為系數(shù)表達式為式中:H3(x),H3(x)為 Hermite 三次插值函數(shù)及其導數(shù);,為插值系數(shù);x,y,m 分別為插值節(jié)點的坐標、函數(shù)和一階導數(shù)。2.2網(wǎng)格生成步驟河勢貼體河道平面二維正交四邊形網(wǎng)格生成方法包括如下三個步驟。2.2.1確定網(wǎng)格控制斷面和節(jié)點選取研究河段的進出口斷面、河段內(nèi)水位/水文站點或測流斷面、河勢控制斷面以及需要重點研究河段的控制斷面等作為網(wǎng)格控制斷面;在所選取的網(wǎng)格控制斷

6、面上確定網(wǎng)格控制節(jié)點,這些節(jié)點可以任意選取在控制斷面的左右岸、深泓點、主流點、中心點等處,所選擇的網(wǎng)格控制節(jié)點即為數(shù)學上 Hermite 三次插值函數(shù)的計算節(jié)點。2.2.2生成河勢擬合曲線利用上述 Hermite 三次插值函數(shù),可以生成一條既通過網(wǎng)格控制節(jié)點,又垂直于網(wǎng)格控制斷面的河道縱向網(wǎng)格控制曲線。通過多次調(diào)整網(wǎng)格控制斷面和節(jié)點,使得所生成的網(wǎng)格控制曲線與研究河段的河勢或主流線相擬合,將最終生成的河道縱向網(wǎng)格控制曲線確定為河勢擬合曲線。這一步是河勢貼體網(wǎng)格生成方法的關鍵和核心。2.2.3構造平面二維網(wǎng)格選取和調(diào)整縱向和橫向網(wǎng)格間距,構造由平行于河勢擬合曲線的曲線簇和垂直于河勢擬合曲線的直線

7、簇(包括網(wǎng)格控制斷面)的河勢貼體平面二維正交四邊形網(wǎng)格。3河道二維數(shù)學模型河勢貼體河道平面二維正交四邊形網(wǎng)格生成 方法 與邊界層坐標系下水深平均流體力學控制方程4以及合適的數(shù)值方法(如 SIMPLER 算法5)之間聯(lián)合使用,可以建立河道二維數(shù)學模型。本文給出了河道二維水流泥沙數(shù)學模型的控制方程和數(shù)值方法。3.1二維水流泥沙控制方程采用邊界層坐標系下簡化的河道二維淺水控制方程和泥沙對流擴散方程4,6式中:x,y,t 為邊界層坐標系下平面二維坐標和時間坐標;R 為 x 軸沿程曲率半徑;U,V 為流速分量;h,為水深和水位;g,n,vt為重力加速度常數(shù)、河道糙率和紊動粘性系數(shù);S 為含沙量;S 為水

8、流挾沙力;s 為泥沙擴散系數(shù);為泥沙恢復飽和系數(shù);為泥沙沉速。上述控制方程中二階導數(shù)項進行了適當?shù)暮喕秃喜?;?R時,上述控制方程可轉換為直角坐標系下水深平均平面二維水流泥沙控制方程;有關定解條件和參、系數(shù)函數(shù)或方程(如水流挾沙力公式等)與直角坐標系下二維水沙數(shù)?;疽恢?,8。3.2數(shù)值求解方法圖 1葛洲壩樞紐至磨盤溪河段河勢圖(1997 年 9 月河道地形)圖 2葛洲壩樞紐至磨盤溪河段河道二維水沙數(shù)模河勢貼體正交四邊形網(wǎng)格圖 3葛洲壩樞紐至磨盤溪河段二維水沙數(shù)模 計算 的流速等值線二維水流控制方程的數(shù)值離散和迭代求解基于 SIMPLER5算法,泥沙對流擴散方程的數(shù)值離散和迭代求解基于有限

9、控制容積法5,離散方程迭代求解方法具體包括以三對角追趕法(TDNA 法)為核心的逐行法和高斯塞德爾點迭代法,并配合塊修正和欠松馳修正技術等。計算過程中采用了動邊界模擬技術,具體處理措施包括:每次迭代根據(jù)二維網(wǎng)格節(jié)點的計算水深值,均要判斷和區(qū)分水域和陸域節(jié)點;對于陸域節(jié)點采用邊界隔墻法5處理,并讓陸域節(jié)點始終保持一個較小的富余水深等。4應用 實例河勢貼體網(wǎng)格生成的主要目的是為河道平面二維數(shù)學模型的 研究 提供二維計算網(wǎng)格及離散節(jié)點。作為應用實例,本文給出了葛洲壩樞紐至磨盤溪近壩河段(圖 1)二維網(wǎng)格生成的具體過程及二維水流泥沙數(shù)學模型研究的部分成果。該項研究的主要目的是 分析 論證葛洲壩下游壅o

10、nclick=g(水工);水工程措施(如胭脂壩左汊布設潛壩)的壅水效果及其對onclick=g(防洪);防洪和航運的 影響 ,而采用壅onclick=g(水工);水工程措施的主要目的是為了解決三峽工程運用初期枯水位時葛洲壩下游引航道通航水深不足問題 。依據(jù)上述河勢貼體網(wǎng)格生成方法及步驟,首先選取0(壩軸線),Y34(廟咀),Y37(宜昌),Y39,Y41(寶塔河),Y44,Y46,Y49 和 Y50(磨盤溪)共 9 個河段河勢控制或水文/水位斷面作為網(wǎng)格控制斷面,選取 9 個斷面的中點作為網(wǎng)格控制節(jié)點;然后,由 Hermite 三次插值函數(shù)生成了本河段的河勢擬合曲線 A(圖 1);最后,通過確

11、定縱向 x,橫向 y 的網(wǎng)格間距,構成本河段河道平面二維計算網(wǎng)格(圖 2)。最終生成的二維計算網(wǎng)格節(jié)點數(shù)為 13741,河勢方向(x 方向)網(wǎng)格間距為100200m,斷面方向(y 方向)網(wǎng)格間距為 50m。采用生成的河勢貼體河道平面二維正交網(wǎng)格,利用上述的河道二維水流泥沙數(shù)學模型,即可進行葛洲壩樞紐至磨盤溪近壩河段的二維水流泥沙數(shù)模計算分析研究。計算程序為自編“HELIU11”程序,有關研究 內(nèi)容 及成果請詳見長江 科學 院“九五”三峽工程泥沙問題研究子題分報告:“葛洲壩樞紐下游近壩段整治二維水流泥沙數(shù)學模型研究”。本文僅給出了宜昌流量 13500m3/s 時該河段二維數(shù)模計算的流速等值線(圖

12、 3)和 1993 年 11 月 1 日至 1997 年 8 月 31日沖淤驗證計算中的沖刷部位等值線(圖 4)。5討論本文提出的河勢貼體河道平面二維正交網(wǎng)格生成方法具有如下幾個主要特點:(1)引入河流動力學的河勢概念,凸現(xiàn)了河道二維數(shù)學模型研究中的主要矛盾,避開了河道岸線擬合這一復雜但相對次要的矛盾;(2)生成的二維正交網(wǎng)格與邊界層坐標系下河道水深平均流體力學控制方程及相應合適的數(shù)值方法(如 SIMPLER 算法)三者之間匹配較好;(3)網(wǎng)格保留了河道一維斷面形式,為長河段河道二維以及河道一、二維嵌套數(shù)學模型的嵌套和聯(lián)合提供了便利的二維網(wǎng)格和一、二維接口條件;(4)生成的二維網(wǎng)格具有河寬方向

13、上分布均勻,河勢方向上可調(diào)整性大,操作上簡便易行等優(yōu)點;(5)該網(wǎng)格生成方法不適用于強彎及鵝頭型分汊等平面形態(tài)的少數(shù)河段。參考 文獻 :1吳江航,韓慶書.計算流體力學的 理論 、方法及應用M.北京:科學出版社,1988.2張瑞瑾.河道水流運動的基本特性A.張瑞瑾onclick=g(論文);論文集C.北京: 中國onclick=g(水利);水利水電出版社,1996.3李慶揚,王能超,易大義.數(shù)值分析M.武漢:華中工學院出版社,1982.4吳望一.流體力學M.北京:北京大學出版社,1982.5美S.V.帕坦卡.傳熱與流體流動的數(shù)值計算M.北京:科學出版社,1980.6DONGYAOHUA. The generat in gmethod of river regime-fitted orriver main stream-fitted or thogonal grid and itsapplicationsA.Lee J.H.,Jayawardena A. W.Wang . Y.Proc

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
  • 4. 未經(jīng)權益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
  • 6. 下載文件中如有侵權或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論