Vortex methods for incompressible flow simulations on the GPU

We present a remeshed vortex particle method for incompressible flow simulations on GPUs. The particles are convected in a Lagrangian frame and are periodically reinitialized on a regular grid. The grid is used in addition to solve for the velocity-vorticity Poisson equation and for the computation of the diffusion operators. In the present GPU implementation of particle methods, the remeshing and the solution of the Poisson equation rely on fast and efficient mesh-particle interpolations. We demonstrate that particle remeshing introduces minimal artificial dissipation, enables a faster computation of differential operators on particles over grid-free techniques and can be efficiently implemented on GPUs. The results demonstrate that, contrary to common practice in particle simulations, it is necessary to remesh the (vortex) particle locations in order to solve accurately the equations they discretize, without compromising the speed of the method. The present method leads to simulations of incompressible vortical flows on GPUs with unprecedented accuracy and efficienc

cuIBM -- A GPU-accelerated Immersed Boundary Method

A projection-based immersed boundary method is dominated by sparse linear algebra routines. Using the open-source Cusp library, we observe a speedup (with respect to a single CPU core) which reflects the constraints of a bandwidth-dominated problem on the GPU. Nevertheless, GPUs offer the capacity to solve large problems on commodity hardware. This work includes validation and a convergence study of the GPU-accelerated IBM, and various optimizations.Comment: Extended paper post-conference, presented at the 23rd International Conference on Parallel Computational Fluid Dynamics (http://www.parcfd.org), ParCFD 2011, Barcelona (unpublished

Density functional theory modeling of vortex shedding in superfluid He-4

Formation of vortex rings around moving spherical objects in superfluid He-4 at 0 K is modeled by time-dependent density functional theory. The simulations provide detailed information of the microscopic events that lead to vortex ring emission through characteristic observables such as liquid current circulation, drag force, and hydrodynamic mass. A series of simulations were performed to determine velocity thresholds for the onset of dissipation as a function of the sphere radius up to 1.8 nm and at external pressures of zero and 1 bar. The threshold was observed to decrease with the sphere radius and increase with pressure thus showing that the onset of dissipation does not involve roton emission events (Landau critical velocity), but rather vortex emission (Feynman critical velocity), which is also confirmed by the observed periodic response of the hydrodynamic observables as well as visualization of the liquid current circulation. An empirical model, which considers the ratio between the boundary layer kinetic and vortex ring formation energies, is presented for extrapolating the current results to larger length scales. The calculated critical velocity value at zero pressure for a sphere that mimics an electron bubble is in good agreement with the previous experimental observations at low temperatures. The stability of the system against symmetry breaking was linked to its ability to excite quantized Kelvin waves around the vortex rings during the vortex shedding process. At high vortex ring emission rates, the downstream dynamics showed complex vortex ring fission and reconnection events that appear similar to those seen in previous Gross-Pitaevskii theory-based calculations, and which mark the onset of turbulent behavior.Comment: 23 pages, 7 figure

Comparison of the finite volume method with lagrangian vortex method for 2d flow simulation around airfoils at intermediate reynolds number

The paper is dedicated to the numerical simulation of two-dimensional viscous incompressible flow. The Viscous Vortex Domains method is considered, which is one of the modifications of Lagrangian vortex methods. This method is implemented in the open-source VM2D code, which developes by the authors. Model problems of exter- nal and internal flow simulation are considered for VM2D testing. For the problem of the flow simulation around two closely spaced circular cylinders, VM2D is compared to OpenFOAM in terms of computational efficiency. As the internal flow problem the flow in a channel with a backward-facing step is considered. For both problems, the results obtained in VM2D are in good agreement with the results of other researchers