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

Combined Inflow and Initial Water Condition

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By IADInocencio

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 4, 2021, 05:03
Default Combined Inflow and Initial Water Condition
  #1
Member
 
Ismael Incoencio
Join Date: Sep 2021
Posts: 45
Rep Power: 7
IADInocencio is on a distinguished road
Is it possible to combine the inflow boundary condition B 60 1 with an initial water condition wherein only a portion of the domain is filled with water? Whenever I try to combine the two I get the "EMERGENCY STOP -- velocities exceeding critical value N 61" error.

So far I've only gotten it to work with either condition, but not both at the same time.
Attached Images
File Type: jpg Capture.jpg (29.2 KB, 38 views)
Attached Files
File Type: txt control.txt (4.6 KB, 10 views)
File Type: txt ctrl.txt (302 Bytes, 12 views)
IADInocencio is offline   Reply With Quote

Old   November 11, 2021, 08:46
Default
  #2
Senior Member
 
Arun Kamath
Join Date: Nov 2014
Location: Trondheim, Norway
Posts: 265
Rep Power: 13
kamath is on a distinguished road
The reason why your simulation crashes is clear: the discharge.
If you set the discharge to 0.0, you will see that the simulation works.

Now as to why this is so, I would guess:
When you set a discharge outlet BC and there is no water at the outlet since you have a "dam break" scenario, the solver the crashes.

Is there a particular reason you stop the water right before the structure?
I also see (may be due to the grid res) that the set water level in the z-direction is lower than the outlet in your object, so there will not be a flow anyway.
What you might also notice is that if you have water filled across the entire x-direction to z=0.6 that you have chosen, the water level to the right of the structure will simply start to reduce to account for the discharge.
On the left hand side, the reinitialisation of the level set function will probably keep restricting the water level to before the structure. Eventually that simulation will crash too.

I hope this gives a different perspective regarding your setup.
__________________
Arun
X years with REEF3D
kamath is offline   Reply With Quote

Old   November 11, 2021, 09:22
Default
  #3
Member
 
Ismael Incoencio
Join Date: Sep 2021
Posts: 45
Rep Power: 7
IADInocencio is on a distinguished road
Hi Arun,

Yes I found that setting W 10 to 0.0 or just turning it off would allow the simulation to continue, however my intent was for water to flow out to simulate a drainage outfall. So that forebay in the picture is to simulate a manhole with water coming from an upstream pipe, and the opening is the pipe leading to a stepped outfall to dissipate the velocity.

Yes, the intent for the initial water to stop just at the structure was just to initialize that volume so I can have a shorter simulation duration than if it was empty. To illustrate, I have a smaller model (as attached) I used to learn and got a working result as shown below, where my initial water level was at A. And with the discharge Q it took some time T for the WSE in the upstream to raise B meters, and also took some time for the water downstream of the outfall to normalize because it was also wet at the start of the simulation.

So in an attempt to shorten the simulation time by T, I was trying to already have the WSE just below the opening and only upstream of it, so it discharges onto a dry stepped outfall.

