-
Notifications
You must be signed in to change notification settings - Fork 12
/
Radiator.m
89 lines (77 loc) · 3.79 KB
/
Radiator.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
81
82
83
84
85
86
87
88
89
classdef Radiator
properties
iteration
id_inlet
id_outlet
id_zone
time_step
matrix_size
matrix_coefficients
right_hand_side_vector
number_of_equations
specific_heat_capacity
mass
specific_heat_capacity_fluid
mass_fluid
temperature_zone
temperature_inlet
temperature_outlet
mass_flow_rate
heat_transfer_coefficient
surface
end
methods
function obj = Radiator(id_inlet, id_outlet, id_zone, solver, specific_heat_capacity, mass, specific_heat_capacity_fluid, mass_fluid, mass_flow_rate,heat_transfer_coefficient,surface)
if nargin > 0
obj.id_inlet = id_inlet;
obj.id_outlet = id_outlet;
obj.id_zone = id_zone;
obj.time_step = solver.time_step;
obj.matrix_size = solver.matrix_size;
obj.specific_heat_capacity = specific_heat_capacity;
obj.mass = mass;
obj.specific_heat_capacity_fluid = specific_heat_capacity_fluid;
obj.mass_fluid = mass_fluid;
obj.mass_flow_rate = mass_flow_rate;
obj.heat_transfer_coefficient = heat_transfer_coefficient;
obj.surface = surface;
obj.iteration = 0;
obj.number_of_equations = 1;
obj.matrix_coefficients = zeros(obj.number_of_equations,solver.matrix_size);
obj.right_hand_side_vector = zeros(obj.number_of_equations,1);
obj.temperature_inlet = solver.temperatures(obj.id_inlet);
obj.temperature_outlet = solver.temperatures(obj.id_outlet);
obj.temperature_zone = solver.temperatures(obj.id_zone);
end
end
% coefficient of inlet temperature of radiator
function c = c_ti(obj)
c = (obj.mass*obj.specific_heat_capacity+obj.mass_fluid*obj.specific_heat_capacity_fluid)/(2*obj.time_step)-obj.mass_flow_rate*obj.specific_heat_capacity_fluid+obj.heat_transfer_coefficient*obj.surface/2;
end
% coefficient of outlet temperature of radiator
function c = c_to(obj)
c = (obj.mass*obj.specific_heat_capacity+obj.mass_fluid*obj.specific_heat_capacity_fluid)/(2*obj.time_step)+obj.mass_flow_rate*obj.specific_heat_capacity_fluid+obj.heat_transfer_coefficient*obj.surface/2;
end
% coefficient of zone temperature of radiator
function c = c_tz(obj)
c = -obj.heat_transfer_coefficient*obj.surface;
end
% right-hand coefficient of radiator
function c = c_r(obj)
c = (obj.mass*obj.specific_heat_capacity+obj.mass_fluid*obj.specific_heat_capacity_fluid)*(obj.temperature_inlet+obj.temperature_outlet)/(2*obj.time_step);
end
% create matrix of coefficients and right-hand side vector
function obj = create(obj, solver)
obj.iteration = obj.iteration + 1;
obj.temperature_inlet = solver.temperatures(obj.id_inlet);
obj.temperature_outlet = solver.temperatures(obj.id_outlet);
obj.temperature_zone = solver.temperatures(obj.id_zone);
obj.matrix_coefficients = zeros(obj.number_of_equations,obj.matrix_size);
obj.right_hand_side_vector = zeros(obj.number_of_equations,1);
obj.right_hand_side_vector = obj.c_r();
obj.matrix_coefficients(obj.id_inlet) = obj.c_ti();
obj.matrix_coefficients(obj.id_outlet) = obj.c_to();
obj.matrix_coefficients(obj.id_zone) = obj.c_tz();
end
end
end