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/)
-   -   [SOWFA] Problem Defining terrainBlockMesherDict, point swl to stl (https://www.cfd-online.com/Forums/openfoam-community-contributions/186122-problem-defining-terrainblockmesherdict-point-swl-stl.html)

chusma April 11, 2017 04:59

Problem Defining terrainBlockMesherDict, point swl to stl
 
Hi to everyone,

I am currently trying to figure how terrainBlockMesher utility works, but I am having some problems with the Dict file. I have defined everything like is told in the tutorial and I have adapted the file to my case, but the same error appears with any modification.

Could somebody help me with this problem? I will appreciate it.

THE ERROR FILE:

PHP Code:

Create time

Reading stl surface
...
...
doneafter 0.01 s.

Running TerrainManager
Selecting transition 
function type 'linear'
   
loading block grading module
   found 5 regions 
for block grading in direction 1.
   found 1 regions 
for block grading in direction 2.
Writing file 
"domainBox.stl"
Writing file "stlBox.stl"

TerrainManagercreating 882 points400 blocks and 840 edges
TerrainManager
Cannot attach point SWL = (0 0 3000to STL.


TerrainManagerError attaching points.

terminate called without an active exception
Aborted 


MY TERRAINBLOCKMESHERDICT:
PHP Code:

stl
{
    
doubleHill.stl
    
{
        
type triSurfaceMesh;
    }
};

writeBlockMeshDict      true;
writePolyMesh           true;

blockManager
{
        
coordinates
        
{
                
baseVectors
                
(
                                ( 
1  0  0// terrain direction 1
                                
0  1  0// terrain direction 2
                                
0  0  1// terrain normal (height direction)
                
);
        }

        
// Specification of the over all box. Height values will be adjusted to stl.
        // p_corner is the lower left box corner, the above base vectors span the box.
        // Dimensions are lengths scales of the box, wrt base vectors given above.
        
p_corner        (0 0 0);

        
dimensions      (396 0.01 85);

        
// a point well above the stl.
        
p_above         (0 0 3000.0);

        
// the numner of blocks in terrain directions 1 and 2, z will be ignored and set to 1
        
blocks          (20 20 1);

        
// the number of cells per block in the three directions
        
cells           (3 3 100);

        
// block grading by regions
        
blockGrading
        
{
           
direction1
           
{
                
region_1
                
{
                        
width           20;
                        
blocks          2;
                        
type            constant;
                }

                
region_2
                
{
                        
width           20;
                        
blocks          4;
                        
type            interpolating;
                }

                
region_3
                
{
                        
width           20;
                        
blocks          8;
                        
type            constant;
                }

                
region_4
                
{
                        
width           20;
                        
blocks          4;
                        
type            interpolating;
                }


                
// the last region is filled automatically
                
region_5
                
{
                        
type            constant;
                }
           }

           
direction2
           
{
                
region_1
                
{
                        
width           0.01;
                        
blocks          20;
                        
type            constant;

           }
        }
}
        
// the maximal distance searched for projection
        
maxDistProj     10000;

        
// the grading command
        
gradingFactors  (1 1 1000);

        
// the patch names
        
patch_name_west         inlet;
        
patch_name_east         outlet;
        
patch_name_north        latLeft;
        
patch_name_south        latRight;
        
patch_name_sky          top;
 
patch_name_ground       terrain;

        
// the patch types
        // (modifyOrography overwrites for makeCyclic-commands)
        
patch_type_west         patch;
        
patch_type_east         patch;
        
patch_type_north        patch;
        
patch_type_south        patch;
        
patch_type_sky          patch;
        
patch_type_ground       wall;



        
// make upward splines orthogonal
        
orthogonalizeUpwardSplines
        
{
                
ignoreBoundary;
                
splineNormalDist        5;
        }


}

checkMesh               false;//true;
checkMeshNoTopology     false;
checkMeshAllGeometry    false;//true;
checkMeshAllTopology    false;//true;

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



All times are GMT -4. The time now is 14:53.