[meshkati]

Hey everyone,
I try to develop a CFD code by solving 3D Navier-Stocks equations in curvilinear coordinates. But I have a problem with coupling of pressure and velocity (pressure correction equation - p’) that I use SIMPLEC method for it.
In initial iterates (that no pressure gradients exists in domain) the momentum equations have no source term and then: ap=sigma(a_neighbors).
Therefore the parameter of [ap - sigma(a_neighbors)] that appears in the denominator of SIMPLEC method equals to zero and then the divides by zero error was occurred.
What can I do?
Thanks for advance.
Hadi
(Some of my code specifications are: Finite volume method, non -staggered grid, incompressible fluid, steady calculation, no free surface calculation, Rhie and Chow interpolation method, laminar or k-e turbulence model, SIMPLEC pressure correction method and the test case is straight channel with uniform velocity boundary condition in its inlet and if any another comments was needed just tell me)