Solidification using smoothed particle hydrodynamics
Solidification using smoothed particle hydrodynamics
复制标题
DOI：
10.1016/j.jcp.2004.11.039
复制
发表时间：
2005-07-01
影响因子：
4.1
通讯作者：
Worster, MG
中科院分区：
文献类型：
--
作者：
Monaghan, JJ;Huppert, HE;Worster, MG
本文介绍了光滑粒子流体动力学（SPH）数值方法如何用于模拟单组分和双组分（二元合金）体系的冻结。首先研究了Ora纯液体的冻结，并将我们的计算结果与一维和二维体系的精确结果进行了比较。包括存在点源和系统边界不规则的情况。与理论的一致性，如果它是可用的，是非常令人满意的。然后，我们考虑了一个两组分系统，它最初是完全液体和模拟它的一组液体SPH粒子与一组虚拟的固体SPH粒子，最初没有质量，在热演化和固化的系统，质量从液体SPH粒子转移到冰（固体）SPH粒子。对于二元熔体，随着固体的体积分数增加，组成或液体中富含不形成固相的合金组分。在盐水的情况下，这一组分是盐。我们发现用SPH计算的温度和液体成分的变化与以前的糊状层理论非常一致，并给出了与实验类似的一致性。在这个初步的研究中，我们简化了计算，假设固体颗粒保持在它们形成的位置：一个很好的近似的情况下，溶液从上面冷却形成冰留下一个相对较轻的残留液体，或从下面冷却留下一个相对较重的液体。我们将我们的结果与硝酸钠水溶液[Nature 314（1985）703]下面的冻结前沿实验进行比较，发现一致性非常令人满意。
We show how the numerical particle method smoothed particle hydrodynamics (SPH) can be used to the simulate the freezing of one and two-component (binary alloy) systems, We first study the freezing Ora pure a liquid, and compare our computations against exact results for one and two dimensional systems. including cases where there are point sources and the boundary of the system is irregular. The agreement with theory, where it is available, is very satisfactory. We then consider a two-component system which is initially entirely liquid and model it by it set of liquid SPH particles together with a set of virtual solid SPH particles which initially have no mass, During the thermal evolution and solidification of the system, mass is transferred from the liquid SPH particles to the ice (solid) SPH particles. For a binary melt, as the volume fraction of the solid increases the composition or the liquid is enriched in the component of the alloy that does not form the solid phase. In the case of salty water, this component is the salt, We find that the variation of temperature and liquid composition calculated with SPH is in close agreement with previous theories of mushy layers and gives similar agreement with experiment. In this initial study we simplify the calculations by assuming the solid particles remain in the position where they are formed: a good approximation for the case where the solution is cooled from above to form ice leaving behind a relatively light residual liquid, or cooled from below leaving behind a relatively heavy liquid. We compare our results with experiments on the freezing front below of an aqueous sodium nitrate solution [Nature 314 (1985) 703], and find that the agreement is very satisfactory, (c) 2005 Published by Elsevier Inc.