Parameterfile Curved Structured Mesh

From HOPR
Jump to: navigation, search












This parameter file can be found in

tutorials/1-05-curved_structured/parameter.ini

See Curved Structured Mesh for a detailed description. The parameter file need not follow a given structure. If mandatory parameters are not given, an error message will be shown.

Parameterfile

!=============================================================================== !
! OUTPUT
!=============================================================================== !
  ProjectName  =CYLINDER_Ngeo5         ! name of the project (used for filenames)
  Debugvisu    =T                        ! Write debug mesh to tecplot file
 
   DebugvisuLevel=1                      ! Level 0: visualization of linear mesh and BC (default), 
                                         ! Level 1: + curved surface visualization (_SplineSurf.dat), if useCurveds
   NVisu        = 8                      ! number of visualization points  per element edge, if useCurved
 
  checkElemJacobians=T                   ! checks the Jacobian and scaled Jacobian for each element
!=============================================================================== !
! MESH
!=============================================================================== !
  Mode         =11                     ! 11 Curved structured block, only hexahedra
  nZones       =1                      ! number of zones (only one possible!)
  nElems       =(/8,8,1/)              ! number of elements in each direction
  BCIndex      =(/2,1,0,4,0,3/)        ! Indices of Boundary Conditions for  six Boundary Faces (z-,y-,x+,y+,x-,z+)
 
  Meshtype     = 3    ! Mesh Types: 1 - Cube     (origin + dimensions)
                      !             2 - Bilinear (8 points CGNS notation)
                      !             3 - Curved (add Whichmapping)
    WhichMapping = 4  ! type of mapping using 6 boundary faces to build the block:
                      !  3...half cylinder, 4...full cylinder
  R_0=0.5             ! inner radius 
  R_INF=20.           ! outer radius of domain
  DZ=2.               ! z dimension: [-DZ,DZ]
  stretchType=(/3,1,0/)
  fac=(/1.5,2.2,10./)
  DXmaxToDXmin=(/6.,100.,1./)
 
  SpaceQuandt  =1.                          ! characteristic length of the mesh
!=============================================================================== !
! CURVED
!=============================================================================== !
  useCurveds   =T                     ! T if curved boundaries defined
  BoundaryOrder= 6                    ! choose order freely!
!=============================================================================== !
! BOUNDARY CONDITIONS
!=============================================================================== !
  BoundaryName=BC_cylinder
  BoundaryType=(/4,0,1,0/)
  BoundaryName=BC_periodicZ-
  BoundaryType=(/1,0,0,1/)
  BoundaryName=BC_periodicZ+
  BoundaryType=(/1,0,0,-1/)
  BoundaryName=BC_farfield
  BoundaryType=(/2,0,0,0/)
 
  vv=(/0.,0.,4./)                    ! vector for periodic BC in z direction (zminus,zplus), index=1
                                       ! here size of domain in z direction is [-2.,2.]

Screen Output

After executing the upper parameter file the following output will be shown:

====================================================================================================================================
           _______     _______    ___________________    ___________________   ___________________                         xX    
          /      /)   /      /)  /                  /)  /                  /) /                  /)     .xXXXXXXXXx.       X     
         /      //   /      //  /      _____       //  /      _____       // /      _____       //    .XXXXXXXXXXXXXXXx  .XXXXx  
        /      //   /      //  /      /)___/      //  /      /)___/      // /      /)___/      //   .XXXXXXXXXXXXXXXXXXXXXXXXXXx 
       /      //___/      //  /      //   /      //  /      //___/      // /      //___/      //  .XXXXXXXXXXXXXXXXXXXXXXXX´     
      /                  //  /      //   /      //  /                  // /                  //  .XX``XXXXXXXXXXXXXXXXXX´        
     /      _____       //  /      //   /      //  /      ____________// /      __      ____//   XX`  `XXXXX`      .X´           
    /      /)___/      //  /      //   /      //  /      /)___________) /      /)_|    |____)   XX     XXX`       .´             
   /      //   /      //  /      //___/      //  /      //             /      //  |    |__     ,X`    XXX´                       
  /      //   /      //  /                  //  /      //             /      //   |      /)   ,X`   .XX´                         
 /______//   /______//  /__________________//  /______//             /______//    |_____//   ,X`   XX´                           
 )______)    )______)   )__________________)   )______)              )______)     )_____)   xX    XXx                            
 
