CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Community Contributions (https://www.cfd-online.com/Forums/openfoam-community-contributions/)
-   -   [cfMesh] pMesh tetMesh fails (https://www.cfd-online.com/Forums/openfoam-community-contributions/221855-pmesh-tetmesh-fails.html)

user007 November 2, 2019 18:56

pMesh tetMesh fails
 
Hi, and thanks for the great tool cfMesh.
I tried some mesh with cartesianMesh since I guess is the fastest one and when I was satisfied with refinement I switch to tetMesh and pMesh, but I encurred in the following error:

Code:

...
 ...
...
 Extracting voronoi mesh
Resizing points!


--> FOAM FATAL ERROR:
Negative size requested.

    From function template<class T, Foam::label Offset>
inline void Foam::LongList<T, Offset>::allocateSize(const label)
    in file lnInclude/OF__LongListI.H at line 73.

  FOAM aborting


If I disable the refinements, the mesh continues without errors, obviously producing a too coarse result, so am I doing something wrong in the definition of refinements?
Has it already been reported that the cartesianMesh is more robust than pMesh and tetMesh?
Is there any way to visualize any areas where geometry can cause problems to the mesher?


I'm using cfMesh 1.1.1 under Windows 10, and in the following the meshDict i'm using:

Code:

surfaceFile "constant/triSurf/Geometry.fms";

enforceGeometryConstraints 0;
keepCellsIntersectingBoundary 0;

maxCellSize 0.128;

localRefinement {

    Main {
        cellSize                    0.012;
        refinementThickness        0.040;
    }

    SubSet_1 {
        cellSize                    0.006;
    }

    SubSet_2 {
        cellSize                    0.006;
    }

    Refinement_1 {
        cellSize                    0.002;
    }
   

    Refinement_2  {
        cellSize                    0.012;
        refinementThickness        0.1;
    }
   
}

objectRefinements {

    Box_Small {
        type        box;
        cellSize    0.016;
        centre      (1.35 0.0 0.0);
        lengthX      5.5;
        lengthY      2.5;
        lengthZ      2.5;
    }

    Box_Medium {
        type        box;
        cellSize    0.032;
        centre      (6.35 0.0 0.0);
        lengthX      10.5;
        lengthY      4.0;
        lengthZ      4.0;
    }

    Box_Large {
        type        box;
        cellSize    0.064;
        centre      (11.35 0.0 0.0);
        lengthX      15.5;
        lengthY      5.0;
        lengthZ      5.0;
    }

}


Krao November 4, 2019 03:50

Hi,

try using additionalRefinementLevels instead of cellSize. If the problem still persists then, more investigation is required

user007 November 9, 2019 04:32

Quote:

Originally Posted by Krao (Post 748784)
Hi,

try using additionalRefinementLevels instead of cellSize. If the problem still persists then, more investigation is required

The error got away recompiling Foam with 64bit labels.


All times are GMT -4. The time now is 22:40.