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

Pipe Junction with interFoam Pressure Bounday Fails

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 16, 2012, 12:58
Default Pipe Junction with interFoam Pressure Bounday Fails
  #1
Senior Member
 
Matthew Denno
Join Date: Feb 2010
Posts: 138
Rep Power: 16
mgdenno is on a distinguished road
Hello,

I am trying to model the flow of water through a pipe junction box that will split the flow of water between two outlet pipes. The run appears to run smoothly for about 7.2 seconds when suddenly the timestep starts to decrease and eventually the simulation stops. Here is an overview of the setup showing the pipe junction and inlet and outlet locations:



What makes this different than anything I have done with OF previously is that the outlets are submerged and as such I believe that I need to specify a pressure (based on p-rgh) and alpha1 = 1, however, this seems like it may be contributing to the problem. If I don't specify alpha1 = 1 the outlet pipes just fill with air, which is not correct. I have tried quite a few different combinations of boundary conditions at both the inlet and outlets, but the most recent is:

Inlet:
alpha1 = fixedValue, value uniform 1
p_rgh = buoyantPressure
U = flowRateInletOutletVelocity, flowRate constant 1.61

Atmosphere:
alpha1 = inletOutlet, inletValue uniform 0
p_rgh = totalPressure, p0 uniform 0
U = pressureInletOutletVelocity, uniform (0 0 0)

Outlet 1 and 2:
alpha1 = fixedValue, value uniform 1
p_rgh = totalPressure, p0 uniform 17294
U = pressureInletOutletVelocity, uniform (0 0 0)

The actual alpha1, p_rgh and U files are linked below. They are kind of a mess because of the number of different setups that I have tried.
alpha1
p_rgh
U

Based on these boundary conditions it runs for about 7.2 seconds then crashes. If I change the Outlet BC as follows it runs fine but of course it is not correct for my case:

Edit1: These BC run without crashing but there is no flow out of the outlets, so I guess it is not ok.

Outlet 1 and 2:
alpha1 = inletOutlet, inletValue uniform 0
p_rgh = inletOutlet, inletValue uniform 0
U = inletOutlet, inletValue uniform (0 0 0)

My fvSchemes and fvSolutions are as follows:
fvSchemes
fvSolution

My mesh was generated using snappyHexMesh based on a couple of STL files generated using salome. CheckMesh says that the mesh is OK.

I am using the interFoam solver version 2.1.0 for this case. I am currently just running laminar flow but am planning to include turbulence at some point after I get it running.

Using setFields I set the initial conditions to have the pipes be full but not have water over the splitter as shown in the t=0 image:



After 7.2 seconds, just before it fails, it looks like this:



Anyone have ideas on how I can maintain the pressure outlet condition but prevent it from crashing?

I have searched the forums extensively; it looks like a number of people are having issues with pressure boundaries in interFoam. It has been suggested that there is an issue with PIMPLE and that the PISO works better. This is one solution that I am pursuing by compiling and using OF 1.7.1, however, I have not been able get it compiled (I am working on this issue in different thread) That said, if there is a solution that doesn't involve PISO I would be preferable.

Thanks,

MD

Last edited by mgdenno; May 16, 2012 at 21:21.
mgdenno is offline   Reply With Quote

Old   May 17, 2012, 00:10
Default
  #2
Senior Member
 
Matthew Denno
Join Date: Feb 2010
Posts: 138
Rep Power: 16
mgdenno is on a distinguished road
A little bit of new information. I ran it again with the following outlet boundaries, however, instead if initializing the outlet pipes as full at t=0 I started them as empty and only the inlet pipe full.

Outlet 1 and 2:
alpha1 = inletOutlet, inletValue uniform 0
p_rgh = inletOutlet, inletValue uniform 0
U = inletOutlet, inletValue uniform (0 0 0)

The result is that as soon as the water reaches the outlet it crashes. It seems like water (alpha1=1) is not being allowed to pass through the outlet boundary.

Picture attached showing outlet2 just before it crashes.
Attached Images
File Type: jpg t=3.5.jpg (36.4 KB, 40 views)
mgdenno is offline   Reply With Quote

Old   May 17, 2012, 03:42
Default
  #3
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19
Phicau is on a distinguished road
Hi

try this:

