|
[Sponsors] |
February 26, 2016, 11:14 |
interFoam, cyclic BC and gravity problem
|
#1 |
New Member
Join Date: Feb 2016
Posts: 2
Rep Power: 0 |
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; } } Code:
internalField uniform 0; boundaryField { inlet { type cyclic; } outlet { type cyclic; } walls { type fixedFluxPressure; value uniform 0; } } Code:
value ( 0.01 -9.81 0 ); 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 Last edited by mhiel; February 29, 2016 at 11:00. |
|
February 29, 2016, 09:17 |
|
#2 |
New Member
Join Date: Feb 2016
Posts: 2
Rep Power: 0 |
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 ) Thank you in advance and kind regards, mhiel Last edited by mhiel; March 1, 2016 at 05:48. |
|
June 18, 2017, 08:14 |
|
#3 |
Member
Amir
Join Date: Jan 2017
Posts: 32
Rep Power: 9 |
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 |
|
June 18, 2017, 11:24 |
|
#4 |
Senior Member
Uwe Pilz
Join Date: Feb 2017
Location: Leipzig, Germany
Posts: 744
Rep Power: 15 |
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) |
|
July 1, 2019, 18:10 |
|
#5 |
Senior Member
Reviewer #2
Join Date: Jul 2015
Location: Knoxville, TN
Posts: 141
Rep Power: 11 |
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. |
|
February 9, 2023, 11:12 |
|
#6 |
New Member
Join Date: Jun 2016
Posts: 7
Rep Power: 10 |
@randolph - Having tried various approaches, your suggestion seems promising. Do you have an example 0/ and boundary file setup? Thanks
Last edited by tias; February 9, 2023 at 11:13. Reason: reply to a member |
|
Tags |
cyclic boundary condition, gravity, interfoam |
Thread Tools | Search this Thread |
Display Modes | |
|
|