K2_SPH方法及其对二维非线性水波的模拟.pdf

第28卷第5期 计 算 物 理 V01．28．No．5 PHYSICS 2011年9月 CH[NESEJOURNALOFCOMPUTATIONAL Sep．，2011 文章编号：1001-246X(2011)05-0659-08 K2一SPH方法及其对二维非线性水波的模拟 郑 兴， 段文洋 (哈尔滨工程大学船舶工程学院，黑龙江哈尔滨150001) 摘要：K2一SPH方法通过泰勒级数展开和联立求解积分方程组的办法得到具有二阶精度的核近似方法．随着核 近似精度的提高，K2一SPH需要对一些关键数值技术进行改进才能成功模拟非线性水波问题，例如自由表面边界和 固壁边界．通过与传统SPH方法计算结果比较，K2一SPH方法在非线性自由表面计算精度和整个粒子系统中有关 变量分布都有显著提高． 关键词：光滑粒子流体动力学；K2一SPH；无网格方法；非线性水波 中图分类号：035 文献标识码：A O 引言 非线性水波是海洋工程与水利工程领域十分普遍的现象，它是模拟波物相互作用的基础．过去几十年， 研究非线性水波的数值方法不断得到更新和发展．最早对水波问题进行研究是基于势流理论的边界元方法 的产生、传播和吸收做了详细研究．贺五州和段文洋¨1采用完全非线性边界元方法对二维摇板造波问题进 行了时域计算．该类方法的优点是可以自然满足质量守恒、不存在能量耗散、计算量较小，但是对粘性效应和 自由表面翻转、破碎现象无法模拟．随着计算机运算能力的不断提高，很多学者发展和应用基于直接求解N- S或者雷诺平均N—S方程的方法卜“，．该类方法最具挑战性的问题是如何处理非稳态非线性自由表面．目 前，常用的方法是VOF方法¨o和LevelSet¨。方法，该类方法都是通过求解每个网格上的一个标量函数来描 述界面的几何特征，标量函数的输运方程形式上完全相同，同属于界面捕获类算法，都有较强的拓扑表达能 力．但是目前该类数值方法对于自由表面，特别是伴随有翻卷、破碎等现象的复杂自由表面问题，所取得的成 就还很有限． 目前，一种被称为无网格的方法引起人们关注，因为在处理大变形的情况下，不需要重新划分网格，该类 方法对于复杂自由表面问题具有较大潜力．本文研究的光滑粒子流体动力学(SPH)方法是无网格方法的一 种．该方法的计算建立在一群离散的粒子上，函数以及导数的求解可以利用周围粒子的相关值积分得到．某 个粒子影响区域的大小称为影响域．影响域内越靠近中心粒子影响越明显，影响域外粒子对计算结果没有影 是对每个粒子的运动轨迹进行模拟，这些粒子能够自动满足自由表面的运动学条件¨2。．目前，SPH方法在自 采用SPH方法对液舱晃荡的相位滞后问题进行研究，之后又实现了SPH方法对液舱晃荡横摇力矩的模 拟¨“．国内对SPH方法的研究已经开展了多年¨¨”o，特别是对强非线性水波问题，采用SPH方法模拟了液 舱晃荡问题【191、推板造波问题Ⅲ】、溃坝问题和楔形体入水问题心1-22J等，但是对一些关键问题的处理，例如 固壁边界处理方法、人工粘性计算方法、密度修正方法等，还处于学习和模仿阶段，需要做进一步的改进和创 收稿日期：2010一10一ll；修回日期：2011—01—09 基金项目：国家自然科学基金及lll计划(1107019)资助项目 yahoo．哪．cn 作者简介：郑兴(1980一)，男，湖北赤壁，讲师，ltl士，主要从事无网格计算方法研究，E-mail：thengxin911060 万方数据 计 算 物 理 第28卷 新．在其他研究领域，国内也开展了许多关于SPH方法的研究旧嗡’． 虽然SPH方法具有巨大潜力，但是仍存在一些关键问题有待解决，例如核近似精度较低、能量守恒性 差、压力计算结果不稳定等．本文介绍的K2一SPH方法通过泰勒级数展开和联立求解积分方程组，得到具有 二阶精度的核近似方法Ⅲ】