WELCOME TO PARABOLIC ELEMENT BEM PROGRAM
I N P U T D A T A
BOUNDARY ELEMENT METHOD APPLIED TO POTENTIAL PROBLEMS
BOUNDARY CONDITIONS ARE ASSIGNED TO ELEMENTS.
IF ELM TYPE(I) = 1, THEN DIRICHLIT.
IF ELM TYPE(I) = 2, THEN NEUMANN.
====> NUMBER OF ELEMENTS= 14
ELM # ---> ELEMENT NUMBER
ELM # I J K ELM TYPE BOUNDARY VALUES
1 1 2 3 2 0.00000 0.00000 0.00000
2 3 4 5 2 0.00000 0.00000 0.00000
3 5 6 7 2 0.00000 0.00000 0.00000
4 7 8 9 2 0.00000 0.00000 0.00000
5 9 10 11 2 0.00000 0.00000 0.00000
6 11 12 13 2 0.00000 0.00000 0.00000
7 13 14 15 1 60.00000 60.00000 60.00000
8 15 16 17 1 60.00000 60.00000 60.00000
9 17 18 19 2 0.00000 0.00000 0.00000
10 19 20 21 2 0.00000 0.00000 0.00000
11 21 22 23 2 0.00000 0.00000 0.00000
12 23 24 25 2 0.00000 0.00000 0.00000
13 25 26 27 1 60.00000 60.00000 60.00000
14 27 28 1 1 60.00000 60.00000 60.00000
====> NUMBER OF BOUNDARY NODES= 28
NODE X-COORDINATE Y-COORDINATE
1 180.00000000 60.00000000
2 150.00000000 60.00000000
3 120.00000000 60.00000000
4 90.00000000 60.00000000
5 60.00000000 60.00000000
6 30.00000000 60.00000000
7 0.00000000 60.00000000
8 -30.00000000 60.00000000
9 -60.00000000 60.00000000
10 -90.00000000 60.00000000
11 -120.00000000 60.00000000
12 -150.00000000 60.00000000
13 -180.00000000 60.00000000
14 -180.00000000 45.00000000
15 -180.00000000 30.00000000
16 -180.00000000 15.00000000
17 -180.00000000 0.00000000
18 -135.00000000 0.00000000
19 -90.00000000 0.00000000
20 -45.00000000 0.00000000
21 0.00000000 0.00000000
22 45.00000000 0.00000000
23 90.00000000 0.00000000
24 135.00000000 0.00000000
25 180.00000000 0.00000000
26 180.00000000 15.00000000
27 180.00000000 30.00000000
28 180.00000000 45.00000000
====> NUMBER OF INTERIOR POINTS= 4
POINT X-COORDINATE Y-COORDINATE
1 0.00000000 20.00000000
2 0.00000000 25.00000000
3 0.00000000 30.00000000
4 0.00000000 35.00000000
====> NUMBER OF POINT SOURCES = 2
X-COORDINATE Y-COORDINATE MAGNITUDE
90.000000 15.000000 -6.0000000
-90.000000 15.000000 -6.0000000
====> NUMBER OF NODES ON FREE SURFACE = 13
NODAL NUMBER
1 2 3 4 5 6 7 8 9 10 11 12 13
END OF INPUT-DATA ECHO PRINT
*** N U M E R I C A L S O L U T I O N ***
-- FLUX ON BOUNDARY -- (QN = HEAT FLUX IN NORMAL DIRECTION)
ELEMENT NODE-I NODE-J NODE-K QN(I) QN(J) QN(K)
1 1 2 3 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
2 3 4 5 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
3 5 6 7 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
4 7 8 9 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
5 9 10 11 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
6 11 12 13 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
7 13 14 15 -0.8018520835E-01 -0.9844336240E-01 -0.1012696858
8 15 16 17 -0.1012696858 -0.1043548332 -0.1070424324
9 17 18 19 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
10 19 20 21 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
11 21 22 23 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
12 23 24 25 0.0000000000E+00 0.0000000000E+00 0.0000000000E+00
13 25 26 27 -0.1070424324 -0.1043548332 -0.1012696858
14 27 28 1 -0.1012696858 -0.9844336240E-01 -0.8018520835E-01
-- FREE TERM AND POTENTIAL VALUES AT NODAL POINTS --
NODE FREE TERM POTENTIAL
1 0.235378 60.000000
2 0.500000 57.135533
3 0.498654 54.072273
4 0.500000 51.450155
5 0.499747 50.434620
6 0.500000 50.253351
7 0.499466 50.226436
8 0.500000 50.253351
9 0.499747 50.434620
10 0.500000 51.450155
11 0.498654 54.072273
12 0.500000 57.135533
13 0.235378 60.000000
14 0.500000 60.000000
15 0.500000 60.000000
16 0.500007 60.000000
17 0.250000 60.000000
18 0.500000 55.148486
19 0.500000 49.945371
20 0.500000 50.133142
21 0.500000 50.194616
22 0.500000 50.133142
23 0.500000 49.945371
24 0.500000 55.148486
25 0.250000 60.000000
26 0.500007 60.000000
27 0.500000 60.000000
28 0.500000 60.000000
-- POTENTIAL AT INTERNAL POINT --
POINT X-COORD Y-COORD POTENTIAL
1 0.00000E+00 20.000 50.208
2 0.00000E+00 25.000 50.212
3 0.00000E+00 30.000 50.216
4 0.00000E+00 35.000 50.219
COORDINATES OF FREESURFACE
X= 180.000000000000 Y= 60.0000000000000
X= 150.000000000000 Y= 57.1355327667767
X= 120.000000000000 Y= 54.0722728009394
X= 90.0000000000000 Y= 51.4501548706977
X= 60.0000000000000 Y= 50.4346197267023
X= 30.0000000000000 Y= 50.2533509052686
X= 0.000000000000000E+000 Y= 50.2264364953581
X= -30.0000000000000 Y= 50.2533509052686
X= -60.0000000000000 Y= 50.4346197267022
X= -90.0000000000000 Y= 51.4501548706976
X= -120.000000000000 Y= 54.0722728009394
X= -150.000000000000 Y= 57.1355327667767
X= -180.000000000000 Y= 60.0000000000000