[Artur]

Hi All,

I'm trying to run a 2D axisymmetric simulation with simpleFoam but not matter what I do it always seems to try to solve for the 3rd component of velocity when I used the "wedge" BC. I have dug into the tutorials and see the same thing happen in interFoam/nozzleFlow2D, for instance. The mesh check also clearly shows that the grid is indeed 2D but for some reason has 3 solution directions:

Checking geometry...
    Overall domain bounding box (-0.65 0 -0.028352602) (1.95 0.65 0.028352602)
    Mesh has 2 geometric (non-empty/wedge) directions (1 1 0)
    Mesh has 3 solution (non-empty) directions (1 1 1)
    Wedge side_pos with angle 2.4976236 degrees
    Wedge side_neg with angle 2.4976236 degrees
    All edges aligned with or perpendicular to non-empty directions.
    Boundary openness (8.0639549e-19 -3.9251831e-15 2.7113778e-14) OK.
    Max cell openness = 2.4730408e-15 OK.
    Max aspect ratio = 125.01216 OK.
    Minimum face area = 2.5901831e-11. Maximum face area = 0.0029599635.  Face area magnitudes OK.
    Min volume = 1.6935059e-14. Max volume = 0.00013206161.  Total volume = 0.0478444.  Cell volumes OK.
    Mesh non-orthogonality Max: 75.256619 average: 4.8024892
   *Number of severely non-orthogonal (> 70 degrees) faces: 1.
    Non-orthogonality check OK.
  <<Writing 1 non-orthogonal faces to set nonOrthoFaces
    Face pyramids OK.
    Max skewness = 0.8198976 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

For this reason I can never achieve good convergence, even though all the other fields converge very well (see attached figure). I will need to run a variant of this set up dozens of times and so stopping the simulation at convergence and not after a fixed amount of iterations is very important to me. Does anyone have any suggestions on how to do it?

Thanks,


Artur

residuals_testrun.png