Implementation of rigid solid boundaries to the Vorticity-Stream function formulation of incompressible Navier-Stokes equations by time dilation

In the present paper, a new penalization method is proposed for implementation of the rigid surfaces on the Navier-Stokes equations in the vorticity-stream function formulation. In this method, a rigid body is considered as a region in the fluid flow, where the time is stopped. Therefore, by stopping the fluid particles, this region plays the role of a rigid body. In this regard, a new transformation is introduced and applied to the governing equations and a set of modified equations are obtained. Then, in the modified equations, the time dilation of the solid region is approached to infinity, while the time dilation of the fluid region remains In the article, the physical and mathematical properties of modified equations are investigated and satisfaction of the no-slip and no-penetration conditions are justified. Then, a suitable numerical algorithm is presented for solving the modified equations. In the proposed algorithm, the modified equation is time integrated via the Crank–Nicolson method, and the spatial discretization with the second-order finite differencing on a uniform Cartesian grid. The method is applied to the fluid flow around a square obstacle placed in a channel, the sudden flow perpendicular to a thin flat plate, and the flow around a circular cylinder. The results show that the no-slip and no-penetration conditions are satisfied accurately, while the flow fields are also high level of accuracy.