[Mat_fr]

Dear Foamers,


I am working with a two-phase flow problem involving relatively high pressure and temperature, and for that I am using compressibleInterFoam



When using adiabatic boundary conditions for the wall of my domain, I can obtain physical results.
However, when using isothermal conditions, the simulation does not converge. I have attached images of the temperature and velocity fields near the wall, in a region of pure gas. The temperature gradient between the cold wall and the warm gas is creating fluid motion towards the wall. Eventually this field becomes too high and deteriorates the other physical variables, generating finally high temperature in the first cell layers.
I understand that the temperature gradient creates fluid motion, but this is unrealistic and especially considering the no-slip wall condition for the velocity.


The thermophysical properties for the gas are:

thermoType
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       sutherland;
    thermo          hConst;
    equationOfState PengRobinsonGas;
    specie          specie;
    energy          sensibleEnthalpy;
}

Any idea about the problem?


I have tried to use an externalWallHeatFluxTemperature boundary condition for the temperature, which is less "strong" than the isothermal one. Although the computation does converge in this case, results are unrealistic as way too much heat is transferred out of the gas.



Regards,

[Mat_fr]

Hello,
No one has an idea regarding this issue ?