The density field can be calculated in SEM or in HOL as follows:

RHO = RHOG * (1-VFOL) + RHOL * VFOL (it's logical)

but

why ? the viscosity field is calculated as follows:

ENUL = ENULG * (1-VFOL) + ENULL * VFOL

Hi, Medhi

In my oppinion it is a simple, but a very nice and strong question? You maybe are "touching" a delicate point inside those methods.

Could you be more specific about why you think it is improper to calculate the viscosity by this expression? In principle, you are averaging the fluids properties.

Kike

Hi, Enrique

If you take a cell and you want to calculate the density, you can write

Mt = Ml+Mg

where

Mt: total mass, Ml: mass of liquid and Mg: mass of gas.

so,

Mt/Vt = Ml/Vt+Mg/Vt

Vt: total volume of the cell.

and finally

Mt/Vt= (Ml/Vl)*(Vl/Vt) + (Mg/Vg)*(Vg/Vl)

Where

Vl: volume of liquid and Vg: volume of gaz

and we have:

RHO = RHOG * (1-VFOL) + RHOL * VFOL

but I don't know why we can calculate the viscosity in this form:

ENUL = ENULG * (1-VFOL) + ENULL * VFOL

Dear Medhi

I am sorry to ansuwer you so late. Effectively there is no physical explanation for this "special average". It is a gross approximation to obtain a transport coefficient for momentum according to the presence of the two phases inside the cell.

In those cells you have the stress balance expresions (normal and tangential to the interface) and you can extract a different relation for the viscosity coefficients ENULG, ENULL and a total viscosity coefficient ENUL from the tangential relation. It is more consistent with the stress continuity at interfaces considering piecewise linear velocity profile in each phase. After some algebraic calculations you will have:

1/ENUL=VFOL*(1/ENULL)+((1-VFOL)*(1/ENULG))

This has been proposed for a professor of CALTEC but I don't remember his name. It was a friend of mine who gave me the hint.

Nevertheless, every body uses the relation you find into PHOENICS and it works for most cases (I mean with respect to experimental results). Try to implement the new relation and see what the differences are. Please tell me what you find after the change in your problem.

Regards

Kike