Hi, I am trying to solve an incompressible navier stokes on an adaptive mesh, and I find that using a petsc fieldsplit preconditioner which works reasonably well on a normal mesh performs very very poorly on an adaptive mesh. My solver options are as follows
-ns_solve_pc_fieldsplit_detect_saddle_point -ns_solve_pc_type fieldsplit -ns_solve_fieldsplit_0_ksp_type gmres -ns_solve_fieldsplit_0_pc_type hypre -ns_solve_fieldsplit_1_pc_type lsc -ns_solve_fieldsplit_1_ksp_type gmres -ns_solve_fieldsplit_1_lsc_pc_type hypre Is there anything in particular about the matrices generated in the navier stokes problem that leads it to fail badly.. ( I usually see the solver converge in 2 steps but with the adaptive mesh, it converges linearly) . I have an associated cahn hilliard solver, and that has no problem when solved after adaptivity. ' Subramanya G Sadasiva, Graduate Research Assistant, Hierarchical Design and Characterization Laboratory, School of Mechanical Engineering, Purdue University. "The art of structure is where to put the holes" Robert Le Ricolais, 1894-1977 ------------------------------------------------------------------------------ Minimize network downtime and maximize team effectiveness. Reduce network management and security costs.Learn how to hire the most talented Cisco Certified professionals. Visit the Employer Resources Portal http://www.cisco.com/web/learning/employer_resources/index.html _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
