2023 Volume 79 Issue 15 Article ID: 22-15022
A multiphase computation method based on the phase-averaging method was applied to water-ice phasechange problems accompanied by the natural convections including density inversion regions in water. The computation method for the Stefan condition was improved in this study to calculate the freezing thickness considering the position of the water-ice interface in the computational cells. After the basic verification of the proposed method using the one-dimensional freezing problem, the experimental results of the unsteady freezing problems and those of the unsteady melting problems in rectangular cavities were calculated with the present method. As a result, it was shown that the changes in the calculated water-ice interfaces are in good agreement with the experimental results and that the natural convections with density inversion are reasonably calculated.
