CFD Online Discussion Forums - FreeFem++

Hi,

I need to solve the first order macroscopic pedestrian flow model in 2D in FreeFem:

rho_t + div(rho*V(rho)*mu) = 0,

where V(rho) is some speed-density relation and mu is constant for my case, but I can't figure out, how to write the code of my problem in FreeFem.

I think that it will be:

Code:

problem pedestrian(rho, rhoh) =

int2d(Th)( rho*rhoh/dt )

+int2d(Th)( ....... )

+on (some boundary condition);

Please, can you help me with it? I'm a beginner in FreeFem.
Thanks.