ECLIPPS: 1. Three-Dimensional CFD Prediction of the Ice Accretion

Recent Developments of the Navier Stokes Multi Block (NSMB) CFD solver

“…91 In NSMB, the supercooled water droplets field and impingement locations are obtained by means of an Eulerian approach similarly to FENSAP-ICE 92 and ONICE3D. 93 The ice thickness is calculated by resolving the Shallow Water Icing Model (SWIM) for three-dimensional and multi-block ice accretion simulations. 92 A Mesh deformation is used to track the ice/air interface through time for multi-steps calculations.…”

A single step ice accretion model using Level-Set method

“…A second flow calculation is performed by imposing a constant and arbitrary wall temperature T wall . The computed heat flux q wall obtained from the RANS solver allows the calculation of the heat transfer coefficient by the following formula (see Montreuil et al (2009) for more details about the calculation) :…”

“…In recent years, icing codes based on the solutions of the Reynolds-Averaged Navier-Stokes (RANS) equations have been developed such as LEWICE3D (NASA) (Bidwell and Potapczuk (1993)),FENSAP-ICE (McGill University) (Beaugendre et al (2006)), CANICE2D-NS (École Polytechnique Montréal) (Hasanzadeh et al (2013)), ONICE3D (ONERA) (Montreuil et al (2009)). Because of a seamless integration to Navier-Stokes solvers, an Eulerian formulation for droplet trajectories initially introduced by (Scott et al (1988)) have been implemented in several icing codes (Bourgault et al (1999); Jung and Myong (2013); Kim et al (2013); Montreuil et al (2009); Zhu et al (2012); Cao et al (2012); Jung et al (2011)). Benefits of an Eulerian approach are numerous, including the direct determination of the impingement efficiency without requiring a statistical averaging process, the computation of the droplet state vector in the entire domain without defining appropriate initial positions of the water droplets and the possibility to compute the droplet field on the same mesh used for the flow.…”

“…In FENSAP-ICE and ONICE3D the heat transfer coefficient is determined with the Navier-Stokes flow solver replacing the boundary layer method historically used in icing simulation. Mesh deformation methods have been employed to track the ice/air interface in an automatic multi-step icing calculations by Beaugendre et al (2003) and Montreuil et al (2009). Grid deformation techniques are usually used because of their easiness of implementation in an automatic process of icing simulation.…”

Development of a three-dimensional icing simulation code in the NSMB flow solver

“…It is evident that a full fluid solution provided by a Navier-Stokes solver in combination with turbulence model would results in more accurate icing simulations. Icing codes such LEWICE3D (NASA) [1],CANICE2D-NS [2] (École Polytechnique Montréal),ONERA3D (ONERA) [3] and FENSAP-ICE (McGill University) [4] were developed for that matter.…”