Numerical method used for the simulation of fluid-solid transitions based on the VOF and SAM models and prediction of snow distributions on large-span roofs

During a snowfall, the surfaces of the snow beds on roofs change dynamically with time and wind conditions. The generation of body-fitted computational fluid dynamics grids for simulating dynamic snow surfaces using the moving-grid method remains challenging. Therefore, a numerical method was developed to simulate the fluid-solid transition based on the volume of the fluid model and solidification and melting models. The snowfall process was divided into n stable steps, and a stable state calculation method was adopted to calculate the air and snow phases. The dynamic snow boundary was simulated at each stage by converting the fluid or solid characteristics of the domain in which the grids resided. Based on the results of the (n-1)th stage, the fluid domain grids buried by the snow were filled with liquid and frozen into the solid domain. Instead, the fluid in the frozen grids, where snow was about to be eroded, was evacuated and converted into fluid-domain grids. Thus, a complex and changeable snow boundary could be adapted without re-meshing. A modified wall friction velocity equation was derived to correct the wall friction velocity at the snow boundary owing to the serrated boundary of the solidified grids, thereby improving calculation accuracy. The numerical simulation results were verified via field observations of snow distribution on the stepped flat and gable roof models. The snow variation on a fullscale, three-span, suspended cable roof under snow-fall conditions was predicted using the proposed method.