Governing Equations
Governing equations are expressed by conservation of mass and momentum. The equations are:
$$
\eqalign{ \frac{D\rho}{Dt}+\rho \nabla \cdot \mathbf{u} =0 }
$$
and
$$
\eqalign{ \frac{D \mathbf{u} }{Dt}=-\frac{1}{\rho}\nabla P +\nu \nabla^2 \mathbf{u} +\mathbf{f} }
$$
where $\rho$ indicates density, $P$ is the pressure, $\mathbf{u}$ the velocity and $\mathbf{f}$ the external force. It should be noted that the first equation is written in the form of compressible flow. In the MPS method, incompressibility is enforced by way of setting $\frac{ D \rho}{D t}=0$ at each particle at each calculation time step.