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

Multiphase interFoam pipeflow One inlet Two outlets

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 20, 2020, 02:22
Default Multiphase interFoam pipeflow One inlet Two outlets
  #1
Member
 
Munaf
Join Date: May 2019
Posts: 38
Rep Power: 7
mwmalkawi is on a distinguished road
I have a multiphase pipe flow case with one inlet (inlet1), two outlets (outlet1 and outlet2) and one vent shaft (atmosphere) as per screen shot below attached with turbulent flow k-e model enabled
When run the solution and simulate results in preview it shows the flow when reaches outlet 1 and 2 it starts to reverse as if there is wall at these outlets that is preventing the flow to exit you can see the simulation video attached.
3D Model.png

flow at outlets.png

reverse flow outlet1 .png

reverse flow outlet2 .jpg
I really appreciate the help in this one as I believe something wrong with the settings of my boundary conditions as I have tried millions of combinations and options for BC for the last 6motnhs.

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object alpha.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 0 0 0 0 0];

internalField uniform 0;

boundaryField
{

wall
{
type zeroGradient;
}


outlet1
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}

outlet2
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}

inlet1
{
type inletOutlet;
inletValue uniform 1;
value uniform 1;
}


}

// ************************************************** *********************** //

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p_rgh;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -2 0 0]; // or [0 2 -2 0 0 0 0];
internalField uniform 0;
boundaryField
{

wall
{
type zeroGradient;
}

outlet1
{
type totalPressure;
p0 uniform 0;
value uniform 0;
}

outlet2
{
type totalPressure;
p0 uniform 0;
value uniform 0;
}

inlet1
{
type fixedFluxPressure;
value uniform 0;
}

}
// ************************************************** *********************** //

FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "0";
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);
boundaryField
{
wall
{
type fixedValue;
value uniform (0 0 0);
}

inlet1
{
type surfaceNormalFixedValue;
refValue uniform -1;
}

outlet1
{
type pressureInletOutletVelocity;
value uniform (0 0 0);
}
outlet2
{
type pressureInletOutletVelocity;
value uniform (0 0 0);
}

}

// ************************************************** *********************** //

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -3 0 0 0 0];

internalField uniform 0.0027;
boundaryField
{
wall
{
type epsilonWallFunction;
value uniform 0.00027;
}

outlet1
{
type zeroGradient;
}
outlet2
{
type zeroGradient;
}


inlet1
{
type fixedValue;
value uniform 0.025;
}


}

// ************************************************** *********************** //

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];

internalField uniform 0.009;
boundaryField
{
wall
{
type kqRWallFunction;
value uniform 0.0052;
}

outlet1
{
type zeroGradient;
}
outlet2
{
type zeroGradient;
}


inlet1
{
type fixedValue;
value uniform 0.002;
}
}

// ************************************************** *********************** //
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -1 0 0 0 0];

