水气二相流光滑粒子流体动力学方法研究进展

从状态方程、连续方程、动量方程和其他数值技术等方面,对光滑粒子流体动力学(SPH)解决水气二相流问题的数值技术和最新进展进行了综述,并对各项技术进行了评价。最后,提出了未来的研究展望。

面向光滑粒子流体动力学(SPH)流体的湍流细节模拟

研究基于光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)方法的恢复涡流和湍流细节的方法,缓解因采用人工黏度而产生的数值耗散.根据由黏性力引起的动能损失率对离散点的涡度进行修正.通过流函数将涡度变化恢复到速度场,避免求解Biot-Savart积分这一耗时过程.该方法基于旋度计算,不仅可以增强现有的涡流,还可在潜在位置产生新的湍流,并且其易于集成到现有SPH方法中.实验证明,该方法可以逼真地模拟湍流细节,并实现不同强度水平的增强效果.

基于光滑粒子流体动力学法的飞机风挡鸟撞失效破坏模拟

运用光滑粒子流体动力学(Smoothed Particle Hydrodynamics,SPH)方法模拟了鸟撞飞机风挡的过程,鸟体模型采用SPH法建立,风挡模型采用Lagrange法定义,鸟体以515 km/h、562 km/h和600 km/h的速度分别撞击风挡对称线上前1/3点。通过仿真模拟,获取了鸟体撞击风挡的应变、位移、应力曲线,风挡失效损伤演化过程,并与试验结果对比,证明所建鸟撞风挡有限元模型可以准确预测风挡的鸟撞动态响应。研究表明:鸟体撞击风挡对称线前1/3点时,风挡最大变形发生在风挡中间部位;风挡在弯曲作用下首先会在内表面发生失效破坏;撞速为600 km/h的风挡发生失效破坏后,在撞击区域形成多条与风挡对称线呈大约50°夹角的裂纹及少部分横向裂纹。

基于PCISPH的流体粒子飞溅改进方法

流体飞溅是自然界中最常见的流体现象。流体模拟是计算机图形学的一个重要研究分支。流体模拟已经广泛应用于电影、游戏和其他工业数值领域。由于飞溅场景中流体的密度和压力的速率变化非常大,因此模拟对离散化解的精度要求很高。对于流体粒子飞溅模拟数值不稳定,且缺乏真实效果,针对经典的预测校正不可压缩SPH(PCISPH),该文提出了基于粒子表面流体飞溅改进的泊松压力解的方法,通过替换压力源项提高了模拟精度。实验结果表明,该方法中流体的压力分布较之经典方法更接近现实的流体粒子飞溅模拟效果,并且改进后的方法经过剧烈运动后帧率仍符合实时性仿真效果的要求,保证了良好的压力稳定性和真实感。该文通过Anaconda3集成Python3和Taichi环境实现流体仿真实验。

光滑粒子动力学方法的发展与应用

光滑粒子动力学(smoothed particle hydrodynamics,SPH)是一种拉格朗日型无网格粒子方法,已经成功地应用到了工程和科学的众多领域.SPH使用粒子离散及代表所模拟的介质,并且基于粒子体系估算和近似介质运动的控制方程.本文分析和综述了SPH模拟方法的发展历程、数值方法与应用进展.介绍了SPH方法的基本思想;从连续性、边界处理、稳定性和计算效率4个方面阐述了SPH方法的研究现状;介绍了SPH方法近年来在可压缩流动、不可压缩流动以及弹塑性材料高速变形与失效方面的一些典型应用;并对SPH方法的发展与应用进行了预测与展望.

光滑质点流体动力学(SPH)及其算法特性

基于Lagrange观点的无网格方法计算是建立在一群可移动的自适应性的质点基础之上,对特大变形区域和自由表面的捕捉有着较高的精度。光滑质点流体动力学(SPH)是无网格方法的一种,本文阐述其理论基础和计算方法,主要包括:SPH的基本概念、核近似和粒子近似的基本原理、核函数的构造方法、SPH形式的流体运动方程、边界处理方法和数值算例等多方面内容。研究工作在了解SPH方法的理论特点和潜在优势的同时,也为相关问题的深入研究和改进奠定了理论和算法基础。

