[Saideep]

Hi guys;

I am stuck with a peculiar problem and would appreciate any sort of advice in this matter.

I am using interfoam solver for the following setup.

I have 4 injection channels {alternate injections of wetting(w), non-wetting(nw) phases} as shown in the figure.
problem1.jpg

An ideal bc at the inlet would be is to fix the velocity and compute the pressure at the ends.

But due to the shape of the domain i am not able to get an exact phase pressure at the inlet and outlet which i need.

So, an alternative is to go for fixing the pressures at the inlets and outlets.
I define the inlet velocity bc as pressureInletOuletVelocity and oulet bc for velocity as zeroGradient.

Interfoam works as expected when the pressure of w phase is less than the nw phase.
But when i try to equal or increase the pressure of the wetting phase the entire flow is reversed in the domain as can be seen from the output.

/*-ve sign here at inlet means flow 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 sign at outlet means flow out of the domain.*/
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
All the variables calculated from thereon are wrong interms of alpha> 4000, p ranging from e-7 to -e-5 Pa and so on.

Is there anything fishy about the bcs that i am using or is there a better bc for the case that i mentoined or a better solution to solve this kind of a case?

Thanks;
Saideep