水库二维分层水流—水质模型的研究和应用.pdf

水库二维分层水流一水质模型的研究和应州 所需空间比较大，很难在实际中推广。目日F人们处理水库污染沿水体垂向混和比 较均匀的问题时，大都采用垂向平均二维方程””和垂向二维方程。7。来描述水体运 动和污染物的浓度分布。 本论文研究了水库的一种分层技术懈’，能够有效的克服水库分层带来的困 难，再导出水流及污染物浓度输移的控制方程，采用交错格点法”…离散、求解水 流的控制方程，建立了水库分层水流和污染物浓度输移的数学模型，并讨论了风 对水流和水质的影响。本论文的研究成果将会为水环境评价、水环境规划与管理 以及水资源保护等提供决策性的依据，对水环境的防治工程有一定的指导意义。 1．3水环境数学模型的计算方法 水流一水质数学模型是将水动力问题和污染物在水体中的迁移转化问题用 数学方程进行描述，并在一定的定解条件下求解这些方程，从而模拟出某个理论 和工程实际问题的水流、水质状况。 数值模拟方法可称为数学模型试验，所有的试验条件都以数字的形式给出， 对水体扰动，不存在缩尺效应，不受试验场地和观测仪器的影响。数值模拟具有 高效、经济、简便的优点。数值模拟也有局限性，它无法模拟数学方程尚不能描 述的物理现象，数值模拟的可靠性和精确度取决于数学方程的离散方法。 l-3．i有限差分法(Finite—DifferenceMethod，FDM) 计算机模拟水流最初用的就是有限差分法“…，至今其应用仍最为广泛。有限 差分法是最早提出的方程离散方法。1928年Courant等人提出了有限差分理论， 但应用于流体力学是在二十世纪六十年代。有限差分法是将微分方程中的各微分 项离散在微小网格上各相邻格点的差商形式，得到一个以各节点上函数数值未知 量的代数方程。并根据原问题的初值和边界值合理地给出离散化代数方程的初 始、边界条件，从而求出控制方程的数值解。主要包括：显、隐式差分、特征线 年，Leendertre出色的将有限差分方法及ADI格式成功地应用于河口与海岸平 泛的研究，取得了很大的进展。程爱杰于1998年给出二维变系数情形的收敛性和 2 水库二维分层水流一农质棱型的研究和应用 稳定性分析，得出了完善的结果““。特征线法是二十世纪五十年代林秉南提出的 一维水流计算的特征线法。破算子法是二十世纪五}‘年代未六十年代初苏联学者 Yananko提出的，并在理论上逐渐加以完善。这些算法的共同点都是基于矩形坐 标网格的有限差分法。FDM建立在经典的数学逼进理论的基础上，简单且为人 们接受，处理效率较高。缺点：控制微分方程表达了质量守恒与动量守恒的物理 定律，而差分方程有时不能保持守恒性质，数值解出现水量、动量不平衡的守恒 误差。 1．3．2有限元泫(Finite—Element Method，FEM) 有限元法…3是1950年有飞机结构师首先提出的一种数值方法，冯康六十年 代初独立创立了有限元法。Ziekiewicz等人1965用有限元法求解势流问题，』F 创了有限元法应用于流体力学，到七十年代有限元法开始应用于平面二维水流模 拟。其原理是分单元对解逼进，使微分方程空间积分的加权残差极小化。其基本 思想可以归纳为：剖分、差值、求权三部分，由此建立FEM方程组给出数值解。 FEM由于计算问题时，对非恒定流每一时间步要求解～个大型线性方程组。FEM 在数学上适于求解椭圆型方程的边值问题，不适于求解对流为主的输运问题。 FEM的最大优点是采用不规则无结构网格，对不规刚的边界拟舍较精确，计算 程序的通用性强；但缺点是程序编制复杂，占有的计算机内存大，计算时间相对 较长，缺乏足够耗散，不适于计算间断。国外水环境工作者研究比较多，如美国 应用于海湾、河口等宽广水域的计算“””…。 1．3．3有限体积法(Finite—Volume Method，FVM) 于气流计算。其计算原理是采用有限元法网格，沿每个控制体各个界面法向，采 用商性能守恒逆风格式计算跨越界面的质量和动量通量，然后对每个控制体将计 算域划分成若干规则或不规则形状的单元或控制体。在计算每个控制体边界沿法 向输入的流量和动量通量后，对每个控制体分别进行水量和动量平衡计算，便得