CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > REEF3D

Flow "sticking" in gaps - REEF3D CFD

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

LinkBack Thread Tools Search this Thread Display Modes
Old   March 22, 2024, 05:06
Default Flow "sticking" in gaps - REEF3D CFD
New Member
Join Date: Mar 2024
Posts: 1
Rep Power: 0
kcoastal is on a distinguished road
Hi all,

I'm a new user to REEF3D and have been very impressed thus far having worked through the tutorials (thanks to the development team).

I was recently tasked with assessing wave attenuation through a "wave screen" on the crest of breakwater structure. Specifically, I need to determine if a structure comprised of vertical members provides better/equivalent/worse attenuation than horizontally placed members.

The wave screen is comprised of 600mm x 300mm x 2000mm structures with a 200mm gap between each. This can be rotated by 90deg to represent vertical / horizontal.

When I run my simulation with the vertical members, water appears to flow in-between the gaps as would be expected. However, when I run the horizontal simulation, the water appears to get "stuck" between the gaps, almost as if the members are "sticky". I found this surprising given everything about the simulations are the same re. domain size and grid spacing (0.006m).

My question: Is there a minimum number of cells that should be used to represent gaps and do the structures have a friction factor? Or could this issue be related to wet/drying of cells?

I have attached a few comparison figures which hopefully explains the issue along with the control and ctrl.txt files.
Attached Images
File Type: jpg Comparision_1.jpg (35.7 KB, 30 views)
File Type: jpg Flow_sticking_n_gaps.jpg (48.4 KB, 24 views)
Attached Files
File Type: txt control.txt (420 Bytes, 5 views)
File Type: txt ctrl.txt (502 Bytes, 6 views)
kcoastal is offline   Reply With Quote

Old   March 24, 2024, 03:26
Felix S.
Join Date: Feb 2021
Location: Germany, Braunschweig
Posts: 85
Rep Power: 6
Felix_Sp is on a distinguished road

this is a problem with the reinitialization of the level-set in REEF3D. It inhibits the propagation on "dry" cells. Due to this problem, the calculation of wave run-up with REEF3D is computationally very costly, as the only effective and physically accurate solution is to increase the grid resolution.

You might try F 40 23 or F 49 0. While F 40 does not help a lot in my experience, F 49 does change free surface propagation quite a lot. However, the overall simulation gets unphysical as the overall mean water level does increase. F 49 0 only worked for me, while staying physical, with very long waves (wave period ~ 60 seconds).

In regard to your second question. If you mean roughness when you say friction factor, then yeah, your Nikuradse roughness is B 50 0.000001 on the solid.

So in conclusion, you might try local or global (if you have the computational power) grid refinement. Also, is your setup/ your wave screen symmetrical? You might save a lot of computational power by having a setup with a symmetry boundary conditions in the mid of the flume. So B 10 0.0 47.0 0.0 2.5 0.0 7.0 and C 12 3 in the control.txt for your setup.

Kind regards
Felix_Sp is offline   Reply With Quote


Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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
Is recirculating flow is free vortex or forced vortex? FluidKo Main CFD Forum 11 July 21, 2022 06:21
PhD in turbulence Hans Main CFD Forum 14 October 8, 2001 03:03
ASME CFD Symposium - Call for Papers Chris Kleijn Main CFD Forum 0 September 25, 2001 10:17
CFD Modeling of Two-phase Flow in Small Dia.Tubes Eric Poindexter Main CFD Forum 2 September 22, 2000 09:21
ASME CFD Symposium, Atlanta, July 2001 Chris R. Kleijn Main CFD Forum 0 September 13, 2000 04:48

All times are GMT -4. The time now is 08:42.