====================================================================================================================================
------------------------------------------------------------------------------------------------------------------------------------
 INIT OUTPUT...
 | Reading from file "parameter.ini":
 |                    projectname |                    CYLINDER_Ngeo5 | *CUSTOM | 
 |                        Logging |                                 F | DEFAULT | 
 |                      DebugVisu |                                 T | *CUSTOM | 
 |                 DebugVisuLevel |                                 1 | *CUSTOM | 
 |                   OutputFormat |                                 0 | DEFAULT | 
 |                       sfc_type |                           hilbert | DEFAULT | 
 |                      doSortIJK |                                 F | DEFAULT | 
 |           useSpaceFillingCurve |                                 T | DEFAULT | 
 INIT OUTPUT DONE!
------------------------------------------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------------------------------------------
 INIT MESH...
 |                     useCurveds |                                 T | *CUSTOM | 
 |                  BoundaryOrder |                                 6 | *CUSTOM | 
 |                           Mode |                                11 | *CUSTOM | 
 |                         nZones |                                 1 | *CUSTOM | 
 |                    stretchType |      Integer array of size (   3) | *CUSTOM |     3    1    0
 |                            fac |         Real array of size (   3) | *CUSTOM |  1.50 2.2010.00
 |                           fac2 |         Real array of size (   3) | DEFAULT |  1.50 2.2010.00
 |                   DxmaxToDxmin |         Real array of size (   3) | *CUSTOM |  6.00***** 1.00
 |                       Meshtype |                                 3 | *CUSTOM | 
 |                        BCIndex |      Integer array of size (   6) | *CUSTOM |     2    1    0    4    0    3
 |                         nElems |      Integer array of size (   3) | *CUSTOM |     8    8    1
 |                   WhichMapping |                                 4 | *CUSTOM | 
 |                            R_0 |                       0.50000E+00 | *CUSTOM | 
 |                          R_INF |                       0.20000E+02 | *CUSTOM | 
 |                             DZ |                       0.20000E+01 | *CUSTOM | 
 |                   preScaleMesh |                                 F | DEFAULT | 
 |                  postScaleMesh |                                 F | DEFAULT | 
 |                      meshScale |                       0.10000E+01 | DEFAULT | 
 |                    SpaceQuandt |                       0.10000E+01 | *CUSTOM | 
 |                 rebuildCurveds |                                 F | DEFAULT | 
 |          doExactSurfProjection |                                 F | DEFAULT | 
 |          nCurvedBoundaryLayers |                                -1 | DEFAULT | 
 |                   BoundaryName |                       BC_cylinder | *CUSTOM | 
 |                   BoundaryType |      Integer array of size (   4) | *CUSTOM |     4    0    1    0
 |                   BoundaryName |                     BC_periodicZ- | *CUSTOM | 
 |                   BoundaryType |      Integer array of size (   4) | *CUSTOM |     1    0    0    1
 |                   BoundaryName |                     BC_periodicZ+ | *CUSTOM | 
 |                   BoundaryType |      Integer array of size (   4) | *CUSTOM |     1    0    0   -1
 |                   BoundaryName |                       BC_farfield | *CUSTOM | 
 |                   BoundaryType |      Integer array of size (   4) | *CUSTOM |     2    0    0    0
 |                             vv |         Real array of size (   3) | *CUSTOM |  0.00 0.00 4.00
 |                  doZcorrection |                                 F | DEFAULT | 
 |                        OrientZ |                                 F | DEFAULT | 
 |                     SplitToHex |                                 F | DEFAULT | 
 |                      nFineHexa |                                 1 | DEFAULT | 
 |                 MeshPostDeform |                                 0 | DEFAULT | 
 |                 ConformConnect |                                 T | DEFAULT | 
 |             checkElemJacobians |                                 T | *CUSTOM | 
 |              jacobianTolerance |                       0.10000E-15 | DEFAULT | 
 INIT MESH DONE!
------------------------------------------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------------------------------------------
 INIT BASIS...
 |                          nVisu |                                 8 | *CUSTOM | 
 |                       nAnalyze |                                 7 | DEFAULT | 
 INIT BASIS DONE!
------------------------------------------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------------------------------------------
 INIT SEARCH...
 |               nElemsNodeSearch |                                25 | DEFAULT | 
 |               RefineSideSearch |                       0.80000E+01 | DEFAULT | 
 INIT SEARCH DONE!
