OM模型中的研究与应用.PDFVIP

动边界技术在POM模型中的研究与应用 申霞，姚琪，王鹏 河海大学环境科学与工程学院 摘 要：POM P（rincetonOceanModel）是基于原始方程组的三维斜压海洋模型，其垂向采用 σ坐标，便 于引入大陆架地形，因此被广泛地应用于河口近岸海域的潮流数值模拟中。但是POM采用固定边界，无 法模拟河口近岸海域由于涨落潮而引起的漫滩流。本文概述了国内外动边界技术的研究进展，分析了各种 方法的利弊，提出了一种适合于POM模型的干湿法动边界处理技术。该方法以单元为干湿判断对象，保 证了水量平衡和计算精度。将该动边界技术引入POM模型中，成功地模拟了南通海域的潮流场，较好地 重现了涨落潮造成的潮间带间歇性出没。验证结果表明计算值与实测值吻合较好，真实地反映了该海域的 水动力特征。 关键词：海洋模型；干湿法；潮间带；南通海域 中图分类号：P731.2 文献标jDWI:A 文章编号：1001-6392(2006)01-0001-07 POM模型是Blumberg和Mellor于1977年建立的，其后该模型不断被改进，并广泛应 用于模拟各种物理海洋问题。POM模型垂向采用o坐标，便于处理复杂地形。但是，其采 用固定边界，无法模拟河口近岸水域常见的漫滩流，所以必须在POM 中引入动边界技术。 动边界技术于 1968年首次由Reid和 Bodine提出，并将其应用于显式有限差分模 型。1970年Leendertse2[首，次将动边界技术应用到ADI隐格式模型中。随后国外发展了多 种动边界技术来处理各种实际问题。Lynch3]运用有限元法实现了岸界的连续运动并满足了 运动学的边界条件，但是由于计算量大而难以推广；Johns4[]通过含时间的代数型坐标变换 达到坐标与岸界的随时拟合，但是其代数型变换是单方向，只能解决局部变边界问题，而且 当岸界弯曲程度大时将导致计算的不稳定；Leclerc等嘲开发了二维有限元模型模拟河口海 岸的滩地出没；Hervouet等嘲提出有限算法 (finitealgorithms)处理湿单元；Defina等川提 出改进的浅水方程来描述边界的干湿变化；Bates[8]在总结了浅水流动动边界问题的基础 上，提出镶嵌网格技术处理动边界问题。 国内应用较多的动边界处理方法有冻结法、窄缝法、开挖法、自适应网格法、线边界法 和干湿法等。冻结法91‘适用于宽浅、底坡较坦的露滩问题，对潮滩相间的河口近岸水域， 会因水量的过分 “冻结”而失真。窄缝法㈣适用于岸边滩的处理，对水域内的浅滩处理则 不适用。开挖法H¨只适于滩地面积占整个计算区域影响不大的情况，而对于水位较浅，坡 高较缓的海岸，这种处理则易失真。自适应网格法12‘1只是针对二维浅海潮波方程组，不能 很好适应三维方程的需要。线边界别131适用于计算域中的水利构筑物尺度远小于空间步长 而无法用网格单元表示时。干湿法的基本思想是先设定一个临界水深，如果单元水深小于临 界水深，判断该单元为干，将其从计算域中去除，反之为湿。本文通过比较单元水深与临界 水深和前一时刻水深的大小，决定单元各边通量是否为零，该方法物理概念更加明确，能够 保证各单元的物质守恒。将其引入POM 模型，建立了基于动边界技术的POM模型，并应 收稿日期收修改稿日期 用于南通海域的潮流数值模拟。 1 基于动边 界技术 的POM模 型 1.1 基本方程 o坐标系下三维潮波的连续性方程、动量方程和物质输运方程如下： ．UD-+业生 +芝-+鲁 =o (1) ar． ay a口 旦．．塑 ar，2D aⅣ∞ au +丝二 u二aD-ap≥Id盯’一 gD塑． 犷 + ax +一『 +警 一佃 e[等 一D 。石； ． +。 aa AMD( l+aa 等) 2‘’ +ic¨。D等 ，+专 [A。Dc等． +熹 ）]+