Hello, I am trying to use ILU preconditioner for Navier stokes problems. None of the standard ILU preconditioners in the SPARSKIT package is working due to the presence of zero's in the pressure (continuity) equation. Does any one has any suggestion about how to implement ILU for Navier stokes problems.

You could try modifying the continuity equation to include a pressure term; e.g. replace div(u) = 0 by something like

k.[p(n) - p(n-1)] + div(u) = 0

where p(n) is the value of p(n) is p at the n'th iteration and k is a constant.

This reminds me of the famous song...

Ye ILU ILU kya hai ? Ye ILU ILU (what is this ILU ILU ?)

ILU ka matlab I love you ! (ILU means I love you...)

I hope having a little fun on this forum is not that bad...

Try to Use Artificial Compressibility Method of

Chorin which include pressure in the continuity equation.

cheer

KK.Khan

Hi: I also want to apply the SPARSKIT to the N-S equations, I think we can use a block implicit method. what is you idea

precence of zeros in diagonal entries is not desirable and is related to poor physical or numerical modeling, but there are some cure (althouth i recommend modifiung modeling, projection method instead probably penalty method), you can use pivoting to change position of zeros to off-diagonals or you can use BILUM package, it is opensource due to Y. Saad and designed for solution of such ill-pose systems, it is not only rubost but also enjoys properties of algebraic multigrids (scalability and fast convergence) free of AMG limitations such as SPD condition.