p_rgh
Code:
    outlet 
    {
        type            totalPressure;
        U               U;
        phi             phi;
        rho             rho;
        psi             none;
        gamma           1;
        p0              uniform 0;
        value           uniform 0;
    }
U
Code:
    outlet
    {
        type            pressureInletOutletVelocity;
        value           uniform (0 0 0);
    }
alpha1
Code:
    outlet
    {
        type            inletOutlet;
        inletValue      uniform 0;
        value           uniform 0;
    }
It should work now.
Phicau is offline   Reply With Quote

Old   May 17, 2012, 12:50
Default
  #4
Senior Member
 
Matthew Denno
Join Date: Feb 2010
Posts: 138
Rep Power: 16
mgdenno is on a distinguished road
Phicau,

Thank you for your suggestion. It seems to work well for the free outflow case.

Any BC suggestions for simulating the condition where the boundary is submerged?

MD
mgdenno is offline   Reply With Quote

Old   May 17, 2012, 13:15
Default
  #5
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19
Phicau is on a distinguished road
Hi

I think those will work too. Did it crush with them?
Phicau is offline   Reply With Quote

Old   May 17, 2012, 13:35
Default
  #6
Senior Member
 
Matthew Denno
Join Date: Feb 2010
Posts: 138
Rep Power: 16
mgdenno is on a distinguished road
I believe that when I increased the totalPressure to account for the depth of water at the outlet it crashed. I think I also had alpha1 set to a fixedValue = 1, which may have been the problem. I will try it again and this time only change the totalPressure, and will report back the results.

Maybe for alpha1 I should set the inletValue to be = 1 such that any inflow would be "water" and not "air". Thoughts?
mgdenno is offline   Reply With Quote

Old   May 17, 2012, 21:59
Default
  #7
Senior Member
 
Matthew Denno
Join Date: Feb 2010
Posts: 138
Rep Power: 16
mgdenno is on a distinguished road
Phicau,

Thanks for your suggestions. It is running and the results look reasonable so far. I really thought that I had tried those BC but I guess I hadn't with the current mesh. I also set the alpha1 inletOutlet inletValue=1 for the outlet1 and outlet2. I am not sure if this was needed, but I think I might have had a problem with the outlet pulling in air before. If I have a chance I will run it with the inletValue set to 0 and see what happens.

MD
mgdenno is offline   Reply With Quote

Old   May 18, 2012, 04:43
Default
  #8
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19
Phicau is on a distinguished road
Sorry that I copied the atmosphere settings without realizing that yours were pipes.

That sounds reasonable, inletValue=1 will introduce water if in-flux velocities appear. I guess this is not the expected behaviour, but can certainly happen specially in the first moments.

Regards
Phicau is offline   Reply With Quote

Old   April 5, 2013, 10:39
Default
  #9
Member
 
Join Date: Mar 2013
Posts: 98
Rep Power: 13
giack is on a distinguished road
Hi to all,
I have a problem with the setting of the boundary condition in a similar problem.
My case is summarized as follow:
at the initial time t=0 the right wall of a closed tube completely filled with liquid is removed, allowing the liquid to exit the domain and, at the same time, allowing the gas to enter.
I set the BC in this way:
left wall: "p" zero gradient "U" fixed value (0,0,0) "alpha1" fixed value 1
pipe wall: "p" zero gradient "U" fixed value (0,0,0) "alpha1" zero gradient
right wall (outlet): "p" total pressure "U" pressureInletOutletVelocity "alpha1" inletOutlet
but I obtain unphysical result

Where is the error?Someone have an idea to set BC for this case?
thank to all
giack is offline   Reply With Quote

Reply

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Pipe Flow - Pressure Drop Daniel L FLOW-3D 2 December 10, 2010 05:23
My Revised "Time Vs Energy" Article For Review Abhi Main CFD Forum 2 July 9, 2002 10:08
Terrible Mistake In Fluid Dynamics History Abhi Main CFD Forum 12 July 8, 2002 10:11
fluid flow fundas ram Main CFD Forum 5 June 17, 2000 22:31
Hydrostatic pressure in 2-phase flow modeling (long) DS & HB Main CFD Forum 0 January 8, 2000 16:00


All times are GMT -4. The time now is 13:37.