径向点插值法在波浪传播数值模拟中的应用 Application of the radial point interpolation method to numerical simulation of waves.pdfVIP

第22卷第2期 水科学进展 V01．22，No．2 2011年3月 ADVANCESINWATERSCIENCE Mar．。20lI 10313．2241．025 DOI：CNKI：32—1309／P．20l 径向点插值法在波浪传播数值模拟中的应用 季小强1，冯卫兵2，张俞2 (1．南京水利科学研究院，江苏南京210024；2．河海大学海岸与海洋工程研究所．江苏南京210098) 摘要：针对波浪数值模拟中基于矩形网格的数值方法在深水到浅水的网格间距选择与复杂边界处理上的缺陷，以 及基于正交曲线网格和无结构网格的数值方法前处理工作复杂的问题。引人最近在计算力学中发展起来的无网格 求解建立近岸波浪传播数学模型。通过椭圆形浅滩地形和环形河道的波浪传播计算验证，表明该无网格方法可较 为有效地模拟近岸波浪的传播变形，且在处理复杂边界时具有较高的精度。 关键词：无网格法；径向点插值法；波浪；双曲型缓坡方程；数学模型 中图分类号：P731．22 文献标志码：A 文章编号：100l-679l(2011)02-0258-08 目前，波浪数值计算中较为常用的方法是基于矩形网格的有限差分法，但从应用情况看，这类求解方法 有两个缺陷：一是当求解区域涉及深水到浅水时，采用一套计算网格会令模型在满足求解精度的同时，使求 解效率降低，而采用多套网格又使模型的求解因网格间的嵌套与耦合变得困难；二是若计算域含有复杂边 界，如港口或海湾，矩形网格只能近似地以阶梯形边界代替实际边界，这往往会使求解问题的精度降低。为 此，一些学者提出采用基于正交曲线网格的有限差分法¨刮或基于无结构网格的有限元法¨。模拟近岸波浪传 播，可较好地解决上述两个问题，然而，这两种数值方法仍属于有网格方法，计算前需花大量的时间处理网 格信息，给这两种方法的应用带来了不便。 Method)于20世纪 为克服有限差分法和有限元法的局限性，一种新的数值方法——无网格法(Meshless 70年代中期迅速发展起来。无网格法容易贴合边界，可根据求解问题的需要设置节点间距，且无需网格单 元，只需要节点信息。具有前处理简单的优点。目前已提出多种无网格方法，其中在波浪数值模拟中得到应 Point 用的有SPH法H1、基于最小二乘原理的无网格法po等。径向点插值法旧’(RadialInterpolationMethod， RPIM)是近几年在计算力学中发展起来的新型无网格方法，其通过径向基和多项式基的线性组合来构造插值 函数．较之以往的无网格法，具有处理本质边界条件简单，形函数及其导数易于确定等优点。本文基于经典 格式o¨建立近岸波浪传播模型，并用于椭圆形浅滩地形和环形河道的波浪传播计算，取得较为满意的结果。 1 无网格波浪数值计算模型 1．1 控制方程 采用Copeland‘。1导出的双曲型缓坡方程作为本文模型的控制方程，该方程形式如下： V蚕=。 (1) 害+瓦C 收稿日期：2010-05-25：网络出版时同：201l-03-13 网络出版地址：http：／／www．cnki．net／kcms／detail／32．1309．P2241．025．html 基金项目：国家自然科学基金资助项目；南京水利科学研究院青年基金资助项目(Y209004) 作者筒介：季小强(1982一)。男．江苏南通人．博士研究生。主要从事波浪计算与海岸防护研究。 E-mall：guyuesk@163．c锄 万方数据 第2期 季小强。等：径向点插值法在波浪传播数值模拟中的应用 259