Improved Algorithm of Boundary Integral Equation Approximation in 2D Vortex Method for Flow Simulation Around Curvilinear Airfoil

The problem of the accuracy improving is considered for vortex method. The general Galerkin-type approach is considered for the numerical solution of the boundary integral equation. It is shown that the airfoil surface line representation as a polygon consisting of rectilinear panels can lead to incorrect behavior of the numerical solution of the boundary integral equation with respect to unknown vortex sheet intensity, especially in case of considerably different lengths of the neighboring panels. However, in this case, the exact analytical formulae are derived for the coefficients of the linear algebraic system, which approximate the integral equation. The approach is developed which makes it possible to take into account explicitly the curvature of the airfoil surface line. The approximate analytical formulae are derived for the coefficients of the linear system, which are expressed through the integrals of the fractions, which denominators can be equal to zero. Numerical example is considered for the elliptic airfoil in the incident flow at its essentially non-uniform discretization. The suggested algorithm has acceptable numerical complexity, however it permits to obtain rather good numerical results.