internalField uniform 0;
boundaryField
{
wall
{
type nutkWallFunction;
value uniform 0;

outlet1
{
type calculated;
value uniform 0;
}
outlet2
{
type calculated;
value uniform 0;
}


inlet1
{
type calculated;
value uniform 0;
}


}

// ************************************************** *********************** //
mwmalkawi is offline   Reply With Quote

Old   January 20, 2020, 04:23
Default
  #2
Senior Member
 
Join Date: Aug 2014
Location: UK
Posts: 213
Rep Power: 12
fresty is on a distinguished road
You may want to link/attach your case folder to evoke more specific support.
fresty is offline   Reply With Quote

Old   January 20, 2020, 04:37
Default
  #3
Member
 
ssa
Join Date: Sep 2018
Posts: 93
Rep Power: 7
ssa_cfd is on a distinguished road
why did you use inletOutlet for the inlet patch in alpha.water and totalPressure for outlet patch in p_rgh.





check the water channel tutorial and modify your boundary conditions.
ssa_cfd is offline   Reply With Quote

Old   January 20, 2020, 07:32
Default
  #4
Member
 
Munaf
Join Date: May 2019
Posts: 38
Rep Power: 7
mwmalkawi is on a distinguished road
Quote:
Originally Posted by ssa_cfd View Post
why did you use inletOutlet for the inlet patch in alpha.water and totalPressure for outlet patch in p_rgh.





check the water channel tutorial and modify your boundary conditions.
I have done all possible cases you would imagine i am attaching to you an excel sheet showing you all cases of BC i have used.
Attached Files
File Type: xlsx interFoam case files BC settings.xlsx (152.3 KB, 38 views)
mwmalkawi is offline   Reply With Quote

Old   January 20, 2020, 08:14
Default
  #5
Member
 
Munaf
Join Date: May 2019
Posts: 38
Rep Power: 7
mwmalkawi is on a distinguished road
Thank you Fresty please see below the zip file for the Case i appreciate anyone who can asist
Attached Files
File Type: zip Case Files.zip (51.2 KB, 25 views)
mwmalkawi is offline   Reply With Quote

Old   January 20, 2020, 11:06
Default Velocity file mistake
  #6
New Member
 
Nope nope nope
Join Date: Nov 2019
Posts: 3
Rep Power: 6
Ferran is on a distinguished road
I'm pretty sure it is because in the U file in the 0 folder, you have set the outlets as "type pressureInletOutletVelocity; value uniform (0 0 0);" so it is acting with velocity=0. Try setting it as a "type zeroGradient;" or another outlet boundary condition.
Let me know if it works, otherwise I could share my boundary conditions settings, I have been running some multiphase simulations.
On the case files, the mesh or the dictionaris to generate the mesh are missing for us to try your case.
Ferran is offline   Reply With Quote

Old   January 20, 2020, 11:30
Default
  #7
Member
 
Munaf
Join Date: May 2019
Posts: 38
Rep Power: 7
mwmalkawi is on a distinguished road
Thank you Ferran will you please share yours if you have a similar pipe flow case with one inlet and two outlets...many thanks
mwmalkawi is offline   Reply With Quote

Old   January 21, 2020, 07:47
Default Example with simple geometry
  #8
New Member
 
Nope nope nope
Join Date: Nov 2019
Posts: 3
Rep Power: 6
Ferran is on a distinguished road
I have run it with a simple geometry, 1 water inlet and 2 pressure outlets. There is a readme listing all the changes I've done with your files as a template.
I do not get the recirculation that you mention, though I have not seen the results of your case.

It took around one hour to calculate with 4 processors, 17.3s of real time.

Without results (~1 MB)
https://drive.google.com/file/d/1IWa...ew?usp=sharing

With results and animation (~1 GB)
https://drive.google.com/file/d/1_D4...ew?usp=sharing

You should be able to download, otherwise tell me.

Hope it is useful.
Ferran is offline   Reply With Quote

Old   January 29, 2020, 18:46
Default
  #9
Senior Member
 
Join Date: Jul 2019
Posts: 148
Rep Power: 6
Bodo1993 is on a distinguished road
Quote:
Originally Posted by Ferran View Post
I have run it with a simple geometry, 1 water inlet and 2 pressure outlets. There is a readme listing all the changes I've done with your files as a template.
I do not get the recirculation that you mention, though I have not seen the results of your case.

It took around one hour to calculate with 4 processors, 17.3s of real time.

Without results (~1 MB)
https://drive.google.com/file/d/1IWa...ew?usp=sharing

With results and animation (~1 GB)
https://drive.google.com/file/d/1_D4...ew?usp=sharing

You should be able to download, otherwise tell me.

Hope it is useful.

Hi Ferran,

Kindly, I would appreciate your assistance. I have simple 2D vertical channel flow with single inlet and single outlet. Fluid 1 initially fills the channel and Fluid 2 is pumped from below (fluids have different densities and viscosities). I have an issue at the outlet - the fluids exists from just a small portion of the outlet boundary.

For the outlet boundary conditions, I use: pressureInletOutletVelocity for U, totalPressure for P and zeroGradient/inletOutlet for alpha.

I am wondering if you can assist me on this case. Thanks and I look forward to hearing from you.
Bodo1993 is offline   Reply With Quote

Old   January 29, 2020, 22:37
Default
  #10
New Member
 
Join Date: Nov 2019
Posts: 24
Rep Power: 6
Nic86 is on a distinguished road
hello everyone,
sorry for the intrusion but could anyone try to answer to my thread?

InterFoam simulation - mesh setup

It's a quite similar case to the one of this post but my problem is that the water reaches the atmosphere patch and i do not know how to solve it.

Thank you
Nic86 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
[IHFOAM] The IHFOAM Thread Phicau OpenFOAM Community Contributions 392 September 8, 2023 18:10
interFoam wave propagation and explosion of Courant number and residuals ChiaraViola OpenFOAM Running, Solving & CFD 1 June 26, 2019 05:36
[waves2Foam] How to set speed at Inlet in waves2Foam, like we set in interfoam, U dictionary Nemo_CFDEngineer OpenFOAM Community Contributions 5 January 15, 2017 12:58
One inlet and two outlets mingzhao OpenFOAM Running, Solving & CFD 1 July 16, 2015 18:03
the inlet pressure about the multiphase flow jasonchang CFX 0 May 25, 2006 08:54


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