<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Oct 23, 2017 at 10:04 AM, Hao Zhang <span dir="ltr"><<a href="mailto:hbcbh1999@gmail.com" target="_blank">hbcbh1999@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">The big picture is I'm solving 3d incompressible Navier-Stokes Equations using staggered/MAC grid with Finite Difference Method. The particular function is poisson pressure solver or Laplacian Mark Adam mentioned. The simulation runs fine for medium size mesh grid. When I try harder to go very fine grid, not DNS level, I'm having some difficulties to have meaning physical results. </div></blockquote><div><br></div><div>I assume that you are getting good solutions when the solver is converging well.<br></div><div><br></div><div>I would not trust that your solution is ever accurate (converged) when you see this (stagnation), unless you get the true residual way, way down. How low do you get the true residual?</div><div><br></div><div><div style="font-size:12.8px">  0 KSP preconditioned resid norm 1.630377897956e+01 true resid norm 7.365742695123e+00 ||r(i)||/||b|| 1.000000000000e+00</div><div style="font-size:12.8px">  1 KSP preconditioned resid norm 3.815819909205e-02 true resid norm 7.592300567353e-01 ||r(i)||/||b|| 1.030758320187e-01</div><div style="font-size:12.8px">  2 KSP preconditioned resid norm 4.312671277701e-04 true resid norm 2.553060965521e-02 ||r(i)||/||b|| 3.466128360975e-03</div><div style="font-size:12.8px">  3 KSP preconditioned resid norm 3.011875330569e-04 true resid norm 6.836208627386e-04 ||r(i)||/||b|| 9.281085303065e-05</div><div style="font-size:12.8px">  4 KSP preconditioned resid norm 3.011783821295e-04 true resid norm 2.504661140204e-04 ||r(i)||/||b|| 3.400418998972e-05</div><div style="font-size:12.8px">  5 KSP preconditioned resid norm 3.011783818372e-04 true resid norm 2.504053673004e-04 ||r(i)||/||b|| 3.399594279422e-05</div><div style="font-size:12.8px">  6 KSP preconditioned resid norm 3.011783818375e-04 true resid norm 2.503984078890e-04 ||r(i)||/||b|| 3.399499795925e-05</div></div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">I would take the first test case (smallest) in which you can see "nonphysical results", use GMRES and ILU with lots of iterations, or a direct solver ideally, and get a good solution. If you can not get a good solution then do not worry about multigrid convergence.</div><div style="font-size:12.8px"><br></div><div style="font-size:12.8px">One possibility is that you solution is going crazy and your density is getting non-smooth and thus the coefficients in the pressure solve are getting large and non-smooth (I assume you have the density term in the pressure solve). This would cause a problem for multigrid.</div><div><br></div><div> <br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">All PETSc solvers converge but come with huge iterations. That's when/how I started using HYPRE.</div><div class="gmail-m_-397246587868680201m_-7915301013834964351m_-3679888773484715993gmail-HOEnZb"><div class="gmail-m_-397246587868680201m_-7915301013834964351m_-3679888773484715993gmail-h5"><div class="gmail_extra"><br></div></div></div></blockquote></div></div></div>