首页/文章/ 详情

同为粒子法的SPH和MPS大PK,谁更胜一筹?

1月前浏览1024



SPH(Smoothed Particle Hydrodynamics)方法和MPS(Moving Particle Semi-Implicit)方法在模拟自由液面的PK:         

1. 自由面形状:

   - SPH方法易于给出更为清晰、光滑的自由面形状。主要是因为SPH方法可以使得粒子移动较为平均,易于形成光滑的自由面。
   - MPS方法得到的粒子分布较为凌乱,尤其在自由面变化剧烈时,MPS的自由面形状不如SPH方法光滑。主要是因为MPS方法中的核函数在粒子接近时导致粒子数密度迅速增加,产生较大的压力,从而导致粒子“飞溅”。
          

2. 收敛性:

   - 随着初始粒子间距的减小,SPH和MPS方法的计算结果均趋于收敛,并收敛于相近的结果。
   - MPS方法的收敛速度比SPH方法快。    
          

3. 时间步长对结果的影响:

   - 对于SPH方法,时间步长对结果的影响较小,不同的时间步长给出的计算结果十分一致。
   - 而对于MPS方法,时间步长的变化会导致计算结果的轻微浮动。。
          

4. 计算效率:

   - SPH方法中压力是由状态方程显式给出的,具有很高的效率,但受声速的限制,需要较小的时间步长以满足CFL条件。
   - MPS方法基于Poisson方程隐式求解压力,计算量较大,且随着粒子数的增加,计算量增加较快。所以,在计算大型问题时,SPH方法更为合适。
   - 由于SPH 的时间步长较MPS方法小10倍,在粒子数较少时,MPS的计算效率较高,但随着粒子数的逐渐增加,MPS的运行时间增加较快,此时SPH 的效率较高。    
   - 对于粘性流动,MPS方法在每个时间步内需要执行两次邻居粒子搜索,而SPH仅需一次,故用时会略少一些
          

5. 数值格式上的差异:

   - SPH方法中,压力是基于状态方程显式计算的,时间积分采用显式格式。
   - MPS方法中,压力是通过压力Poisson方程隐式求解的,时间积分采用半隐式格式。
   - SPH方法的核函数不仅作为权重函数,其导数直接参与梯度等微分算子的离散格式中。
   - MPS方法的核函数仅作为权重函数,梯度的计算还需借助径向函数来完成。    
          

6. 边界条件处理:

   - 两种方法都需要对固壁边界条件进行特殊处理,以避免粒子穿透和计算失真。
          

7. 适用性:

   - 对于处理带有自由面问题的流动,尤其是自由面发生大变形的情况,SPH和MPS方法都具有良好的表现。
   - SPH方法因其较高的计算效率和能够给出光滑自由面形状的特性,可能更适合于大型复杂流动问题的模拟。
   - MPS方法在每个时间步内需要执行两次邻居粒子搜索,而SPH仅需一次,因此在粘性流动问题中,MPS的计算量较大。
   - MPS方法虽然收敛速度快,但由于计算量大,可能更适合于中等规模的问题,或者在可以接受较长计算时间的情况下使用。    
          
SPH方法和MPS方法在处理溃坝问题这类大变形自由面问题时都表现出了较好的可靠性。然而,SPH方法在计算效率和自由面形状的光滑度方面具有优势,而MPS方法在收敛速度上表现更好。SPH方法适合求解大型、中等、小规模复杂流动问题,MPS方法更适合于中等和小规模的流动问题。    






来源:CFD饭圈
通用
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-09-08
最近编辑:1月前
CFD饭圈
硕士 分享CFD文章,感谢关注
获赞 22粉丝 23文章 378课程 0
点赞
收藏
作者推荐

极度深入,探讨下涡流形成和雷诺系数的关系

