2023 Volume 15 Issue 1 Pages 27-35
This study proposes a numerical method to treat water freezing and ice melting simultaneously in order to apply the previously proposed numerical method for water-ice phase change on an orthogonal structured grids to a wider range of temperature conditions assuming realistic phenomena. In the proposed method, fluid flow, heat conduction in solids, and solid-liquid phase change are all calculated on orthogonal structured grids. In the calculation of the solid-liquid phase change, the approximate location of the interface in each calculation cell is taken into account, and the Stefan condition is calculated under the assumption that the interface temperature is always equal to the freezing temperature. After a basic check of the reproducibility of natural convection including density inversion region, the proposed method is applied to a problem in which freezing and melting occur simultaneously. From the calculation results, it was confirmed that the natural convection of water including the density inversion region has a significant effect on the freezing of water and melting of ice, and that the temperature of the water-ice interface agrees well with the freezing temperature of water.
