HydrostaticPressureAlgorithm
The hydrostatic pressure algorithm solves the poisson equation derived in
DerivePoissonEquationForPressure .
The numerical representation of the equation is
\begin{align} -\frac{1}{\Delta t^{2}}\frac{1}{\rho }\frac{\partial \rho }{\partial p}\cdot \tilde{p}_{hyd,i}^{n+1}+\sum\limits_{j=1}^{N}{c_{ij}^{\nabla ^{T}\left( \frac{1}{\rho }\nabla \right)}}\cdot \tilde{p}_{hyd,j}^{n+1}=-\frac{1}{\Delta t^{2}}\frac{1}{\rho }\frac{\partial \rho }{\partial p}p_{hyd,i}^{n}+\nabla ^{T}\left( \frac{1}{\rho }\nabla \mathbf{S}_{s}^{n} \right)_{i}+\nabla ^{T}\mathbf{g}_{i}\end{align}