基于改进核函数的SPH方法及其数值模拟.pdfVIP

基于改进核函数的SPH方法及其数值模拟

摘要

光滑粒子流体动力学（SmoothedParticleHydrodynamics，SPH）方法是一种无网格

计算方法，在模拟溃坝、涌浪和水下爆炸等大变形问题时具有一定优势，因此越来越多

的学者们使用SPH方法来研究各类大变形问题。但是，目前SPH方法还不够成熟，还

存在粒子间的应力不稳定、计算效率低、边界条件施加困难等问题。其中，应力不稳定

在实际模拟的过程中造成的后果较为严重，会导致粒子的非物理聚集，甚至导致计算崩

溃。

由于流体只能承受压力，不能承受拉力，所以流体中的应力不稳定就变成了压力不

稳定。为解决该问题，本文提出了一种改进的SPH方法，可以改善流体的压力不稳定现

象，并将改进方法应用于波浪水池以及浮体与波浪相互作用的模拟中，模拟结果表明，

本文改进的SPH方法在理论上是可行的。

首先，针对压力不稳定现象，本文提出了一种新型核函数，使得该核函数的一阶导

数的变化趋势与压力的变化趋势相符，在理论上，改进核函数可以改善压力不稳定现象。

其次，为验证改进核函数的效果，用不同核函数对波浪进行了数值模拟。模拟结果

表明，与传统核函数相比，改进核函数使粒子分布更规则，验证了改进核函数对于改善

压力不稳定现象的有效性。

最后，用改进的SPH方法模拟了浮体在波浪中的运动情况，对模拟结果进行水动力

分析，验证了改进核函数的SPH方法应用于浮体与波浪相互作用问题的可行性。

关键词：光滑粒子流体动力学；无网格法；核函数；浮体与波浪相互作用

基于改进核函数的SPH方法及其数值模拟

ABSTRACT

SmoothedParticleHydrodynamics(SPH)methodisameshlessmethod,whichhascertain

advantagesinthesimulationoflargedeformationproblemssuchasdambreaks,surgewaves

andunderwaterexplosions.Forthisreason,anever-increasingnumberofscholarsstudyvarious

largedeformationproblemsusingtheSPHmethod.However,thecurrentSPHmethodisnot

matureenough,andtherearestillproblemssuchasstressinstabilitybetweenparticles,low

computationalefficiency,anddifficultyinapplyingboundaryconditions.Amongthem,the

consequencesofstressinstabilityintheactualsimulationarerelativelyserious,whichcanbring

aboutthenon-physicalaggregationoftheparticles,andevenbringaboutthecomputationto

crash.

Duetothefactthatfluidscanonlywithstandpressionandcannotwithstandtension,stress

instabilityinfluidsbecomespressioninstability.Tosolvethisproblem,animprovedSPH

methodisproposedinthisarticle,whichcanimprovethepressioninstabilityoffluids.The

improvedmethodisappliedtothesimulationofthew