Body force in fvOptions in interFoam
 

Hello Foamers,
I am thinking to model a dam break using interFoam solver. The dam has a rectangular shape as most of the dam break tutorials, but with a solid rectangular block(constant density) on the top of it. As per my understanding the impact of the block on cells beneath it is a pressure force of magnitude density*gravity* depth of the block acting on the vertical surface area(xy plane) of the cells.

Now my question is how can add this force as a coded source.

My first thought was: providing the force as P = (0 0 rho*g*h) and to be included into the momentum equation it should be force per unit volume. I was lost after this. The actual dimension of the force P is dimension of pressure. I have to change it to weight (force) by: F = P &mesh.sf(). And dividing F by mesh.V() provides the dimension required. However, I am not sure if this is possible.

Any recommendation please?

Hi

Have you got an answer for this problem?