基于有限体积法的二维水流.ppt

基于有限体积法的二维水流水质模拟及其可视化研究 疯狂猫 目录 绪论 浅水动力学的控制方程 有限体积法原理 浅水模型的程序实现 数据可视化 验证算例 绪论 1，本课题的主要工作 使用有限体积法（Finite Volume Methods）在无结构网格上离散二维浅水方程，求得水流速度场，并把求得的流场作为输入条件进一步求解对流扩散方程得出污染物浓度场。 同时，计算过程将会产生大量的数据，本课题使用Delphi语言＋DirectX技术实现流场和浓度场的可视化。 绪论 2，课题的现实意义 本课题的浅水模型可以用来实现河流、河口、湖泊、近海的水流水位模拟；而对流扩散方程是水质模型的基本控制方程，可以实现污染物（BOD，COD）、泥沙、热、盐的数值模拟 ，具有重大的现实意义。 应用举例：河道流量及洪水预报、溃坝决堤、洪水漫滩、河口潮汐、盐水入侵、都市排水、海上油膜扩散等。 第一章 浅水动力学的控制方程 1，何谓浅水 这里所说的浅水（Shallow Water）并非是指水深比较浅的水，而是具有它自身特定的含义。我们把满足以下四个方面的水称为浅水： 有自由表面； 以重力为主要驱动力，同时考虑风应力 和地转柯氏力； 水平流速沿垂线近似均匀分布； 水平运动尺度远大于垂直运动尺度； 满足浅水假设的包括：河流、湖泊、河口和近海。 第一章 浅水动力学的控制方程 2，一维浅水控制方程 浅水方程又称圣维南方程（Saint－Venant），是由圣维南1871年在法国科学院汇刊73卷中提出的，它的一维形式为： 第一章 浅水动力学的控制方程 下面我们来推导以一维浅水方程，考虑以下的一维单元模型： 3，从一维到两维 我们上述的浅水方程是在只考虑横向尺度（X方向）的情况下得到的，对上述方程进行纵向尺度（Y方向）扩展，可以很容易地得到二维的浅水方程： 第一章 浅水动力学的控制方程 或写成向量形式： 式中： 第一章 浅水动力学的控制方程 4，求解浅水方程常用的数值方法 第一章 浅水动力学的控制方程 有限体积法和有限元一样将计算域分成若干单元。在计算出每个控制体边界沿法向输入（出）的流量和动量通量之后，对每个控制体分别进行水量和动量平衡计算，得到计算时段末各控制体的平均水深和流速。因此，FVM正是对于推导原始微分方程所用控制体途径的回归。 本课题使用FVM进行浅水方程的求解。 上述三种方法在具体应用中并没有绝对的好，和绝对的坏，根据具体问题的不同，可以用不同的数值方法。比如，DHI的二维浅水动力学软件Mike 21采用的就是差分中的ADI－Quick格式，而著名的SMS（Surface Water Modeling System）使用的是有限元法，英国Wallingford的InfoWorks中的RS（River Simulation）模块使用的是最为传统的四点差分格式。 第二章 有限体积法原理 1，有限体积法离散原理 第二章 有限体积法原理 2,从二维到一维 第二章 有限体积法原理 4，FVS格式求解黎曼问题 第二章 有限体积法原理 5，边界条件 第三章 浅水模型的程序实现 1，软件构架设计 第三章 浅水模型的程序实现 2，预处理－网格的生成 第三章 浅水模型的程序实现 3，模型数据结构设计 第三章 浅水模型的程序实现 4，模型求解过程的封装－类的使用 第三章 浅水模型的程序实现 5，模型求解的流程 第四章 数据可视化 1，数据可视化的需求分析 第四章 数据可视化 2，流场的绘制 第四章 数据可视化 3，等值云图的绘制 第五章 验证算例 1,溃坝算例-水流的验证 第五章 验证算例 2,均匀流场中污染物的迁移-水质的验证 第五章 验证算例 3,某河道的水流水质模拟 结尾 报告结束,多谢各位光临!? {TNode} //节点的数据结构 TNode=Record X:Real; //x坐标 Y:Real; //y坐标 q:Array[0..3] of Real; NeibElement:TIntArray;//相邻单元 End; 第三章 浅水模型的程序实现