伴随数据同化法反演涌潮河口开边界-大气科学学院.PDFVIP

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