Solid-liquid phase change algorithm with Euler-Lagrange iteration
During the solid-liquid phase change process,the external force can cause relative motion of the solid phase change material within the liquid phase change material,which can seriously affect the flow and heat transfer.An Euler-Lagrange iterative solid-liquid phase change algorithm is proposed to numerically solve the above problem.The Lagrange iteration for predicting the solid relative motion is externally coupled to the Euler iteration for calculating the phase change flow and heat transfer,which can stably and accurately simulate the flow and heat transfer of the solid-liquid phase change and the relative motion of the solid phase change material.Based on the present algorithm,the close-contact melting processes in the square cavity are investigated under gravity.Furthermore,the influence of different mushy zone coefficients and gravity accelerations on this algorithm is explored.The results show that the average error of the algorithm in predicting the liquid phase volume fraction is 4.93%,and the numerical oscillation of the solid relative motion prediction is reduced by 51.42%.For paraffin materials,the mushy zone coefficient of 1010 is recommended for this algorithm.The increase of the gravity acceleration can improve the melting rate and accelerate the relative downward motion,yet has less impact on the overall trends of liquid phase fraction and solid phase relative motion.The results can be used as the theoretical guideline for the design and optimization of the efficient solid-liquid phase change energy storage devices.