黄河三角洲海域的潮汐流动数值模拟 Numerical modeling of the Yellow River delta sea field tidal flow.pdfVIP

第32卷第1期 水道 港 口 Vbl．32No．1 Journalof andHarbor Feb．20ll 2011年2月 Waterway 黄河三角洲海域的潮汐流动数值模拟 张宏伟t，袁仲杰2 (1．大连花园口经济区规划建设局，大连116023；2．国家海洋环境监测中心，大连116023) 摘要：文章将一种有效处理浅水流动计算中的Roe型二阶精度迎风型FVM格式应用到平面二维浅水 流动数学模型中。为保证Roe格式的稳定性，应用Ropers格式以解决通量梯度项与源项的平衡问题。为 避免源项引起的计算不稳定，对源项采用分步法及半隐式求解。在此基础上，应用该模型对地形复杂的黄 河三角洲海域潮汐进行了实际数值模拟，并对其水动力特性进行了详细解析。结果表明该计算模型能够 很好地应用于具有复杂岸线条件以及计算域的近岸潮汐数值模拟中。 关键词：Roe迎风离散格式；有限体积法；黄河三角洲 242．1 中国分类号：P731．23：O 文献标识码：A 文章编号：1005—8443(2011)仇—0012—06 本文把三角结构网格上Roe型二阶精度迎风型FVM格式应用于对具有复杂水深地形和计算域的黄河 三角洲海域的潮汐数值模拟中。该格式相对其他FVM离散格式而言，具有相容性较好、物理意义清晰、计算 精度较高、编程容易等优点。采用水位插值代替水深插值的方法来解决通量梯度项与源项的平衡，从而避免 了因通量梯度项与源项的不平衡而引发的静水体时产生的虚假流动问题，使Roe格式能够适用于复杂地形 条件下的计算。为增加计算格式的稳定性，对摩擦力源项采用分步半隐式[1]求解。 1基本方程 为保证计算格式的和谐性，使得Roe格式能够应用于复杂地形条件下的计算，本文采用Ropers格式处 理底坡源项，以解决通量梯度项与源项的不平衡问题。 Ropers格式的二维浅水方程的守恒形式可表达为 以√(U)一g(U)部 (1) 『f 1 『曲 1『t，D 1 f0 1 式中：皓luDI为守恒向量；产l(u2D+g(铲+2hg)／2)1；g：luvDI；日_I劬一积潲)1为源向量；s。．，， 【vD J j 【uvD 【(勘饥g(弭2^f沈)JI一彩i嘲炙s—勘)J ‰分别为海床在x,y方向的底坡源项和摩擦溉其值分另『J为鼢(器憔]_【畿霈】； 7．29xl旷rad／s为地球自转角速度；咖为地球纬度。 2方程离散 Centered 采用CC(Cell Scheme)方式的有限体积方法，把变量存在单元的中心，单元的边界为控制体。 第i号控制单元记为d继，在d以控制单元内，对式(1)两边同时积分。得到二维浅水方程的积分离散形式 收稿日期：2010-08—16；修回日期：2010-lo．19 基金项目：国家973项目(2002CB412410) 作者简介：张宏伟(1980-)，男。辽宁省凌源市人，博士，工程师，主要从事海域管理以及水动力数学模型研究工作。 Biography：ZHANGHong-wei(1980-)，male，engin眦 万方数据 201 1年2月 张宏伟，等黄河三角洲海域的潮汐流动数值模拟 13 击Lm以Fi投L日d以 (2) 式中：d以为面积元积分微元；邙为线积分微元；-+ni=(nh，凡打)=(cos0，sin口)， n曲，％分别为第i号控制单元边界的单位外法向分量。图1为非结构网格控 制单元及其外法向示意图。 · 击琚丽1荟脚瓦