[论文解读] The pear-shaped fate of an ice melting front

James N. Hewett, Mathieu Sellier|arXiv (Cornell University)|May 7, 2017

Icing and De-icing Technologies参考文献 9被引用 3

一句话总结

本研究通过动态网格技术结合节点重排算法，数值模拟了接近水密度反常点（4.0 °C）的加热水平圆柱体周围冰的熔化过程，以追踪移动的固液界面。结果表明，由于热羽流的影响，熔化前沿呈现梨形：当壁温高于4.0 °C时，热羽流向上；当壁温低于4.0 °C时，热羽流向下。这证实自然对流显著影响熔化速率和界面形态，且与实验数据在形状和整体熔化进程上具有良好的一致性。

ABSTRACT

A fluid-structure interaction problem with the melting of water around a heated horizontal circular cylinder is analysed with numerical simulations. Dynamic meshing was used for evolving the flow domain in time as the melting front extended radially outward from the cylinder; a node shuffle algorithm was used to retain mesh quality across the significant mesh deformation. We simulated one case above the density inversion point of water and one case below, yielding pear-shaped melting fronts due to thermal plumes either rising or falling from the cylinder, respectively. Results were compared with previous experimental studies and the melting front profiles matched reasonably well and melting rates were in agreement. We confirm that natural convection plays a significant role in the transport of energy as the melt zone increases, and needs to be considered for accurately modelling phase change under these conditions.

研究动机与目标

• 研究水在约4.0 °C附近的密度反常对加热水平圆柱体周围相变过程中熔化前沿形态与动力学的影响。
• 验证在自然对流效应下，通过显式界面追踪的斯蒂芬问题数值模拟的可靠性。
• 将模拟得到的熔化前沿形状与速率与实验数据进行对比，涵盖反常点以上和以下温度两种情况。
• 评估自然对流在熔融区域扩展过程中对增强传热的作用。
• 评估采用节点重排的动态网格方法在处理移动边界问题中大变形时的准确性。

提出的方法

• 采用二维、不可压缩、层流模型，并应用Boussinesq近似以考虑浮力驱动的对流。
• 通过动态网格技术显式追踪移动边界（熔化前沿），并利用节点重排算法保持网格质量。
• 在每个时间步长，基于更新后的网格在准稳态假设下求解流场。
• 水的热物理性质除密度外均假设为常数，密度随温度线性变化以捕捉密度反常效应。
• 界面运动由局部热通量和斯蒂芬条件确定，未引入人工相变源项。
• 模拟针对两种壁温进行：2.3 °C（低于反常点）和14.1 °C（高于反常点），与实验条件一致。

实验结果

研究问题

• RQ1 水在约4.0 °C附近的密度反常如何影响加热水平圆柱体周围熔化前沿的形状与演化？
• RQ2 随着熔融区域扩大，自然对流在多大程度上主导传热过程？其影响在反常点以上与以下有何差异？
• RQ3 采用节点重排的动态网格能否准确捕捉实验中观察到的复杂梨形界面形态？
• RQ4 为何模拟在从导热主导到对流主导的过渡阶段（τ ≈ 0.3）与实验相比出现熔化速率偏差？
• RQ5 二维稳态假设如何影响对三维实验中观测到的瞬态流动不稳定性与增强传热的预测？

主要发现

• 两种情况下均因热羽流形成梨形熔化前沿：T_w = 14.1 °C时向上，T_w = 2.3 °C时向下，与实验观测一致。
• 对应的瑞利数分别为6700（倒置）和7400（正置），表明尽管流动结构不同，对流强度相近。
• 随着熔融体积增加，自然对流成为主导传热机制，显著影响界面演化。
• 在整个模拟周期内，界面形状与实验剖面吻合良好，证实了模型的有效性。
• 在过渡阶段（τ ≈ 0.3）熔化速率的偏差归因于二维稳态模型中缺失的非定常三维效应与瞬态不稳定性。
• 正置梨形界面表现出更复杂的流动特征，包括圆柱体上方的反向旋转涡旋，这归因于密度反常点的引入。