hyperion/mesh_dry_run.py at main · festim-dev/hyperion · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
import gmsh


def generate_mesh(mesh_size=2e-4, fname="mesh_solid_only.msh"):
    # Initialize the GMSH API
    gmsh.initialize()
    gmsh.model.add("HYPERION_mesh_solid_only")

    # ##### key coordinates (m) #####
    x_in = 0.039  # inner vertical Ni wall (3.90 cm)
    x_out = 0.041  # outer vertical Ni wall (4.10 cm)

    y0 = 0.000  # bottom
    y_bT = 0.002  # 0.20 cm (bottom Ni top)
    y_mB = 0.022  # 2.20 cm (middle thin Ni bottom)
    y_mT = 0.024  # 2.40 cm (middle thin Ni top)
    # y_fT = 0.0319  # 3.19 cm (fluid top)
    y_tIn = 0.1091  # 10.91 cm (top Ni bottom / inner top)
    y_tOut = 0.1111  # 11.11 cm (top Ni outer top)

    # make rectangles
    gmsh.model.occ.addRectangle(0, y0, 0, x_in, y_bT - y0, tag=1)
    gmsh.model.occ.addRectangle(0, y_mB, 0, x_in, y_mT - y_mB, tag=2)
    gmsh.model.occ.addRectangle(0, y_tIn, 0, x_in, y_tOut - y_tIn, tag=3)
    gmsh.model.occ.addRectangle(x_in, y0, 0, x_out - x_in, y_tOut - y0, tag=4)
    # salt_rectangle = gmsh.model.occ.addRectangle(0, y_mT, 0, x_in, y_fT - y_mT, tag=5)

    # fuse all solid parts
    gmsh.model.occ.fuse([(2, 1), (2, 2), (2, 3)], [(2, 4)], tag=10)

    # synchronize CAD kernel with GMSH model
    gmsh.model.occ.synchronize()
    print(gmsh.model.getEntities(2))

    # mark volumes

    solid = gmsh.model.addPhysicalGroup(2, [10], tag=2)
    gmsh.model.setPhysicalName(2, solid, "solid")

    # mark boundaries
    print(gmsh.model.getEntities(1))
    out = gmsh.model.addPhysicalGroup(1, [1, 2, 3], tag=3)
    gmsh.model.setPhysicalName(1, out, "out")

    left_bc_top_Ni = gmsh.model.addPhysicalGroup(1, [4], tag=42)
    gmsh.model.setPhysicalName(1, left_bc_top_Ni, "left_bc_top_Ni")

    left_bc_middle_Ni = gmsh.model.addPhysicalGroup(1, [8], tag=43)
    gmsh.model.setPhysicalName(1, left_bc_middle_Ni, "left_bc_middle_Ni")

    left_bc_bottom_Ni = gmsh.model.addPhysicalGroup(1, [12], tag=44)
    gmsh.model.setPhysicalName(1, left_bc_bottom_Ni, "left_bc_bottom_Ni")

    top_cap_Ni = gmsh.model.addPhysicalGroup(1, [5], tag=5)
    gmsh.model.setPhysicalName(1, top_cap_Ni, "top_cap_Ni")

    top_sidewall_Ni = gmsh.model.addPhysicalGroup(1, [6], tag=6)
    gmsh.model.setPhysicalName(1, top_sidewall_Ni, "top_sidewall_Ni")

    bottom_sidewall_Ni = gmsh.model.addPhysicalGroup(1, [10], tag=7)
    gmsh.model.setPhysicalName(1, bottom_sidewall_Ni, "bottom_sidewall_Ni")

    mem_Ni_top = gmsh.model.addPhysicalGroup(1, [7], tag=8)
    gmsh.model.setPhysicalName(1, mem_Ni_top, "mem_Ni_top")

    mem_Ni_bottom = gmsh.model.addPhysicalGroup(1, [9], tag=9)
    gmsh.model.setPhysicalName(1, mem_Ni_bottom, "mem_Ni_bottom")

    bottom_cap_Ni = gmsh.model.addPhysicalGroup(1, [11], tag=10)
    gmsh.model.setPhysicalName(1, bottom_cap_Ni, "bottom_cap_Ni")

    set(gmsh.model.getBoundary([(2, 10)], oriented=False, recursive=False))

    # set CharacteristicLengthMax
    gmsh.option.setNumber("Mesh.CharacteristicLengthMax", mesh_size)

    gmsh.model.occ.synchronize()

    gmsh.model.mesh.generate(2)

    # write to file
    gmsh.write(fname)

    gmsh.finalize()


if __name__ == "__main__":
    generate_mesh()
    # gmsh.fltk.run()