function [dQ_ztraty,dT_ztraty] = tepelne_ztraty(hustota,viskozita,Reynolds,dw_cat,rho_cat,area_reactor,cp_SUM,dp_cat,dr,temp_K,lr,Mw_str,m0_SUM,lambda_spaliny,w_cat_max,dtempdw,lreactor)
% pocita tepelne ztraty delkoveho elementu reaktoru dl

% zadani konstant
Re = Reynolds; % [-] 
g = 9.81; % [m.s-2]
deltaT = 10; % [K] rozdil mezi teplotou steny reaktoru a plynu - hruby odhad 10 K
T = temp_K; % [K]
Tokoli = 273.15+25; % [K]

    % vedeni a prestup tepla materialem
    lambda_ocel = 22; % [W.m-1.K-1]  priblizna hodnota pro ocel 1.4541 (Cr 18% Ni 10,5%)
    d_ocel = 0.002; % [m] tloustka plechu

    lambda_sibral = 0.095; % [W.m-1.K-1]  pouzita izolace Sibral 0.095
    d_sibral = 0.03; % [m] tloustka izolace

    alfa_vzduch = 5; % [W.m-2.K-1]  z literatury v rozmezi 1-60 W.m-2.K-1

    % uz nepouzivane - lambda_spaliny = 0.052;  [W.m-1.K-1] z tabulek pro spaliny 400st C
    
    

% vypocet delkoveho elementu reaktoru dl a elementu plochy dS 
dV = dw_cat/rho_cat; % [m^3] objemovy element katalytickeho loze
dl = dV/area_reactor; % [m] delkoveho elementu katalytickeho loze
dS = 4*dl*dr; % [m^2] plosny element katalytickeho loze
dlreactor = lreactor/(w_cat_max/dw_cat); % [m] delkovy element celeho reaktoru
dSreactor = 4*dlreactor*dr;  % [m^2] plosny element celeho reaktoru

% prepocet cp na [J.kg-1.K-1] 
cp_kg = cp_SUM*1000/Mw_str; % [J.kg-1.K-1]

% vypocet bezrozmernych kriterii

Pr = cp_kg*viskozita/lambda_spaliny; % [-] Prandtlovo cislo
Gr = dr^3*hustota^2*g*deltaT/(viskozita^2*T); % [-] Grashofovo cislo 
% - je pouzit charakteristicky rozmer reaktoru (dr)
GrPr =Gr*Pr; % [-]

    % vypocet Nusseltova cisla (plati pouze pro laminarni proudeni Re<2300)
    if GrPr > 2*10^5; % volna konvekce se vyznamne uplatnuje, Nu se vypocte podle vztahu:
    Nu = 0.74*(Re*Pr)^0.2*GrPr^0.1; % [-]
    
    else    % volna konvekce se neuplatnuje, Nu se vypocte podle vztahu:
    dviskozita = 0.9; % 0.04*10^(-5) visk/viskwall - odhad   
    Nu = 1.86*(Re*Pr*dr/lr)^0.333*(dviskozita)^0.1; % [-]
    end

alfa_spaliny = Nu*lambda_spaliny/dr; % [W.m-2.K-1]     

k = 1/(1/alfa_spaliny+d_ocel/lambda_ocel+d_sibral/lambda_sibral+1/alfa_vzduch); % [W.m-2.K-1] soucinitel prestupu tepla
 
dQ_ztraty_okoli = (T-Tokoli)*k*dSreactor;  % [W] tepelne ztraty do okoli na delkovem elementu dl

cp_ocel = 460; % [J.kg-1.K-1] tepelna kapacita oceli
m_vestavba_ocel = 0; % [kg] hmotnost konstrukce vestavby (dle výkresu Solid) 30.9

% m0_SUM [kg/s]
Tout = -(dQ_ztraty_okoli/(m0_SUM*cp_kg)-T); % [K] vypocet teploty smesi po tepelnych ztratach
dT_ztraty = T-Tout; % pokles teploty diky tepelnym ztratam na krokovem elementu dw resp. dl a dS

dQ_akumulace_konstrukce = (dtempdw*dw_cat-dT_ztraty)*cp_ocel*m_vestavba_ocel/3600; % [W] tepelne ztraty do konstrukce na delkovem elementu dl

dQ_ztraty = dQ_akumulace_konstrukce + dQ_ztraty_okoli; % [W]

Tout = -(dQ_ztraty/(m0_SUM*cp_kg)-T); % [K] vypocet teploty smesi po tepelnych ztratach
dT_ztraty = T-Tout; % pokles teploty diky tepelnym ztratam na krokovem elementu dw resp. dl a dS
end