------------------------------------------------------------------------------------------------------------------------------------
------------------------------------------------------------------------------------------------------------------------------------
 THE FOLLOWING INI-FILE PARAMETERS WERE IGNORED:
------------------------------------------------------------------------------------------------------------------------------------
====================================================================================================================================
 
                                              Entering fillMesh 
 
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 BUILDING CURVED CARTESIAN MESH...
DONE!   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   [ time: .038s ]
 Number of Elements:           64
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
WRITING THE DEBUGMESH...
   #Elements           64
   WRITE DATA TO VTX XML BINARY (VTU) FILE... CYLINDER_Ngeo5_Debugmesh.vtu   DONE
    Mesh visualized for debug purposes in file : CYLINDER_Ngeo5_Debugmesh
DONE!   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   [ time: .000s ]
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
WRITING THE BC MESH...
   #BCSides          144
   WRITE DATA TO VTX XML BINARY (VTU) FILE... CYLINDER_Ngeo5_Debugmesh_BC.vtu   DONE
    Boundary mesh visualized for debug purposes in file : CYLINDER_Ngeo5_Debugmesh_BC
DONE!   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   [ time: .000s ]
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 Mesh connect starts
Count sides, delete Null BCs and check periodic BCs...
-----------------------------------
number of sides          :          384
number of Inner sides    :          240
number of BC sides       :          144
number of periodic sides :          128
-----------------------------------
Insert periodic sides...
Eliminate multiple Nodes...
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
GLOBAL UNIQUE NODES ...
   All Nodes sorted...
   Number of nodes to check:        10342
    size of tolerance box:                    1
  Number of deleted nodes         502
  Number of unique nodes        9840
DONE!   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   [ time: .011s ]
Connect Conforming inner and periodic sides...
    -->          368  sides of          368   sides connected.
Mesh Connect completed with success.  
DONE!   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   [ time: .011s ]
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
BUILD EDGES ...
DONE!   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   [ time: .000s ]
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
GLOBAL UNIQUE NODES ...
   All Nodes sorted...
   Number of nodes to check:         9840
    size of tolerance box:                    1
  Number of deleted nodes           0
  Number of unique nodes        9840
DONE!   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   [ time: .010s ]
 =======================CHECK NODE CONNECTIVITY =========================================
###### CHECK CORNER NODES
###### CHECK EDGE CURVED NODES
###### CHECK SURFACE CURVED NODES
SORT ELEMENTS ON SPACE FILLING CURVE, TYPE hilbert ...
... DONE
 write space filling curve to sfc.dat
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
WRITING THE DEBUGMESH...
   #Elements           64
   WRITE DATA TO VTX XML BINARY (VTU) FILE... CYLINDER_Ngeo5_Debugmesh.vtu   DONE
    Mesh visualized for debug purposes in file : CYLINDER_Ngeo5_Debugmesh
DONE!   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   [ time: .000s ]
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
WRITING THE BC MESH...
   #BCSides          144
   WRITE DATA TO VTX XML BINARY (VTU) FILE... CYLINDER_Ngeo5_Debugmesh_BC.vtu   DONE
    Boundary mesh visualized for debug purposes in file : CYLINDER_Ngeo5_Debugmesh_BC
DONE!   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   [ time: .000s ]
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
WRITE CURVED SURFACE VISUALIZATION...
Nothing to visualize: No Boundaries with CurveIndex>0 found.
DONE!   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   [ time: .000s ]
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Check Elements Jacobian...
 Number of element with scaled Jacobians ranging between:
   <  0.0  <  0.1  <  0.2  <  0.3  <  0.4  <  0.5  <  0.6  <  0.7  <  0.8  <  0.9  <  1.0 
     0 |     0 |     0 |     0 |     0 |    16 |    16 |     8 |     8 |     8 |     8 | 
DONE!   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   [ time: .041s ]
====================================================================================================================================
 
 
                                        GOT mesh (incl. rasterfahndung) 
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
 WRITE DATA TO HDF5 FILE...
 Mesh statistics:
 Element Type | number of elements
 104        |        0
 204        |        0
 105        |        0
 115        |        0
 205        |        0
 106        |        0
 116        |        0
 206        |        0
 108        |        0
 118        |        0
 208        |       64
DONE!   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   [ time: .003s ]
DONE!   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   .   [ time: .105s ]
====================================================================================================================================
 HOPR successfully finished: Mesh "CYLINDER_Ngeo5_mesh.h5" written to HDF5 file.
====================================================================================================================================