伴随数据同化法反演涌潮河口开边界-水科学进展.PDFVIP

第 18 卷 第 6 期 水 科 学 进 展 Vol18 ,No6 2007 年 11 月 　　 ADVANCES IN WATER SCIENCE 　 Nov. , 2007 　 伴随数据同化法反演涌潮河口开边界 1 ,2 1 1 ,2 1 ,2 陈耀登 , 王如云 , 李 　未 , 卢长娜 ( 1 河海大学水文水资源与水利工程科学国家重点实验室 , 江苏 南京 　210098 ; 2 河海大学海洋学院 , 江苏 南京 　210098) 摘要 : 构造了非线性守恒型浅水方程的伴随方程 , 并给出无结构三角网格的开边界校正表达式 , 利用 Godunov 格式 和 Riemann 间断解 Roe 通量格式的设计思想 , 建立了无结构三角网格有限体积法的伴随开边界反演模型 。进行了概 化涌潮河口和钱塘江涌潮河 口M2 分潮开边界反演的数值实验 , 经过同化得到开边界上振幅的平均误差分别为 0000 6 m 和 0053 3 m , 证实了本文构造有限体积伴随模型的可行性 , 也表明本模型能够适应间断解问题 。 关 　键 　词 : 开边界 ; 伴随数据同化 ; 涌潮河口; 有限体积法 ; 无结构网格 ( ) 中图分类号 : P7312 ; P73123 　　　文献标识码 : A 　　　文章编号 : 2007 对于有限区域水动力数值模拟 , 准确给定开边界条件是一件较困难而又重要的工作 。不少学者基于有限差 ( ) ( ) [ 1～4 ] 分法 FDM 和有限单元法 FEM , 进行了伴随数据同化方法反演开边界条件的研究 。FDM 虽处理灵活计算 效率高 , 却不能较好拟合边界 , 而 FEM 虽可较好拟合边界 , 却带来了相当大的计算量 。实际工作中 , 往往会 碰到区域形状复杂和物理变量在较小区域内变化较大 , 甚至出现涌潮等间断现象 。针对复杂的边界形状 , 人们 [5 ] [6 ] 提出了曲线网格 , 无结构网格 等多种解决方法 ; 针对物理变量在区域内存在变化较大的问题 , 利用了许多 高分辨格式[7 ] , 如 TVD 格式[8 ] , Roe 格式[9 ] , ENO 格式[ 10 ] 等 。至今尚未见到能够处理复杂计算域中带有大梯 ( ) 度水流现象的伴随模型 。鉴于有限体积法 FVM 集中了 FDM 计算简便和 FEM 边界拟合精度高等优点 , 本文利 用有限体积法的设计思想 , 建立了适应间断解问题的无结构网格有限体积法的伴随数据同化模型 ,