Step02 Geometry Create

 

The data file for generation of the geometry is in: Geol_002b\Case01\02_Geometry_Create\Data. It consists of Geol_002b_Case01_abaq.inp and the Geol_002b_Case01_prop_IJK.spat files output during the Zmap import in the previous step. For convenience an edit in the data has been performed as shown below.

 

 

Data File


 

*PARAGEO CONVERSION, TYPE=GROUP, PARAGEO FILE=Geol_002b_Case01_geom

*CONVERSION OUTPUT, LEVEL=2, GEOMETRY=YES, MESH=YES

**

*NODE

  1,   33000.0000    ,   65000.0000    ,   5830.05847

  2,   34000.0000    ,   65000.0000    ,   5810.06664

     (...)

  9016,   74000.0000    ,   80000.0000    ,   6334.25155  

     (...)

**

*spatial_grid, INPUT=Geol_002b_Case01_prop_IJK.spat, OUTPUT=Geol_002b_Case01_prop_IJK_ReNum.spat

**

The ouput file name defined in PARAGEO FILE has been edited for convenience to "Geol_002b_Case01_geom".

 

Note the *SPATIAL_GRID keyword written at the end of the .inp file. The INPUT= command defines the .spat file created during Zmap import which contains the Spatial grid with the corresponding element properties. By including the spatial grid in the .inp file during geometry creation the re-numbering being applied to the mesh elements will also be applied to the elements defined in the Spatial_grid. The OUTPUT= command writes the new grid in a new file.

 

The Geol_002b_Case01_abaq.inp is run in ParaGeo to perform the conversion and obtain the results.

 

 

 

Results

 

The results for the present workflow steps are provided in Geol_002b\Case01\02_Geometry_Create\Results

 

The key results comprise:

 

1.Geol_002b_Case01_geom.geo file containing the geometry and mesh in ParaGeo format.

2.Geol_002b_Case01_geom.geometry containing the generated Group_data for the 10 formations and 26 Geometry_set data structures.

3.Geol_002b_Case01_prop_IJK_ReNum.spat containing the I, J, K spatial grid with the elements being renumbered consistently according to the generated geometry element numbers.

 

Data File


 

! Initial grid

* Spatial_grid                     NUM=1

! ----------------------------------------

 Name                     "grid_IJK"

 File_name                "Geol_002b_Case02_prop_IJK.spat"

 Type                     "Element"

 Element_variables  IDM=3

    "I"

    "J"

    "K"

 Cell_values  IDM=3  JDM=7661

   1.00000        1.00000        1.00000  

   2.00000        1.00000        1.00000

    (...)

   41.0000        15.0000        20.0000

 Element_numbers  IDM=7661

    1   2   3  ...   14093    14094

 

 

 

! Converted grid

* Spatial_grid                     NUM=2

! ----------------------------------------

 Name                     "grid_IJK"

 File_name                "Geol_002b_Case02_prop_IJK_ReNum.spat"

 Type                     "Element"

 Element_variables  IDM=3

    "I"

    "J"

    "K"

 Cell_values  IDM=3  JDM=7661

   1.00000        1.00000        1.00000

   2.00000        1.00000        1.00000

    (...)

   41.0000        15.0000        20.0000

 Element_numbers  IDM=7661

    1   2   3  ...   7660   7661

 

The differences in the IJK spatial grid before and after geometry creation (element renumbering) are highlighted in yellow.

 

 

The considered mesh contains 7661 elements and 9016 nodes and may be visualised by loading the file Geol_002b_Case01_abaq_mesh.xmf in ParaView.