Dear vonboett,
Thanks for your replay. In fact, I was running my case in the opposite way where I impose the inlet BC from the bottom at z=0 and flow should go out at the top where z=0.017 m. So when I shrink the domain length to half so less, liquid comes out smoothly not like when z=0.017. Any how the case now is working for the long domain when I swamped the BC upside down which I do not why???. If you have any explanation you can help me with because it is not logical... Could you please show me what was the inlet and walls BC you used for the case you ran and sent me?
best wishes,
Sandy,

Hi guys;

I presently have some sort of b.c issues that you had previously. Hope anyone can help me out.:)

I am using interfoam solver and trying to inject two fluids simulataneously through alternate channels from the inlet {water, oil, water, oil} and expect a mixed flow{water+oil} at the outlet boundary as shown in figure.

I fix pressure[Pa] at the inlet(x) and at the outlet(0).

For velocity i use a "zeroGradient" at outlet and "pressureInletOutletVelocity" for inlet.

When i have lower pressure injection for water relative to oil, everything seems to be fine as expected. But when i try to make all the inlet channels have the same pressure values, then
initially everything looks fine but after some time of run, i have a backward flow at one of the inlet{that is the flow is going out of the domain at the inlet} and alpha gets out of range. Also the flow at the outlet changes it direction of flow. [Can be seen in the posted data. Negative sign at inlet is flow into domain, positive sign at outlet is outflow through the domain. Above one is something at initial time and the bottom data is at a later time].

Expression totalInflow_0 : sum=-3.92639e-15
Expression totalInflow_1 : sum=-7.33046e-16
Expression totalInflow_2 : sum=-1.85515e-17
Expression totalInflow_3 : sum=-1.06606e-15
Expression inflow_0_Water : sum=-3.92639e-15
Expression inflow_1_Water : sum=0
Expression inflow_2_Water : sum=-1.85515e-17
Expression inflow_3_Water : sum=0
Expression inflow_0_Oil : sum=0
Expression inflow_1_Oil : sum=-7.33046e-16
Expression inflow_2_Oil : sum=0
Expression inflow_3_Oil : sum=-1.06606e-15
Expression outflow_0_Water : sum=2.93572e-15
Expression outflow_1_water : sum=2.28004e-16
Expression outflow_2_Water : sum=1.86927e-17
Expression outflow_3_water : sum=4.26874e-22
Expression totalOutflow_0 : sum=2.96816e-15
Expression totalOutflow_1 : sum=5.51016e-16
Expression totalOutflow_2 : sum=1.05719e-15
Expression totalOutflow_3 : sum=1.16767e-15
Courant Number mean: 0.0292251 max: 0.300384
Interface Courant Number mean: 0.0116064 max: 0.23988


---------------------------
Expression totalInflow_0 : sum=-1.11676e-13
Expression totalInflow_1 : sum=2.21668e-13
Expression totalInflow_2 : sum=-5.45042e-14
Expression totalInflow_3 : sum=-1.20261e-14
Expression inflow_0_Water : sum=-1.11676e-13
Expression inflow_1_Water : sum=0
Expression inflow_2_Water : sum=-5.45042e-14
Expression inflow_3_Water : sum=0
Expression inflow_0_Oil : sum=0
Expression inflow_1_Oil : sum=2.21668e-13
Expression inflow_2_Oil : sum=0
Expression inflow_3_Oil : sum=-1.20261e-14
Expression outflow_0_Water : sum=-1.9006e-14
Expression outflow_1_water : sum=-9.30047e-15
Expression outflow_2_Water : sum=-4.94652e-16
Expression outflow_3_water : sum=-2.04992e-19
Expression totalOutflow_0 : sum=-1.90947e-14
Expression totalOutflow_1 : sum=-1.34784e-14
Expression totalOutflow_2 : sum=-6.41545e-15
Expression totalOutflow_3 : sum=-4.47371e-15
Courant Number mean: 0.00731049 max: 0.301158
Interface Courant Number mean: 0.00405459 max: 0.195399


Though i guess symmetry of my structure{spreads out the pressure} and the viscosity values(oil>water) have a role to play in the displacement i am not able to understand why there is an outflow through the inlet with positive pressure(gradient) and inflow at the region of outlet.

Thanks;
Saideep

Do you have gravity acting in -y direction? could you show the pressure distribution along a center slice cut through the domain?

Hi vonboett;

Thanks for your reply.

The gravity has been turned off.

I taught rather than plots the attached pics make sense in following order. Just let me know if you find it the other way round.

I was using an injection pressure: 10000Pa {all inlets} and outlet Pressure: 0Pa.

You can notice a negative Ux that starts to occur at the 1st inlet causing an outflux. [numbering 0 to 3 from bottom to top] see fig1 and fig2 which are successive time steps data.
Attachment 42532 Attachment 42533

Also attached is a plot at the outlet end where the velocities turn negative which start to cause an inflow at the outlet region[fig3].
Attachment 42534

Coming to pressure, it looks quite physical initially[fig4] {last time step before it starts to behave bad}:
Attachment 42535
and at the time where it bursts out[fig5] In the last figure please look at the legend values.
Attachment 42536

Overall it makes sense of the flow behavior but why is there a change in the direction of flow and sudden pressure burst throughout the domain?
I notice that this happens when the water injection pressure >= oil injection pressure but not the other way round.:mad:

Sorry for the long post;
Thanks;
Saideep

