It is common for an aircraft to encounter icing weather conditions, which would be dangerous to the flight. Thus, there is a need to study the detail of icing effect and the process of ice accretion on the aircraft. In this paper, considering three different icing models according to weather conditions, i.e., sharp-angled ice, blunt-nosed ice and double horn ice, the Reynolds-averaged N-S equations and the S-A turbulence model are used to analyze the flow field for an iced wing/body configuration with a multi-block strategy and structured grid technique. The numerical result is compared with the experimental data. A flow solver is developed based on the Euler equations to investigate the ice accretion process. The droplets are tracked by using the Lagrangian method. In addition, a revised Messinger model is proposed to simulate the ice accretion. This numerical simulation is conducted for the ice accretion on an M6 wing and a wing/body/tail configuration. The presented results preliminarily show that the numerical methods are feasible and effective.