I am now working on Navier-Stokes equations with PETSc. I am trying to solve the equations in a
125*100*100 domain, and build this cube domain mesh by Gmsh. However, when I try to solve it with a higher resolution (finer mesh), the simulation will return to a Nan error. The solver I used is:
set(NS, sparams = “-pc_type lu”);
I want to know why this error occurs and if there is any way to solve it.
Thanks in advance!