DIVEMesh :: User Guide
DIVEMesh 18.12
Contact
Hans Bihs
Associate Professor
Marine Civil Engineering
NTNU Trondheim
[Link]@[Link]
Team
Arun Kamath
Mayilvahanan Alagan Chella
Nadeem Ahmad
Ankit Aggarwal
Weizhi Wang
Tobias Martin
Petter Fornes
ii
Contents
1 The ‘[Link]‘ file 1
1.1 B :: Boundary . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2 C :: Channel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.3 D :: Data Interpolation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.4 F :: Freecoor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.5 G :: Geodat . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.6 M :: MPI . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1.7 O :: Objects . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
1.8 S :: Solid . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
iii
iv
1
Chapter 1
The ‘[Link]‘ file
1.1 B :: Boundary
B 1 double cell size dx
default: 0.0
B 2 int number of cells in x-, y- and z-direction
When activated, it will overwrite B 1.
default: 0 0 0
B 10 * double rectangular domain; xstar t , xend , ystar t , yend , zstar t , zend
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
B 22 * double sphere; xcenter , ycenter , zcenter , r adi us
default: 0.0 ; 0.0 ; 0.0 ; 0.0
B 31 * double straight pipe in x-direction; xstar t , xend , ycenter , zcenter , r adi us
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
B 32 * double straight pipe in y-direction; ystar t , y xend , xcenter , zcenter , r adi us
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
B 33 * double straight pipe in z-direction; zstar t , zend , xcenter , ycenter , r adi us
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
B 101 int type of grid stretching function in x-direction
CHAPTER 1. THE ‘[Link]‘ FILE 2
0 OFF
1 center focus (sinh, using B113)
2 wall focus (tanh, using B113)
5 point focus using sinh (requires B 114)
6 point focus using exponents (requires B 114)
8 three zones with constant dx (requires B 121)
9 three zones with stretching (requires B 124)
default: 0
B 102 int type of grid stretching function in y-direction
0 OFF
1 center focus (sinh, using B113)
2 wall focus (tanh, using B113)
5 point focus using sinh (requires B 115)
6 point focus using exponents (requires B 115)
8 three zones with constant dx (requires B 122)
9 three zones with stretching (requires B 125)
default: 0
B 103 int type of grid stretching function in z-direction
0 OFF
1 center focus (sinh, using B113)
2 wall focus (tanh, using B113)
3 lid focus (expontential, using B113)
4 lid focus (sinusoidal)
5 point focus using sinh (requires B 116)
6 point focus using exponents (requires B 116)
8 three zones with constant dx (requires B 123)
9 three zones with stretching (requires B 126)
default: 0
B 111 double stretching factor in x-direction
default: 1.0
B 112 double stretching factor in y-direction
default: 1.0
CHAPTER 1. THE ‘[Link]‘ FILE 3
B 113 double stretching factor in z-direction
default: 1.0
B 114 double focus point for x-direction stretching B 101 5
default: 0.0
B 115 double focus point for y-direction stretching B 102 5
default: 0.0
B 116 double focus point for z-direction stretching B 103 5
default: 0.0
B 121 input for the three zones of B 101 8: int N1, double x1, int N2, double x2, int N3
The sum of the given N1,N2,N3 result in the number of elements in x-direction.
default: na
B 122 input for the three zones of B 101 8: int N1, double y1, int N2, double y2, int N3
The sum of the given N1,N2,N3 result in the number of elements in y-direction.
default: na
B 123 input for the three zones of B 101 8: int N1, double z1, int N2, double z2, int N3
The sum of the given N1,N2,N3 result in the number of elements in z-direction.
default: na
B 124 input for the three zones of B 101 9: int N1, double x1, double f1, int N2, double xf, double
f2, int N3, double x2, double f3
The sum of the given N1,N2,N3 result in the number of elements in x-direction. The fac-
tors f1 and f3 are the linear stretching factors for the layers. The parameter x1 gives the border
between the first and second layer, x2 between second and third. The parameter xf gives the
location of the stretching focus in the second layer, where sinh stretching based on the factor f2
is used.
default: na
CHAPTER 1. THE ‘[Link]‘ FILE 4
B 125 input for the three zones of B 102 9: int N1, double y1, double f1, int N2, double yf, double
f2, int N3, double y2, double f3
see B 124.
default: na
B 126 input for the three zones of B 103 9: int N1, double z1, double f1, int N2, double zf, double
f2, int N3, double z2, double f3
see B 124.
default: na
1.2 C :: Channel
2
1 4
3
5
z
Figure 1.1: Definition of cell sides within DIVEMesh and REEF3D.
C 11 int Boundary Condition on Surfside 1
1 inflow
3 symmetry plane
6 wave generation
7 numerical beach
21 wall
default: 21
C 12 int Boundary Condition on Surfside 2
3 symmetry plane
CHAPTER 1. THE ‘[Link]‘ FILE 5
6 wave generation
7 numerical beach
8 active beach
21 wall
default: 21
C 13 int Boundary Condition on Surfside 3
3 symmetry plane
6 wave generation
7 numerical beach
8 active beach
21 wall
default: 21
C 14 int Boundary Condition on Surfside 4
2 outflow
3 symmetry plane
6 wave generation
7 numerical beach
8 active beach
21 wall
default: 21
C 15 int Boundary Condition on Surfside 5
3 symmetry plane
21 wall
default: 21
C 16 int Boundary Condition on Surfside 6
3 symmetry plane
21 wall
default: 21
CHAPTER 1. THE ‘[Link]‘ FILE 6
1.3 D :: Data Interpolation
D 10 int turn data interpolation on/off
0 OFF
1 ON
default: 0
D 11 double ∆x, ∆y, ∆x
default: 0.0 ; 0.0 ; 0.0
D 12 double factor x-coordinate, factor y-coordinate, factor z-coordinate,
default: 0.0 ; 0.0 ; 0.0
D 13 int Read data every i th iteration
default: 1
D 14 int Type of interpolation
1 inverse distance
2 kriging
default: 1
D 15 int Number of smoothing iterations
default: 0
D 16 double Factor for data smoothing
default: 0.5
D 17 double Factor for data inverse distance
default: 16.0
D 18 double Range factor for kriging
default: 0.3
D 19 int Read dummy letter in front of 3D coordinates
0 OFF
1 ON
CHAPTER 1. THE ‘[Link]‘ FILE 7
default: 0
D 23 int reverse sign of data variable in [Link]
0 OFF
1 ON
default: 0
D 24 double distance limiter for inverse distance
default: na
1.4 F :: Freecoor
F 3 double zstart, zend
default: 0.0
F 11 * double left side wall LINE; xstar t , xend , ystar t , yend
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ;
F 21 * double right side wall LINE; xstar t , xend , ystar t , yend
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ;
1.5 G :: Geodat
G 10 int turn geodat on/off
0 OFF
1 ON
default: 0
G 11 double ∆x, ∆y, ∆x
default: 0.0 ; 0.0 ; 0.0
G 12 double factor x-coordinate, factor y-coordinate, factor z-coordinate,
default: 0.0 ; 0.0 ; 0.0
G 13 double rotation angle of geo coordinates around vertical axis
CHAPTER 1. THE ‘[Link]‘ FILE 8
default: 0.0
G 14 double x-coordinate and y-coordinate of origin for the rotation angle of geo coordinates around
vertical axis
default: 0.0 ; 0.0
G 15 int interpolation scheme
1 global inverse distance interpolation
2 local inverse distance interpolation
default: 2
G 17 int dij lower limit for local inverse distance interpolation
default:
G 19 int read a letter in front of the coordinates in the [Link] file
0 OFF
1 ON
default: 0
G 20 int use automatic grid size
0 OFF
1 ON
default: 0
G 21 double Margins for automatic grid size xstar t , ystar t , zstar t
default: 0.0 ; 0.0 ; 0.0
G 22 double Margins for automatic grid size xend , yend , zend
default: 0.0 ; 0.0 ; 0.0
G 23 int reverse sign of vertical coordinate in [Link] file
0 OFF
1 ON
default: 0
G 24 double raise topography above the level h by dz
CHAPTER 1. THE ‘[Link]‘ FILE 9
default: 0.0 ; 0.0
G 25 double multiply topography above the level h by factor f z
default: 0.0 ; 0.0
G 31 int Number of smoothing iterations
default: 0
G 32 double Factor for topography data smoothing
default: 0.5
G 35 double Factor for data inverse distance
default: 16.0
G 41 int print SWAN bottom file from interpolated geo points
0 OFF
1 ON
default: 0
1.6 M :: MPI
M 10 int Number of processes
default: 1
M 11 int Partition in x-direction
0 OFF
1 ON
default: 1
M 12 int Partition in y-direction
0 OFF
1 ON
default: 1
M 13 int Partition in z-direction
CHAPTER 1. THE ‘[Link]‘ FILE 10
0 OFF
1 ON
default: 0
M 20 int Decomposition method
1 Standard rectangular base domain
2 Improved rectangular base domain
default: 1
M 31 double Variance parameter for decomposition method 2
default: 0.1
1.7 O :: Objects
O 10 * double rectangular object; xstar t , xend , ystar t , yend , zstar t , zend
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
O 32 * double cylinder in y-direction; xcenter , zcenter , r adi us
default: 0.0 ; 0.0 ; 0.0
O 33 * double cylinder in z-direction; xcenter , ycenter , r adi us
default: 0.0 ; 0.0 ; 0.0
O 61 * double wedge object; xstar t , xend , ystar t , yend , zstar t , zend
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
1.8 S :: Solid
S 1 int read STL file and generate solid
IMPORTANT: the STL file needs to be in ASCII format! Most CAD programs export to a
binary STL file. It is possible to use e.g. Paraview for STL conversion from binary to ASCII.
0 OFF
1 ON
default: 0
CHAPTER 1. THE ‘[Link]‘ FILE 11
S 2 int auto-generate mesh extend based on STL max/min coordinates
0 OFF
1 ON
default: 0
S 3 double margins for auto-generated mesh xmstar t , xmend , y mstar t , y mend , zmstar t , zmend
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
S 4 double Scale STL geometry
default: 1.0
S 5 double 3D rotation; xor igin , yor igin , zor igin , φ, θ, ψ
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
S 6 int print transformed STL model to ”REEF3D_Solid.stl”
0 OFF
1 ON
default: 0
S 7 double change origin of STL model dx, dy , dz
default: 0.0 ; 0.0 ; 0.0
S 8 double Turn STL geometry in horizontal xy-plane in degree o
default: 0.0
S 9 int Invert inside/outside for STL geometry
1 regular
2 invert
default: 1
S 10 * double rectangular object; xstar t , xend , ystar t , yend , zstar t , zend
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
S 11 * double rectangular object array; xor igin , yor igin , zor igin , box length L, gap G, number of objects
in each direction ni , nj , nk
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
CHAPTER 1. THE ‘[Link]‘ FILE 12
S 32 * double cylinder in y-direction; xcenter , zcenter , r adi us
default: 0.0 ; 0.0 ; 0.0
S 33 * double cylinder in z-direction; xcenter , ycenter , r adi us
default: 0.0 ; 0.0 ; 0.0
S 34 * double cylinder with flexible orientation and front face orthogonal to x-plane
xstar t , ystar t , zstar t , r adi usstar t , xend , yend , zend , r adi usend
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
S 35 * double cylinder with flexible orientation and front face orthogonal to y-plane
xstar t , ystar t , zstar t , r adi usstar t , xend , yend , zend , r adi usend
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
S 36 * double cylinder with flexible orientation and front face orthogonal to z-plane
xstar t , ystar t , zstar t , r adi usstar t , xend , yend , zend , r adi usend
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
S 37 * double cylinder with flexible orientation and front face orthogonal to the cylinder axis
xstar t , ystar t , zstar t , r adi usstar t , xend , yend , zend , r adi usend
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
S 38 * double cylinder with flexible orientation and a vertical front face in line with the cylinder axis
xstar t , ystar t , zstar t , r adi usstar t , xend , yend , zend , r adi usend
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
S 41 * double cone in x-direction; ycenter , zcenter , xstar t , xend , r adi usstar t , r adi usend
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
S 42 * double cone in y-direction; xcenter , zcenter , ystar t , yend , r adi usstar t , r adi usend
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
S 43 * double cone in z-direction; xcenter , ycenter , zstar t , zend , r adi usstar t , r adi usend
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
S 51 * double sphere; xcenter , ycenter , zcenter , r adi us
CHAPTER 1. THE ‘[Link]‘ FILE 13
default: 0.0 ; 0.0 ; 0.0 ; 0.0
S 52 * double ellipsoid; xcenter , ycenter , zcenter , aaxis , baxis , caxis
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
S 53 * double semi ellipsoid with vertical base; xcenter , ycenter , zcenter , aaxis , baxis , caxis , hbase
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
S 54 * double semi ellipsoid with vertical base with rotation around the center;
xcenter , ycenter , zcenter , aaxis , baxis , caxis , hbase , φ, θ, ψ
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
S 61 * double wedge object in x-direction; xstar t , xend , ystar t , yend , zstar t , zend
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
S 62 * double wedge object in y-direction; xstar t , xend , ystar t , yend , zstar t , zend
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
S 63 * double wedge object in z-direction; xstar t , xend , ystar t , yend , zstar t , zend
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
S 81 * double tetrahedon object, each of the 4 points is given by the coordinates x1 , y1 , z1 , x2 , y2 , z2 , x3 , y3 , z3 , x4 , y4 , z4
default: [4x] 0.0 ; 0.0 ; 0.0
1 2
Figure 1.2: Definition of the tetrahedon points.
CHAPTER 1. THE ‘[Link]‘ FILE 14
S 82 * double pyramid object, each of the 5 points is given by the coordinates x1 , y1 , z1 , x2 , y2 , z2 , x3 , y3 , z3 , x4 , y4 , z4 ,
x5 , y5 , z5
default: [5x] 0.0 ; 0.0 ; 0.0
4
3
1 2
Figure 1.3: Definition of the pyramid points.
S 83 * double wedge object, each of the 6 points is given by the coordinates x1 , y1 , z1 , x2 , y2 , z2 , x3 , y3 , z3 , x4 , y4 , z4 ,
x5 , y5 , z5 , x6 , y6 , z6
default: [6x] 0.0 ; 0.0 ; 0.0
4
5
1 2
Figure 1.4: Definition of the wedge points.
S 84 * double hexahedron object, each of the 8 points is given by the coordinates x1 , y1 , z1 , x2 , y2 , z2 , x3 , y3 , z3 , x4 , y4 , z4 ,
x5 , y5 , z5 , x6 , y6 , z6 , x7 , y7 , z7 , x8 , y8 , z8
default: [8x] 0.0 ; 0.0 ; 0.0
CHAPTER 1. THE ‘[Link]‘ FILE 15
8
7
5
6
4
3
1 2
Figure 1.5: Definition of the hexahedron points.
S 121 * double vertical ogee weir, coordinates of upstream bottom corner, width, downstream height
and hydraulic head x, y , z, b, Pd , H0
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
S 122 * double vertical ogee weir, K, n, xc and yc
default: 0.5 ; 1.85 ; 0.22 ; 0.075
S 123 * double vertical ogee weir, R1 and R2
default: 0.45 ; 0.2
S 131 * double semi-circular cylinder in y-direction, xc , zc , ystar t , yend , r adi usstar t , r adi usend
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
S 141 * double arch bridge; xstar t , xend , ystar t , yend , zstar t , zend , r adi us
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
S 201 * double thin plate x-normal; xlocation , ystar t , yend , zstar t , zend
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
S 202 * double thin plate y-normal; ylocation , xstar t , xend , zstar t , zend
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
S 203 * double thin plate z-normal; zlocation , xstar t , xend , ystar t , yend
default: 0.0 ; 0.0 ; 0.0 ; 0.0 ; 0.0
