CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Meshing & Mesh Conversion

[blockMesh] Mesh defining with parameters

Register Blogs Community New Posts Updated Threads Search

Like Tree4Likes
  • 2 Post By HPE
  • 1 Post By GerhardHolzinger
  • 1 Post By HPE

 
 
LinkBack Thread Tools Search this Thread Display Modes
Prev Previous Post   Next Post Next
Old   February 9, 2021, 14:06
Default Mesh defining with parameters
  #1
Member
 
Wowbagger's Avatar
 
Michael S.
Join Date: May 2020
Location: Hamburg, Germany
Posts: 50
Rep Power: 6
Wowbagger is on a distinguished road
Send a message via Skype™ to Wowbagger
Hi everyone,

I want to simulate a kind of tubular heat exchanger and prepared the mesh via blockmesh.
Because of its symmetry I made the model as a wedge.

Maybe this is irrelevant depending my question but I get the same warning for a change in the geometry setup.
For some bigger angles I get the known (?) error below.

Quote:
--> FOAM Warning : From function virtual void Foam::wedgePolyPatch::calcGeometry(Foam::PstreamBu ffers&) in file meshes/polyMesh/polyPatches/constraint/wedge/wedgePolyPatch.C at line 70 Wedge patch 'symmetrie_back' is not planar.At local face at (2.925 119.75 0.102082) the normal (-0.0348784 -1.37971e-17 0.999392) differs from the average normal (-0.0348783 -3.67866e-20 0.999392) by 7.88202e-15Either correct the patch or split it into planar parts
So far so good. This should not really a problem as I read but to avoid this I use to go to really tiny angles.
All constructive data are set as parameters in the first part of the blockMeshDict file.
There comes the strange thing, I recognised that it is not a problem to do some calculation with the parametric setup but for the grading it is not possible.

I did the grading as shown below according to this file: http://www.wolfdynamics.com/wiki/mes..._blockmesh.pdf

Code:
    hex (5 7 8 9 23 25 26 27) ($ZR3 $ZL2 1) simpleGrading           // Block 3 Fluid 
    (
            (
            //x-Achse
            ($BB $DD $EE) 
            ($AA $CC 1) 
            ($BB $DD $FF) 
            )            
        (
        //y-Achse
        1;
        )
        (
        //z-Achse
        1
        )
    )
But this is only possible if the FF variable is not calculated.
At this time my blockMeshDict file looks like this and it works fine.
But I would prefer to calculate the FF dependent on EE as shown in the commented lines:


Code:
    // Grading der Zellen, Ringspalt
    AA 0.2;                  // Anteil Kernbereich
    BB #calc "(1-$AA)/2";    // Bereich Grenzschicht
    CC 0.1;                  // Anteil Zellen Kernbereich
    DD #calc "(1-$CC)/2";    // Anteil Zellen Grenzschicht
    EE 4;                    // Grading Da
    FF 0.25;
//    FF #calc "(5-$EE)";    // Grading Di
//    FF #calc "pow ($EE, -1)";        // Grading Di - Berechnen geht nicht
//    FF #calc "1/($EE)";              // Grading Di - Berechnen geht nicht
If I start the meshing with one of the comment lines I get the same Warning as mentioned above (--> FOAM Warning ..... Either correct the patch or split it into planar parts) and paraFoam will close after a short time telling me "Segmentation fault (core dumped)".

So far I realised in the last few attemts that the meshing works with values for EE<2. I have no problem setting one more parameter for FF but I would like to know why this calculation is not working in that intented use. Does anyone have a clue?

Best regards
Michael
Wowbagger is offline   Reply With Quote

 


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[snappyHexMesh] Creating multiple multiple cell zones with snappyHexMesh - a newbie in deep water! divergence OpenFOAM Meshing & Mesh Conversion 0 January 23, 2019 04:17
[snappyHexMesh] SnappyHexMesh no layers and no decent mesh for complex geometry pizzaspinate OpenFOAM Meshing & Mesh Conversion 1 February 25, 2015 07:05
[Gmsh] 2D Mesh Generation Tutorial for GMSH aeroslacker OpenFOAM Meshing & Mesh Conversion 12 January 19, 2012 03:52
Convergence moving mesh lr103476 OpenFOAM Running, Solving & CFD 30 November 19, 2007 14:09
Icemcfd 11: Loss of mesh from surface mesh option? Joe CFX 2 March 26, 2007 18:10


All times are GMT -4. The time now is 04:01.