CFD Online Discussion Forums - Implementation of a very simple solver

Dear All,

I am trying to figure out how to manage a very simpe solver. I need only the continuity equation to be solved.

What I want to do is the following:

1. my domain is very simple: an empty cube with 5 walls and 1 inlet;
2. I insert a certain flow rate from the inlet and I want to fill all the domain;
3. the idea is to evaluate the density of each cell;
4. in the cells with density < 1, everything is fine; where the density is bigger, I do the following: Delta_rho = rho - 1. Then I "the" the Delta_rho in the neibouhring cells.

Does anyone have an idea about how to solve this?

Thanks a lot,
Samuele