Commit faf089e7 authored by Steffen Schotthöfer's avatar Steffen Schotthöfer
Browse files

restructured test case folder, restructured validation tests

parent b3a2a352
import pygmsh as pg
import numpy as np
import itertools
import os
def add_block(x0,y0,length,char_length,geom):
coords = np.array([
[x0, y0, 0.0],
[x0+length, y0, 0.0],
[x0+length, y0+length, 0.0],
[x0, y0+length, 0.0]
])
return geom.add_polygon(coords, char_length)
char_length = 0.4
geom = pg.opencascade.Geometry()
domain = add_block(0, 0, 7, char_length, geom)
xpos = ypos = [1, 2, 3, 4, 5]
pattern = list(itertools.product(xpos, ypos))[::2]
pattern.pop(7)
boxes = [domain]
for pos in pattern:
boxes.append(add_block(pos[0], pos[1], 1, char_length, geom))
geom.boolean_fragments(boxes,[])
geom.add_physical(domain.lines, label="void")
mesh_code = geom.get_code()
with open("checkerboard.geo","w") as mesh_file:
mesh_file.write(mesh_code)
os.system('gmsh checkerboard.geo -2 -format su2 -save_all')
os.system('rm checkerboard.geo')
This diff is collapsed.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Checkerboard Benchmarking File MN %
% Author <Steffen Schotthöfer> %
% Date 10.11.2020 %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% ---- File specifications ----
%
OUTPUT_DIR = ../../tests/result
OUTPUT_FILE = rtsn_test_checkerboard_MN
LOG_DIR = ../../tests/result/logs
MESH_FILE = checkerboard.su2
%
% ---- Problem specifications ----
%
PROBLEM = CHECKERBOARD
%
% ---- Solver specifications ----
%
CFL_NUMBER = 0.5
TIME_FINAL = 3.2
SOLVER = MN_SOLVER
MAX_MOMENT_SOLVER = 2
ENTROPY_FUNCTIONAL = MAXWELL_BOLTZMANN
ENTROPY_OPTIMIZER = NEWTON
%
% ----- Newton Solver Specifications ----
%
NEWTON_FAST_MODE = NO
NEWTON_ITER = 1000
NEWTON_EPSILON = 0.01
NEWTON_STEP_SIZE = 0.7
NEWTON_LINE_SEARCH_ITER = 1000
%
% ---- Boundary Conditions ----
%
BC_DIRICHLET = ( void )
%
% ---- Quadrature ----
%
QUAD_TYPE = GAUSS_LEGENDRE_TENSORIZED
QUAD_ORDER = 8
%
% ----- Output ----
%
VOLUME_OUTPUT = (MINIMAL)
This diff is collapsed.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Checkerboard Benchmarking File PN %
% Author <Steffen Schotthöfer> %
% Date 10.11.2020 %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% ---- File specifications ----
%
OUTPUT_DIR = ../../tests/result
OUTPUT_FILE = rtsn_test_checkerboard_PN
LOG_DIR = ../../tests/result/logs
MESH_FILE = checkerboard.su2
%
% ---- Problem specifications ----
%
PROBLEM = CHECKERBOARD
%
% ---- Solver specifications ----
%
CFL_NUMBER = 0.5
TIME_FINAL = 3.2
SOLVER = PN_SOLVER
MAX_MOMENT_SOLVER = 2
%
% ---- Boundary Conditions ----
%
BC_NEUMANN = ( void )
%
% ----- Output ----
%
VOLUME_OUTPUT = (MINIMAL)
This diff is collapsed.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Checkerboard Benchmarking File SN %
% Author <Steffen Schotthöfer> %
% Date 10.11.2020 %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% ---- File specifications ----
%
OUTPUT_DIR = ../../tests/result
OUTPUT_FILE = rtsn_test_checkerboard_SN
LOG_DIR = ../../tests/result/logs
MESH_FILE = checkerboard.su2
%
% ---- Problem specifications ----
%
PROBLEM = CHECKERBOARD
%
% ---- Solver specifications ----
%
CFL_NUMBER = 0.5
TIME_FINAL = 3.2
SOLVER = SN_SOLVER
%
% ---- Boundary Conditions ----
%
BC_NEUMANN = ( void )
%
% ---- Quadrature ----
%
QUAD_TYPE = LEBEDEV
QUAD_ORDER = 15
%
% ----- Output ----
%
VOLUME_OUTPUT = (MINIMAL)
This diff is collapsed.
NDIME= 2
NELEM= 96
5 38 48 32 0
5 34 47 30 1
5 42 45 31 2
5 36 46 27 3
5 44 51 50 4
5 38 60 48 5
5 34 59 47 6
5 39 52 51 7
5 27 58 36 8
5 31 57 42 9
5 50 51 4 10
5 41 50 49 11
5 51 52 4 12
5 49 52 35 13
5 6 33 30 14
5 33 35 34 15
5 8 28 7 16
5 11 40 27 17
5 13 29 12 18
5 18 25 17 19
5 23 26 22 20
5 37 38 32 21
5 16 43 31 22
5 21 37 32 23
5 33 34 30 24
5 37 39 38 25
5 11 27 10 26
5 6 30 5 27
5 7 33 6 28
5 16 31 15 29
5 21 32 20 30
5 25 44 43 31
5 26 39 37 32
5 28 33 7 33
5 25 43 17 34
5 26 37 22 35
5 28 35 33 36
5 29 41 40 37
5 12 40 11 38
5 29 40 12 39
5 17 43 16 40
5 22 37 21 41
5 28 36 35 42
5 29 42 41 43
5 49 50 4 44
5 4 52 49 45
5 14 45 13 46
5 9 46 8 47
5 19 48 18 48
5 24 47 23 49
5 28 46 36 50
5 29 45 42 51
5 13 45 29 52
5 23 47 26 53
5 8 46 28 54
5 18 48 25 55
5 5 54 0 56
5 10 53 1 57
5 15 55 2 58
5 20 56 3 59
5 0 54 24 60
5 1 53 9 61
5 2 55 14 62
5 3 56 19 63
5 40 58 27 64
5 43 57 31 65
5 35 52 34 66
5 39 51 38 67
5 48 60 25 68
5 47 59 26 69
5 50 57 44 70
5 49 58 41 71
5 36 49 35 72
5 42 50 41 73
5 14 55 45 74
5 9 53 46 75
5 19 56 48 76
5 24 54 47 77
5 30 54 5 78
5 27 53 10 79
5 31 55 15 80
5 32 56 20 81
5 25 60 44 82
5 26 59 39 83
5 41 58 40 84
5 44 57 43 85
5 46 53 27 86
5 47 54 30 87
5 45 55 31 88
5 48 56 32 89
5 42 57 50 90
5 36 58 49 91
5 44 60 51 92
5 39 59 52 93
5 52 59 34 94
5 51 60 38 95
NPOIN= 61
-0.6 -0.6 0
0.6 -0.6 1
0.6 0.6 2
-0.6 0.6 3
0 0 4
-0.3999999999999999 -0.6 5
-0.1999999999999995 -0.6 6
6.661338147750939e-16 -0.6 7
0.2000000000000013 -0.6 8
0.4000000000000014 -0.6 9
0.6 -0.3999999999999999 10
0.6 -0.1999999999999995 11
0.6 6.661338147750939e-16 12
0.6 0.2000000000000013 13
0.6 0.4000000000000014 14
0.3999999999999999 0.6 15
0.1999999999999995 0.6 16
-6.661338147750939e-16 0.6 17
-0.2000000000000013 0.6 18
-0.4000000000000014 0.6 19
-0.6 0.3999999999999999 20
-0.6 0.1999999999999995 21
-0.6 -6.661338147750939e-16 22
-0.6 -0.2000000000000013 23
-0.6 -0.4000000000000014 24
-0.0883323991063687 0.4335312117598478 25
-0.4335312117598478 -0.08833239910636875 26
0.408417471068293 -0.2827770063845272 27
0.09836539647744563 -0.4277386580268777 28
0.4277386580268777 0.0983653964774456 29
-0.2991525261438244 -0.4306455909888335 30
0.2827770063845272 0.408417471068293 31
-0.4306455909888335 0.2991525261438244 32
-0.1001311882777293 -0.4275939876646933 33
-0.2078135244157358 -0.2745810407467589 34
-0.006062372576521488 -0.2693268128117099 35
0.1965430057863219 -0.279816127364245 36
-0.4285594132868549 0.1018033545062421 37
-0.2747419450171192 0.2080922187871546 38
-0.2762725138155953 0.01714894082566772 39
0.4296818547442399 -0.08736893761343628 40
0.2535898384862244 9.351322400109502e-16 41
0.2771932983099974 0.1975534012157423 42
0.0890411038419491 0.4306472803664015 43
-9.361119987333181e-16 0.2535898384862243 44
0.4402164668679376 0.2929876845412847 45
0.2928192853030479 -0.4406536050436456 46
-0.4176849844018014 -0.2856621985635512 47
-0.2856786173666491 0.4177134226029665 48
0.1000000000000012 -0.1267949192431118 49
0.1267949192431118 0.1000000000000011 50
-0.1038469482946126 0.122107485866346 51
-0.1228819694783086 -0.1031192711093941 52
0.4535898384862245 -0.4535898384862241 53
-0.4535898384862241 -0.4535898384862245 54
0.4535898384862241 0.4535898384862245 55
-0.4535898384862245 0.4535898384862241 56
0.1476383677033232 0.2830882411612322 57
0.2828866267682136 -0.1477353741141124 58
-0.2883451608840447 -0.1398020157735396 59
-0.1398020157735396 0.2883451608840447 60
NMARK= 1
MARKER_TAG= void
MARKER_ELEMS= 24
3 0 5
3 5 6
3 6 7
3 7 8
3 8 9
3 9 1
3 1 10
3 10 11
3 11 12
3 12 13
3 13 14
3 14 2
3 2 15
3 15 16
3 16 17
3 17 18
3 18 19
3 19 3
3 3 20
3 20 21
3 21 22
3 22 23
3 23 24
3 24 0
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Linesource Benchmarking File MN %
% Author <Steffen Schotthöfer> %
% Date 01.07.2020 %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
% ---- File specifications ----
%
OUTPUT_DIR = ../../tests/result
OUTPUT_FILE = rtsn_test_linesource_MN_MB
LOG_DIR = ../../tests/result/logs
MESH_FILE = linesource.su2
%
% ---- Problem specifications ----
%
PROBLEM = LINESOURCE
SCATTER_COEFF = 1
%
% ---- Solver specifications ----
%
CFL_NUMBER = 0.7
TIME_FINAL = 0.2
SOLVER = MN_SOLVER
MAX_MOMENT_SOLVER = 2
ENTROPY_FUNCTIONAL = MAXWELL_BOLTZMANN
ENTROPY_OPTIMIZER = NEWTON
%
% ----- Newton Solver Specifications ----
%
NEWTON_FAST_MODE = NO
NEWTON_ITER = 1000
NEWTON_EPSILON = 0.01
NEWTON_STEP_SIZE = 0.7
NEWTON_LINE_SEARCH_ITER = 1000
%
% ---- Boundary Conditions ----
%
BC_DIRICHLET = ( void )
%
% ---- Quadrature ----
%
QUAD_TYPE = GAUSS_LEGENDRE_TENSORIZED
QUAD_ORDER = 8
%
% ----- Output ----
%
VOLUME_OUTPUT = (MINIMAL)
# vtk DataFile Version 4.2
vtk output
ASCII
DATASET UNSTRUCTURED_GRID
POINTS 61 double
-0.6 -0.6 0 0.6 -0.6 0 0.6 0.6 0
-0.6 0.6 0 0 0 0 -0.4 -0.6 0
-0.2 -0.6 0 6.6613381478e-16 -0.6 0 0.2 -0.6 0
0.4 -0.6 0 0.6 -0.4 0 0.6 -0.2 0
0.6 6.6613381478e-16 0 0.6 0.2 0 0.6 0.4 0
0.4 0.6 0 0.2 0.6 0 -6.6613381478e-16 0.6 0
-0.2 0.6 0 -0.4 0.6 0 -0.6 0.4 0
-0.6 0.2 0 -0.6 -6.6613381478e-16 0 -0.6 -0.2 0
-0.6 -0.4 0 -0.088332399106 0.43353121176 0 -0.43353121176 -0.088332399106 0
0.40841747107 -0.28277700638 0 0.098365396477 -0.42773865803 0 0.42773865803 0.098365396477 0
-0.29915252614 -0.43064559099 0 0.28277700638 0.40841747107 0 -0.43064559099 0.29915252614 0
-0.10013118828 -0.42759398766 0 -0.20781352442 -0.27458104075 0 -0.0060623725765 -0.26932681281 0
0.19654300579 -0.27981612736 0 -0.42855941329 0.10180335451 0 -0.27474194502 0.20809221879 0
-0.27627251382 0.017148940826 0 0.42968185474 -0.087368937613 0 0.25358983849 9.3513224001e-16 0
0.27719329831 0.19755340122 0 0.089041103842 0.43064728037 0 -9.3611199873e-16 0.25358983849 0
0.44021646687 0.29298768454 0 0.2928192853 -0.44065360504 0 -0.4176849844 -0.28566219856 0
-0.28567861737 0.4177134226 0 0.1 -0.12679491924 0 0.12679491924 0.1 0
-0.10384694829 0.12210748587 0 -0.12288196948 -0.10311927111 0 0.45358983849 -0.45358983849 0
-0.45358983849 -0.45358983849 0 0.45358983849 0.45358983849 0 -0.45358983849 0.45358983849 0
0.1476383677 0.28308824116 0 0.28288662677 -0.14773537411 0 -0.28834516088 -0.13980201577 0
-0.13980201577 0.28834516088 0
CELLS 96 384
3 38 48 32
3 34 47 30
3 42 45 31
3 36 46 27
3 44 51 50
3 38 60 48
3 34 59 47
3 39 52 51
3 27 58 36
3 31 57 42
3 50 51 4
3 41 50 49
3 51 52 4
3 49 52 35
3 6 33 30
3 33 35 34
3 8 28 7
3 11 40 27
3 13 29 12
3 18 25 17
3 23 26 22
3 37 38 32
3 16 43 31
3 21 37 32
3 33 34 30
3 37 39 38
3 11 27 10
3 6 30 5
3 7 33 6
3 16 31 15
3 21 32 20
3 25 44 43
3 26 39 37
3 28 33 7
3 25 43 17
3 26 37 22
3 28 35 33
3 29 41 40
3 12 40 11
3 29 40 12
3 17 43 16
3 22 37 21
3 28 36 35
3 29 42 41
3 49 50 4
3 4 52 49
3 14 45 13
3 9 46 8
3 19 48 18
3 24 47 23
3 28 46 36
3 29 45 42
3 13 45 29
3 23 47 26
3 8 46 28
3 18 48 25
3 5 54 0
3 10 53 1
3 15 55 2
3 20 56 3
3 0 54 24
3 1 53 9
3 2 55 14
3 3 56 19
3 40 58 27
3 43 57 31
3 35 52 34
3 39 51 38
3 48 60 25
3 47 59 26
3 50 57 44
3 49 58 41
3 36 49 35
3 42 50 41
3 14 55 45
3 9 53 46
3 19 56 48
3 24 54 47
3 30 54 5
3 27 53 10
3 31 55 15
3 32 56 20
3 25 60 44
3 26 59 39
3 41 58 40
3 44 57 43
3 46 53 27
3 47 54 30
3 45 55 31
3 48 56 32
3 42 57 50
3 36 58 49
3 44 60 51
3 39 59 52
3 52 59 34
3 51 60 38
CELL_TYPES 96
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
5
CELL_DATA 96
FIELD FieldData 1
radiation%20flux%20density 1 96 double
0.023075686389 0.023075686389 0.023075686389 0.023075686389 3.2567884905 0.023075686389 0.023075686389 2.9000325329 0.023075686389
0.023075686389 4.0517913501 3.1542357631 3.9630162161 2.7957825828 0.023075686389 0.023075686389 5.0580418059e-101 0.023075686389
5.0580418059e-101 1.6357402065e-101 1.6357402065e-101 0.023075686389 0.023075686389 0.023075686389 0.023075686389 0.023075686389 1.0925488099e-124
1.67533772e-128 4.807237643e-101 1.0925488099e-124 1.67533772e-128 0.023075686389 0.023075686389 0.023075686389 0.023075686389 0.023075686389
0.023075686389 0.023075686389 5.1070999977e-101 0.023075686389 3.5741746143e-101 3.353049222e-101 0.023075686389 0.023075686389 3.9616384051
3.9219732906 2.8869513693e-129 2.6164279107e-129 5.1807883453e-126 5.2355333002e-126 0.023075686389 0.023075686389 0.023075686389 0.023075686389
0.023075686389 0.023075686389 1.609494535e-180 1.609494535e-180 1.609494535e-180 1.609494535e-180 1.609494535e-180 1.609494535e-180 1.609494535e-180
1.609494535e-180 0.023075686389 0.023075686389 0.76456256015 0.79934508806 0.023075686389 0.023075686389 1.48065084 1.2673622482
0.65568641859 0.81185421685 0.023075686389 0.023075686389 0.023075686389 0.023075686389 0.023075686389 0.023075686389 0.023075686389
0.023075686389 0.023075686389 0.023075686389 0.023075686389 0.023075686389 0.023075686389 0.023075686389 0.023075686389 0.023075686389
0.023075686389 0.023075686389 1.2625328905 0.94584263198 0.023075686389 0.023075686389
POINT_DATA 61
FIELD FieldData 1
radiation%20flux%20density 1 61 double
1.609494535e-180 1.609494535e-180 1.609494535e-180 1.609494535e-180 3.9746048155 0.007691895463 0.007691895463 0.007691895463 0.007691895463
0.007691895463 0.007691895463 0.007691895463 0.007691895463 0.007691895463 0.007691895463 0.007691895463 0.007691895463 0.007691895463
0.007691895463 0.007691895463 0.007691895463 0.007691895463 0.007691895463 0.007691895463 0.007691895463 0.019229738658 0.019229738658