|
[Sponsors] |
September 4, 2017, 06:41 |
Volume fraction at outlet not known
|
#1 |
Member
Kamal Bisht
Join Date: Jun 2015
Location: Germany
Posts: 57
Rep Power: 11 |
Hello,
I am simulating a two phase (Gas & liquid) problem (see figure). Problem statement: Liquid is flowing in a rectangular channel which has an opening at the top. The flow is driven by a motor at the outlet. As outlet flow rate increases, pressure in the channel drops, subsequently the free surface of the liquid at the opening collapsed and ingested air in the from of bubbles. The frequency of the ingested bubbles depends on the outlet flow rate. The known boundary conditions are the pressure at the inlet recorded with sensor and the outlet mass flow rate. I am using pressure inlet and velocity outlet BC as mass flow outlet is not supported with VOF. At velocity boundary, I need to define the Volume Fraction of the two phases which is not known in advance for air (bubble). When I defined vf of liquid as 1.0 at outlet, the simulation works fine until the bubble reached to the outlet boundary, where it couldn't pass through the outlet due to volume fraction of air defined as zero. I would like to know how to proceed in this case, the vf of air is unknown. Am I using the correct BC for the problem? |
|
September 4, 2017, 06:48 |
|
#2 |
Senior Member
Cees Haringa
Join Date: May 2013
Location: Delft
Posts: 607
Rep Power: 0 |
why aren't you using a pressure outlet b.c.? That doesn't discriminate between phases and allows both gas and liquid to flow out.
Edit: I see the rationale between your choice of b.c.; however, if you do know the inlet volume fractions, then it may be easier to use velocity inlet and pressure outlet - the net flow is the same at in and outlet due to conservation of mass anyway. |
|
September 4, 2017, 07:45 |
|
#3 |
Member
Kamal Bisht
Join Date: Jun 2015
Location: Germany
Posts: 57
Rep Power: 11 |
Hi Cees,
At inlet the volume fraction of liquid is always 1.0, the pressure is also always constant regulated through some other means. The mass flow rate at the outlet is known which is controlled by a pump placed in the downstream. The problem is that air ingested through an opening at the top when the flow rate is higher than a specific value. As you can see I don't know the bubble size or the frequency of the bubble ingestion beforehand in this case. So, at inlet its always liquid and at outlet its liquid with intermittent bubbles. |
|
September 4, 2017, 08:04 |
|
#4 |
Senior Member
Cees Haringa
Join Date: May 2013
Location: Delft
Posts: 607
Rep Power: 0 |
But currently there is no opening at the top, it's a symmetry. Should't there be a pressure inlet condition or so at the top to draw air from?
|
|
September 4, 2017, 08:24 |
|
#5 |
Member
Kamal Bisht
Join Date: Jun 2015
Location: Germany
Posts: 57
Rep Power: 11 |
Ah sorry for the vague representation of the image. While initializing the flow field I patched the top box with VF of air as 1 (image attached)
|
|
September 4, 2017, 08:55 |
|
#6 |
Senior Member
Cees Haringa
Join Date: May 2013
Location: Delft
Posts: 607
Rep Power: 0 |
I see, but how is the air being drawn from the top chamber? Is it being replaced by water somehow then? Or being replenished by air from another location - in the latter case, you will need an inlet at the top too, right?
|
|
September 4, 2017, 09:08 |
|
#7 |
Member
Kamal Bisht
Join Date: Jun 2015
Location: Germany
Posts: 57
Rep Power: 11 |
That is my other concern, which I am currently investigating. But yeah you are right the top should be a Pressure inlet to replicate atmosphere.
However, how to handle the bubble at outlet? The mass flow rate should be satisfied at outlet, which can be done by Velocity or Mass flow Boundaries. The current situation is that velocity boundary can handle only one phase, Mass flow outlet is not available with VOF. Can it be done by using mass flow inlet at outlet? |
|
September 4, 2017, 10:43 |
|
#8 |
Senior Member
Cees Haringa
Join Date: May 2013
Location: Delft
Posts: 607
Rep Power: 0 |
You are dealing with a rather unique system in your conditions; I don't suppose you can consider the liquid inflow to be constant, right - it is really affected by what happens at the outlet? From what I gather, it seems you probably need to look into a custom UDF boundary condition, that maintains a constant total flux through the boundary with a variable composition, based on the volume fractions of the neighboring cells.
|
|
September 5, 2017, 03:38 |
|
#9 |
Member
Kamal Bisht
Join Date: Jun 2015
Location: Germany
Posts: 57
Rep Power: 11 |
Yeah inflow depends on the outlet condition. I will look into the UDF manual to deal with this issue.
Thanks for your help. |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
How to use a UDF to set the volume fraction in the cells next to a wall? | DF15 | Fluent UDF and Scheme Programming | 33 | August 20, 2020 14:36 |
Is it possible that multi-phase flow shows gradation of volume fraction? | kjh9537 | Fluent UDF and Scheme Programming | 0 | June 23, 2017 15:16 |
[blockMesh] BlockMesh FOAM warning | gaottino | OpenFOAM Meshing & Mesh Conversion | 7 | July 19, 2010 15:11 |
interDyMFoam - change in volume fraction | gopala | OpenFOAM Running, Solving & CFD | 0 | April 27, 2009 11:46 |
the Eulerian model and the lower volume fraction | hx li | FLUENT | 2 | October 27, 2005 07:17 |