I have to add a source term in the mass conservation equation. This source term is a function of local 'y' co-ordinate and also local 'Uy' (y-component of U)
Lets say, for example....
ddt(rho) + div(phi) == -(rho*Uy)/y

For the correct addition of a massSource have a look at the dieselFoam-sources.

__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request