%% Zustandsberechnung mittels Kühlturmanalogie
clear

%% Eingabewerte 
CPr = 0.516;
g = 9.81;
% Massenströme
mL = 0.02;
mW = 0.1274;
% Stoffgrößen
cpWg = 1.86;
cpWf = 4.18;
cpL = 1;
lambda = 26.425*10^(-6);
eta = 1.875*10^(-5);
Pr = 0.7135;
hv0 = 2500;
phiein = 1;
rhoL = 1.15;
%Temperaturen und Druck
p = 101.325; %kPa
Tumg = 20;
Twein = 60;
TLein = 40;

    Tm = (Twein+TLein)/2;
    % Geometriegrößen
    Aaussen = 0.72;
    L = 0.7;
    a = 180;
    V = 0.032;

    %% anfängliche Berechnungen
    % Stoffaustauschkoeff., Temperaturausdehnungskoeff., spez. Eingangsenthalpie
    % Wasser
    K = (mW/(a*V))*0.52*(mL/mW)^(-0.16); 
    volkoeff = 1/(273.15+Tumg);
    hwein = cpWf*Twein;

    % Luftzustand
    Xein = 2.19*10^(-6)*TLein^3-1.85*10^(-4)*TLein^2+7.06*10^(-3)*TLein-0.077;
    hLein = 0.00585*TLein^3-0.497*TLein^2+19.87*TLein-207.61;

    

    x0 = [TLein+10,Twein-10]; %geratener Startwert
    my_fun = @(x)parameterfun(x,hwein,cpWf,mL,CPr,lambda,volkoeff,g,rhoL,Tm,Tumg,L,eta,Pr,mW,hLein,K,a,V,Aaussen);
    y = fsolve(my_fun,x0);
    y1 = real(y(1));
    y2 = real(y(2));

    Twaus = y2;
    TLaus = y1;

    Xaus = 2.19*10^(-6)*TLaus^3-1.85*10^(-4)*TLaus^2+7.06*10^(-3)*TLaus-0.077;
   
    mKond = (Xaus-Xein)*mL;
    Wasser = mKond*3600*18;