Energetische Nutzung von feuchter Biomasse in ... - tuprints

Anhang
mmGluc = 0.180 ; (* molare Masse Glucose in kg mol - 1 *)
(* Kinetik Braunprodukte *)
eaBP[temp_] := 88*10^3 /; temp (325 + 273); (* Aktivierungsenergie in J mol - 1 *)
k0BP[temp_] := 423 /; temp (325 + 273);
nBP[temp_] := 3 /; temp (325 + 273);
(* ab 420°C bleibt C - Bilanz konstant *)
fBP[temp_] := 1 /; temp < (420 + 273);
fBP[temp_] := 0 /; temp >= (420 + 273);
(* 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 Energie- und Komponentenbilanzen *)
(* Startwerte bei x=0: Anfangsgrädigkeit 55 °C, Feed tritt mit ca. 20 °C ein (hier:27 °C), 10 % (g g-1) Glucose
oder Cellulose, 100 % C-Bilanz bzw. keine Braunprodukte *)
sol1 = NDSolve[{
temp1'[x] == ((k*2*rInn/(rOut^2 - rInn^2 ))/(m/(Pi*(rOut^2 - rInn^2))*cp[temp1[x]]))
*(temp1[x] - temp2[x]),
temp2'[x] == ((k*2/rInn)/(m/(Pi*rInn^2)*cp[temp2[x]]))*(temp1[x] - temp2[x]),
yGluc'[x] == -k0*Exp[-ea/(rGas*temp2[x])]*yGluc[x]/uloc2[temp2[x]],
yCell'[x] == -k0/10*Exp[-ea/(rGas*temp2[x])]*yCell[x]/uloc2[temp2[x]],
yBP2'[x] == -fBP[temp2[x]]*k0BP[temp2[x]]*Exp[-eaBP[temp2[x]]/(rGas*temp2[x])]
*yBP2[x]^nBP[temp2[x]]/uloc2[temp2[x]],
temp1[0] == 300 + 55, temp2[0] == 300, yGluc[0] == 0.10, yCell[0] == 0.10, yBP2[0] == 100},
{temp1, temp2, yGluc, yCell, yBP2}, {x, 0, lnge0}];
lsg = %[[1]];
ntemp1[x_] = temp1[x] /. lsg;
ntemp2[x_] = temp2[x] /. lsg;
(* Graphen ausgeben *)
Plot[Evaluate[uloc1[ntemp1[x]]], {x, 0, lnge0}, AspectRatio -> 0.5, PlotStyle -> {Thickness[0.002]},
AxesLabel -> {"Weg / m", "uloc 1 / m s-1"}];
Plot[Evaluate[uloc2[ntemp2[x]]], {x, 0, lnge0}, AspectRatio -> 0.5, PlotStyle -> {Thickness[0.002]},
AxesLabel -> {"Weg / m", "uloc 2 / m s-1"}];
Plot[Evaluate[NIntegrate [(1/uloc1[ntemp1[x]]), {x,0, lngex}]], {lngex, 0,lnge0}, AspectRatio -> 0.5,
PlotStyle -> {Thickness[0.002]}, AxesLabel -> {"Weg / m", "VWZ1 / s"}];
Plot[Evaluate[NIntegrate [(1/uloc2[ntemp2[x]]), {x, 0, lngex}]], {lngex, 0, lnge0}, AspectRatio -> 0.5,
PlotStyle -> {Thickness[0.002]}, AxesLabel -> {"Weg / m", "VWZ2 / s"}];
Plot[{(Evaluate[temp1[x] /. sol1] - 273), (Evaluate[temp2[x] /. sol1] - 273)}, {x, 0,lnge0},
AspectRatio -> 0.5, PlotStyle -> { {RGBColor[1, 0, 0]}, {Thickness[0.005]}}, AxesLabel ->
{"Weg / m", "Temp. / °C"}, PlotRange -> {0, 600}];
Plot[{Evaluate[yGluc[x] /. sol1]} , {x, 0, lnge0}, AspectRatio -> 0.5, PlotStyle -> {Thickness[0.005]},
AxesLabel -> {"Weg / m", "yGlucose / kg kg-1"}, PlotRange -> {0, 0.1}];
Plot[{Evaluate[yCell[x] /. sol1]} , {x, 0, lnge0}, AspectRatio -> 0.5, PlotStyle -> {Thickness[0.005]},
AxesLabel -> {"Weg / m", "yCellulose / kg kg-1"}, PlotRange -> {0, 0.1}];
Plot[{Evaluate[yBP2[x] /. sol1]} , {x, 0, lnge0}, AspectRatio -> 0.5, PlotStyle -> {Thickness[0.005]},
AxesLabel -> {"Weg / m", "C-Bilanz / % (mol mol-1)"}, PlotRange -> {0, 100}];
312