基于cfd的芙蓉湖流场的二维数值模拟-环境工程专业论文.docxVIP

西安建筑科技大学硕士学位论文 西安建筑科技大学硕士学位论文 西安建筑科技大学硕士学位论文 西安建筑科技大学硕士学位论文 湖泊资源，保护湖泊水环境，维持湖泊生态平衡。开展湖泊水环境系统综合研究 的重要方法之一就是建立确定性数学模型 ，精细而全面地模拟水体流动。该方法 具有投资少、不受时间限制、 可完全自由地改变和控制流体性质等优点 。70 年代 以来，随着计算机技术的发展和普及，数值模拟方法开始应用到湖泊水环境的研 究领域中，主要任务是根据流域湖泊水体的数量与质量的历史及现时数据、资料 信息等，按其运动 、变化与发展规律，很好地展示目前的状况并预测其未来的状 况，从而为湖泊水环境的合理利用和管理提供辅助决策的依据。 1.2 湖泊水流数值模拟研究概况 研究流体流动的主要方法有三种: 1.理论分析方法，建立力 学模型(通常是 对实际流动问题的合理简化)，运用物理学的基本定律结合流动的特点推导相应的 数学方程，用数学分析方法结合初始与边界条件求解方程，检验和解释求解结果。 正确的理论分析结果可揭示流体运动的本质特性和规律，因此具有普遍的适用性。 这些研究成果推动了流体力学的发展，奠定了今天数值模拟的基础，很多 方法仍 是目前解决实际问题是常采用的方法 181 2.实验研究方法，在相似理论的指导下， 在实验室内建立模型实验装置，用流体测量技术测量模型实验中的流动参数，处 理和分析实验数据并将它归纳为经验公式。其缺点是从实验室中得到的经验公式 的普适性较差。 3.数值研究方法，随着计算机技术突飞猛进的发展，过去无法解析 的流体力学微分方程，现在可以用计算机数值方法求解。数值研究的一般过程是: 对流体力学数学方程作简化和数值离散化，编制程序作数值计算，将计算结果与 实验或理论解析结果比较。其优点是能计算解析方法无法求解的流动问题，能模 拟多种工况的流动问题，比实验方法省时省钱。随着数值模拟技术的发展和计算 机性能的不断提高，以及己有理论和试验成果的积累，将使得数值模型得到更多 的发展和应用。 湖水运动是湖泊的重要物理特征之一，它直接关系到湖泊中物质、能量的输 送与转换，以及湖泊与外界的物质和能量的交换和转移，在很大程度上影响、决 定着水质的宏观变化过程。因此，湖水运动是湖泊中 一切物理因子变化的基本动 力条件，也是 引起湖泊环境变化的一个极其重要的驱动 力闷。迄今为止，国内外在 湖泊水流方面做了大量的工作，运用数值模拟的方法取得了很多具有重要意义的 结果。 1.2.1 湖泊数值模拟的数值计算方法 湖泊水流模型数值计算方法目前 的离散方法主要有 :有限差分法 (FDM) 、有 限单元法 (FEM) 、有限体和、法 (FVM)、边界单元法(BEM) 等。 ( 1)有限差分法 (FDM) 有限差分法是最早采用，至今仍被广泛运用的 一种数值方法。基本思想是把 计算区域划分为矩形或正交曲线网格，在网格线交点(即节点)上，将控制方程中的 每一个微商用差商 来代替，从而将连续函数的微分方程离散为网格节点上定义的 差分方程，每个方程 中包含了本节点及其附近一 些节点上的待求函数值，通过求 解这些代数方程就可获得所 需的数值解。常用的有迎风格式、 ADI 法、蛙跳格式、 Rusanov 格式及特征格式等。 有限差分法的优点是它建立在经典的数学逼近理论的基础上，容易被人们理 解和接受，求解直观迅速，数学概念清晰 ，表达简单，编程及资料的整理比较简 单，占用计算机内存少，计算的收敛性和稳定性理论也比较成熟:有限差分法的 主要缺点是对于复杂流体区域的边界形状处理不方便 ，网格布置不灵活，难以概 化和贴近不规则的边界，处理得不好将影响计算精度。为了克服其局限性，许多 学者都对如何更好地处理不规则边界问题进行了研究，如美国的Th ompson .J卫9J 等人提出的边界拟合坐标系方法，谢利曼 [1创所提出的任意网格差分法，华祖林 [11] 等也对边界拟合中的调节因子进行了有益的探索。 (2) 有限单元法 (FEM) 有限元法是一种高效能、常用的计算方法。基本思想是把适定的微分问题的 计算区域进行离散化，将其剖分成有限个相连结又互不重叠且具有 一定规则几何 形状的有限个子区域(如 :二维问题可以划分为三角形或囚边形:三维问题可以划 分为四面体或六面体等)，这些子区域称之为单 元或有限元，单元之间以节 点联结。 函数值被定义在节点上，对所有 单元进行独立的处理和运算，以便对 一个整体问 题进行局部化处理，在单元中选择基函数(又称插值函数)，使每个单元的计算都变 得非常简单，以节点函数值与函数乘积的线性组合构成单元的近似解，逼近单元 中的真解。整个计算区域上总体的基函数可以看为由每个单 元基函数组成的，则 整个计算区域内的解可以看作是由所有单元上的