27 February 2022 01:30:17 PM
FEM2D_HEAT
C++ version:
Compiled on Feb 27 2022 at 13:29:53.
Solution of the time dependent heat equation
on an arbitrary triangulated region D in 2 dimensions.
Ut - Uxx - Uyy + K(x,y,t) * U = F(x,y,t) in D
U = G(x,y,t) on the boundary.
U = H(x,y,t) at initial time.
The finite element method is used with
6 node quadratic triangular elements ("T6").
The time derivative is approximated using the
backward Euler method.
Node file is "square_nodes.txt".
Element file is "square_elements.txt".
Number of nodes = 81
First 10 nodes
Row: 0 1
Col
0: 0 0
1: 0.125 0
2: 0.25 0
3: 0.375 0
4: 0.5 0
5: 0.625 0
6: 0.75 0
7: 0.875 0
8: 1 0
9: 0 0.125
Element order = 6
Number of elements = 32
First 10 elements
Row: 1 2 3 4 5 6
Col
1 1 19 3 10 11 2
2 21 3 19 12 11 20
3 3 21 5 12 13 4
4 23 5 21 14 13 22
5 5 23 7 14 15 6
6 25 7 23 16 15 24
7 7 25 9 16 17 8
8 27 9 25 18 17 26
9 19 37 21 28 29 20
10 39 21 37 30 29 38
Quadrature order = 7
The matrix half bandwidth is 18
The matrix bandwidth is 37
The storage bandwidth is 55
Initial time = 0
Final time = 0.5
Step size = 0.05
Number of steps = 10
FEM2D_HEAT:
Normal end of execution.
27 February 2022 01:30:17 PM