[Tetragramm]

I'm having problems with my mesh causing divergence. I'm trying to replicate the results in this PDF here but the mesh I've created has issues are causing simpleFoam to diverge. For some reason, the solution is only diverging at higher angles of attack.

I've put together a set of pictures showing the problem and my geometry. They are located here. For some reason, all along the join between the four hex blocks, the solution diverges. The most extreme example is the portion pictured in the album.

I'm not sure if the non-ortho faces that checkMesh finds are the source of the problem. There are some right at the wingtip, where the most extreme divergence is, but the divergence continues all the way out, but there are no non-ortho faces that far out.

I'm hoping that one of you can tell me how to fix the problem without having to re-do the whole blockMesh file. Or, if you say that's not possible, to give me recommendations on how to make this geometry without causing bad cells and faces in the mesh. I'd like any suggestions to keep in mind that I will be changing things like airfoils and aspect ratio, so it should be automation friendly.

Thank you for you help.

I've attached my blockMesh file and this is my checkMesh file

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.1.1-221db2718bbb
Exec   : checkMesh
Date   : Sep 22 2012
Time   : 13:54:28
Host   : "Heterodyne"
PID    : 2781
Case   : /home/tetragramm/OpenFOAM/0012
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
    points:           4040101
    faces:            12040000
    internal faces:   11920000
    cells:            4000000
    boundary patches: 4
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     3960000
    prisms:        40000
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    0
    polyhedra:     0

Checking topology...
    Boundary definition OK.
    Cell to face addressing OK.
    Point usage OK.
    Upper triangular ordering OK.
    Face vertices OK.
    Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces ...
    Patch               Faces    Points   Surface topology                  
    movingWall          40000    40001    ok (non-closed singly connected)  
    fixedWalls          40000    40001    ok (non-closed singly connected)  
    symmetricWall       40000    40400    ok (non-closed singly connected)  
    allAround           0        0        ok (empty)                        

Checking geometry...
    Overall domain bounding box (-2.59808 -2.84808 0) (2.59808 2.34808 3)
    Mesh (non-empty, non-wedge) directions (1 1 1)
    Mesh (non-empty) directions (1 1 1)
    Boundary openness (-2.07313e-15 -5.50671e-15 5.79448e-15) OK.
    Max cell openness = 1.74778e-15 OK.
    Max aspect ratio = 149.528 OK.
    Minumum face area = 4.19645e-07. Maximum face area = 0.00338268.  Face area magnitudes OK.
    Min volume = 1.46671e-09. Max volume = 7.5958e-05.  Total volume = 30.7017.  Cell volumes OK.
    Mesh non-orthogonality Max: 86.9482 average: 28.9597
   *Number of severely non-orthogonal faces: 24795.
    Non-orthogonality check OK.
  <<Writing 24795 non-orthogonal faces to set nonOrthoFaces
    Face pyramids OK.
 ***Max skewness = 11.086, 9003 highly skew faces detected which may impair the quality of the results
  <<Writing 9003 skew faces to set skewFaces
    Coupled point location match (average 0) OK.

Failed 1 mesh checks.

End