|
[Sponsors] |
Boundary Condition for an Open Channel in InterFoam |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
January 16, 2014, 04:52 |
Boundary Condition for an Open Channel in InterFoam
|
#1 |
New Member
Join Date: Oct 2012
Posts: 19
Rep Power: 13 |
Hi,
I'm trying to setup an open channel case in Interfoam but I found some problem with the boundary condition. I couldn't find a BC for the inlet that allows water and air to enter in the domain (with a defined free surface level) in a quantity that has to be calculated from the flow rate imposed at outlet. To make it more clear: Domain: you can imagine it as a sink Inlet: water and air with imposed free surface level (I prefer to don't split the boundary in water_Inlet an air_Inlet to make it easier to change the free surface level) Outlet: imposed massflowrate (mainly water with some air ingestion). I tried to setup the inlet as InletOutlet and outlet as a negative flowRateInletVelocity but I had serious problem with convergence. I'll really appreciate any advice. Probably the solution will be straight forward and I don't know it because I'm new to OF. Thank you in advance. |
|
January 16, 2014, 05:24 |
|
#2 |
Senior Member
|
Hi,
To fix water level at the inlet you can set calculated BC for alpha1 (so it will use values of alpha1 set initially with setFields), then you can set for example fixedValue (calculated from inflow rate) for velocity. For outlet you can try variableHeightFlowRate, which will calculate velocity from water level (using alpha1) and discharge rate. Though AFAIR I had problems with turbulence model equations divergence with this BC. |
|
January 16, 2014, 05:39 |
|
#3 |
New Member
Join Date: Oct 2012
Posts: 19
Rep Power: 13 |
Thank you for your reply.
I forgot to write it but I've already set the Inlet alpha1 as calculated. The problem is related to how to impose the flow rate for the air because it will depends on the system solution if there will be air ingestion or not. Another important think that i forgot to clarify is that the inlet has a x normal while the outlet has a z normal and is completely submerged. (as in a sink) |
|
January 20, 2014, 04:06 |
|
#4 |
New Member
Join Date: Oct 2012
Posts: 19
Rep Power: 13 |
don't you have any suggestion?
|
|
January 20, 2014, 07:56 |
|
#5 |
Senior Member
|
It's rather difficult to have suggestions without knowing the whole problem.
From what you've said I can guess that the inlet is something with two phases and you'd like to set the speed of these phases independently. I can not propose a way to do this without whether using two different patches with different BCs, or implementing your own BC, or using codedFixedValue. For the outlet I again can try to guess something though I can't yet imagine what's going on with outlet from your description. So, in short: no, I don't have any |
|
January 20, 2014, 09:06 |
|
#6 | |
New Member
Join Date: Oct 2012
Posts: 19
Rep Power: 13 |
Quote:
|
||
January 20, 2014, 09:18 |
|
#7 |
New Member
Join Date: Oct 2012
Posts: 19
Rep Power: 13 |
I start to think that it can be also a problem related to time step and fvschemes.
Because I have some problem of divergence and because the k start to be bounded since few iterations after the beginning. I'm using this setting for adjustable time step Code:
adjustTimeStep yes; maxCo 0.5;// [max courant number] maxAlphaCo 0.5; // [ maxDeltaT 0.05; // [s] Code:
divSchemes { div(rho*phi,U) Gauss upwind; div(phi,alpha) Gauss vanLeer; div(phirb,alpha) Gauss interfaceCompression; div(phi,k) bounded Gauss upwind; div(phi,omega) bounded Gauss upwind; div((muEff*dev(T(grad(U))))) Gauss linear; } |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Wind turbine simulation | Saturn | CFX | 60 | July 17, 2024 05:45 |
Trouble compiling utilities using source-built OpenFOAM | Artur | OpenFOAM Programming & Development | 14 | October 29, 2013 10:59 |
how to use Up Stream boundary condition in open channel Flow | ms.jafarinik | FLUENT | 0 | April 18, 2010 14:38 |
Open Channel Boundary Conditions via journal | Matteo | FLUENT | 0 | January 21, 2008 11:05 |
Convective Heat Transfer - Heat Exchanger | Mark | CFX | 6 | November 15, 2004 15:55 |