%clear all
%close all
%clc

%%  Fyzikalni konstanty

c_Voda = 4200;                      %   Merna tepelna kapacita voda
c_Vzduch = 1010*1.4;                %   Merna tepelna kapacita sucheho vzduchu
c_pp = 1840;                        %   Merna tepelna kapacita vodni pary
c_chladici = 2331*0.35+4200*0.65;   %   Merna tepelna kapacita chladici kapaliny
c_Med = 383;                        %   Merna tepelna kapacita medi

kapa_Vzduch = 1.4;
kapa_Vodni_para = 1.33;

lambda_med = 386;                   %   soucinitel tepelne vodivosti medi
hustota_med = 8960;                 %   hustota medi

l_voda23 = 2500000;     %   Merne vyparne teplo vody

Rm = 8.3144598;         %   Univerzalni plynova konstanta J.K-1.mol-1
M = 0.02896;            %	Molova hmotnost vzduchu kg.mol-1
M_voda = 0.01802;       %	Molova hmotnost vody

%%  Globalni parametry

T_Okoli = 273.15+20;      %   Teplota okolniho vzduchu
p_Okoli = 100000;   %   Tlak vzduchu v okoli
fi_Okoli = 0.6;     %   Relativni vlhkost vzduchu v okoli

%   parametry tepelneho vymeniku
C_vymenik = 152138;            %   Tepelna kapacita vymeniku
C_Topne_teleso = 1834;
R_vymenik_vzduch = 0.006555605;    %   Tepelny odpor vymeniku na prechodu vzduch - vymenik
R_vymenik_kapalina = 0.0000588347131321683;  %   Tepelny odpor vymeniku na prechodu kapalina - vymenik
R_Vymenik_teleso = 0.018223256;
T_Vymenik_Pocatecni = T_Okoli;  %   Pocatecni teplota vymeniku
m_chladici = 42.53817762;
delay = 206.0581943;
K_topne_teleso = 1;
T_Cirk1_Chlazeni_zasobnik = 8 + 273.15;

%   parametry chladiciho potrubi
d_potrubi = 0.028;  %   prumer potrubi
S_potrubi = 84.823e-6; %   plocha prurezu potrubi
o_potrubi = pi*(d_potrubi-2);   %   vnitrni obvod potrubi

%%  Hlavní komora
T_Komora_pocatecni = T_Okoli;

V_Komora = 1.19*1.06*1.90;   %   Objem komory

p_Komora_pocatecni = 100000;    %   Pocatecni tlak vzduchu v komore

fi_Komora_pocatecni = 0.28;      %   Pocatecni relativni vlhkost v komore

C_veci_v_komore = 12677;        %   tepelna kapacita veci, ktere se nachazi v komore

R_veci_v_komore = 0.116145332;          %   prumernytepelny odpor veci, ktere se nachazi v komore

Par_Komora_prostup_okoli =[T_Okoli; %   Pocatecni teplota
                        0.025;      %   Merna tepelna vodivost
                        50;         %   Hustota izolacniho materialu
                        1500;       %   Merna tepelna kapacita izolacniho materialu
                        1.9*1.4 + 2*(1.9*2.25 + 2.25*1.4);     %   Plocha prostupu tepla  1.9 * 1.4 * 2.25 vnejsi rozmery
                        0.08];	%	Tloustka izolace
                    
Par_Komora_prostup_okoli2 =[T_Okoli; %   Pocatecni teplota
                        0.00166568724341481;    %   Merna tepelna vodivost
                        50;    %   Hustota izolacniho materialu
                        1500;    %   Merna tepelna kapacita izolacniho materialu
                        1.9*1.4;     %   Plocha prostupu tepla  1.9 * 1.4 * 2.25 vnejsi rozmery
                        0.01];	%	Tloustka izolace
                    
S0 = 1e-4;  %   Plocha  otvoru v komore
t_netesnost = 5;    %   dopravni zpozdeni netesnosti
                    
Parc_tlak_vodni_para = exp(23.58-(4044.2/(T_Komora_pocatecni-37.55)))*fi_Komora_pocatecni;
Parc_tlak_suchy_vzduch = p_Komora_pocatecni-Parc_tlak_vodni_para;
M0_suchy_vzduch = (Parc_tlak_suchy_vzduch*V_Komora*M)/(Rm*T_Komora_pocatecni);
M0_vodni_para = (Parc_tlak_vodni_para*V_Komora*M_voda)/(Rm*T_Komora_pocatecni);
x0 = M0_vodni_para/M0_suchy_vzduch;
%%  Cirkulacni ohrev a chlazeni 1

Qm_Cirk1_Chlazeni = 2;
Qm_Parameter = 1;
Qv_Cirk1 = 0.091626145;%0.1;    %   Objemovy prutok cirkulace
Pz_Cirk1_vent = 128.8;    %   Ztratovy vykon ventilatoru v okruhu
T_chladici_pocatecni = 273.15 + 20.2;

%   Parametry trubky 1
t_Cirk1_1 = 1;  %   Dopravni zpozdeni mezi komorou a topnym telesem

