An improved characteristic based volume penalization method for the Euler equations towards icing applications. (30th May 2021)

• Record Type: Journal Article
• Authors: Lavoie, Pierre; Radenac, Emmanuel; Blanchard, Ghislain; Laurendeau, Éric; Villedieu, Philippe
• Is Part Of: Computers & fluids. Volume 222(2021)
• Publication Date: 2021-05-30
• Subjects: CFD -- Euler equations -- Volume penalization -- Compressible flow -- Inviscid flow -- Immersed boundary method -- Icing -- Ice accretion; Fluid dynamics -- Data processing -- Periodicals; 532.050285
• DOI: 10.1016/j.compfluid.2021.104917
• Languages: English
• ISSNs: 0045-7930

Abstract

Highlights: A new penalization method is developed specifically for the Euler equations. A curvature-based pressure extrapolation is used at the immersed boundary. Entropy and total enthalpy are conserved in the wall normal direction. The method is more accurate on high curvature geometries typical of ice accretion. The method is suitable for subsonic and transonic flows, and on general meshes. Abstract: Immersed boundary methods (IBMs) are an interesting alternative to the usual body-fitted mesh approach when dealing with complex geometries as they allow simpler mesh generation. The volume penalization method (an IBM) is commonly used for incompressible and compressible viscous flows but only one application to compressible inviscid flows can be found, which uses the characteristic-based volume penalization (CBVP) method. This approach penalizes the Euler equations to enforce a no-penetration velocity and an adiabatic wall while accounting for wall curvature. A new penalization method based on the CBVP 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 than the CBVP on coarser meshes and better at retrieving attached flows for curved geometries. The paper concludes that the proposed method is suitable for general aerospace applications and beneficial for icing simulations which can exhibit highly curved geometries.