COMSOL 软件建模教程:如何模拟自由液面 (二)

在之前的文章中,我们讨论了使用基于场的方法(水平集和相场)为自由液面建模。动网格是另一种可选方法,它能够处理没有发生拓扑变化的自由液面。本文中我们将展示如何使用动网格方法为自由液面建模,并将结果与基于场的方法进行对比。

自由液面问题

为了演示动网格功能,我们继续以介绍相场和水平集方法的博客中的问题为例。实心杆模型一半浸没到小型管道内的水中。值得注意的是,实际上在水平集与相场方法中,动网格功能还可用于指定小矩形杆在液面上来回移动,在动网格方法中亦是如此。

COMSOL 软件建模教程:如何模拟自由液面 (二)的图1
示例问题的几何结构和定义。

自由液面的动网格方法限定了矩形杆的位移并追踪液面的位移。我们通过在动量方程中添加一个源来考虑模型内的重力。为了将动网格方法得到的结果与相场和水平集方法得到的进行对比,我们对壁使用了 Navier 滑移边界条件。滑移长度等于单元长度。

使用动网格为自由液面建模

在 COMSOL Multiphysics® 软件中,对于相同的问题,动网格自由液面建模功能与上一节讨论的水平集与相场法属于完全不同的方法。利用动网格方法,自由液面被建模为分隔两个域的几何面。表面张力及其他表面力被直接施加在自由表面,作为边界条件。

利用自由表面 特征,我们可以计算出任意时间点上的液面流体速度,将之作为自由液面的位移速度。动网格方程的解可以平滑地取代流体域中的网格节点。我们随即计算了移动坐标系中的纳维-斯托克斯方程,并同时求解动网格方程,从而获得移动量。

通过该方法,我们完成了液体的建模,与此同时,我们以表面张力和压力效应的形式对自由液面上方的气体域内的流场进行了处理。因此,我们不能使用当前方程获得气相流场(当然,你可以通过手动添加第二个流体流动接口,或者通过使用“微流体模块”的两相流,动网格 接口中的流体-流体界面 特征改变这种情况)。此外,COMSOL Multiphysics 中的动网格方程无法处理自由液面的拓扑变化,比如破碎波。

在某种程度上,相比于水平集和相场法,利用动网格方法为自由液面建模显得更加简单干脆,因为如前文所述,我们可以直接将表面张力及其他表面力用作边界条件。不求解自由液面上方空气域内的流体流动有利于大大提高计算速度,因为纳维-斯托克斯系统的自由度数量几乎减少到基于场的方法的一半。在这种情况下,我们之所以忽略空气域的影响,是因为水和空气的密度与动力粘度比值很大。所造成的差异将在下一节中详述。

比较动网格与相场法的结果

下图比较了分别使用动网格与相场法计算的自由液面。我们可以看到两种方法的结果非常一致,自由液面的形状和速度场的流线都很相似。

不过模型并非完全相同。在相场法的案例中,自由液面上方的空气域在液面上产生了微小的阻尼效应,而动网格案例中不存在空气域,并且液面只“看见”流体表面上气压恒定不变。换句话说,动网格案例中的自由液面不必移动空气,并且可以利用这种能量使水波更高,表面波动更大。

COMSOL 软件建模教程:如何模拟自由液面 (二)的图2

使用两相流动网格接口(左)和相场方法(右)计算得到的不同时间下的自由液面形状和速度场。

下方动画展示了利用动网格方法求得的动态自由液面,我们可以将它与上一篇中只用相场法生成的动画作比较。可以清晰地看到,与相场动画相比,自由液面的波动幅度更大,反应也更快。这可能是因为动网格中没有空气域,而水平集和相场方法中的空气域会阻碍自由液面运动。

利用动网格方法获得的自由液面动画。

我们还可以将默认的自由液面动网格功能与两相流动网格 接口进行比较,后者能够分析液相和气相的流场。可以看到,不管是包含两相的动网格,还是相场法,流场和速度矢量的大小都非常相似。对于三种情况(动网格、相场及包含两相的动网格),自由液面的形状均相似,但是在此例中,两个动网格案例所对应的形状更加相似。根据仿真结果,结论是空气域对流体的速度场具有一定的阻尼效应。此外,动网格和相场法在处理相边界的壁面上的差异似乎是造成自由表面形状不同(尽管差异很小)的原因。

COMSOL 软件建模教程:如何模拟自由液面 (二)的图3

COMSOL 软件建模教程:如何模拟自由液面 (二)的图4

COMSOL 软件建模教程:如何模拟自由液面 (二)的图5

使用两相流动网格接口得到的结果。

自动重新划分网格与自由液面

针对相场、水平集和动网格方法,我们还考虑到使用自动划分网格功能,以防单元质量太差而不能满足质量要求。质量与单元的最大角有关,还与最大和最小边的关系有关。如果角度很大,或单元被压缩得很厉害,就意味着单元质量很差,若勾选了自动重新划分网格功能,便会触发软件重新划分网格。下图显示了对示例问题应用自动重新划分网格的效果。0.35 秒后单元质量太差,从而触发了软件重新划分网格,这将大幅改善网格的质量。

