黎曼近似解模型在太湖藻类浓度场求解中应用-3卷第1期.PDF

JournalofHohaiUniversity(NaturalSciences) Jan．2005 黎曼近似解模型在太湖藻类浓度场求解中的应用 洪晓瑜，逢 勇 (河海大学环境科学与工程学院，江苏 南京 210098) g 摘要：利用太湖常规监测点的监测资料，在黎曼近似解模型的基础上建立了太湖二维藻类生长模 r 型 ．该模型采用有限体积法及黎曼近似解对二维水流水质方程组逐时段、逐单元进行求解，将每个 时段各个单元 p(TN)Z~p(TP)的模拟值传递到藻类的生长项，从而模拟 出藻类的生长过程．率定验 o 证结果表明，该模型能较好地模拟藻类的生长过程 ，并能对太湖水体富营养化进行预测． . 关键词 ：太湖；黎曼近似解模型；藻类生长；浓度场 r 中图分类号 ：X143 文献标识码 ：A 文章编号：1000—1980(2005)01—0041—04 te 随着沿湖地区工农业的迅猛发展，水体富营养化以及由此引发的相关问题成为太湖 目前面临的主要环 a 境问题 ．水体富营养化的重要指标是叶绿素a的浓度 ，而叶绿素 a的产生与 TN和 11P的浓度有密切关系．因 此，本文应用二维藻类生长模型对太湖水质指标TN和 11P以及富营养化指标藻类叶绿素a进行模拟 ，并对该 w 网 模型进行率定验证 ，以期对水体富营养化进行预测． 1 控制微分方程及其积分离散ity水 1．1 控制方程 c 镇 二维浅水方程和对流一扩散方程的守恒形式为 ’ a nat+十 a +十 a1， 一：l， )， (1，) 城 q= (h，hu，hv，hp)T ，(口)= ( ， + ／2，huv，hup)T g(q)= (hv，huv，hv+ ／2，hvp)r i 6(口)= (o，gh(S0一 厅)，gh(S0y一 )， -(Di (^P))+S／A—r,hp)T h 国 式中：口—— 守恒物理量；，(口)—— 向通量 ；g(口)—— ，，向通量 ；b(口)—— 源汇项；^—— 水深 ； ，l，—— c 向和Y向垂线平均水平流速分量 ；JD—— 污染物质量浓度． 1．2 方程离散 . 中 定义矩阵F(口)=[，(q) g(q)]，在任意形状的单元 上对式(1)进行积分．利用散度定理 ，可得有 限 w 体积法的基本公式…1 w w lIqtd8：一如IF(q)·ndl+ (q)d6 (2) 式中：Jl—— 单元边界a 的单位外法向量；d叫，d卜一 面元积分和线元积分 ；F(q)·Jl——法向数值通量． 离散方程 (2)并利用 ，(口)和 g(q)的旋转不变性n]，得有限体积法的基本方程为 A ：一∑ ()一，()+A6(口) (3) 』#1 式中：，()叫 (口)的增广矩阵；——单元第 ．『边的长度 ；A—