Dear developers,

For the Navier-Stokes equations, there is a null space in the block matrix 
(right-bottom block)
--                         --
|  -laplace     grad  |
|    div          0       |
--                         --

I use bjacobi+ilu as the preconditioner for GMRES in PETSc to solve the system. 
Sometimes it works, but most of time it does not converges.
My code is just like the example: /systems_of_equations_ex2.

I wonder if I try the way in example: /systems_of_equations_ex3, would it be 
better?
Or should I use other solver/perconditioner ?

I am appreciated for your suggestions.


--
Li Luo



------------------------------------------------------------------------------
What NetFlow Analyzer can do for you? Monitors network bandwidth and traffic
patterns at an interface-level. Reveals which users, apps, and protocols are 
consuming the most bandwidth. Provides multi-vendor support for NetFlow, 
J-Flow, sFlow and other flows. Make informed decisions using capacity planning
reports.http://sdm.link/zohodev2dev
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to