Dear all,
I would like to simulate flow in a "bathtube", espcially the outflow.
Therefore I wrote a grid in ucd format which I read using GridIn. A
problem arises when the channel of outflow has a cylindrical shape. I
applied the CylinderBoundary description to the grid. This doesn't work
properly where the cylinder is connected with the cube. Can someone
explain what is happening there and if I can fix this problem?
What I do looks like this:
GridIn<3> grid_in;
grid_in.attach_triangulation(triangulation);
std::ifstream input_file_3d ("tube.inp");
grid_in.read_ucd (input_file_3d);
const double radius = std::sqrt(2)/5;
const Point<3> direction = Point<3> (0,0,1);
const Point<3> point_on_axis = Point<3> (5,1,0);
static const CylinderBoundary<3> boundary_description(radius,
direction, point_on_axis);
triangulation.set_boundary(9, boundary_description);
Cheers,
Bärbel
64 82 0 0 0
1 0 0 0
2 2 0 0
3 4 0 0
4 6 0 0
5 0 2 0
6 2 2 0
7 4 2 0
8 6 2 0
9 0 0 1
10 2 0 1
11 4 0 1
12 6 0 1
13 0 2 1
14 2 2 1
15 4 2 1
16 6 2 1
17 0 0 2
18 2 0 2
19 4 0 2
20 6 0 2
21 0 2 2
22 2 2 2
23 4 2 2
24 6 2 2
25 4.4 0.4 1.5
26 5.6 0.4 1.5
27 4.4 1.6 1.5
28 5.6 1.6 1.5
29 4.8 0.8 0.5
30 5.2 0.8 0.5
31 4.8 1.2 0.5
32 5.2 1.2 0.5
33 4.4 0.4 1
34 5.6 0.4 1
35 4.4 1.6 1
36 5.6 1.6 1
37 4.8 0.8 0
38 5.2 0.8 0
39 4.8 1.2 0
40 5.2 1.2 0
41 4.4 0.4 0
42 5.6 0.4 0
43 4.4 1.6 0
44 5.6 1.6 0
45 4.9 0.9 0
46 5.1 0.9 0
47 4.9 1.1 0
48 5.1 1.1 0
49 4.9 0.9 -0.5
50 5.1 0.9 -0.5
51 4.9 1.1 -0.5
52 5.1 1.1 -0.5
53 4.9 0.9 -1
54 5.1 0.9 -1
55 4.9 1.1 -1
56 5.1 1.1 -1
57 4.8 0.8 -0.5
58 5.2 0.8 -0.5
59 4.8 1.2 -0.5
60 5.2 1.2 -0.5
61 4.8 0.8 -1
62 5.2 0.8 -1
63 4.8 1.2 -1
64 5.2 1.2 -1
1 7 hex 9 10 14 13 1 2 6 5
2 7 hex 10 11 15 14 2 3 7 6
3 7 hex 17 18 22 21 9 10 14 13
4 7 hex 18 19 23 22 10 11 15 14
5 7 hex 29 30 32 31 37 38 40 39
6 7 hex 33 34 30 29 41 42 38 37
7 7 hex 34 36 32 30 42 44 40 38
8 7 hex 35 31 32 36 43 39 40 44
9 7 hex 33 29 31 35 41 37 39 43
10 7 hex 33 34 36 35 29 30 32 31
11 7 hex 25 26 28 27 33 34 36 35
12 7 hex 19 20 26 25 11 12 34 33
13 7 hex 20 24 28 26 12 16 36 34
14 7 hex 23 27 28 24 15 35 36 16
15 7 hex 19 25 27 23 11 33 35 15
16 7 hex 19 20 24 23 25 26 28 27
17 7 hex 15 35 36 16 7 43 44 8
18 7 hex 12 16 36 34 4 8 44 42
19 7 hex 11 33 35 15 3 41 43 7
20 7 hex 11 12 34 33 3 4 42 41
21 7 hex 37 38 46 45 57 58 50 49
22 7 hex 38 40 48 46 58 60 52 50
23 7 hex 39 47 48 40 59 51 52 60
24 7 hex 37 45 47 39 57 49 51 59
25 7 hex 57 58 50 49 61 62 54 53
26 7 hex 58 60 52 50 62 64 56 54
27 7 hex 59 51 52 60 63 55 56 64
28 7 hex 57 49 51 59 61 53 55 63
29 7 hex 45 46 48 47 49 50 52 51
30 7 hex 49 50 52 51 53 54 56 55
31 4 quad 1 5 6 2
32 4 quad 2 6 7 3
33 4 quad 3 41 42 4
34 4 quad 3 7 43 41
35 4 quad 4 42 44 8
36 4 quad 7 8 44 43
37 2 quad 9 10 2 1
38 2 quad 10 11 3 2
39 2 quad 11 12 4 3
40 1 quad 12 16 8 4
41 3 quad 13 14 6 5
42 3 quad 14 15 7 6
43 3 quad 15 7 8 16
44 2 quad 17 18 10 9
45 5 quad 17 21 22 18
46 2 quad 18 19 11 10
47 5 quad 18 22 23 19
48 2 quad 19 20 12 11
49 5 quad 19 23 24 20
50 1 quad 20 24 16 12
51 3 quad 21 22 14 13
52 3 quad 22 23 15 14
53 3 quad 23 15 16 24
54 9 quad 37 57 59 39
55 9 quad 37 38 58 57
56 9 quad 38 40 60 58
57 9 quad 39 59 60 40
58 4 quad 41 43 39 37
59 4 quad 41 37 38 42
60 4 quad 42 38 40 44
61 4 quad 43 44 40 39
62 6 quad 53 55 56 54
63 9 quad 57 61 63 59
64 9 quad 57 58 62 61
65 9 quad 58 60 64 62
66 9 quad 59 63 64 60
67 6 quad 61 63 55 53
68 6 quad 61 53 54 62
69 6 quad 62 54 56 64
70 6 quad 63 64 56 55
71 9 line 61 62
72 9 line 62 64
73 9 line 64 63
74 9 line 63 61
75 9 line 57 58
76 9 line 58 60
77 9 line 60 59
78 9 line 59 57
79 9 line 37 38
80 9 line 38 40
81 9 line 40 38
82 9 line 38 37
_______________________________________________
