# Problem modeling hydraulic jump in an open channel with interFoam (VoF)

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

November 5, 2013, 14:42
Problem modeling hydraulic jump in an open channel with interFoam (VoF)
#1
New Member

Arnau
Join Date: Jan 2012
Posts: 17
Rep Power: 13
Dear FOAMers,

I am trying to model a hydraulic jump in a rectangular open channel using interFoam but so far I have not managed to define successfully the boundary conditions. I want to force a supercritical inlet (imposing a given velocity and the water depth) and a subcritical outlet (imposing an experimentally-determined water depth and a hydrostatic pressure profile). Please see the attached figure for further details.

I have tried thousands of boundary conditions using groovyBC but none of them worked, apparently, due to a bad definition of the outlet boundary condition. I modeled hydraulic jumps in the past by means of downstream steps, slopes and so on. But this time I have to do it directly imposing the water depth at the outlet. The core of the problem seems to be how to make the outlet boundary condition raise the water level.

Has anybody done something similar before with OpenFOAM? Any advice? I have not found any solution in other threads.

Thank you very much,

Arnau.
Attached Files
 Hydraulic Jump Scheme.pdf (27.9 KB, 214 views)

 November 6, 2013, 06:24 #2 Senior Member   Join Date: Aug 2010 Location: Groningen, The Netherlands Posts: 216 Rep Power: 17 Dear Arnau, I don't think that your problem is related to the inlet or outlet boundary conditions. The problem is that for a hydraulic jump you need a high velocity in the beginning which is then slowed down and the water level increases. Usually this is no problem when using sloped geometries like a backwards facing slope. However dealing only with a plain ground this is slightly more difficult. What you can do now is in my opinion to increase the length of the domain. With now having more surface which causes friction maybe an hydraulic jump will occur. If this does not help I would also increase the speed of the water at the inlet. I hope my explanations help kind regards Colin

 November 6, 2013, 06:39 #3 Senior Member   Albrecht vBoetticher Join Date: Aug 2010 Location: Zürich, Swizerland Posts: 234 Rep Power: 15 See the trick of kflora in http://www.cfd-online.com/Forums/ope...tml#post391836 at 5th November 2012. Was aswell applied by pythagOra5 and compared to flume experiments with satisfying results.

 November 6, 2013, 23:00 #4 Senior Member   Matthew Denno Join Date: Feb 2010 Posts: 138 Rep Power: 15 vonboett, Could you provide a link or reference for the flume experiments comparison? Thanks, Matt

 November 7, 2013, 14:51 #5 Senior Member   Albrecht vBoetticher Join Date: Aug 2010 Location: Zürich, Swizerland Posts: 234 Rep Power: 15 well Matthias send me the comparison to his experiments by mail so I can't tell if its O.K. to post them here, maybe you could directly contact him: http://www.cfd-online.com/Forums/mem...ythag0ra5.html His results deviate from the experiments due to turbulence at an obstacle in the center of the domain, but his flow height at the outlet is stable.

 November 11, 2013, 04:20 Trying Arnau case with kflora suggested bc #6 New Member   Antonio Join Date: Jan 2013 Posts: 11 Rep Power: 12 I was simulating a lab channel in which among other things, there is a hydraulic jump. I started simulating the hydraulic jump isolated, trying for that purpose of creating it from a step and establishing zeroGradient conditions in the outlet over the step. Nonetheless, after reading Vonboett, I decided to give a try to kflora suggested BC, as they are closer to typical input data in hydraulic problems. At the moment, I have implemented Arnau sketch, but establishing buoyant pressure in either inlet and outlet. A fixed speed in the inlet (with alpha1 = 1), and a fixed speed in the whole outlet (in order to get y2) for air and/or water. Over the inlet I have try two different bc, the atmospheric one and also a wall (nutkWallFunction) with 0 speed. On the channel floor I fixed a zero speed (U) and zeroGradient (alpha1) with a nutkRoughWallFunction. In both cases the hydraulic jump starts moving towards upstream but instead of stopping in the channel it locates exactly in the inlet. If wall over the inlet is chaged with an atmospheric bc, there is even some small negative flow rate over the inlet. I have try enlarging or shortening the channel, but hydraulic jump position remains the same. Even changing lightly y1 or U1. I am using the k-epsilon model for the turbulence (and interFoam as solver) I upload here the case and some plots of the commented results Animation where Ux is plotted during 100s until a steady state is reached. The volume plotted corresponds to alpha1>0.50 https://drive.google.com/file/d/0Bzj...it?usp=sharing In the steady-state, this figure shows U: https://drive.google.com/file/d/0Bzj...it?usp=sharing Basic case files: https://drive.google.com/file/d/0Bzj...it?usp=sharing

 November 14, 2013, 05:51 #7 Senior Member   Albrecht vBoetticher Join Date: Aug 2010 Location: Zürich, Swizerland Posts: 234 Rep Power: 15 to get a solution comparable to experiments you should go for three dimensions and maybe use LES

