Hydrodynamics on a moving Voronoi mesh - Obstacle modeling

Hydrodynamics on a moving Voronoi mesh. Obstacles can be modeled conveniently by placing non-moving mesh generating points along both sides of the surface of the obstacle. Furthermore, reflective boundary conditions are applied to the edges composing the surface. The code is publicly available: