Diagonalized alternating-direction implicit (DADI) method is implemented in the Eulerian hyperbolic droplet solver, ICEPAC, for efficient high-order accurate analysis of aircraft icing. Detailed techniques for implementing the DADI method considering hyperbolicity characteristics are discussed. For the Eulerian droplet equation system to be strictly hyperbolic, additional source terms regarding artificial droplet pressure are included. Validations of the present implicit solver are conducted using two- and three-dimensional steady benchmark tests: NACA0012 airfoil, NACA23012 airfoil, and a swept wing. Also, the oscillating airfoil SC2110 case was analyzed to verify the robustness and efficiency of the proposed solver. In addition, the computational cost of the current implicit solver is considerably lower than that of the explicit multi-stage solver.
