clc; clear; 

% Strömungsgeschwindigkeit

d_Schlauch= 25.4e-3;                     % Innendurchmesser PFA-Schlauch [m]
A_Schlauch=d_Schlauch^2*pi/4;            % Fläche PFA-Schlauch [m^2]
Q=5.83e-4;                               % Volumenstromvpumpe [m^3/s]

c=4*Q/(d_Schlauch^2*pi);                 % Strömungsgeschwindigkeit [m/s] S.45 (3.2b )

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Kinematische Viskosität

roh=1250;                                % Dichte 40%HNO3 bei 20°C & 1,013bar [1,25g/cm^3 = 1250kg/m^3]
n=1.55e-3;                               % Viskosität [1,55mPa*s = 1,55e-3kg/(m*s)] 

v=n/roh;                                 % Kinematische Viskosität [m^2/s]

%http://www.schweizer-fn.de/stoff/begriffe/v2_begriffe.htm
%http://www.convertworld.com/en/dynamic-viscosity/Centipoise.html
%http://www.cactus2000.de/de/unit/massvis.shtml

% Reynoldszahl

Re=c*d_Schlauch/v                       % Strömung Re<=2300 laminar / Re>2300 turbolent 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Rohrwiderstandsbeiwert 
%http://www.schweizer-fn.de/stroemung/druckverlust/v2_druckverlust.htm#druckverlustzeta

Ra_Schlauch=0.0002;                     % Oberflächenrauhigkeit PFA-Schlauch [mm]

Lambda=0.3164/(Re^0.25)                 % Rohrwiderstandsbeiwert für hydr. glattes Rohr k=0 / S.49/Bild3-1

%http://www.win-ing.de/lambdaturbulent.php

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Druckverlust delta_P

l_Schlauch_waag=40;                     % Summe Schlauchlänge waagrecht verlegt  [m]
l_Schlauch_steig=3;                     % Summe Schlauchlänge steigend in Strömungsrichtung verlegt [m]
l_Schlauch_fall=2;                      % Summe Schlauchlänge fallend in Strömungsrichtung verlegt [m]

g=9.81;                                 % Erdbeschleunigung [m/s^2]

h1=0.5;                                 % Höhe Pumpe  [m]
h2=3.5;                                 % Höhe Rohrtrasse  [m]
h3=1.5;                                 % Höhe Abfüllkopf  [m]

Zeta_Bogen_Fitting=0.14;               % r/d=36/25,4=1.4 entspricht bei 90° "0.14"  S.52/Tab.3.2
Zeta_Bogen_eigengefertigt=0.09;        % r/d=700/25,4=27.6 entspricht bei 90° "0.09"  S.52/Tab.3.2

Zeta_sum_Bogen_Fitting=0.14*5;              
Zeta_sum_Bogen_eigengefertigt=0.09*4;

% Widerstandsbeiwert Zeta=Lambda*l_Schlauch/(d_Schlauch*10^3)

delta_p_waag=(Lambda*l_Schlauch_waag/(d_Schlauch))*(roh*c^2)/2;                  % Druckverlust für waagrecht verlegten Schlauch [Pa]

delta_p_steig=(Lambda*l_Schlauch_steig/(d_Schlauch))*(roh*c^2)/2+roh*g*(h2-h1);  % Druckverlust für steigend verlegten Schlauch [Pa]

delta_p_fall=(Lambda*l_Schlauch_fall/(d_Schlauch))*(roh*c^2)/2+roh*g*(h2-h3);    % Druckverlust für steigend verlegten Schlauch [Pa]

delta_p_Bogen_Fitting=Zeta_sum_Bogen_Fitting*(roh*c^2)/2;                        % Druckverlust für Summe Entegris-Bogen 90° [Pa]

delta_p_Bogen_eigengefertig=Zeta_sum_Bogen_eigengefertigt*(roh*c^2)/2;           % Druckverlust für Summe selbstgefertigte Bögen 90° [Pa]

delta_p_sum=delta_p_waag+delta_p_steig+delta_p_fall+delta_p_Bogen_Fitting+delta_p_Bogen_eigengefertig

%????deltap3=Re*Q^2

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Anlagenkennlinie (aktuell ohne statischem Druck)

Zeta_Schlauch_waag=Lambda*l_Schlauch_waag/(d_Schlauch*10^3);                                         % Widerstandsbeiwert 
Zeta_Schlauch_steig=Lambda*l_Schlauch_steig/(d_Schlauch*10^3);                                       % Widerstandsbeiwert 
Zeta_Schlauch_fall=Lambda*l_Schlauch_fall/(d_Schlauch*10^3);                                         % Widerstandsbeiwert 

Zeta_sum=Zeta_Schlauch_waag+Zeta_Schlauch_steig+Zeta_Schlauch_fall+Zeta_sum_Bogen_Fitting+Zeta_sum_Bogen_eigengefertigt

Q=0:0.001:0.0010

%delta_h=0.5;                                                        % Höhendifferens Tank Pumpe [m]
%p_statisch=roh*g*delta_h
%p_dynamisch=Zeta_sum*(roh/2)*(Q/A_Schlauch).^2;

%delta_p_Anlage=p_statisch+p_dynamisch;

delta_p_Anlage=(1+Zeta_sum)*(roh/2)*c.^2                             % [(kg/m^3)*(m^3/s)^2*(1/m^2)^2 = kg/(m*s^2) = Pa]

plot(Q,delta_p_Anlage);grid;
title('Anlagenkennlinie');
xlabel('Volumenstrom Q [m³/s]'); 
ylabel('Druck p [Pa]'); 

%Kv-Werte Ventil und Drucktransmitter sowie Filter fehlen noch!