The modification of the shallow water icing model to handle de-icing phenomenon is the main focus of this study. As stated in the original model [1], the runback water is modeled utilizing a lubrication assumption for the water film velocity profile. A constant film temperature Tf(t,x) is then calculated under the thin-film hypothesis. Unlike the simplified icing model, the temperature field within the ice layer Tice(t,x,z) is no longer assumed to be constant. Instead, a temperature profile is utilized, enabling the generation of a static film on the wall when a heat conduction source term from a resistance is present [2]. A Temperature profile Ts(t,x,z) is also used in the static film layer if the model predicts the occurrence of this state. In the energy equation for both the solid ice and liquid portion of the static water film, transverse transfers are not considered, a 1D heat equation is then resolved. An integral approach and proper boundary conditions are used to close the problem. The validity of the integral method deteriorates as the thickness over which vertical integration is performed increases. To avoid this problem, a multi-layer approach is proposed. The thickness of the ice block is then divided into three layers of identical thickness. The purpose of this study is to offer a straightforward and robust method suitable for conducting industrial test cases. The model will first be introduced, followed by a description of the numerical approach. Subsequently, validation test cases will be conducted. Realistic de-icing scenarios will then be designed to evaluate the model [3]. Additionally, non uniform roughness effects will be examined.
Published on 23/10/24
Submitted on 23/10/24
Volume Advances in Numerical Methods for Solution Of PDEs, 2024
DOI: 10.23967/eccomas.2024.031
Licence: CC BY-NC-SA license
Are you one of the authors of this document?