Numerical investigation of the thermal state of a layer of anisotropically scattering ice

Using the methods of mathematical modeling in a single-phase formulation of the Stefan problem, the problem of melting of an ice layer, which scatters radiation, was formulated and solved. To solve the radiation part of the problem, a modified method of mean fluxes was used. Anisotropic scattering was taken into account using the method of expansion of the scattering indicatrix in terms of Legendre polynomials. The influence of the scattering direction on the rate of ice layer melting is shown. Satisfactory agreement of the calculation results with the experimental data and with the data of calculation by the transport method has been obtained.