具有复杂计算域和地形的潮汐流动数值模拟.pdfVIP

第20卷第1期 水科学进展 V01．20．No．1 2009年1月 ADvANCESINWATER Jan．．2009 sClENCE 具有复杂计算域和地形的潮汐流动数值模拟 张宏伟，康海贵 (大连理工大学海岸和近海国家重点实验室，辽宁大连116023) 摘要：采用无结构网格上的Roe型二阶精度迎风型FVM格式的有限体积方法对具有复杂计算域和地形的胜利油田 海域的潮汐流动进行数值模拟。为保证Roe格式能够应用于复杂地形条件下的计算．采用Ropers格式来解决通量梯 度项与源项的平衡问题，对摩擦力源项采用分步法求解以增加格式的稳定性。应用此方法对地形复杂的胜利油田海 域进行了实际模拟，数值计算结果和实测结果吻合较好。 关键词：浅水方程；Roe迎风型离散；有限体积法；复杂地形；潮汐流动；数值模拟 中图分类号：P731．23 文献标识码：A 文章编号：1001．6791(a】0910l-oll8．06 大多数学者已将空气动力学中成熟的FVM格式直接移植到不考虑底坡和摩阻作用的齐次浅水波方程，且 利用一些经典的间断流算例对格式进行验证u“J，但并未进一步研究格式在多变地形下的实际应用。随着 中【5娟J。本文尝试采用无结构网格上Roe型二阶精度迎风型FVM格式对具有复杂计算域和水深地形的胜利油 田海域进行潮汐数值模拟，该格式相对其他FfvM离散格式而言，具有相容性好，物理意义更为清晰，编程易 于实现等优点。为保证该格式的和谐性，本文采用将通量梯度项与源项的平衡的方法，对于该方法一些学者也 做了大量的研究工作旧一2|，本文采用Ropers格式来解决因通量梯度项与源项的不平衡而引起的静止的水体时 产生的虚假流动问题，使用水位插值代替水深插值，从而使得Roe格式能够应用于复杂地形条件下的计算。对 摩擦力源项采用半隐式¨3j求解以增加格式的稳定性。 1 基本方程 主￡中u=[三孑]为号f铂亘向量量；，=[：uDi孑+gc r- o ] H=I跏一g善(S。+瓢)I为源向量。Sox，S町，％，勘分别为海床在算，Y方向的底坡源项和摩擦源项，其值 收稿日期：2008-01．04 基金项目：国家重点基础研究发展计划(973)资助项目(2002CB412410) 作者简介：张宏伟(1980一)，男，辽宁凌源人，博士研究生，主要从事海洋水文及其工程应用。 E-mail：hwr2maghem@163．{30111 第1期 张宏伟．等：具有复杂计算域和地形的潮汐流动数值模拟 1 19 速度；声是地球纬度。 2方程离散 采用CC(CellCentered scheme)方式的有限体积方法，把变量存在单元的中心，单元的边界为控制体。将第i 元取銎竺要塑：。．翌苎兰毫登墨毒擘控掣方程0． 图J非结构网格控制元及外法向示意图 非结构网格控制元及其外法向示意图，如图1所示。 Sketch Fig．1 0f哪叽二‘m=。j：：二 删 私=南善E庐‘+玩 (3) 式中 E，，为通过第i单元第J条边上的法向数值通量。△Z，元为第．『条边的边长。 由式(3)可见，问题归结为确定法向通量R。，可通过解局部一维黎曼问题求得。 2．1 Roe型数值通量的构造 U+A(U)以=0 (4) 式中 跳跃量的关系得