Just a simple question : why incomplete factorization ? Why don't you consider all the terms in your factorization and calculate the complete LU ?

Moreover the laplacian operator involved in your equation is symmetric, thus you may use the LLt factorization (Cholevski) which is faster (half the time) since you only have to perform the L matrix calculation instead of L and U

As Paul mentioned, you are likely to find all the algorithms you need on www.netlib.org

Sir i am working on research project "revemping of Francis Turbine applying CFD". I have got book "An introduction to computational fluid dynamics (finite volume) by malalasakra" it is very good book and i have understand SIMPLE, SIMPLER AND SIMPLEC method of solving Navier Stock equation but in this book just presentation that how Navier Stock equation can be solved but there is no solve example. Do you know any book or article in which this type of example is present. I mean that for any given problem(good if turbine rotor)firstly assume pressure field and then find velocity field then check pressure correction equation and repeat this procedure still convergence occure. waiting for your prompt reply thanks