[论文解读] Three-Dimensional Smoothed Particle Hydrodynamics Simulation for Liquid Metal Solidification Process

Raden Ahnaf Faqih S, Christian Fredy Naa|arXiv (Cornell University)|Sep 17, 2013

Fluid Dynamics Simulations and Interactions参考文献 6被引用 4

一句话总结

本研究采用基于能量法的热传导模型与非等温相变公式，对快速冷却条件下液态金属的三维光滑粒子流体动力学（SPH）模拟进行了研究。该方法成功捕捉了温度演化、相变过程及缩孔的形成，但因液-固相互作用建模不足，裂纹预测能力仍有限。

ABSTRACT

The solidification behavior of liquid metal in a container under rapid cooling process is one of the major concerns to be analyzed. In order to analyze its fundamental behavior, a three- dimensional (3D) fluid dynamics simulation was developed using a particle-based method, known as the smoothed particle hydrodynamics (SPH) method. Governing equations that determine the fluid motion and heat transfer involving phase change process are solved by discretizing their gradient and Laplacian term with the moving particles and calculating the interaction with its neighboring particles. The results demonstrate that the SPH mehod can successfully reproduce the behavior and defect prediction of liquid metal solidification process.

研究动机与目标

• 采用无网格粒子方法模拟快速冷却下液态金属的三维凝固行为。
• 开发并验证一种基于SPH的热传导模型，考虑潜热释放与非等温相变，以准确模拟温度演化。
• 研究凝固过程中缺陷的形成，特别是缩孔与潜在裂纹。
• 评估液-固相互作用模型在三维SPH模拟中的性能。
• 识别当前建模在预测快速凝固中复杂缺陷（如裂纹）方面的局限性。

提出的方法

• 通过核加权粒子相互作用近似梯度与拉普拉斯算子，对纳维-斯托克斯方程与能量方程进行SPH离散化。
• 采用立方样条核函数与平滑长度 $ h $，在 $ 2h $ 半径范围内计算粒子间相互作用。
• 应用焓法建模热传导过程中的潜热释放，实现非等温凝固模拟。
• 采用与温度相关的粘度模型，模拟冷却与相变过程中的流体行为。
• 施加边界力（重力与排斥力）以模拟容器壁面与热量提取过程。
• 将固相与糊状相视为高粘度流体，实现无网格依赖的凝固动力学模拟。

实验结果

研究问题

• RQ1 SPH方法能否准确模拟快速冷却下三维液态金属凝固过程中的温度分布与相变？
• RQ2 基于焓法的热传导模型在捕捉潜热释放与非等温相变方面表现如何？
• RQ3 凝固过程中会形成哪些缺陷（如缩孔或裂纹），其形成受粒子密度与相间相互作用的影响如何？
• RQ4 尽管处于快速冷却条件，为何当前模拟中未形成裂纹？液-固相互作用模型中缺失了什么？
• RQ5 不同层的密度演化如何影响缺陷形成，特别是在最终凝固阶段？

主要发现

• SPH模型成功再现了温度分布与相变过程：液相温度高于 589.7 °C（红色），糊状相温度为 536.1 °C（蓝色），固相温度低于 536.1 °C（银色）。
• 所有液相粒子在 2.7272 秒内完全凝固，归因于高热导率与初始金属（600 °C）和容器壁（200 °C）之间显著的温度梯度。
• 由于固相密度高于液相，导致体积收缩，最终阶段观察到缩孔的形成。
• 温度历史显示糊状相与固相均经历快速冷却，相变后温度出现急剧下降。
• 尽管冷却速率快，但因糊状相与固相界面力模型不足，未形成裂纹。
• 模拟全过程保持稳定，表明当前SPH格式与边界处理具有良好的鲁棒性。