%kontroludregning for energi ind og ud i systemet ved steady state
format short;
clc, clear;

stoker = 9;     %duty cycle (%) stoker i steady state (SS)
blower = 56.36; %duty cycle (%) blæser SS

tp_procent = 95*10^-6;      % masse tp der fremføres/sek/%stoker kg/s
blow_procent = 0.259*10^-3; % masse luft der indblæses/sek/%blower kg/s

H = 17.8*10^6;      %J/kg, brændværdi  træpiller
c_vand = 4.18*10^3; %J/(kg*K), specifik varmekapacitet vand
c_luft = 1007;      %J/(kg*K), specifik varmekapacitet luft @ 20C
c_rg = 1012;        %J/(kg*K), specifik varmekapacitet luft @ 100C

m_tp = tp_procent*stoker;       %kg/s, masseflow træpiller
m_luft = blower*blow_procent;   %kg/s, masseflow luft
m_rg = (m_tp*0.85 + m_luft);    %kg/s, masseflow røggas

T_omg = 20;         %omgivelsestemperatur grader C
T_retur = 50.00;    %returvandstemperatur grader C
T_frem = 57.65;     %temperatur i fremløbsvand grader C
T_rg = 106.13;      %temperatur i røggas ud i skorsten grader C


S = solve('T_omg*c_luft*m_luft + H*m_tp + m_vand*c_vand*(T_retur) - m_rg*c_rg*(T_rg)-m_vand*c_vand*(T_frem)','m_vand');
m_vand = eval(S);    %Udregn masseflow af vand ud fra Formel 3.3

E_luft = T_omg*c_luft*m_luft;       %effekt ind med luft
E_tp = H*m_tp ;                     %effekt ved afbrænding af træpiller
E_retur = m_vand*c_vand*(T_retur);  %effekt ind med returvand
E_rg = m_rg*c_rg*(T_rg);            %effekt ud med røggas
E_frem = m_vand*c_vand*(T_frem);    %effekt ud med fremløbsvand

%% Udskriver resultater i command window

m_vand
E_luft
E_tp
E_retur
E_rg
E_frem

% disp('sum af input')
% E_ind = E_tp+E_retur + E_luft
% 
% disp('Sum af output')
% E_ud = E_frem + E_rg
% 
% disp('forskel E_ind - E_ud')
% E_ind-E_ud