An Improved Characteristic Based Volume Penalization Method for the Euler Equations Towards Icing Applications

Pierre Lavoie (1, 2) , Emmanuel Radenac (2) , Ghislain Blanchard (2) , Éric Laurendeau (1) , Philippe Villedieu (2)

1 EPM - École Polytechnique de Montréal
2 DMPE, ONERA, Université de Toulouse [Toulouse]

The numerical prediction of in-flight ice accretion generally involves geometry updates and re-meshing as the ice builds up. However, the generation of body-fitted meshes around complex ice shapes is not trivial and needs to be repeated several times to obtain the final ice shape. The use of an immersed boundary method can simplify the mesh generation and help in the automation of ice accretion simulations. In this paper, the development of an immersed boundary method for Euler flows is detailed. A new penalization method is proposed to impose the conservation of entropy and total enthalpy in the normal direction to the wall instead of the classical adiabatic condition. The two approaches are compared and numerically tested on several cases: weakly compressible flow around a circular cylinder, subsonic flow around a NACA0012 airfoil and flow around a challenging high curvature ice horn. The new method is found to be more accurate on coarser meshes and better at retrieving attached flows for curved geometries. The paper concludes that the proposed method is suitable for icing simulations while being still appropriate for general aerospace applications.

Mots clés: CFD; Euler Equations; Volume Penalization; Compressible Flow; Inviscid Flow; Immersed Boundary Method; Icing; Ice Accretion

Domaines: Sciences de l'ingénieur [physics]; Physique [physics]

Pierre Lavoie, Emmanuel Radenac, Ghislain Blanchard, Éric Laurendeau, Philippe Villedieu. An Improved Characteristic Based Volume Penalization Method for the Euler Equations Towards Icing Applications. Computers & Fluids, 2021, 222, pp.104917. ⟨10.1016/j.compfluid.2021.104917⟩. ⟨hal-02872954v2⟩