<div dir="ltr">Hi Folks,<div><br></div><div>I am trying to solve the momentum equation in a projection preconditioner using GAMG or Hypre solver. The equation looks like for velocity variable <i><span><span style="font-style:italic">v</span></span></i> looks like:</div><div><br></div><br><div> <img src="cid:ii_m2tn1esj1" alt="Screenshot 2024-10-28 at 4.15.17 PM.png" width="542" height="72"><br><div><br></div><div>Here, <span><span style="font-style:italic">μ</span></span> is spatially varying dynamic viscosity and <span><span style="font-style:italic">λ</span></span> is spatially varying bulk viscosity. I understand that I need to specify rigid body nullspace modes to the multigrid solver in order to accelerate its convergence. Looking into this routine MatNullSpaceCreateRigidBody() (<a href="https://urldefense.us/v3/__https://petsc.org/release/manualpages/Mat/MatNullSpaceCreateRigidBody/__;!!G_uCfscf7eWS!eKqgIJjCdMzIU76f7X65AmGxrU_-lC7W02BMWafJ77DNf_IuQk6O1X3qU1x9Ez8NJ20vZEL-mF6T1yNmDnwv0eWa2w$">https://petsc.org/release/manualpages/Mat/MatNullSpaceCreateRigidBody/</a>), I see that I need to provide the coordinates of each node. I am using staggered grid discretization. Do I need to provide coordinates of staggered grid locations?</div><div><br></div><div>Thanks,  </div><span class="gmail_signature_prefix">-- </span><br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div>--Amneet <br><br></div><div><br></div><div><br></div></div></div></div></div>