SPectral Elements in Elastodynamics with Discontinuous Galerkin
Input files generation

In this section we described the flow chart used for the generation of input files FileName.mesh and LS.input.

Mesh file

To obtain the FileName.mesh in the format described in Mesh File we perform the following steps:

  • Generate with CUBIT the geometry of the problem (FileName.cub) and export it in the Exodus format (FileName.e)
  • Convert FileName.e into an ASCII file (FileName.txt)
  • Rewrite FileName.txt in the format read by SPEED (FileName.mesh)

How to generate the geometry of the problem

Here, we report the CUBIT commands needed to generate the computational domain for the tutorials Test SEm, Test SEm not-honoring and Test DG, respectively. To generate the mesh and export it is sufficiet to play the journal files TEST1_SEM.jou and TEST1_DG.jou in CUBIT. The journal file are contained in the folder TUTORIALS.

1_TEST_SEM.jou

     reset
     brick x 100 y 100 z 100
     move Volume 1  x 50 y 50 z -350 include_merged 
     brick x 100 y 100 z 200
     move Volume 2 x 50 y 50 z -200 include_merged 
     brick x 100 y 100 z 100
     move Volume 3 x 50 y 50 z -50 include_merged 
     merge all
     vol 1 2 3 size 50
     mesh vol 1 2 3 
     block 1 vol 3 
     block 2 vol 2
     block 3 vol 1
     block 4 surf 2
     block 5 surf 18 15 9 12 6 3 17 11 6 16 10 4  5
     draw block 1 2 3 4 5
     set large exodus file off
     export mesh "1_TEST_SEM.e" dimension 3 block 1 2 3 4 5 overwrite        

3_TEST_SEM_NH.jou

     reset
     brick x 100 y 100 z 100
     move Volume 1  x 50 y 50 z -350 include_merged 
     brick x 100 y 100 z 200
     move Volume 2 x 50 y 50 z -200 include_merged 
     brick x 100 y 100 z 100
     move Volume 3 x 50 y 50 z -50 include_merged 
     merge all
     vol 1 2 3 size 50
     mesh vol 1 2 3 
     block 1 vol 3 
     block 2 vol 2
     block 3 vol 1
     block 4 surf 2
     block 5 surf 18 15 9 12 6 3 17 11 6 16 10 4  5
     draw block 1 2 3 4 5
     set large exodus file off
     export mesh "3_TEST_SEM_NH.e" dimension 3 block 1 2 3 4 5 overwrite 
     create surf from surf 7 13
     surf 19 20 scheme trimesh size 25
     mesh surf 19 20
     set large exodus file off
     block 19 surf 19
     block 20 surf 20
     set large exodus file off
     export mesh "ALL.e" dimension 3 block 19 overwrite
     export mesh "XYZ.e" dimension 3 block 20 overwrite
     draw block 19 20

4_TEST_DG.jou

     reset
     brick x 100 y 100 z 100
     move Volume 1  x 50 y 50 z -350 include_merged 
     brick x 100 y 100 z 200
     move Volume 2 x 50 y 50 z -200 include_merged 
     brick x 100 y 100 z 100
     move Volume 3 x 50 y 50 z -50 include_merged 
     vol 1 2 size 50
     vol 3 size 25
     mesh vol 1 2 3 
     block 1 vol 3 
     block 2 vol 2
     block 3 vol 1
     block 4 surf 2
     block 5 surf 18 15 9 12 6 3 17 11 6 16 10 4  5
     block 6 surf 1
     block 7 surface 8
     block 8 surf 14
     block 9 surf 7
     draw block 1 2 3 4 5 6 7 8 9
     set large exodus file off
     export mesh "4_TEST_DG.e" dimension 3 block 1 2 3 4 5 6 7 8 9 overwrite    

FileName.e to FileName.txt

To convert all Exodus files into ASCII file it is sufficient to digit on a terminal

ncdump  FileName.e > FileName.txt 

on Unix/Linux OS or

ncdump.exe  FileName.e > FileName.txt 

on Windows OS.

Note
The executable ncdump.exe is provided along with the tutorials.

FileName.txt to FileName.mesh

To convert the ASCII file just obtained you have to run the MATLAB script MAKE_MESHinput.m contained in the folder TUTORIALS.

Note
Fill the part TO BE FILLED BY USER

ALL(XYZ).txt to ALL(XYZ).out

To convert the ASCII files ALL.txt and XYZ.txt needed by the not-honoring technique run the MATLAB scripts MAKE_ALLinput.m and MAKE_XYZinput.m, repsectively, contained in the folder TUTORIALS.

Note
Fill the part TO BE FILLED BY USER

LS.input file

To obtain the LS.input in the format described in Monitors File we perform the following steps:

  • Generate with CUBIT a surface of monitored points (LS.cub) and export it in the Exodus format (LS.e)
  • Convert LS.e into an ASCII file (LS.txt)
  • Rewrite LS.txt in the format read by SPEED

LS.input generation

Here, we report the CUBIT commands needed to generate a surface of montitored points to be used for the tutorials Test SEm, Test SEm not-honoring and Test DG, respectively. To generate and export the surface of monitored points it is sufficiet to play the journal file LS.jou in CUBIT. The journal file is contained in the folder TUTORIALS.

LS.jou

      create vertex 0 0 0
      create vertex 0 100 0
      create vertex 100 100 0
      create vertex 100 0 0
      create surf vertex 25 to 28
      draw surf 19
      surf 19 scheme trimesh size 25
      mesh surf 19
      block 10 surf 19
      undo group begin
      set large exodus file off
      export mesh "LS.e" dimension 3 block 10 overwrite
      undo group end    

From LS.e to LS.txt

To convert Exodus file into ASCII file il is sufficient to digit on a terminal

ncdump  LS.e > LS.txt 

on Unix/Linux OS or

ncdump.exe  LS.e > LS.txt 

on Windows OS.

Note
The executable ncdump.exe is provided along with the tutorials.

From LS.txt to LS.mesh

To convert the ASCII file run the MATLAB script MAKE_LSinput.m contained in the folder TUTORIALS.

Note
Fill the part TO BE FILLED BY USER