光滑粒子动力学方法及其应用

光滑粒子动力学(Smoothed Partic le Hydrodynam ics,SPH)方法是近年来得到广泛发展和应用的无网格方法的一个重要分支,它是一种纯Lagrangian方法。本文对现有的光滑粒子动力学方法进行了综述,介绍了该方法的理论基础、连续介质守恒方程、方法稳定性的改善等,重点阐述了边界条件的处理,并给出了SPH方法的算例。最后,介绍了SPH方法的最新进展状况。

光滑粒子流体动力学方法研究

光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)方法在模拟大变形问题时具有明显的优势,但是由于粒子的不连续性,致使其计算精度较低。文中对光滑粒子流体动力学方法中函数及其一阶导数的核估计进行详细研究。讨论传统的SPH方法和改进的CSPH(corrective smoothed particle hydrodynamics)方法的离散思想,在Taylor展开的基础上引入修正的MSPH(modified smoothed particle hydrodynamics)方法,并推导一维和二维情况下函数的核估计和函数的一阶导数核估计的离散形式。最后通过数值算例,对三种不同的SPH方法的计算精度进行详细比较,结果表明,CSPH和MSPH方法可以极大提高边界粒子的计算精度,在二维情况下,MSPH方法的计算精度要优于CSPH方法。

基于光滑粒子流体动力学方法的溃坝水流模拟

受全球气候变暖的影响,由极端天气引发的类似溃坝等问题发生的概率大大增加,深入研究溃坝水流的水动力特性势在必行。在分析光滑粒子流体动力学基本原理的基础上,提出了一种改进的边界处理方法,即将接近壁面的流体视为层流,在耦合动力边界附近引入层流黏性近似边界层理论。采用该方法对溃坝水流进行数值模拟,将SPH数值模拟得到的外轮廓、自由液面高度以及压力与实验结果进行了比较和分析。结果表明:改进的边界处理方法较完整地得到了水流与壁面相互作用而产生的多种复杂的物理现象,其外部轮廓与实验非常吻合;自由表面融合过程中液面间冲击的能量耗散会导致融合后的液面高度存在一些差异;不同监测点处压力随时间的变化基本落在置信区间之内。数值模拟结果与实验结果吻合度较高,验证了改进方案的可靠性和计算结果的准确性。

基于修正剑桥模型的光滑粒子流体动力学的土柱崩塌模拟

土柱崩塌是一种被广泛研究的岩土大变形模型。考虑非关联流动法则的修正剑桥(MCC)模型描述土体行为;并将其导入光滑粒子流体动力学(SPH)算法来模拟土柱崩塌。在算法设计中,利用有限元方法中常用的中心回映算法来计算应力积分,以提高算法的精度。应用所设计的SPH算法模拟初始高宽比分别为1.5、3和5.5的土柱崩塌模型;并与实验结果比较,表明SPH算法在研究土柱崩塌问题的可靠性。

交互式角色运动特效的光滑粒子流体动力学仿真

长期以来,由于流体仿真和物体变形计算都具有相当程度的复杂性,使得流体与刚体的交互模拟,特别是和带有复杂动画的角色交互的效果,难以达到实时计算和渲染。在此,笔者提出了一个新的方法,用于生成沿着角色运动而产生交互的流体特效。为了实现这类效果的生成,控制流体特效与运动角色的交互,首先针对角色运动轨迹进行跟踪,根据轨迹的几何性质而生成初始状态的流体特效;然后借助光滑流体动力学(SPH)对流体粒子进行仿真。其中针对基于SPH技术的复杂性,流体仿真的过程借助GPU并行计算的能力,采用了一种新的高效粒子搜索算法,最终实现普通用户级个人计算机上实时渲染具有流体运动特征的角色运动特效。

