三维修正SPH方法实现的海水运动效果模拟 .pdfVIP

三维修正SPH方法实现的海水运动效果模拟

欧训勇;陈美伊;鲍海琴;康小平

【摘要】围绕SPH的标准控制方程,对SPH的标准核函数使用了所谓的线性再生

核近似方法进行修正,这种方式所取得的梯度效果令人满意.修正的SPH方法模型用

来模拟流体自由表面运动.仿真效果展现了水平面上波的形成和传播过程.最后在

OpenGL环境下实现了的逼真模拟效果.

【期刊名称】《微型机与应用》

【年(卷),期】2017(036)012

【总页数】4页(P71-74)

【关键词】SPH;自由表面流体;流体运动

【作者】欧训勇;陈美伊;鲍海琴;康小平

【作者单位】海南热带海洋学院海洋信息工程学院,海南三亚572022;海南热带海

洋学院海洋信息工程学院,海南三亚572022;海南热带海洋学院海洋信息工程学

院,海南三亚572022;海南热带海洋学院海洋信息工程学院,海南三亚572022

【正文语种】中文

【中图分类】TP311.1

光滑流体动力学方法(以下称为SPH)是1977年由Lucy、Gingold和Monaghan

率先提出的，至今该方法已经得到了广泛的研究及应用，并成为通用性和有效性最

强的应用于物理、力学、工程等领域计算离散粒子的方法之一。SPH方法成为有

重要意义的计算方法可以归因于它几个突出的优势：概念简单、容易编程实现、处

理大变形问题。作为一种无网格的纯拉格朗日方法，SPH方法已成功应用于一系

列现象的数值模拟，如动态冲击、金属损伤和断裂、海浪运动、传热过程、水下爆

炸等。

当然，SPH方法也被用来解决各种流体力学问题，包括有重要的水文现象和离岸

工程等。例如，文献[1]中运用这一方法研究水上的重力流和孤立波；文献[2]中分

析了在一个垂直墙面和倾斜的海滩上的孤立波运动过程；文献[3]研究了经典的溃

坝问题，包括水撞击墙体的分析；文献[4-5]使用SPH方法研究了物体入水等自由

表面流动问题；文献[6]研究利用SPH方法模拟了二维和三维波浪在倾斜的海滩上

的破碎问题；文献[7]利用SPH研究了弹性流体结构相互作用的问题；文献[8]利用

不可压缩流的SPH方法研究流体自由表面运动的模拟过程。

以上所提到关于对于流体力学问题的研究，SPH方法的各种变体已经应用到了很

多领域之中。早期的研究成果中，有代表性的，例如文献[1]中，传统的SPH方法

是以水作为一个弱可压缩流体[1]，和其中所有的场变量利用一个单一的插值函数

(称为SPH光滑核)，在其作用域内近似离散成粒子。在后来的研究文献中，不断

提出各种改进的方法，以实现更好的真实流体建模，并加强了计算的准确性[4-5]。

本文提出一种修正的SPH方法模拟海水自由表面运动。

本文中，把水视为一个可压缩的粘性液体，并假设它的运动发生在等温条件下，完

全是由重力引发的运动。流体质量守恒的连续方程如下：

式中,D/Dt是时间微分量，δ是流体密度，t是时间量，υ是流体的速度矢量是微

分算子。动量守恒的欧拉方程如下：

式中，p是流体压力，b指的是外力，本文所指的外力在没有其他作用力时其实就

是重力。如果考虑其他力，那么b就是这些力和重力之和。

假设流体受正压，在这种情况下，压力是由流体密度唯一确定。根据传统的SPH

方法，压力p和密度之间的变化规律如下：

p(δ)=P0

式中，P0和δ0压力和密度的参考值，参数τ的取值如果为水流体是7，为空气

则1.4。相应地，声音的速度表示为：

Cs=

方程(3)的突出作用是表明了压力和密度的直接关系，这就使得没有必要解决任何

额外的压力方程(这是不可压缩流的泊松方程)。因此，快速显示的时间步长算法可

以在计算中采用，避免在泊松方程的影响下使用隐式方法。

在SPH方法中，任何系统都认为是由一套随意离散分布的粒子构成的，每个粒子

都承载相关的物理信息，如：质量、密度、动量、温度等。由于粒子分布是离散的，

独立运动，因此该方法具有完全无网格的特征。场变量是通过特定的插值函数来近

似的，这就是光滑核。这些函数表述了给定空间或质点的场量值。通常核函数由限

定域的非零函数构成，这个限定域称为支持域。设x为位置矢量，a和b为两个粒

子，场函数f(x)在粒子a处的取值由核函数W以粒子a为中心近似估算取得。由

于核支持域的紧支性，只有位于这个域内的粒子才有被插值累加。其方程如下：

式中，fb=f(xb)是粒子b的离散值，Vb是粒子b的体积，N是