[HiF]

Hello everyone,

I'm trying to calculate a simple rotating cylinder with temperature gradient on the wall by dsmcFoam. As you see in the attached file, I think I correctly made the mesh and boundary condition, but I don't know why the pressure distribution seems strange (wall side should be higher)

Does anyone have some idea or comment about this problem?

Thank you in advance and sorry for bad English.

(In the attached file, boudary condition and blockMeshdict are contained)

convertToMeters 1;

vertices
(
( 0 0 0)
( 0.665 0 0)
( 0 0.0921480096 -0.0080619062)
( 0.665 0.0921480096 -0.0080619062)

( 0 0.0921480096 0.0080619062)
( 0.665 0.0921480096 0.0080619062)

);

blocks
(
hex (0 1 3 2 0 1 5 4) (100 40 1) simpleGrading (1 0.2 1)
);

edges
(
);

patches
(
symmetryPlane rotAxis1
(
(0 1 1 0)
)
wall UrotatingWall
(
(0 0 4 2)
)
wall BrotatingWall
(
(1 1 3 5)
)

wall wallShroudLateral
(
(2 4 5 3)
)
wedge front
(
(0 1 5 4)
)
wedge back
(
(0 2 3 1)
)
);

mergePatchPairs
(
);

// ************************************************** *********************** //