This paper presents a two-dimensional implementation of the high-order penalized vortex in cell method applied to solve the flow past an airfoil with a vortex trapping cavity operating under moderate Reynolds number. The purpose of this article is to investigate the fundamentals of the vortex trapping cavity. The first part of the paper treats with the numerical implementation of the method and high-order schemes incorporated into the algorithm. Poisson, stream-velocity, advection, and diffusion equations were solved. The derivation, finite difference formulation, Lagrangian particle remeshing procedure, and accuracy tests were shown. Flow past complex geometries was possible through the penalization method. A procedure description for preparing geometry data was included. The entire methodology was tested with flow past impulsively started cylinder for three Reynolds numbers: 550, 3000, 9500. Drag coefficient, streamlines, and vorticity contours were checked against results obtained by other authors. Afterwards, simulations and experimental results are presented for a standard airfoil and those equipped with a trapping vortex cavity. Airfoil with an optimized cavity shape was tested under three angles of attack: 3°, 6°, 9°. The Reynolds number is equal to Re = 2 × 10⁴. Apart from performing flow analysis, drag and lift coefficients for different shapes were measured to assess the effect of vortex trapping cavity on aerodynamic performance. Flow patterns were compared against ultraviolet dye visualizations obtained from the water tunnel experiment.

In this paper, the motion of the n-vortex system as it collapses to a point in finite time is studied. The motion of vortices is described by the set of ordinary differential equations that we are able to solve analytically. The explicit formula for the solution demands the initial location of collapsing vortices. To find the collapsing locations of vortices, the algebraic, nonlinear system of equations was built. The solution of that algebraic system was obtained using Newton’s procedure. A good initial iterate needs to be provided to succeed in the application of Newton’s procedure. An unconstrained Leverber–Marquart optimization procedure was used to find such a good initial iterate. The numerical studies were conducted, and numerical evidence was presented that if in a collapsing system $\mathit{n}=50$ point vortices include a few vortices with much greater intensities than the others in the set, the vortices with weaker intensities organize themselves onto the vortex sheet. The collapsing locations depend on the value of the Hamiltonian. By changing the Hamiltonian values in a specific interval, the collapsing curves can be obtained. All points on the collapse curves with the same Hamiltonian value represent one collapsing system of vortices. To show the properties of vortex sheets created by vortices, the passive tracers were used. Advection of tracers by the velocity induced by vortices was calculated by solving the proper differential equations. The vortex sheets are an impermeable barrier to inward and outward fluxes of tracers. Arising vortex structures are able to transport the passive tracers. In this paper, several examples showing the diversity of collapsing structures with the vortex sheet are presented. The collapsing phenomenon of many vortices, their ability to self organize and the transportation of the passive tracers are novelties in the context of point vortex dynamics.