CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Pre-Processing (http://www.cfd-online.com/Forums/openfoam-pre-processing/)
-   -   Pressure driven flow in vertical pipe (http://www.cfd-online.com/Forums/openfoam-pre-processing/123028-pressure-driven-flow-vertical-pipe.html)

Hale September 3, 2013 04:00

Pressure driven flow in vertical pipe
 
1 Attachment(s)
Hi to all,

I have a tank connected to a vertical pipe. The tank is full of water and the sides of the tank have the following boundary condition

Inlet
U: zeroGradient;
p: zeroGradient;
alpha1: fixedValue uniform 1;

The water in the tank would then fall into the vertical pipe. The outlet of the pipe has the following boundary condition

outlet
U: zeroGradient;
p: fixedValue uniform -250;
alpha1: zeroGradient;

The upper side of the tank is then open to the atmosphere with following boundary condition

atmosphere
U: inletOutlet;
p: totalPressure (p0=0);
alpha1: zeroGradient;

In this way I want the water to fall down based on the pressure difference between outlet and inlet. But with the mentioned boundary conditions the water moves upward and the simulation crashes after some seconds. I guess that the hydrostatic pressure at the outlet causes problems, since there is a positiv hydrostatic pressure corresponding to the water in the pipe.

Am I choosing wrong boundary conditions?

I really appreciate any help and idea.

NB. I'm using interFoam.

nimasam September 4, 2013 13:47

well whats your OpenFOAM version ?

Hale September 4, 2013 17:44

It's version 2.0.1

nimasam September 5, 2013 00:53

then you should assign p_rgh not p! p will be calculated then ;) through solving problem whats your p_rgh configurations?

Hale September 5, 2013 04:15

Quote:

Originally Posted by nimasam (Post 449980)
then you should assign p_rgh not p! p will be calculated then ;) through solving problem whats your p_rgh configurations?

Sorry I didn't write that in the discribtion but I am actually assigning p_rgh and not p.

By the way do you know how to change the values for Sp and Su. In the page below, post #7 decribes that the problem with unbounded alpha1 can be solved by changing Sp and Su but I have no idea how and where to change them. Do you know that?

http://www.cfd-online.com/Forums/ope...tml#post449990


Thanks for your time

nimasam September 5, 2013 09:17

Quote:

Originally Posted by Hale (Post 450007)
Sorry I didn't write that in the discribtion but I am actually assigning p_rgh and not p.

why do u assign p_rgh a negative value?
Quote:

By the way do you know how to change the values for Sp and Su. In the page below, post #7 decribes that the problem with unbounded alpha1 can be solved by changing Sp and Su but I have no idea how and where to change them. Do you know that?

http://www.cfd-online.com/Forums/ope...tml#post449990


Thanks for your time
well, you should change the source code, look at alphaEqn.H, in
Code:

solvers/multiphase/interFoam
and
solvers/multiphase/compressibleInterFoam

then compile your code ;), however i think atleast now you should get sure about your BC ;) it can be future step

Hale September 13, 2013 08:16

Thanks a lot Nima. I have solved the problem.


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