# Body force in fvOptions in interFoam

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 LinkBack Thread Tools Search this Thread Display Modes
 January 27, 2022, 04:55 Body force in fvOptions in interFoam #1 New Member   Join Date: Mar 2020 Posts: 16 Rep Power: 4 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, 14:39
#2
Member

Mahmoud
Join Date: Nov 2020
Location: United Kingdom
Posts: 37
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

 Thread Tools Search this Thread Search this Thread: Advanced Search Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are Off Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post ds4719 Main CFD Forum 14 February 18, 2022 19:05 hxax STAR-CCM+ 3 January 16, 2020 11:08 be_inspired OpenFOAM Programming & Development 10 September 14, 2018 12:12 traens OpenFOAM Programming & Development 2 July 17, 2017 05:14 selvam R.P CFX 4 November 25, 2005 05:01

All times are GMT -4. The time now is 17:57.

 Contact Us - CFD Online - Privacy Statement - Top