<div>Hi,</div>
<div>&nbsp;</div>
<div>I have been using a few different Krylov&nbsp;linear solver package such as nspcg,sparsekit and now petsc to solve the linear eqns for my NS solver momentum and poisson eqn. For moving bodies simulation, while solving the poisson eqn, sometimes the result suddenly changes drastically, although the iterations seems to have converged&nbsp;
e.g. at the previous time step, everything is fine, but at the next time step, the velocity suddenly explodes from 2-3 to 20-30+. </div>
<div>&nbsp;</div>
<div>The strange thing is that this deviation happens at different time step for different linear solver package ie for nspcg, it&nbsp;happens at time=20 but&nbsp;for petsc, it happens&nbsp;at time=40. It seems &quot;random&quot; in this sense. Changing the type of solver such as gmres or bcgs does not change the time of deviation. 
</div>
<div>&nbsp;</div>
<div>Does anyone know why this is happening? Is it a characteristic of krylov solver to happen once in a while? It does not happen in the stationary body simulation though.</div>