基于光滑粒子流体动力学方法与TANH本构方程的钛合金切屑形态预测

引入TANH本构方程,建立了钛合金切削过程的光滑粒子流体动力学(SPH)方法的切削模型.新模型解决了基于有限元方法(FEM)的传统切削模型经常出现的网格畸变问题,另外利用控制变量法标定出TANH本构方程的修正系数.相较于传统的Johnson-Cook本构模型,该模型将材料的应变软化现象考虑在内,更加准确描述出钛合金在大应变和动态再结晶情况下材料的动态力学性能.同时,新模型很好解释了钛合金切削过程中锯齿形切屑的形成过程与形成机理.实验结果与模拟结果对比显示,基于SPH方法与TANH本构方程的切屑模型可以准确可靠地预测钛合金切屑形态与切削力.

基于SPH方法的电场下液滴动力学数值模拟

建立了适用于电场作用下液滴动力学的光滑粒子流体动力学(SPH)模型,并验证了模型的准确性。运用此模型研究了电场强度对液滴聚并时间的影响。结果表明:电场强度增加,液滴接触时的极化变形程度增加,不利于液滴融合,当电场强度过高时,液滴会发生破碎;当电场强度达到2860 kV/m时,液滴粒子迅速向边壁移动,液滴在其端部向四周扩张之前,已经发生破碎,液滴的破碎时间趋于稳定;当电场强度达到3580 kV/m时,液滴两端受到了更大的电场力,液滴被迅速拉伸变形,液滴从中间部位直接破碎为左右两个子液滴;当液滴直径达到6 mm时,固壁边界抑制液滴变形,液滴的破碎时间迅速上升。

基于SPH模型的尾矿库溃坝风险模拟研究

作为一种具有高势能的危险源,尾矿库可能发生溃坝事故,进而产生泥石流,造成严重的人员伤亡和经济损失。传统的深度积分方法在模拟泥石流演进过程时存在局限性,而光滑粒子流体动力学(Smoothed Particle Hydrodynamics,SPH)是一种新型的数值模拟方法,可以弥补传统方法的缺陷。以攀枝花市某尾矿库为例,构建溃坝泥石流的SPH模型,模拟不同工况下泥石流的演进过程,获得泥石流的流量数据,研究溃口宽度与深度对泥石流演进特征的影响。模拟结果与实际情况较为吻合,因此SPH模型可以有效模拟尾矿库溃坝风险。研究结果为尾矿库安全管理和风险防控提供新方法,有利于保障人民群众的生命财产安全,促进矿业可持续发展。

基于光滑粒子流体动力学方法的液滴蒸发问题数值模拟研究

基于光滑粒子流体动力学方法(Smoothed Particle Hydrodynamics,SPH),开展了SPH新算法在蒸发燃烧领域的研究。建立了适用于SPH方法的蒸发数值模型,推导了基于傅立叶热传导公式和菲克扩散定律的SPH离散方程;借鉴VOF方法(Volume of Fluid)的思想,提出了SPH粒子的液相质量分数的概念,以有效表征蒸发过程中的相变问题。采用SPH方法对高温环境中单个液滴的蒸发过程进行数值模拟,结果符合D2定律,与理论模型相一致;在强迫对流环境中,液滴的蒸发过程受到对流作用及表面张力的影响,蒸发速率加快;进一步对双液滴在静止、对流环境中的蒸发过程进行数值模拟研究。结果表明,液滴的间距、滴径对多个液滴的蒸发过程影响至关重要,液滴间距至少在两倍的液滴直径以上,相互之间的影响才可以近似忽略。通过本文研究,拓宽了SPH方法在蒸发相变领域的应用范围,研究结果也能够为进一步的燃烧问题研究奠定基础。

三维光滑粒子流体动力学程序的并行计算

