A high-order accuracy method for calculating the initial icing stage of civil aircraft structural elements

An effective approach based on the Discontinuous Galerkin Method (DG) of a high-order accuracy for calculating the initial stage of aircraft wing icing is presented. The problem is solved in the Euler approximation for small-sized water droplets that do not affect the main flow. Navier-Stokes equation system, Euler model equation system for the Liquid Water Content and some relations of the ice growth thermodynamics equations are written. Initial and boundary conditions are formulated. A supercomputer DG implementation is proposed to solve these equation systems. The efficiency of the parallel version for the code is investigated. Comments are given on the peculiarities of the calculation procedure organization. The accuracy of the calculation using DG schemes of different accuracy order is investigated. Test cases on the finely dispersed flow of supercooled droplets around the cylinder and the NACA0012 profile are presented. A comparison of numerical and experimental data is performed. The conclusion is made about the possibility of applying the developed methodology in practice.