A numerical scheme is developed to simulate a flow around airfoils by using vortex methods. For this scheme, a numerical algorithm is constructed and exact analytical expressions are obtained for the...

The possibilities of applying the pure Lagrangian vortex methods of computational fluid dynamics to viscous incompressible flow simulations are considered in relation to various problem formulations. The modification of vortex methods—the Viscous Vortex Domain method—is used which is implemented in the VM2D code developed by the authors. Problems of flow simulation around airfoils with different shapes at various Reynolds numbers are considered: the Blasius problem, the flow around circular cylinders at different Reynolds numbers, the flow around a wing airfoil at the Reynolds numbers ${10}^4$ and ${10}^5$, the flow around two closely spaced circular cylinders and the flow around rectangular airfoils with a different chord to the thickness ratio. In addition, the problem of the internal flow modeling in the channel with a backward-facing step is considered. To store the results of the calculations, the POD technique is used, which, in addition, allows one to investigate the structure of the flow and obtain some additional information about the properties of flow regimes.

A high-accuracy numerical scheme is proposed for vortex methods of flow simulation around airfoils of arbitrary shape including airfoils with sharp edges, because it does not require the solution...

This article describes the open-source C++ code VM2D for the simulation of two-dimensional viscous incompressible flows and solving fluid-structure interaction problems. The code is based on the Viscous Vortex Domains (VVD) method developed by Prof. G. Ya. Dynnikova, where the viscosity influence is taken into account by introducing the diffusive velocity. The original VVD method was supplemented by the author’s algorithms for boundary condition satisfaction, which made it possible to increase the accuracy of flow simulation near the airfoil’s surface line and reduce oscillations when calculating hydrodynamic loads. This paper is aimed primarily at assessing the efficiency of the parallelization of the algorithm. OpenMP, MPI, and Nvidia CUDA parallel programming technologies are used in VM2D, which allow performing simulations on computer systems of various architectures, including those equipped with graphics accelerators. Since the VVD method belongs to the particle methods, the efficiency of parallelization with the usage of graphics accelerators turns out to be quite high. It is shown that in a real simulation, one graphics card can replace about 80 nodes, each of which is equipped with 28 CPU cores. The source code of VM2D is available on GitHub under GNU GPL license.

When simulating incompressible flow around airfoil using vortex methods of computational hydrodynamics, it is necessary to solve the boundary integral equation with respect to the intensity of the vortex sheet at the airfoil surface line. For its approximate solution, the Galerkin method with a piecewise-constant and piecewise-linear representation of the solution on the airfoil surface line can be applied. The coefficients of the system of linear equations, arising from the procedure of discretization of the integral equation, e xpress the influence of the vortex sheet located on the panels of the airfoil surface line on other panels; the components of the right-hand side vector express the influence of point vortices that simulate a vortex wake in the flow domain. For the case of the airfoil surface line approximation with a polyline consisting of rectilinear panels, exact analytical expressions are obtained for the coefficients of the matrix and the right-hand side vector; for the case of taking into account the curvature of the panels, an approach for their approximate calculation is proposed and all the necessary formulae are derived.