%*************************************************************
% Modelo de un calentador con volumen constante.
% Se considera un temperatura UNICA de la camisa y una UNICA del tanque
% Se considera densidad constante
% M. Rodriguez 2003
%****************************************************************

function dx_dt=heater(t,x);

F=10; %Caudal de entrada/salida del tanque
Fj=15; %Caudal de la camisa
Ti=50; %Temperatura de entrada al tanque
Tji=150; %Temperatura de entrada a la camisa
V=25; %Volumen del tanque
Vj=10; %Volumen de la camisa
rho=1; %Densidad fluido en el tanque
cp=1.25; %Calor especifico fluido en el tanque
rhoj=1.1; %Densidad del fluido en la camisa
cpj=1.3; % Calor especifico fluido en la camisa
U=5; %Coeficiente global de transmision
A=31; %Area de intercambio

T=x(1);
Tj=x(2);

if t>30
    Fj=20;
end
if t>60
        Fj=15;
    end
    if t>90
            Tji=200;
end 

%***Ecuaciones del modelo***
dT_dt=(F/V)*(Ti-T)+U*A*(Tj-T)/(V*rho*cp);
dTj_dt=(Fj/Vj)*(Tji-Tj)-U*A*(Tj-T)/(Vj*rhoj*cpj);
dx_dt(1)=dT_dt;
dx_dt(2)=dTj_dt;
dx_dt=dx_dt';
