CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > REEF3D

Dam Break process problems

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 25, 2021, 21:18
Default Dam Break process problems
  #1
New Member
 
NSW
Join Date: May 2020
Posts: 20
Rep Power: 7
pipi is on a distinguished road
Dear Hans,

Good morning!

I use REEF3D to simulate the dam break process with dry bed and compare the result with Ritter (1982) analytical solution. There are some problems I got.

1. I set the water level at 0.12m, but when I extract data from REEF3D result, the water level has slightly lower than the setting value. And the vertical surface of the dam is slightly inclined. (Figure 1).

2. I already set the slip bed, with no friction, but the result shows the bed is not slipped (Figure 2).

Do you have some suggestions for these?

Kind regards
Dongfang
Attached Images
File Type: jpg Figure 1.jpg (31.7 KB, 24 views)
File Type: png Figure 2.png (64.0 KB, 25 views)
Attached Files
File Type: txt control.txt (103 Bytes, 9 views)
File Type: txt ctrl.txt (194 Bytes, 5 views)
pipi is offline   Reply With Quote

Old   October 28, 2021, 17:01
Default
  #2
Super Moderator
 
Hans Bihs
Join Date: Jun 2009
Location: Trondheim, Norway
Posts: 403
Rep Power: 19
valgrinda is on a distinguished road
Hi Dongfang,

I suggest the following changes to the setup:

- in the control.txt, C 12 and C13 need to be 3, because you run in 2D and you need symmetry plane boundary conditions there

in the ctrl.txt:
- B 20 2 works better

- N 40 3 (TVD RK3) is much better than N 40 1 (AB 2)

- you can add F 42 1.0 to speed up the initialization process

- N 47 (CFL number) should be fine at 0.1

- there must be a problem in how you extract the data, Paraview shows the free surface to be located at exactly 0.12 m - so everything is ok
__________________
Hans Bihs
Team REEF3D
www.reef3d.com
valgrinda is offline   Reply With Quote

Old   November 3, 2021, 03:52
Default
  #3
New Member
 
NSW
Join Date: May 2020
Posts: 20
Rep Power: 7
pipi is on a distinguished road
Dear Hans,

Appreciate your suggestions!

The N 40 3 (TVD RK3) and N42 1.0 works well!

The reason I set B 20 1 (slip) is that if I set B 20 2 (no-slip bed), the result from REEF3D does not match the Ritters analytical solution very well near the bed. However, even I set slip bed (B 20 1) and no friction, no viscosity, it still shows something resists the water propagates near the bed in REEF3D. The result still cannot create a sharp front as Ritter's analytical solution.

Do you have some ideas for this point?

Thanks in advance.

Kind regards
Dongfang
Attached Images
File Type: jpg Figure.jpg (38.4 KB, 28 views)
pipi is offline   Reply With Quote

Old   November 3, 2021, 06:39
Default
  #4
Super Moderator
 
Hans Bihs
Join Date: Jun 2009
Location: Trondheim, Norway
Posts: 403
Rep Power: 19
valgrinda is on a distinguished road
Hi Dongfang,

I think the results look pretty good ! Differences like you show should be acceptable for a numerical model, anything else would be surprising.
__________________
Hans Bihs
Team REEF3D
www.reef3d.com
valgrinda is offline   Reply With Quote

Old   November 3, 2021, 06:45
Default
  #5
Super Moderator
 
Hans Bihs
Join Date: Jun 2009
Location: Trondheim, Norway
Posts: 403
Rep Power: 19
valgrinda is on a distinguished road
Hi Dongfang,

sorry for the previous post, I somehow read your figure differently. I think you are still running a 2D simulation with walls in y-direction, see my previous comments. You need to set C 12 and C 13 to 3. Then B 20 2 will be better as well.

Also, I think you are running a domain that is a bit too large for your phenomenon at hand. A tank with smaller length would also do.
__________________
Hans Bihs
Team REEF3D
www.reef3d.com
valgrinda is offline   Reply With Quote

Old   November 7, 2021, 20:59
Default
  #6
New Member
 
NSW
Join Date: May 2020
Posts: 20
Rep Power: 7
pipi is on a distinguished road
Thanks, Hans,

The results are the same as the figures I show before after I change the settings.

But they are already pretty good enough.

Thanks for your help!

Kind regards
Dongfang
pipi is offline   Reply With Quote

Old   November 9, 2021, 10:45
Default
  #7
Super Moderator
 
Hans Bihs
Join Date: Jun 2009
Location: Trondheim, Norway
Posts: 403
Rep Power: 19
valgrinda is on a distinguished road
Hi Dongfang,

could you please upload your most recent control files once more? I am pretty sure that the REEF3D results should be better than depicted in your figure.
__________________
Hans Bihs
Team REEF3D
www.reef3d.com
valgrinda is offline   Reply With Quote

Old   November 9, 2021, 21:30
Default
  #8
New Member
 
NSW
Join Date: May 2020
Posts: 20
Rep Power: 7
pipi is on a distinguished road
Dear Hans,

Thanks!

The attachment is my control files and results (Figure 1).

The initial water level change has been fixed because I use the line probe in the model directly instead of extracting data from Paraview.

Paraview defines the bottom of our model domain at the middle of the bottom cell (Figure 2), which means if we want to use Paraview to read water surface level, we lose half of our mesh-size height, that's the reason why the previous figures show some water level loss (Figure 3).

If you have any suggestions for this case, please contact me, very much appreciate it!

Kind regards
Dongfang
Attached Images
File Type: png Figure 1.png (62.0 KB, 20 views)
File Type: jpg Figure 2.jpg (119.7 KB, 21 views)
File Type: png Figure 3.png (68.6 KB, 19 views)
Attached Files
File Type: txt control.txt (137 Bytes, 5 views)
File Type: txt ctrl.txt (202 Bytes, 7 views)
pipi is offline   Reply With Quote

Reply


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 Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
OpenFOAM Tutorial: Simulation of Dam break using interFoam foamiste OpenFOAM Running, Solving & CFD 0 August 16, 2021 06:12
The fl process could not be started because of UDF majid_kamyab Fluent UDF and Scheme Programming 6 December 15, 2015 09:42
compressible dam break simulation smsanth OpenFOAM Running, Solving & CFD 0 September 2, 2013 03:51
help- free surface level- dam break TUL FLUENT 1 December 6, 2012 00:29
Dam break problem sega OpenFOAM Running, Solving & CFD 3 April 20, 2008 11:03


All times are GMT -4. The time now is 18:58.