A newly developed two-dimensional ice accretion and antiicing simulation code, CANICE2D-NS, is presented. The method is used to predict iced airfoil shapes and performance degradation with a multistep approach. A multiblock Navier-Stokes code, NSMB, has been coupled with the CANICE2D icing framework, supplementing the existing panel method-based flow solver. Attention is paid to the roughness iniplementation within the turbulence model and to the convergence of the steady and quasi-steady iterative procedure. The new conpling allows fully automated multilayer icing simulation, whereas also permitting flow analysis and performance prediction of iced airfoils. Effects of uniform surface roughness in quasi-steady ice accretion simulation are analyzed through different validation test cases. The results demonstrates the benefits and robustness of the new framework in predicting ice shapes and aerodynamic performance parameters, as well as iced airfoil surface pressure coefficients. Finally, the convergence of the quasi-steady algorithm is verified and identifies the need for an order of magnitude increase in the number of multitime steps in icing simulations.
Nomenclaturethe nearest wall, m í/() = offset in the wall distance to account for wall roughness, m /(,., = turbulent convective heat transfer coefficient, W/m-•K K¡ = equivalent sand grain roughness height normalized hy chord k = von Kármán constant k^ = equivalent sand grain roughness height, m LWC = liquid water content, kg/m' MED = mean droplet diameter, m Pr, = turbulent Prandtl number 5? = Stanton number Stf. = roughness Stanton number S = transformed vorticity u = longitudinal velocity component, m/s Ug = boundary-layer edge velocity, m/s «J = friction velocity, m/s 1/ = kinematic viscosity, m^/s t/, = turbulent eddy viscosity, m-/s D = transported quantity in the S-A model w = wall value y = distance along the wall normal, m y+ = normalized wall distance •p = density, kg/m-