# Body force in fvOptions in interFoam

 Register Blogs Members List Search Today's Posts Mark Forums Read

 January 27, 2022, 03:55 Body force in fvOptions in interFoam #1 New Member   Join Date: Mar 2020 Posts: 16 Rep Power: 5 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?

June 23, 2022, 13:39
#2
Member

Mahmoud
Join Date: Nov 2020
Location: United Kingdom
Posts: 43
Rep Power: 4
Quote:
 Originally Posted by Geb1313 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?

 Tags fvoptions, interfoam