%% Parameter zur Simulation der zeitabhängigen Raketengleichung

%Globale Variablen
global rakete
global druck

%Allgemeine Variablen
ml=0.223;                   %Masse der leeren Rakete [kg]
d2=0.009;                   %m; Durchmesser an der Austrittsstelle (Gardena-Düse)[m]
d1=0.0923;                  %m; Durchmeesser des Flugkörpers an der Stelle 
                            %der Spiegelhöhe zur Zeit t=0 [m]
x=0.355;                     %Eingabe erfolgt in Litern [l]
Vw=x*10^(-3);               %Volumen des Brennstoffs [m³]
%m=0.58;                    %kg; Gesamtmasse Rakete zur Zeit t=0 [kg]
%hf=0.34;                   %m; gesamte Flugkörperhöhe [m]

A1=pi*d1^2/4;               %Querschnittsfläche des Flugkörpers [m²]
A2=pi*d2^2/4;               %Querschnittsfläche an der Austrittsstelle [m²]

h1=Vw/A1;                   %Höhe des Wasserspiegels [m]
Vges=0.0015;                %Fassungsvolumen des Flugkörpers [m³]
Vl=Vges-Vw;                 %Luftvolumen [m³]
hf=Vges/A1;                 %Flugkörperhöhe [m]


p0=1.013*10^5;              %Atmosphärendruck [Pa]
rholuft=1.225;              %Dichte Luft bei Normathmosphäre nach DIN ISO 2533 [kg/m3]
pa=8*10^5;                  %Anfangsdruck im Flugkörper [Pa]
pa2=druck(max(Schritt_1),1);%Anfangsdruck Phase 2 [Pa]

%Konstanten
g=9.81;                     %Gravitationskonstante [m/s²]
rho=998.2;                  %Dichte des Wassers [kg/m3]


mb=Vw*rho;                  %Brennstoffmasse [kg]

m=ml+mb;                    %Anfangsmasse der Rakete [kg]

cw=0.34;                    %cw-Wert Halbkugel konvexe Seite
alpha=pi/4;                 %Abschusswinkel

%Phase2 Flug mit Treibstoff Luft

psimax=0.484;               %maximale Auslfussfunktion
mue=0.82;                   %Ausflusszahl (abgerundete, glatte Kurzdüse; 
                            %zylindrisches Mündungsstück)
RLuft=287.05;               %Gaskonstante Luft [J/(kg*K)]
T=293.15;                   %Temperatur [K]
mLufta2=pa2*Vges/(RLuft*T); %Masse Luft nach Brennschluss Wasser [kg]
rhoa2=mLufta2/Vges;         %Dichte Luft Anfang Phase 2 [kg/m3]
kappaLuft=1.402;            %Isentropenexponent Luft
pk=(2/(kappaLuft+1))^(kappaLuft/(kappaLuft-1));
pa2u=p0/pk;                 %Druck nach überkritischem Ausströmen [Pa]
                            %Dauer überkritisches Ausströmen [s]
format long                           
t2=(-1)*(Vges*log(pa2u/pa2)/(mue*A2*sqrt(2*RLuft*T)*(2/(kappaLuft+1))^(1/(kappaLuft-1))*sqrt(kappaLuft/(kappaLuft+1))));
rhokrit=0.634;              %kritisches Dichteverhältnis
mLufta2b=pa2u*Vges/(RLuft*T);%Masse Luft nach überkritischem Ausströmen

mlu=p0*Vges/(RLuft*T);