I am trying to implement a cahn-hilliard navier stokes solver using libMesh. I got the whole thing working without adaptivity. When I add adaptivity to the mix though - the cahn hilliard part of the system solves without any problems, however, when I get to the navier stokes part the solver fails for even very few elements (<1000) . The solver options that I am using are as follows
Are there limitations on the discretizations or linear solvers that people know of while solving the navier stokes equation with adaptive mesh refinement?
"But memory eventually fades. Turbulences damp out, internal strains yield to plastic flow, concentration inhomogeneities diffuse to uniformity. Systems tend to subside to very simple states,independent of their specific history."
Herbert Callen . Thermodynamics and an Introduction to Thermostatics.