%   Parametry trubky 2
t_Cirk1_2 = 1;  %   Dopravni zpozdeni mezi topnym telesem a komorou

%   Parametry chladici trubky
t_zpozdeni_chladici_in = 45.8159622;
t_zpozdeni_chladici_out = 13.13854088;
t_zpozdeni_chladici_nadrz = 1.11934155;
%   Vypocet pocatecnich podminek pro dopravni zpozdeni

Qm0_Cirk1_vzduch = M*Qv_Cirk1*p_Komora_pocatecni/(T_Komora_pocatecni*Rm);
P0_Cirk1_vzduch = 7/2*p_Komora_pocatecni*Qv_Cirk1 + x0*Qm0_Cirk1_vzduch*(l_voda23+c_pp*T_Komora_pocatecni);
T0_Cirk1 = T_Komora_pocatecni;
Qm0_Cirk1_voda = x0*Qm0_Cirk1_vzduch;

%%  Komora se svetly

P_LED = 0;      %   Prikon LED osvetleni

R_sklo = 0.0109759;    %	Tepelny odpor skla mezi hlavni komorou a LED komorou
C_sklo = 11440;    %	Tepelna kapacita skla mezi hlavni komorou a LED komorou

C_LED = 34156;  %   Tepelna kapacita panelu se svetly
R_LED = 0.0623549;    %   Tepelny odpor prostupu tepla panelu se svetly

C_strop = 179498;  %   Tepelna kapacita stropu fytotronu
R_strop = 0.059736947;   %   Tepelny odpor stropu fytotronu

LED_pomer_ztrat = 0.1;  %   Pomer ztratoveho vykonu svetel ku prikonu svetel

P_LED_max = 925.376;   %   Maximalni prikon svetel
P_LED_offset = 221.41;

P_Zareni = 2.515741854;
T_LED_max = 310;

Qv_LED = 0.01070430;



V_Komora_LED = 1.9*1.4*0.05;

Par_LED_prostup_okoli =[T_Okoli; %   Pocatecni teplota
                        0.03645;    %   Merna tepelna vodivost
                        50;    %   Hustota izolacniho materialu
                        1500;    %   Merna tepelna kapacita izolacniho materialu
                        1.9*1.4;     %   Plocha prostupu tepla  1.9 * 1.4 * 2.25 vnejsi rozmery
                        0.01];	%	Tloustka izolace


%   Vypocet pocatecnich podminek pro dopravni zpozdeni

% Qm0_LED_Sani_vzduch = M*Qv_LED_sani*p_Komora_pocatecni/(T_Okoli*Rm);
% P0_LED_Sani_vzduch = 7/2*p_Komora_pocatecni*Qv_LED_sani + x0*Qm0_LED_Sani_vzduch*(l_voda23+c_pp*T_Okoli);
% T0_LED_Sani = T_Okoli;
% Qm0_LED_Sani_voda = x0*Qm0_LED_Sani_vzduch;

%%  Odvlhcovac vzduchu

Qv_odvl_suchy = 290/3600/10;   %   Objemovy prutok odvlhcovaneho vzduchu
Qv_odvl_vlhky = 80/3600/10;    %   Objemovy prutok regeneracniho vzduchu

t_Odvl_1 = 1;     %   Dopravni zpozdeni mezi komorou a vstupem odvlhcovace
t_Odvl_2 = 1;     %   Dopravni zpozdeni mezi vystupem odvlhcovace a komorou

deltaT1 = 2.514;
deltaT2 = 20.36;

%   Vypocet pocatecnich podminek pro dopravni zpozdeni
x0 = 0.0073;
Qm0_Odvl_vzduch = M*Qv_odvl_suchy*p_Komora_pocatecni/(T_Komora_pocatecni*Rm);
P0_Odvl_vzduch = 7/2*p_Komora_pocatecni*Qv_odvl_suchy + x0*Qm0_Odvl_vzduch*(l_voda23+c_pp*T_Komora_pocatecni);
T0_Odvl = T_Komora_pocatecni;
Qm0_Odvl_voda = x0*Qm0_Odvl_vzduch;
P_odvlh_od = 100;

%%  Zvlhcovac
Qm_Zvlhcovac_k = -4e-5;%-2.615e-5;    %   Hmotnostni prutok vody ve zvlhcovaci je dana vztahem Qm = k*fi + c
Qm_Zvlhcovac_c = 40e-6;%22.346e-6;
Qm_Zvlhcovac = 10e-5;
t_Zvlh = 6;     %   dopravni zpozdeni vody ve zvlhcovaci
T_Zvlh = 273.15+20; %   teplota vody ve zvlhcovaci

% %%  Cidla
% PT100delay = 10;
% %[PT100num,PT100den] = pade(PT100delay,1);
% PT100num = 1/PT100delay;
% PT100den = [1 1/PT100delay];
t_Otvor=10;

%%  Hlavni regulator
Kr = 0.0382;
Ti = 2660;
Kr_chl = 2.6;
T_zadana = 273.15 + 20;

%%  Regulator poruchy osvetleni
K_osvetleni = 4800;
C_osvetleni = 1880;
Td_osvetleni = 10;
T_osvetleni = 5;

