[CFD2D!]

Good afternoon,
I am new here in this forum. I hope that there is someone who can help me eventually.

I simulate the flow of a liquid (oil) on a capillary surface.My model is 2D (axisimetric), Eulerian multiphase, VOF, laminar, Phase Interactions (surface tension force), Phase 1: Liquid (prime) Phase 2: Gas (Air, secondary).

I am rebuilding a publication and have exactly the same domain (meshoptions, base size, number of cells).

My problem:
The system is under atmospheric pressure and there are 1 inlet (top) and 2 pressure outlets (one at the top next to the inlet and one at the bottom).
The inlet is for the fluid and I have added a field function to the initial conditions in the physics:

((${RadialCoordinate} > 0.000475) && (${RadialCoordinate} < 0.001))? 0:1

This field function says for the liquid phase (composite N-1) that at the beginning of the simulation my wire is partly covered with liquid and the rest is automatically lust (N-1). This function is for the Volumefractions (Initial conditions).

My Inlet from which only the liquid comes has a Volumefraction of 1natural and air is present in the system through the Initial Condition. But how do I define the volume fractions of both outlets? I get backflows all the time when I define the outlets, no matter how. No liquid should come out of the top of the outlet, but air can come out by displacement and air can come out of the bottom. Is there any physical explanation for this and maybe you can develop a field function for it.


And of course I have taken gravity into account. Because I have an axis-symmetric model around the x-axis, my gravity is [9.81,0.0] m/s. I have also included a piece of the capillary (imprint) and defined the capillary center as axis of rotation. Where is my error?
And is there a possibility to display each phase in a different colour?

Greeting

CFD2D!