Three-Dimensional Numerical Simulation of Ice Crystal Melting in Jet Engine

A method was presented for calculating the melting rate of ice crystals in different axial positions in a compressor. The volume of fluid model and enthalpy-porosity model were used to simulate the melting process of single ice crystal. The result was validated by experimental data from literature. The change of the liquid fraction of the ice crystals with the axial position was obtained when the temperature distribution inside the compressor was introduced into the melting process. The liquid fraction of the spherical ice crystals was fitted with the ice crystal diameter as a variable, and the general distribution function for computing liquid fraction in any axial position in the compressor was obtained. The function of axial position as a variable for non-spherical ice crystal was also got.