一种改进的SPH固壁边界处理方法.pdfVIP

第二十一届全国水动力学研讨会暨第八届全国水动力学学术会议 暨两岸船舶与海洋工程水动力学研讨会文集 一种改进的一种改进的 SPHSPH 固壁边界处理方法固壁边界处理方法∗ 一种改进的一种改进的 SPHSPH 固壁边界处理方法固壁边界处理方法 龚 凯 1 刘 桦 1 (1. 上海交通大学船舶海洋与建筑工程学院，上海，200242) 摘要摘要: 提出一种适用于光滑质点水动力学(SPH)方法的改进的边界处理方法。在这种方法 摘要摘要 中，边界粒子的压力可通过其周围的流体粒子的压力插值得到 ，从而改进了耦合边界法在边 界上压力不准的问题。运用这种改进的边界处理方法模拟了二维方形水箱 中的非线性晃荡问 题。模拟结果与实验结果吻合较好，证明了此改进的边界处理方法是有效的。 关键词关键词: 光滑质点水动力学(SPH);边界处理;晃荡 关键词关键词 1 引言 SPH 方法是一种纯拉格朗日质点方法[1-2] 。尽管已经发展多年，非常适合处理复杂的大变 [3-4] [5] 形边界问题 ，但固壁边界的处理仍是 SPH 方法的主要难点之一 。目前还没有一种公认较 [3] 好的通用边界处理方法，就以往的文献来看，比较常用的有三种处理方法：边界力法 ，镜像 [4-5] [6] 法 ，耦合边界法 。这三种方法都各有优缺点。 在边界力法中，固壁边界上布置一系列的边界点，假定边界点对靠近它的流体质点施加 一个大小适当的中心排斥力，以阻止流体质点穿过固壁边界。其优点在于不受固壁边界形状 的影响，无论固壁边界形状如何，都可以很方便地形成边界点，缺点在于容易形成初始扰动 且排斥力是人为给定的，目前仍没有最好的排斥力模型。 镜像法是在固壁边界之外布置和内部真实粒子关于边界对称的虚拟粒子，他们和真实粒 子密度相同，速度沿边界镜面对称，从而保证无穿透边界条件。这种方法守恒性比较好，精 度高，但是生成自由表面的镜像粒子时过于复杂，计算时间大大增加，固壁边界形状复杂时 更是如此。 耦合边界法是把固壁离散成一系列的边界点，令这些边界点具有流体质点的部分属性(质 ∗博士点基金项目(20060248046)和国家自然科学基金项目资助. - 347 - 第二十一届全国水动力学研讨会暨第八届全国水动力学学术会议 暨两岸船舶与海洋工程水动力学研讨会文集 量，密度，压力等) ，在求解方程时，在边界点上像在流体质点上一样求解连续方程和状态方 程，从而将边界问题耦合到方程的求解过程中，而不需要提出显式的边界条件。其优点在于 无需对边界进行特殊处理，缺点在于通过求解连续方程和状态方程所得到的边界点上的压力 并不是真实压力。 针对耦合边界法的不足，改进了边界粒子压力的求解方法，其压力通过其周围的流体粒 [5] 子的压力插值得到。边界压力插值方法是从 Oger 等 给出的插值方法改进而来，主要不同有 三点：首先，Oger 等的文章中，其边界插值方法只是用来求边界上某点的压力，其边界处理 方法采用的是镜像法，而本研究中则是用于求边界粒子上的压力，并将此压力用于固壁边界 处理。其次，在Oger 等的插值方法中，当有几个流体粒子在边界上的投影相同时，则其中的 部分点的压力贡献会被忽略，而本研究的插值方法则可以避免这一缺陷。最后，在Oger 等的 插值方法中，每个投影点的压力贡献在每个时刻都是不同的，因此会增加很多计算量，而本 研究的插值方法则没有这一问题。 2 SPH 形式的控制方程