[daev]

Hello,

I would like to solve flow in complex duct where the flow conditions are determined by pressure at the inlet and at the outlet.

Firtsly I am trying to solve flow in a simple rectangular channel using simpleFoam solver, no turbulence, where the 'system' directory is taken from tutorial simpleFoam/pitzDaily (everything only 2D). The channel is made by blockMesh, 0.1 m long and 0.02 m height (I will upload the case in next post)

straightMesh.jpg

At the inlet and at the outlet pressure is specified while velocity is prescribed as pressureInletVelocity and pressureInletOutletVelocity, respectively for inlet and outlet. This runs without problems, the results match perfectly with analyticaly calculated laminar velocity profile. (see results in attached straightU.jpg)

straightU.jpg

Secondly I tried the same scenario but now the geometry (same lenght and channel height) is slightly skewed, so the channel is tilted (again I will upload the case in next post)

tiltedMesh.jpg

This solution does not converge and eventually crashes as the 'time step continuity errors' explodes to astronomical values. Mesh and results are in attached tiltedMesh.jpg, tiltedU.jpg and tiltedP.jpg

tiltedU.jpg

tiltedP.jpg

I also tried to engage nOrthogonalCorrections but without any change.

Does anybody have idea what's wrong?

Thanks and Marry X-mas

[daev]

Here are the cases.
I also forgot to mention, both results are after 1000 iterations (not converged). Straight case was ran up to 20,000 iterations, when solution was compared with analytical, while tilted one started to boost 'time step continuity errors' after approx 2000 iterations and crashed after 13,000 iterations.
checkMesh is in both cases fine
OpenFoam 1.7.1 was operated.

straight.zip

tilted.zip