非結(jié)構(gòu)三角形網(wǎng)格下三維淺水流動(dòng)的數(shù)值模擬.pdf

1、在模擬具有間斷或解梯度較大物理過(guò)程時(shí)必須采用高精度的計(jì)算格式.建立了一種基于有限體積法的三維的淺水流動(dòng)及物質(zhì)輸移數(shù)學(xué)模型.在理論上,模型通過(guò)求解深度積分方程計(jì)算水位,作為外模式；通過(guò)求解三維淺水運(yùn)動(dòng)控制方程,計(jì)算流速垂向分布；采用κ-ε紊動(dòng)模型封閉模型.在計(jì)算方法上,利用隱格式求解三維流速場(chǎng),以增強(qiáng)數(shù)值穩(wěn)定性；平面上采用無(wú)結(jié)構(gòu)三角形網(wǎng)格,以適應(yīng)不規(guī)則地形；垂向上采用σ坐標(biāo)系,以簡(jiǎn)化表面和底部邊界條件；采用、Osher格式計(jì)算界面數(shù)值通量

2、,使水流模型具有很好的守恒性、逆風(fēng)性、單調(diào)性和高分辨率地捕捉間斷的能力,并且在邊界的處理上不需引入冗余的數(shù)值邊界條件.輸運(yùn)方程采用TVD格式進(jìn)行離散,并引入高階數(shù)值黏性項(xiàng),使模型既可高精度捕捉間斷,又有良好的穩(wěn)定性. 利用模型模擬了一個(gè)圓錐形濃度場(chǎng)在旋轉(zhuǎn)流場(chǎng)中的運(yùn)動(dòng).分別采用迎風(fēng)格式、中心格式、TVD格式計(jì)算對(duì)流流通量.結(jié)果顯示TWD格式在保持濃度的形狀和精度方面表現(xiàn)最好. 利用模型模擬了因突擴(kuò)引起的明渠回流區(qū)水流運(yùn)動(dòng),

3、重點(diǎn)模擬了由剪切層不穩(wěn)定引起的渦的產(chǎn)生、輸移及其被主渦吞并的過(guò)程以及對(duì)應(yīng)一個(gè)共振周期內(nèi)回流區(qū)流場(chǎng)、水位變化過(guò)程；將斷面流速分布的計(jì)算結(jié)果與相應(yīng)的實(shí)測(cè)結(jié)果進(jìn)行了比較,結(jié)果表明二者是一致的. 將模型模擬應(yīng)用于一個(gè)平直理想河道,模擬了河道中鹽水與淡水混合的過(guò)程.計(jì)算了鹽淡水混合過(guò)程中鹽度平面分布以及沿河道中心縱剖面分布；同時(shí)還給出了水位分布和縱剖面速度場(chǎng).結(jié)果表明模型計(jì)算結(jié)果正確的反映了鹽淡水混合過(guò)程中的鹽度分布、水位分布和速度場(chǎng)的情