Energetische Nutzung von feuchter Biomasse in ... - tuprints

(* Tabellenausgabe *)
xdata = Flatten[Table[x, {x, 0, lnge0, (lnge0/30)}]]
ydata1 = Flatten[Table[(Evaluate[temp2[x] /. sol1]), {x, 0, lnge0, (lnge0/30)}]]
ydata2 = Flatten[Table[(Evaluate[yBP2[x] /. sol1]), {x, 0, lnge0, (lnge0/30)}]]
Null
Quellcode 2:
(* Heatric/Rohrbündel-WT: Berechnung von Eintrittsseite heissem Reaktoraustrag aus *)
(* mit rho = f(T) und cp = f(T) bei 50 MPa *)
(* Stoffdaten - Fits ohne Unstetigkeit *)
(* Strom 1 im Aussenbereich = Reaktoreffluent , Strom 2 in Innenrohren (Rohrbündel) = Biomassefeed *)
Anhang
Remove["Global`*"]; (* Löschen aller Variablen und Funktionen *)
lnge0 = 82 ; (* Reaktorlänge in m; variiert bis Feedeintrittstemperatur ca. 20 °C *)
rInn = 0.0075; (* Radius Innenrohr in m; hier Beispiel konventioneller RB-WT *)
rOut = rInn*2^0.5 (* Radius virtuelles Aussenrohr pro Rohrbündel-Rohr in m, so gewählt,
dass Fläche A1 = A2 *)
usoll = 0.2; (* angestrebte Mindest-Strömungsgeschwindigkeit in WT-Rohren in m s - 1 *)
mges = 10; (* Massenstrom gesamt, in t h - 1 *)
nRohre = mges/(3.6*1020*3.14*rInn^2*usoll) (* Rohranzahl berechnen... *)
(* nRohre = 55; *) (* ...oder direkt angeben ! *)
m = mges/3.6/nRohre; (* Massenstrom pro Rohr mges/N(Rohre) in kg s - 1 *)
k = 1000; (* Wärmedurchgangskoeffizient in W m - 2 K - 1; hier Beispiel konventioneller RB-WT *)
(* Wärmekapazität Wasser (abhängig von temp in K) in J kg - 1 K - 1 *)
... (siehe oben)
(* Dichte (abhängig von temp in K) in kg m - 3*)
... (siehe oben)
(* Geschwindigkeit im Aussenrohr in m s - 1 *)
uloc1[temp_] := m/(rho[temp]*Pi*(rOut^2 - rInn^2));
(* Geschwindigkeit in Innenrohren in m s - 1 *)
uloc2[temp_] := m/(rho[temp]*Pi*rInn^2);
(* DGL der Energiebilanzen *)
(* Startwerte bei x=0: Endgrädigkeit 50 °C, Feed tritt mit 550 °C aus *)
sol1 = NDSolve[{
temp1'[x] == ((k*2*rInn/(rOut^2 - rInn^2 ))/(m/(Pi*(rOut^2 – rInn^2))*cp[temp1[x]]))
*(temp2[x] - temp1[x]),
temp2'[x] == ((k*2/rInn)/(m/(Pi*rInn^2)*cp[temp2[x]]))*(temp2[x] - temp1[x]),
temp1[0] == 873, temp2[0] == 873 - 50},
{temp1, temp2}, {x, 0, lnge0}];
lsg = %[[1]];
ntemp1[x_] = temp1[x] /. lsg;
313