点击上方蓝下图是Nature Communications 期刊计算了不同雷诺系数下的涡流形成。涡流(vortex)是流体流动中一种典型的旋涡状结构,其形成与流体内部的流动不稳定性密切相关。涡流的产生、演化及其影响因素众多,其中一个重要的关联参数就是雷诺数(Reynolds number)。雷诺数是流体动力学中衡量惯性力与粘性力相对重要性的无量纲数,其定义为:其中:·ρ 是流体的密度,·v 是特征流速(如平均流速或最大流速),通常取于所考虑的流动区域或特征长度处,·L 是特征长度(如管道直径、物体特征尺寸等),代表流动的几何尺度,·μ 是流体的动力粘度。 雷诺数与涡流形成的关系可以从以下几个方面来详细阐述: 一、转捩雷诺数雷诺数是判断流体流动状态(层流与湍流)的重要指标。在低雷诺数(Re << 1)条件下,流体内部的粘性力显著,足以迅速抑制任何微小的流动扰动,使得流体分子保持有序、平行的层状流动,即层流。然而,随着雷诺数增大,惯性力逐渐增强,对流动扰动的影响加大。当雷诺数超过某个临界值(通常称为转捩雷诺数,记作 RecritRecrit)时,流体内部原有的层流结构开始不稳定,扰动得以放大并传播,最终形成涡旋结构,流动进入湍流状态。 二、涡旋生成机制在雷诺数较高的湍流流动中,涡旋的形成主要有以下几种机制:·剪切层不稳定:在层流与湍流过渡区域,流体速度梯度较大,导致剪切应力增加。当雷诺数超过阈值时,剪切层内的小扰动因惯性力的作用而不断增长,形成涡旋。 ·卡门涡街:在流体绕过圆柱、圆管或其他障碍物时,由于边界层分离和逆压梯度,会在物体后方形成交替排列的涡旋对,即卡门涡街。雷诺数对此现象的发生有直接影响,高雷诺数下涡街更易形成且更为明显。·旋涡脱落:在特定的几何结构(如弯头、扩缩管等)中,流体流动受到强烈扰动,边界层内的涡旋可能从壁面剥离并卷入主流,形成大规模的涡旋结构。雷诺数高时,这种脱落过程更加频繁且强烈。 三、雷诺数与层流-湍流转变:·层流:当雷诺数远低于某个临界值时,流体流动表现为层流(Laminar flow)。在层流状态下,流体质点沿平行的流线有序流动,彼此间几乎没有横向混合,流场平滑且可预测。此时,粘性力占主导地位,能够有效地抑制流场中的微小扰动,防止它们发展成涡旋。·转捩雷诺数(RecritRecrit):在特定流动条件下,存在一个特定的雷诺数,当流体的雷诺数接近或超过这个值时,层流开始变得不稳定,流动开始向湍流过渡。这个临界雷诺数被称为转捩雷诺数。对于简单管道流动或平板边界层,转捩雷诺数的经典经验值约为2300,但实际流动中,由于流道几何形状、流体性质、壁面粗糙度等因素的影响,转捩雷诺数可能会有所不同,范围大致在2000至2300之间,甚至更高或更低。 ·湍流:当雷诺数超过转捩雷诺数后,惯性力逐渐占据主导地位,使得原本被粘性力抑制的小扰动得以放大并持续传播。这些扰动在流场中发展成涡旋结构,并进一步相互作用、合并、分裂,形成复杂的涡旋网络,这就是湍流(Turbulence)。湍流中流体质点的轨迹变得高度不规则,流场呈现出强烈的时空变化和混沌特性。 四、涡流形成的具体雷诺数范围:尽管转捩雷诺数提供了从层流到湍流转变的粗略界限,但在实际应用中,涡流的形成可能在更高的雷诺数下更为显著:·过渡区间:在雷诺数略高于转捩雷诺数(如2000至4000之间)时,流动可能处于层流与湍流的过渡状态,涡旋开始出现但并不十分明显,局部可能存在涡旋结构,整体流动仍然部分保持层流特性。·典型湍流区间:当雷诺数进一步增大至4000以上时,流动通常被认为是充分发展的湍流。在这个范围内,涡旋生成、演化和相互作用非常活跃,形成了典型的湍流特征,如强烈的涡旋混合、能量耗散和速度脉动。 需要注意的是,上述雷诺数阈值是针对特定流动配置(如圆管流动、平板边界层等)的一般经验准则。对于其他复杂几何形状(如弯管、扩缩管、螺旋桨等)或特定流体(如非牛顿流体、多相流等),涡流的形成可能在不同的雷诺数下发生。此外,壁面粗糙度、温度效应、外部扰动等因素也会影响涡流形成的雷诺数阈值。 五、涡旋的演化与相互作用在高雷诺数湍流中,涡旋不仅易于生成,还表现出强烈的非线性行为。涡旋可以相互吸引、排斥、合并、分裂,形成复杂的涡旋结构网络。雷诺数越大,涡旋的生成、演化速度越快,涡旋尺度分布更广,能量传递跨越多个尺度,从而导致流动的复杂性和不可预测性增加。 六、雷诺数对涡流影响的实验观察实验研究表明,通过改变流速、管径或流体性质(从而调整雷诺数),可以直观地观察到流场从层流向湍流的转变以及涡旋的出现与发展。例如,在管道流实验中,随着雷诺数增大,通过可视化技术(如粒子成像测速PIV、激光诱导荧光LIF等)可以看到流场从规则的层流条纹逐渐转变为混乱的涡旋结构。 七、雷诺数在工程应用中的意义理解雷诺数与涡流形成的关系对于许多工程应用具有重要意义。例如,在设计管道、叶片、船体等流体设备时,通过合理控制雷诺数,可以预测和避免不利的涡流现象(如旋涡脱落引起的振动、阻力增加),或者利用涡流实现特定功能(如混合、传热增强)。此外,在环境流、生物流等领域,了解雷诺数对涡流形成的影响有助于解释自然现象和优化相关技术。 总结来说,雷诺数是描述流体流动状态转变、涡旋生成机制、涡旋演化过程以及工程应用中涡流控制的关键参数。它量化了惯性力与粘性力之间的相对强度,高雷诺数倾向于引发和维持涡旋结构,而低雷诺数则有利于保持层流流动。通过调控雷诺数,可以有效地预测和管理流体系统中的涡流现象。 ‍‍‍来源:CFD饭圈

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