[doodek]

Hello,

I am trying to solve flow in convergent-divergent "nozzle" (2D) using rhoSimplecFoam solver. My computational domain you can see below. However, after ~150 iterations residual values increased and they are oscillating. I am interested in temperature field, but it seems that something bad happened and maximum temperature raised after 800 iterations, as you can see in the attached plot. Decreasing relaxation factors does not help.

Do you have any idea how to solve this problem?

My settings are the following:

RASModel realizableKE;

thermoType hPsiThermo<pureMixture<constTransport<specieThermo <hConstThermo<perfectGas>>>>>;

relaxationFactors
{
U 0.5;
nuTilda 0.8;
e 0.5;
p 0.3;
epsilon 0.8;
rho 0.5;
k 0.8;
omega 0.8;
h 0.5;
}

rhoMin rhoMin [1 -3 0 0 0 0 0] 50;
rhoMax rhoMax [1 -3 0 0 0 0 0] 100;

BC's:

inlet:
p, totalPressure: 3.1e6
u, pressureInletOutletVelocity
T, inletOutlet: 270

outlet:
p, fixedValue: 3e6
u, pressureInletOutletVelocity
T, inletOutlet: 270

wall:
fixedValues and wallFunction for k and epsilon

frontAndBackPlanes: empty