第二章溃坝水流数值模拟.ppt

* 第二章 溃坝水流数值模拟 M. Quecedo， and M. Pastor. A reappraisal of Taylor–Galerkin algorithm for drying–wetting areas in shallow water computations，Int. J. Numer. Meth. Fluids 2002; 38:515–531 Application of the Level Set Method to the Finite Element Solution of Two-Phase Flows Comparison of two mathematical models for solving the dam break problem using the FEM method The Finite Element Method for Fluid Dynamics, Sixth Edition 一 传统经验算法 1 溃坝坝址峰顶流量计算 瞬间全溃 瞬间部分溃 逐渐溃 2 坝址流量过程线的计算 平底无阻力瞬间全溃 斜底无阻力瞬间全溃 斜底无阻力瞬间部分溃 平底有阻力瞬间全溃 3 洪水演进的近似算法 洪峰展平法 线性河道法 概化过程线法 统一公式： 二 流体力学基本方程 几何方程 平衡方程 本构方程 （不可压缩流体力学控制方程组） 1 N-S方程 2 浅水方程 非守恒连续方程： 非守恒动量方程： 张量形式表示： 浅水方程组： 其中 可以化为： 3 守恒式浅水方程 三 浅水方程数值离散 将 对时间在 附近进行二阶Taylor展开： 四 有限元实现 五 程序编制 六 数值算例 七 自由面问题的LEVEL SET方法 1 Level Set方法基本思路 Osher等人提出用来求解两相流问题的一种方法引入标记函数，是一个距离函数，定义其为某空间点到交界面的符号距离，在交界面以上为正，以下为负，而零值处便代表了交界面。随水流输运，从而可以达到捕捉界面的目的 2 控制方程 3 实施步骤 4 离散方法 （1）CBS格式 对多维情形有： （2）指标符号表示的CBS格式离散过程 应用加权余量法： A：右边第一项： B：右边第二项： C：右边第三项： D：右边第四项： E：右边第五项： 若记： （3）张量表示的CBS格式离散过程 采用伽辽金加权余量法时，取权函数为形函数时即为 A：右边第一项： B：右边第二项： C：右边第三项： D：右边第四项： E：右边第五项： （4）边界条件 对于纯对流问题，微分方程是一阶的，只需在进口施加边界条件；进口、出口的定义是根据流向来的确定的 （5）时间步长 5 N-S方程的离散 Step 1: Step 2: Step 3: 6 LEVEL SET方程数值离散 沿特征线展开得到 7 LEVEL SET函数重新初始化的数值离散 沿特征线展开得到 8 LEVEL SET算例 *