COMSOL 软件建模教程:如何模拟自由液面 (二)的图6

COMSOL 软件建模教程:如何模拟自由液面 (二)的图7

自动重新划分网格之前(上图)及之后(下图)的网格。我们可以看到小矩形杆右侧的单元在重新划分网格前被拉长了,但在重新划分网格后变得更加各向同性。

流-固耦合与自由液面

在示例问题中,我们规定了矩形杆在液面上的位移。我们可以在矩形杆上施加作用力,并计算矩形杆在流体的反作用力下产生的位移,从而直接地对问题进行扩展。这就得到了所谓的流-固耦合(FSI)问题。

COMSOL Multiphysics 软件 5.3a 版本引入了新的流-固耦合接口,方便用户定义包含流-固耦合的自由液面问题(或两相流问题)。下图展示了经典的流体壁塌陷问题,但是现在一个小小的障碍物挡在流动路径上。我们将表面张力的影响、液体界面的位置及流动与被流体冲击的晶须的结构位移、应力和应变完全耦合。

根据示例问题得到的一些结论

我们在本系列博客的第一篇中得出结论:当对受表面张力影响的系统的自由液面建模时,相场法在性能与精确度平衡方面优于水平集方法。示例问题的解表明采用自由表面特征的动网格方法在性能与精确度平衡方面甚至更出色。不过,它也存在两个缺点:

  • 不能处理拓扑变化。

  • 未默认分析自由液面上方的空气(或其他气体)。

我们可以得出结论:如果没有预见拓扑变化,动网格方法是自由液面建模的首选。如果存在表面张力和拓扑变化,那么相场方法是最优选择。或许在之后的博客文章中,我们可以选择表面张力可忽略不计,且有或无拓扑变化的问题,全面比较这三种方法。

来源:COMSOL

COMSOL

COMSOL 软件建模教程:如何模拟自由液面 (二)的评论0条

    暂无评论

    COMSOL 软件建模教程:如何模拟自由液面 (二)的相关视频课程

    COMSOL 软件建模教程:如何模拟自由液面 (二)的相关案例教程

    业界领先的多物理场仿真软件COMSOL Multiphysics® 最新版本推出用于仿真数据管理的平台通用功能,以及用于不确定性量化分析的新模块,并为已有产品带来全面的功能更新和性能提升。 美国马萨诸塞州,伯灵顿(2021年12月14日)—— 业界领先的多物理场仿真解决方案提供商 COMSOL 公司发布了全新的COMSOL Multiphysics® 6.0版本。 新版本推出的“模型管理器”为 C
    3 总结和展望 在锂离子电池的研究中,仍存在许多科学问题尚未解决,这些问题严重影响着锂离子电池的安全性能和使用寿命。例如,锂枝晶的生长演化、SEI膜的形成和破裂演化、正极颗粒在循环中的破裂、电池寿命预测、热失控、以及电池组的电池状态实施监测和管理等问题。这些问题涉及到电场、浓度场、力场和温度场等多个物理场之间的耦合,很难通过单一的实验表征手段对各个驱动力进行分别观测,更难以给出多场耦合的综合结果。
    在 COMSOL Multiphysics 中可以使用 AC/DC 模块中的非线性磁性材料数据库中的非线性磁饱和曲线进行频域仿真。您也可以使用有效非线性磁曲线计算器仿真 App 将关联的 B-H 或 H-B 曲线(以前仅支持稳态和瞬态研究)转换为有效的 B-H 或 H-B 曲线。这篇文章我们将讨论如何在频域仿真中使用这个仿真 App。 频域中的非线性磁性材料 一个常见的建模假设是在本构关系中指定线
    如果你抬头看看夜空,特别是在远离城市灯光的地方,就会发现很多闪烁的星星。古往今来,美丽的星空是诗人和浪漫主义者的灵感来源,但也一直是研究夜空和宇宙的天文学家们面临的一个挑战。它使自适应光学得以发展。如今,自适应光学被用来提高光学系统的功率,以消除光学介质所带来的不良影响。 自适应光学:消除光学介质带来的不必要影响 当一束来自天文物体(例如恒星)的光,穿过地球大气层时,大气层的湍流会导致这些恒星得图
    正/逆压电效应与材料本身的各向异性程度紧密相关,反过来又与压电材料的晶体结构存在关联,而各向异性的程度同时又受到极化过程的影响。下面,我们将介绍如何在 COMSOL 软件中正确地模拟压电材料的晶体取向和极化方向。 压电效应简介 让我们快速回顾一下压电效应的概念:正压电效应指材料受到机械力的作用时,其电极化会发生改变;而逆压电效应指对材料施加外部电场后,材料会发生变形。 压电效应源自晶体结构 在 3
    影响力
    粉丝
    内容
    获赞
    收藏
      项目客服
      培训客服
      0 1