CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

interFoam, cyclic BC and gravity problem

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 26, 2016, 11:14
Exclamation interFoam, cyclic BC and gravity problem
  #1
New Member
 
Join Date: Feb 2016
Posts: 2
Rep Power: 0
mhiel is on a distinguished road
Hi everyone,

I'm trying to simulate turbulent water flow through a pipe (the pipe is half-filled with water, the upper half is air) using interFoam and cyclic boundary conditions with the goal to examine the velocity profile when the flow is fully developed.

The g file is set to:

value ( 0.01 -9.81 0)

to simulate a slight inclination of the pipe (about 1 per mille) which induces a constant potential on the water.
In ParaView I expected to see an acceleration of the water (increasing U) until the flow reaches developed state.
However the results are quite irritating as the U values in the water region are constantly fluctuating on low magnitudes (as can be seen in the attached image) and never reach fixed values. Moreover the velocities are far too slow for the water phase and there is no acceleration of the water phase as anticipated.

I assume I did something wrong with the cyclic boundary conditions or other settings but I'm new to OF and would appreciate any help/ideas.

Here's a bit more information:

my U file:

Code:
internalField   uniform ( 0 0 0 );

boundaryField
{
    inlet
    {
        type            cyclic;
    }
    walls
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    outlet
    {
        type            cyclic;
    }
}
p_rgh:

Code:
internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            cyclic;
    }
    outlet
    {
        type            cyclic;
    }
    walls
    {
        type            fixedFluxPressure;
        value           uniform 0;
    }
}
g:
Code:
value           ( 0.01 -9.81 0 );
I attached the whole case if you'd like to run it yourself.
Note that in the case the orientation of the pipe axes differs from what I've described above (which i did to avoid confusion). Thus in the case g reads: ( 0 0.01 -9.81 )

mhiel
Attached Images
File Type: jpg Screenshot_2016-02-26_15-42-08.jpg (122.1 KB, 84 views)
Attached Files
File Type: zip cyclicWaterAirPipe.zip (11.6 KB, 10 views)

Last edited by mhiel; February 29, 2016 at 11:00.
mhiel is offline   Reply With Quote

Old   February 29, 2016, 09:17
Default
  #2
New Member
 
Join Date: Feb 2016
Posts: 2
Rep Power: 0
mhiel is on a distinguished road
I have already searched the forum for possible solutions.
I found this thread

http://www.cfd-online.com/Forums/ope...ion-issue.html

which discusses a problem with cyclic BC related to g.
There it is suggested to change the code of the interFoam solver, but unfortunately this didn't help in my case.

Edit: I also tested the case with a real inclination of the pipe and g set to:

Code:
value           ( 0 0 -9.81 )
The results are pretty much the same. I'd be glad if anybody could provide some advice on this problem.

Thank you in advance and kind regards,
mhiel
Attached Files
File Type: zip cyclicPipeInclined.zip (11.7 KB, 7 views)

Last edited by mhiel; March 1, 2016 at 05:48.
mhiel is offline   Reply With Quote

Old   June 18, 2017, 08:14
Default
  #3
Member
 
Amir
Join Date: Jan 2017
Posts: 30
Rep Power: 5
albet is on a distinguished road
Hello Foamers,

Does anybody have any solution for this problem?
I tried cyclic boundary condition with interFoam but it doesn't work as mhiel mentioned.
there is a kind of wave moving in the flow.
I am thinking about a long pipe to pass the entrance length for fully developed Turbulent flow, but it doesn't sound reasonable.

I hope somebody can help,
regards,
Amir
albet is offline   Reply With Quote

Old   June 18, 2017, 11:24
Default
  #4
Senior Member
 
piu58's Avatar
 
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 590
Rep Power: 7
piu58 is on a distinguished road
It is physical nonsense to have a cyclic b.c. AND to expect an accellearation. Cyclic means the inlet is equal the outlet which at best simulates a piece of the pipe where the flow is quasi static.
Even in this case one may expect a velocity which fits to the problem. Probably you have to set this velocity as starting condition.

I would experiment with other boundary conditions. If you want to see the full developed flow it may be easier making the pipe longer.
__________________
Uwe Pilz
--
Die der Hauptbewegung überlagerte Schwankungsbewegung ist in ihren Einzelheiten so hoffnungslos kompliziert, daß ihre theoretische Berechnung aussichtslos erscheint. (Hermann Schlichting, 1950)
piu58 is offline   Reply With Quote

Old   July 1, 2019, 18:10
Default
  #5
Member
 
Join Date: Jul 2015
Location: Gainesville,FL
Posts: 93
Rep Power: 6
randolph is on a distinguished road
Here is how I get around this problem.

Centrally what you need is a hybrid between "variableHeightFlowRateInletVelocity" and "mapped" boundary condition.

You want to use the "variableHeightFlowRateInletVelocity" to scale the velocity BC to keep the right flowrate, while the the velocity profile that scaled by the "variableHeightFlowRateInletVelocity" should sample from some distance downstream (take a look at the "mapped" BC).

And run the model after a long time, the depth and velocity should adjust themselves to the equilibrium condition.
randolph is offline   Reply With Quote

Reply

Tags
cyclic boundary condition, gravity, interfoam

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



All times are GMT -4. The time now is 10:30.