-
Notifications
You must be signed in to change notification settings - Fork 1
/
main.m~
80 lines (70 loc) · 1.52 KB
/
main.m~
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
clc
clear
close all
%Define Constants
k = 0.02605; %W/m^2*K @ 25 degrees C
T_ign = 632 + 273.15; %K
q_gen = 522.43 * 10^6; %W/m^3
q_out = 706554; %W/m^2
delta_x = 20.0 * 10^-3; %m
%Number of nodes
N = 30;
%Balance matrix
filename1 = 'balance_matrix.csv';
B = csvread(filename1);
%Position matrix
filename2 = 'position_matrix.csv';
p = csvread(filename2);
%Coefficient matrix
filename3 = 'coefficient_matrix.csv';
A = csvread(filename3);
%Constants matrix
filename4 = 'constants_vector.csv';
C = csvread(filename4);
%Solve
T = zeros(N,1);
T(1:27,1) = A(1:27,1:27)\C(1:27,1);
for i = 28:30
T(i) = T_ign;
end
Temp_grid = zeros(5,8);
Temp_grid(1,1) = T(1);
Temp_grid(1,2) = T(2);
Temp_grid(1,3) = T(3);
Temp_grid(1,4) = T(4);
Temp_grid(1,5) = T(5);
Temp_grid(1,6) = T(6);
Temp_grid(1,7) = T(7);
Temp_grid(1,8) = T(8);
Temp_grid(2,2) = T(9);
Temp_grid(2,3) = T(10);
Temp_grid(2,4) = T(11);
Temp_grid(2,5) = T(12);
Temp_grid(2,6) = T(13);
Temp_grid(2,7) = T(14);
Temp_grid(2,8) = T(15);
Temp_grid(3,3) = T(16);
Temp_grid(3,4) = T(17);
Temp_grid(3,5) = T(18);
Temp_grid(3,6) = T(19);
Temp_grid(3,7) = T(20);
Temp_grid(3,8) = T(21);
Temp_grid(4,4) = T(22);
Temp_grid(4,5) = T(23);
Temp_grid(4,6) = T(24);
Temp_grid(4,7) = T(25);
Temp_grid(4,8) = T(26);
Temp_grid(5,5) = T(27);
Temp_grid(5,6) = T(28);
Temp_grid(5,7) = T(29);
Temp_grid(5,8) = T(30);
%Visualization of solution
x = 0:-0.04:-0.28;
y = 0.02:0.04:0.18;
[X, Y] = meshgrid(x,y);
alpha = transpose(X);
beta = transpose(Y);
gamma = transpose(
figure
C = contourf(X,Y,Temp_grid);
clabel(C)