# 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, 15:42
Problem modeling hydraulic jump in an open channel with interFoam (VoF)
#1
New Member

Arnau
Join Date: Jan 2012
Posts: 17
Rep Power: 5
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, 91 views)

 November 6, 2013, 07:24 #2 Senior Member   Join Date: Aug 2010 Location: Groningen, The Netherlands Posts: 216 Rep Power: 9 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, 07:39 #3 Senior Member   Albrecht vBoetticher Join Date: Aug 2010 Location: Zürich, Swizerland Posts: 178 Rep Power: 6 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 7, 2013, 00:00 #4 Senior Member   Matthew Denno Join Date: Feb 2010 Posts: 137 Rep Power: 7 vonboett, Could you provide a link or reference for the flume experiments comparison? Thanks, Matt

 November 7, 2013, 15:51 #5 Senior Member   Albrecht vBoetticher Join Date: Aug 2010 Location: Zürich, Swizerland Posts: 178 Rep Power: 6 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, 05:20 Trying Arnau case with kflora suggested bc #6 New Member   Antonio Join Date: Jan 2013 Posts: 8 Rep Power: 4 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, 06:51 #7 Senior Member   Albrecht vBoetticher Join Date: Aug 2010 Location: Zürich, Swizerland Posts: 178 Rep Power: 6 to get a solution comparable to experiments you should go for three dimensions and maybe use LES

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

Antonio
Join Date: Jan 2013
Posts: 8
Rep Power: 4
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.

 Thread Tools 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 On 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 19:13 Gildeh OpenFOAM Running, Solving & CFD 0 February 6, 2012 16:39 shoumo_30 FLUENT 1 December 1, 2011 03: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 18:03.

 Contact Us - CFD Online - Top