[marango]

Hello everybody,

I'm a newbie in OpenFOAM and have some problems with the implementation in OpenFOAM. As you can see in the title I'm examing flow through porous media using a new approach.

Similar to the Darcy-Forchheimer-Equation I want to add a new term to rhs of the incompressible Navier-Stokes-Equation. As a basis I use the simple-solver in OpenFOAM. Now I'd like to change the implemented differential equation. The rhs looks in x-direction as follows:

... = dp/dx - a* u_x - b* (u_x)² .

( a,b - scalar values)
I tried a lot to implement it. But certainly I don't know the right commands. Can anybody help me?

Many thanks in advance!