The motion of a viscous fluid flow is described by the well-known Navier-Stokes equations. The Navier-Stokes equations contain the conservation laws of mass and momentum, and describe the spatial-temporal change of the fluid velocity field. This paper aims to investigate numerical solvers for the incompressible Navier-Stokes equations in two and three space dimensions. In particular, we focus on the second-order projection method introduced by Kim and Moin, which was extended from Chorin’s first-order projection method. We impose periodic boundary conditions and apply Fourier-Spectral methods for the periodic boundary condition. Numerically, we discretize the system using centered differences scheme on Marker and Cell (MAC) grid spatially and the Crank-Nicolson scheme temporally. We then apply the fast Fourier transform to solve the resulting Poisson equations as sub-steps in the projection method. We will verify numerical accuracy and perform von Neumann stability analysis. In addition, we will simulate the particles’ motion in the 2D and 3D fluid flow.
Received: August 27, 2018.
Published online: November 10, 2018.
Vorleak YekDepartment of Mathematics and Statistics, California State University, Long Beach, 1250 Bellflower Blvd, Long Beach, CA 90840, USA.