Tr_AMTD.m

function Tr = Tr_AMTD(q,Ts,Ti,qo,Tso,Tro,Tio,n)
%%  Calculates the return temperature from the radiator unit
%   emprical equation based on AMTD by Soumerai H. 
%   Ref: Phetteplace - Optimal Design of Piping Systems for District Heating
%   prepared by Hakan ibrahim Tol, PhD

%% INPUT&OUTPUT
%   Input
%   q   : Heat demand - actual condition            [kW]
%   Ts  : Supply temperature - actual condition     [degC]
%   Ti  : Set degree of indoor temperature - actual [degC]

%   qo  : Heat demand - design condition            [kW]
%   Tso : Supply temperature - design condition     [degC]
%   Tro : Return temperature - design condition     [degC]
%   Tio : Set degree of indoor temperature - design [degC]
%   n   : Emprical radiator constant                [-]

%   Output
%   Tr  : Radiator return temperature - actual      [degC]

%% Calculation of AMTD for the design condition
%   AMTD: Arithmetic Mean Temperature Difference    [degC]

AMTDo=tAMTD(Tso,Tro,Tio);

%% Calculation

Tr=2*(Ti+(AMTDo*(q/qo)^(1/n)))-Ts;

if or(Tr<=20,Tr>=Ts)
    Tr=NaN;
    return
end

%% Evaluation of error (AF: Approach Factor)
AF=(Tr-Ti)/(Ts-Ti);

if AF>=0.5
    fprintf('User Warning: Approach factor is %f - Error less than 0.04',AF);
else
    fprintf('User Warning: Approach factor is %f - Error larger than 0.04',AF);
end

end