I did some other tests (attached .txt's) with an even smaller domain and simpler structures and found water can inflow with a PatchBC, but I can't seem to get it to work correctly. After the initial pulse of water due to the numerical warm up causing some overtopping, it seems the forebay never fills up. https://youtu.be/7qS-H0HbP5A
Attached Images
File Type: jpg Capture.jpg (40.5 KB, 24 views)
Attached Files
File Type: txt control.txt (207 Bytes, 4 views)
File Type: txt ctrl.txt (304 Bytes, 7 views)
valgrinda likes this.
__________________
Still learning
IADInocencio is offline   Reply With Quote

Old   November 11, 2021, 11:52
Default
  #4
Member
 
Ismael Incoencio
Join Date: Sep 2021
Posts: 45
Rep Power: 7
IADInocencio is on a distinguished road
Here is a better screenshot of what I was trying to do, where the outfall is discharging but the downstream should initially be dry.
Attached Images
File Type: jpg dry.jpg (36.0 KB, 38 views)
__________________
Still learning
IADInocencio is offline   Reply With Quote

Old   November 15, 2021, 09:10
Default
  #5
Senior Member
 
Arun Kamath
Join Date: Nov 2014
Location: Trondheim, Norway
Posts: 265
Rep Power: 13
kamath is on a distinguished road
Interesting problem you are dealing with. My question is:
How essential is it to start from a dry condition?

The reason I ask this is:
If you look at literature on, say, dam break cases you will see that there is a certain (though very small) amount of water that exists downstream of the break on start up. eg. here and here

The reason being something that you have already experienced. The beginning of a "wet film" due to overflow of the fluid is hindered by the reinitialisation of the LS function. An extended BC on this would probably help somewhat, but I am unsure to what extent it would be good.

So, unless it is of utmost importance to have a dry start, I would recommend a wet start with an initial water level of about 1-2 cells on the downstream side.
__________________
Arun
X years with REEF3D
kamath is offline   Reply With Quote

Old   November 15, 2021, 21:29
Default
  #6
Member
 
Ismael Incoencio
Join Date: Sep 2021
Posts: 45
Rep Power: 7
IADInocencio is on a distinguished road
Hi Arun, it doesn't have to be 100% dry, having an initial downstream depth of around 0.2-0.5 is alright for me actually. My intent was to just have the steps I defined as dry. I realized from your previous reply too that the water should at least reach the downstream BC.

Thank you for the links, it's actually exactly like the initial water setup in Figure 2 of the 2nd article that I was trying to recreate but could not figure out the correct way to implement it using the F 51-F 56, F 70, and F 71 parameters. I was always ending up with just a flat initial water surface or just in the forebay; not both.

Again thank you for your time! There are lots engineering problems that I have been wanting to analyze using CFD and I am learning quite a lot on how to go about it.
__________________
Still learning
IADInocencio is offline   Reply With Quote

Old   November 16, 2021, 05:07
Default
  #7
Senior Member
 
Arun Kamath
Join Date: Nov 2014
Location: Trondheim, Norway
Posts: 265
Rep Power: 13
kamath is on a distinguished road
Glad I was able to help.
Good luck ahead!
__________________
Arun
X years with REEF3D
kamath is offline   Reply With Quote

Old   November 16, 2021, 09:28
Default
  #8
Super Moderator
 
Hans Bihs
Join Date: Jun 2009
Location: Trondheim, Norway
Posts: 403
Rep Power: 19
valgrinda is on a distinguished road
Hi Ismael,

maybe check B 77 2?
__________________
Hans Bihs
Team REEF3D
www.reef3d.com
valgrinda is offline   Reply With Quote

Old   November 16, 2021, 20:56
Default
  #9
Member
 
Ismael Incoencio
Join Date: Sep 2021
Posts: 45
Rep Power: 7
IADInocencio is on a distinguished road
Good morning Prof Hans, thanks for suggesting B77 2! I activated this and started off with a high initial water surface instead and the downstream water level is become lower as expected. (Side note, the entry in the manual for B 77 refers to a B 76 but is not listed?)

Thanks Arun; it looks like using the B 77 is another simple approach to my solution. Rather than my initial approach of letting the water in the forebay rise and spill to the steps, I'll let it work backward (in a sense) and fill the entire domain with a high initial water level and, hopefully, using the B77 will allow the downstream to normalize to a shallower depth.

Will let this run for a while and provide updates for science
__________________
Still learning
IADInocencio is offline   Reply With Quote

Old   November 30, 2021, 22:03
Default
  #10
Member
 
Ismael Incoencio
Join Date: Sep 2021
Posts: 45
Rep Power: 7
IADInocencio is on a distinguished road
Quote:
Originally Posted by kamath View Post
Interesting problem you are dealing with. My question is:
How essential is it to start from a dry condition?

The reason I ask this is:
If you look at literature on, say, dam break cases you will see that there is a certain (though very small) amount of water that exists downstream of the break on start up. eg. here and here

The reason being something that you have already experienced. The beginning of a "wet film" due to overflow of the fluid is hindered by the reinitialisation of the LS function. An extended BC on this would probably help somewhat, but I am unsure to what extent it would be good.

So, unless it is of utmost importance to have a dry start, I would recommend a wet start with an initial water level of about 1-2 cells on the downstream side.
Hi Arun, apologies but this scenario has become an itch that I must scratch and have gone back to attempting simulating a domain where the water is higher on one side of a barrier with a discharge coming in. I have been able to successfully implement two different water levels as you recommended using the F 72 parameter but flows were still not going over the obstacle. So to compare I made two different simulations:
1) Scenario 1
- Inflow set by B 411 and B 442
- Constant water level for the whole domain set by F 72
2) Scenario 2
- Inflow set by W 10
- Initial water level set by F 60

Scenario 2 works as expected (although looks like accuracy could be better, likely due to a 0.02 grid size and sudden influx of water), but in Scenario 1 water just seems to disappear. The screenshot below is at 3 seconds, with the red surface set to Scenario 1, and the blue surface as Scenario 2.

I've also attached the ctrl files.

Which of these parameters am I implementing incorrectly?
Attached Images
File Type: jpg Capture.jpg (34.1 KB, 27 views)
Attached Files
File Type: txt 1 control.txt (205 Bytes, 6 views)
File Type: txt 2 control.txt (204 Bytes, 11 views)
File Type: txt 2 ctrl.txt (378 Bytes, 18 views)
File Type: txt 1 ctrl.txt (410 Bytes, 13 views)
__________________
Still learning
IADInocencio 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
Correct BCs for known outlet conditions and unknown inlet conditions cfdcheckers OpenFOAM Running, Solving & CFD 20 June 17, 2021 03:56
High number iteration of pressure equation mgg OpenFOAM Running, Solving & CFD 4 June 3, 2019 00:33
HronTurekFsi3 foam-extend3.1 Maimouna OpenFOAM Running, Solving & CFD 112 April 20, 2019 22:11
Proper boundary & initial value for a water turbine Yohan OpenFOAM Running, Solving & CFD 0 January 31, 2018 02:05
Final residual greater than initial residual? jsumner OpenFOAM Running, Solving & CFD 2 March 2, 2011 09:19


All times are GMT -4. The time now is 19:40.