Lava flows represent complex thermofluid phenomena in which surface cooling leads to the formation of a solidified surface layer. Understanding the influence of such a surface layer on fluid flow is an important issue in lava flow modeling. It also shares essential characteristics with a wide range of engineering problems involving surface solidification. However, the role of plastic surface skin in controlling flow deceleration and stopping behavior has not been sufficiently clarified in existing models. In this study, two-dimensional smoothed particle hydrodynamics (SPH) simulations were conducted to investigate the influence of surface skin formation on lava flow dynamics. The temperature dependence of viscosity was introduced to reproduce a plastic surface skin. The skin was represented as a low-temperature, high-viscosity region. Comparisons with simulations without surface skin formation demonstrated that the surface skin exhibits a suppressive effect on the flow. This behavior was consistent with qualitative observations of flowing lava. It was also found that this surface skin caused the successive deceleration characteristic in Bingham fluids. As a result, both the flow velocity and the flowing distance are affected. These results suggest that accurate lava flow simulations require models that incorporate both surface skin effects and non-Newtonian behavior.
TOMITA et al. (2026) studied this question.