有限体积法和激波捕捉格式-黎曼问题解析解.pptx

浅水流动的特性与数值模拟 朱德军 旧水利馆 206,zhudejun@tsinghua.edu.cn 2 黎曼问题 广义黎曼不变量 稀疏波 间断波 接触间断 内容回顾 3 第五章 激波捕捉格式求解浅水方程组 有限体积法和Godunov格式 对流方程几种常见的数值格式 双曲方程组守恒格式 黎曼问题 黎曼问题解析解 黎曼问题近似解 Godunov格式的二阶修正 4 其中: 5.5 黎曼问题解析解 5 时， ，λk空间为线性退化空间 时， ，λk空间为完全非线性空间 6 对浅水方程，第一、三簇波既可能是稀疏波，也可能是间断波；第二簇波是接触间断。 所以 稀疏波特征线dx/dt=λi两侧，变量满足如下关系: 间断波两侧，变量满足如下R-H关系: 接触间断同时满足稀疏波和间断波的条件 7 穿过第一簇特征线广义黎曼不变量： 扇形区内部： 8 穿过第三簇特征线广义黎曼不变量： 扇形区内部： 9 浅水方程间断两侧变量的关系： 10 11 从(1)和(3)两式可以得到： 12 13 从上两式中消去未知量s，得到： 只有 满足 的条件 14 整理得到： 15 设 当α趋于0时， 16 相空间，雨奥尼高(Hugoniot)轨迹和稀疏波积分曲线 17 穿过第二簇特征线广义黎曼不变量： 18 穿过第二簇特征线，利用R-H条件，有如下关系： 变形后得到： 19 对浅水方程，第一、三簇波既可能是稀疏波，也可能是间断波；第二簇波是接触间断。 20 如果第一簇波为稀疏波： 如果第一簇波为间断波： 21 如果第三簇波为稀疏波： 如果第三簇波为间断波： 22 解的型式 可以暂时不考虑v 23 根据第一簇特征线两侧关系： 根据第三簇特征线两侧关系： 24 h*是下面关于h的方程的解 (2)-(1)=0， 25 26 27 为了保证h0，必须f(0)0，即： 28 29 30 例1：渠道长20 km，坝位于10 km处，瞬时全溃 uL=uR=0, hL=10 m, hR=5 m，计算溃坝100 s后的水面线 31 所以不会出现干涸，以下进行波形判断 32 ,左侧第一簇波为稀疏波，右侧为间断波 33 根据第一簇特征线两侧关系： 根据第三簇特征线两侧关系： (2)-(1)=0， 34 扇形区内部： 得到： 35 根据R-H条件： 得到： 36 最终得到水面线如下： 37 如果h*=0，得到： 间断波不可能和河床直接相连 不失一般性，设 h0，根据间断波R-H关系，得到： 与 h0矛盾，证明了间断波不能直接和河床相连 有干湿变化情况的解析解 38 间断波不可能和河床直接相连 39 40 两种出现干湿变化的情况 41 如果uL=0,则为干床上的溃坝问题 42 43 例2：渠道长20 km，坝位于10 km处，瞬时全溃 uL=-25 m/s，uR=25 m/s，hL=hR=10 m, 计算溃坝100 s后的水面线 44 所以出现干涸，第一簇波和第二簇波均为稀疏波 45 左侧扇形区内部： 得到： 46 右侧扇形区内部： 得到： 47 48 相空间 雨奥尼高(Hugoniot)轨迹 稀疏波积分曲线 浅水方程黎曼问题、溃坝问题 解的结构 稀疏波-稀疏波、激波-激波、稀疏波-激波、激波-稀疏波 有干湿变化的情况 无干湿变化的情况 小结 49 溃坝问题Godunov型格式模拟 假定一矩形河道长4km，宽10m，一堤坝位于2km处，坝上游水深10m，下游水深5m。从t=0时刻起，堤坝突然溃决，挡水能力完全消失，形成向上下游传播的溃坝波。河床水平，曼宁系数n=0.0001。 空间步长ds=40m，上游边界条件为水位保持10m不变，下游边界条件为水位保持5m不变。初始条件为：t=0s时，从断面1到断面51水位z=10m，流量Q=0，从断面52到断面101水位z=5m，流量Q=0。要求采用Godunov型格式求解，界面通量采用解析法计算，试选择计算时间步长,并 计算t=10s、30s、60s、120s时的水面线 计算坝址处的流量和水位过程线 作业 50 重构步：构造黎曼问题 演进步：求解黎曼问题 平衡步：根据黎曼问题的解，获得新时刻的节点变量 Godunov型格式计算步骤 思考题 如果在演进步采用黎曼问题的解析解，是否说明，采用Godunov格式求得的就是原问题的解析解？为什么？ 51