icofoam doubt
 

Hello all,
I am new to openfoam and was just trying to understand icofoam code.
in the follwing portion of the code
set up the pressure equation
fvScalarMatrix pEqn
(
fvm::laplacian(rUA, p) == fvc::div(phi)
);
for setting pr eqn divergence of phi is used. if phi is scalar field(inner product of face velocity and surface area normal) then how can its divergence be defined ?
thanks in advance.

phi is actually a surfaceScalarField, which is a scalar (as you have noted) defined at each face of the particular cell. The divergence of this is then the summation of this surfaceScalar over the owner and neighbor faces. If you want to learn more about this, have a look at Dr. Jasak's Thesis here .

Thanks for ur kind response.