ONERA Three-Dimensional Icing Model

A three-dimensional icing model has been developed at ONERA to calculate ice accretion shapes for aerodynamic components that can not be predicted using conventional two-dimensional codes. It is described, emphasizing the original parts with respect to the two-dimensional existing models. The model includes Euler inviscid flow calculation. Droplet trajectories are calculated in a three-dimensional grid. The remesh on the leading edge is adapted to follow aerodynamics singularities. The boundary layer is calculated using a mixing length formulation to model the wall roughness influence on convective heat transfer. Runback paths are integrated. The heat balance is calculated in a grid created along the runback paths. The domain of validity of the three-dimensional icing code is described; compared with the two-dimensional model this domain is wider, especially for high speeds. The three-dimensional model is shown to simulate well a uniform ice deposit on a three-dimensional rotor blade tip. Then, a comparison of the three- and two-dimensional codes on an infinite swept wing shows that the corrected two-dimensional code predicts the catch efficiency but not the ice shape. Finally, it is shown that the continuum flux hypothesis prevents the three-dimensional model from simulating correctly the "lobster tail" ice shape (nonuniform ice deposit).