[EvgenyMogilevskiy]

Hello,

I'm trying to simulate a laminar axissymnnetric flow of argon through a convergent-divergent nozzle with openFOAM. I use rhoCentralFoam solver.

BC's are the following:
inlet:
pressure timeVaryingTotalPressure
velocity pressureInletOutletVelocity
temperature totalTemperature

wall:
pressure zeroGradient
velocity 0 0 0
temperature zeroGradient

outlet
pressure, temperature waveTrassmissive
velocity zeroGradient

Sutherland transport properties are used with Pr=1.

To start the simulation, I set pressure about 1000 Pa in divergent part and 5000 Pa in convergent one. Then total pressure in the inlet starts to grow up to desired value of 50 000 Pa. Initial temperature is 300 K, total temperature at the inlet is also 300 K

I see very strange behavior of T field. First, when the inlet pressure starts to grow I see something like a shock wave with quite a big temperature (up to 400 K), but pressure goes up slowly. Then near the wall, hot boundary layer is formed. T reaches 500 K there. I tried to refine mesh in this layer to resolve temperature gradient but it seems this result is mesh-independent.

Does anyone have any suggestions?