光滑粒子流体动力学方法，是近20多年来发展起来的一种无网格Lagrange流体动力学算法，它既保留了拉氏计算的描述物质界面准确的优势，又具备欧拉方法的长处，因此适宜计算大变形问题。SPH的基本思想是插值。在SPH方法中，物质被离散为一系列“粒子”，各种宏观物理量（密度、压力速度、内能等）被定义在粒子中心，粒子的物理量及其空间导数可以通过邻近的相互作用粒子的物理量插值得到，这样拉氏流体力学偏微分方程组就变得容易求解，这也是SPH方法实质性优点所在。

基于HBP本构模型的泥石流动力过程SPH数值模拟

本构模型是描述泥石流流变特性的关键,也是决定其动力过程数值模拟准确性的核心问题之一。泥石流流体属多相混合物,现有的研究已证实其存在剪切增稠或剪切变稀的现象,传统基于Bingham及Cross线性本构关系的数值模型难以准确描述泥石流流变特性。文中探讨了Bingham模型在低剪应变率下的数值发散问题,在光滑粒子流体动力学(SPH)方法框架上建立了整合Herschel-Bulkley-Papanastasiou(HBP)本构关系的稀性泥石流动力过程三维数值模型。相比传统基于浅水波假设的二维数值模型,所述方法从三维尺度建立SPH形式下的泥石流浆体纳维-斯托克斯方程并进行数值求解,可获取泥石流速度场时空分布及堆积形态,同时采用HBP本构关系描述泥石流流变特性,能在确保数值收敛的前提下反映泥石流流体在塑性屈服过渡段及大变形状态下应力-应变的非线性变化。为验证提出方法的合理性,结合小型模型槽实验观测进行了对比,结果表明数值模拟与实测结果基本吻合。

基于GIS空间数据的滑坡SPH粒子模型研究

针对光滑粒子流体动力学方法(SPH)在滑坡模拟中建立粒子模型的难题,提出了基于地理信息系统(GIS)栅格数据的粒子排列与插入方法。根据该方法,建立了滑坡SPH粒子模型及相关粒子生成程序,进一步以结合摩尔-库仑破坏准则的SPH宾汉流体模型为核心,实现了运用SPH方法模拟滑坡破坏后三维运动的过程。该SPH模型在对唐家山滑坡的模拟中得到了验证,并预测了金坪子滑坡破坏后的影响范围。结果表明:基于GIS空间数据的滑坡SPH粒子模型具有可行性与良好的适用性。以GIS数据库为基础,开展滑坡灾害的模拟研究,将大大提高对滑坡等地质灾害的仿真分析,为滑坡灾害的预测与防治提供参考。

基于光滑粒子法边坡失稳影响范围的模拟

基于光滑粒子流体动力学(SPH,smoothed particle hydrodynamics),以经典弹塑性本构模型为核心,结合一系列数值处理方法,探讨精确和便捷地对边坡失稳破坏影响范围的模拟方法与技术。文中借助地理信息系统(geographic information system,GIS)平台与C#编程建模,将边坡地形细节与特征创建成为SPH粒子,对模型部分边界粒子排列方式进行改进,最后应用到边坡模拟的实例中。结合GIS平台与Matlab开发了滑坡SPH粒子自动排列程序,将此模型与方法应用在美国华盛顿州Oso滑坡实例中,经过参照对比,数值模拟结果与滑坡现场收集数据较为吻合,验证了光滑粒子法及运动模型在模拟边坡失稳影响范围的适用性。

SPH的核近似和粒子近似

光滑粒子流体动力学(SPH,Smoothed Particle Hydrodynamic)方法是一种具有纯La-grangian性质的无网格方法,在众多学科和领域得到了广泛发展和重要应用。文中主要对用SPH方法近似偏微分方程(PDEs)的核近似法和粒子近似法进行了描述分析,并说明了两个的利弊。除了一些简单情况,PDEs的解析解是很难求出的,因此描述分析这两种方法对寻得近似解就具有十分重要的意义。