clear all;
close all;
clc;

%Variable-Def
d = 0.001; %Durchmesser Messlanze [m]
L = 0.3; %Länge Messlanze [m]
T_0 = 28+273.15; %Lanzentempertur [C]
T_F = -30+273.15; % Fluidtemperatur [C]
Q_v = 0.003; %Verlustwärmestrom [W]

%%
%Messlanze (Al)
ro_Al = 2700; %Dichte [kg/m3]
cp_Al = 1030; %spez.Wärmekapazität [J/kgK]
l_Al = 236; %Lambda [kg/m^3]

%%
%Messlanze (Cu)
ro_Cu = 8960; %Dichte [kg/m3]
cp_Cu = 401; %spez.Wärmekapazität [J/kgK]
l_Cu = 400; %Lambda [kg/m^3]

%%
%Strömung
a = 35070; %Wärmeübergangskoeffizient [W/m^2K]
t_tot = 15*60; %Simulationszeit [s]

%%
%Formeln
k_Al = (a^(-1) + L / (10 * l_Al))^(-1); % Wärmeduchrgangkoef Al
k_Cu = (a^(-1) + L / (10 * l_Cu))^(-1); % Wärmeduchrgangkoef Cu

A = d * pi *L; % Oberfläche im Fluid
V = pi * (d/2)^2 * L; % Volumen Messlanze

m_Al = ro_Al * V; % Masse Messlanze AL [kg]
m_Cu = ro_Cu * V; % Masse Messlanze Cu [kg]

%%
%Initialisation
n=10000;
dt=t_tot/n; %Zeitschrittweite
t(1)=0;
T(1)=T_0;

%%
%Main-AWL-Loop

for i=2:n
    t(i)=dt*(i-1);
    dT(i)=[(((k_Al*A)/(m_Al*cp_Al))*T_F-(((k_Al*A)/(m_Al*cp_Al))*T)-Q_v)]*(t(i)-t(i-1));
end
    
    
%%
%FUP - MainLoop
simOPT=simset('Solver','ode5','FixedStep',dt);
sim('EXP03_SimDD',Tsim,simOPT);

%%
%Plot
figure(1)
hold on
box off
plot(t/60,dT,'r-')
plot(tout./60,Tsim,'b--')
%axis([0 -11 5 30])
xlabel('Zeit in Minuten')
ylabel('Temperatur in °C')
set(gca,'FontSize',14)