November 18, 2013, 03:52
#8
New Member

Antonio
Join Date: Jan 2013
Posts: 11
Rep Power: 12
Quote:
 Originally Posted by vonboett to get a solution comparable to experiments you should go for three dimensions and maybe use LES
Thanks for your answer... I was trying first to make it comparable to Berlanguer equation... now I have been trying to increase momentum, keeping q constant, and at the same time following your suggestion I have re-run simulation in 3D with a symmetry plane (but still with k-e). I will comment here results later on.

 September 2, 2021, 04:40 Hydraulic jump does not occur where it should #9 Member   Join Date: Apr 2015 Posts: 32 Rep Power: 9 hello everyone, I'm simulating a free hydraulic jump with InterFoam. In the laboratory model, the jump occurs at a distance of one meter from the gate, but in the simulation, despite changes in various parameters such as mesh and inlet velocity, this does not happen and the jump occurs immediately after the gate. Can anyone help me? I will appreciate that,

 September 28, 2021, 06:52 #10 Senior Member   Albrecht vBoetticher Join Date: Aug 2010 Location: Zürich, Swizerland Posts: 234 Rep Power: 15 Hi What turbulence model do you use, what is the Froude number of your flow? Hydraulic jumps in mountain torrents were captured by interFoam quite well using hybrid URANS-LES...

 September 28, 2021, 09:46 #11 Member   Join Date: Apr 2015 Posts: 32 Rep Power: 9 hi my Fr=4.25 and turbulence model is k-e RNG

 September 29, 2021, 05:46 #12 Member   Michael Sukham Join Date: Mar 2020 Location: India Posts: 44 Rep Power: 5 Btw, Froude no. is used in determining whether the flow is sub critical or supercritical to allow fr the hydraulic jump. If instead there is a moving wall instead of inlet, the movement of the wall should induce the critical Froude no so that hydraulic jump occurs. I did simulations using interDyMFoam but the interface height is not what is shown by shallow water considerations of the NS eq. Any tips. I though the topic actually matches a lot.

December 26, 2022, 05:54
#13
Member

Join Date: Apr 2015
Posts: 32
Rep Power: 9
Hi, I am trying to simulate a hydraulic jump in an open channel with a Froude number of 13.4 and a Reynolds number of 10^5 with interisofoam with k-e turbulence model and high turbulence intensity of 12%, these are my initial and boundary conditions in zip file, but the jump does not form properly and failed but the simulation didn't diverge as it is in the movie. the mesh quality is good, I also reduced the Courant number and changed B.c conditions and solver to interfome but didn't help me either. Do you have any suggestions for me?
thank you
Attached Images
 Untitled1.jpg (30.7 KB, 7 views) Untitled3.jpg (32.8 KB, 8 views) Untitled2.jpg (34.3 KB, 5 views) Untitled4.jpg (31.6 KB, 4 views)
Attached Files
 New folder.zip (2.6 KB, 9 views)

Last edited by luccy; December 26, 2022 at 17:47.

 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 arshiya4 Fluent UDF and Scheme Programming 3 March 6, 2012 18:13 Gildeh OpenFOAM Running, Solving & CFD 0 February 6, 2012 15:39 shoumo_30 FLUENT 1 December 1, 2011 02:48 Matthew Roberts FLUENT 6 July 31, 2009 12:52 yan FLUENT 1 May 21, 2005 00:21

All times are GMT -4. The time now is 11:47.

 Contact Us - CFD Online - Privacy Statement - Top