If you have multiple fluids at the inlet i guess that setting a fixed value boundary condition for pressure in both fluids is not a good idea. The pressure inside the two fluids is, in theory, different due to capillary pressure but you are trying to fix the same value. I think this is the problem. I have observed similar behavior at the outlet when both fluids leave the domain at the same time.
I would switch to a velocity fix boundary condition and leave the pressure free to adjust depending on the fluid at the inlet. Or fix the pressure only at the inlet where you have fluid 1 (or fluid 2).

Hope this help

Andrea

what is your alpha boundary condition at the inlet and outlet? Could you take a look if you get negative concntrations or values grater 1 for alpha?

Try out
outlet
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}
for alpha1 together with a pressureInletOutletVelocity and zeroGradient for U at inlet?

Thanks Andrea and vonboett;

@Andrea:
Actually i started with fixed velocity values for the inlets which is an ideal case. Due to the shape of the model{throats arent symmetric at the ends of the model} I am not able to maintain a constant pressure at inlets of specific phase injection channels{which i am interested in, to calculate the grad(P)_oil/water}, so i switched to fixed pressure boundaries. May be an alternative could be is to extend my model{i.e. include throats at the top and bottom also so that it could symmetric at all inlet channels}.

@vonboett:
I am using the following bc for alpha.
left_0
{
type fixedValue;
value uniform 1; /*phase1= water*/
}

left_1
{
type fixedValue;
value uniform 0; /*phase0 = oil*/
} /*continue to other inlets also*/
and outlet channels with zeroGradient and equilibrium contact angle over the surface of the solid.

Do you mean inletOutlet and outletInlet for U at inlet and P at oulet respectively or for alpha?

And yes once this happens nothing remains to be computed accurately in terms of volume fraction, alpha[increases every step of calculation c.a. 4000], pressure etc.

Thanks a lot Adrea and vonboett for your insight on this;:)
Saideep

I am curious about what happends with pressureInletOutletVelocity for U at the outlets
and
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}
for alpha at the outlets...

Hi and sorry for my late reply.

It practically made not difference.
Data written as output:

/*-ve sign is inflow into the domain*/
Expression totalInflow_0 : sum=-1.92084e-13
Expression totalInflow_1 : sum=4.19385e-13
Expression totalInflow_2 : sum=-1.78747e-13
Expression totalInflow_3 : sum=-1.6061e-14

/*+ve shall be the outflow through the outlet*/
Expression totalOutflow_0 : sum=-6.52957e-15
Expression totalOutflow_1 : sum=-1.24967e-14
Expression totalOutflow_2 : sum=-7.9281e-15
Expression totalOutflow_3 : sum=-5.53935e-15

Max(alpha1) = 176.701

I used a fixed P at inlet for water at 10500Pa and for oil at 10000 Pa.
Inlet and outlet velocities mentioned as pressureInletOutletVelocity.
Also changed the alpha outlet bc to inletOutlet{I still dont understand why you meant to use this over here but made no difference}!!

Any further ideas.
Thanks!!
Saideep

Hi Saideep,

I work with turbulent open channel flows and there it can happen that you get some inflow at the outlet in some parts for some moment, and therefore I like to use inletOutlet for alpha to make sure that it sucks in atmospheric air and not water from the downstream. If you use
U
{ type inletOutlet
inlet value uniform (0 0 0)
value uniform (0 0 0)
}
at the outlet you should at least get rid of inflow at the outlet. Just as a test to see that everything else is right, how about specifying the inflow velocity at the inlets and use p_rgh zeroGradient at the outlet (as well as alpha zeroGradient)?

Cheers,

Albrecht

Hi, Albrecht

I was trying for the past couple months to figure up how to control my surface water level , and I think I manage to do that. I am simulating flow in river past plate positioned inside the flow domain. My goal is to do LES but for know I am using RANS with interFoam with Re=50000. I got the result , however the flow fields don't looks ok, and also when I compare the results with fluent the drag signature is not correct. The openfoam results shown a constant increase in the drag coefficient! while it should be settled after awhile!

I attached the case files and a pictures of alpha and velocity field. Any advice are appreciated.

Bashar

Hello foamers,
I am using OF 2.3.0 and I have a problem with BC for open-channel flow as well.
I hope somebody can help.
my geometry is a 2D channel with an inlet, outlet, top and bottom.

inlet:
p-rgh: fixedFluxPressure
U: variableHeightFlowRateInletVelocity
alpha.water: variableHeightFlowRate

bottom:
p-rgh: fixedFluxPressure
U: fixedValue (0 0 0)
alpha.water: zeroGradient

outlet, top:
P-rgh: totalPressure (po=0)
U: pressureInletOutletVelocity
alpha.water: inletOutlet (inletValue uniform 0;value uniform 0; )

in this case, the velocity of the water near the outlet increase and the water level decrease. (I think based on what I read about it in the forum the water lever should stay constant until the end!!)

on the other case,
I changed the outlet to
outlet:
p-rgh: outletInlet
U: inletOutlet
alpha.water: zeroGradient

the result is the same as the first case.

in the other case:

oulet:
p-rgh: zeroGradient
U: inletOutlet
alpha.water: zeroGradient

in this case the water level increase and the velocity decrease in outlet.

the only method I can keep the water level constant is what kflora suggested, but I am not sure that it gives the correct velocity profile. because I used very long channel but it seems that it will affect the results in some way, finally.
I am not very experienced in the simulation. I hope somebody can help.

Regards,
Amir