[klilla]

Dear All,

I try to develop a solver: incompressible 2phase, temperature dependent VOF flow solver embedded in the cht multiregion environment.

What I have: vof+temperature without cht working fine
cht + single phase solver working fine

Where I fail is the definition of alpha (void fraction) in the multiregion environment. My alpha do not get updated if I try:

creatFluidFields.H : PtrList<volScalarField> alpha1Fluid(fluidRegions.size());
alpha1Fluid.set
(
i,
new volScalarField
(
thermoFluid[i].alpha1()
)
);

setRegionFluidFields.H: volScalarField& alpha1 = alpha1Fluid[i];

I thought to do it in an other way:

creatFluidFields.H : nothing of alpha1

setRegionFluidFields.H: volScalarField& alpha1 = thermo.alpha1();

with thermo = thermoFluid[i]

This is updating the field, but just close to the walls and gives very wierd results, though the equations I'm solving are the same as I use for the solver without cht.

Can anyone help?
Thanks,
Lilla

[kmooney]

You may need to adapt the CHT boundary conditions to take into account the variable density and viscosity (and thus the local heat transfer coeffs)..

Just a guess!

Kyle

[klilla]

Dear Kyle,

thanks for your answer. In fact, that I did already. In any case, right now I run an isotherm case, my problem is that the free-surface is not moving...

Thanks again,
Lilla