3 Katalytische Performance der Mo/V(/W)-Mischoxide - tuprints

188 Anhang
(MOX, MOXO, 18MOXO, Bulk, Bulk-O, Bulk-18O) nicht, d. h. ihre Verweilzeit im
Reaktor ist unbegrenzt.
Die Gesamtdauer des Experiments beträgt 1200 s. Der sprunghafte Tausch (Dirac-
Funktion) von 16 O2 gegen 18 O2 erfolgt nach 600 s und wird durch eine Gaussfunktion am
Eingang des Leitungssystems (R0) approximiert. Dafür wird je eine Randbedingung vom
Typ Danckwerts für die Sauerstoffkonzentrationen am Eingang als Funktion der Zeit
aufgestellt (kein fester Wert). Die Implementierung in PRESTO ® erfolgt in Form von
sogenannten Interpreterfunktionen. Für 16 O2 lautet die entsprechende Danckwerts-
Randbedingung „springen16Od.fun“:
w = getkp("w_v")
D = getkp("Dax_v")
beta = -w/D
result1 = beta
result2 = eval("Sprung16O2", acttime(),0)
Diese ruft wiederum die Interpreterfunktion „Sprung16O2“ auf, die in einer Bibliothek
eingetragen ist:
// arg1 und damit x ist die aktuelle Zeit in der Simulation
x = arg1
pos = 600
sigma = 0.2
result1 = 1
//
if (x > pos)
result1 = eval("Gauss", pos-x, sigma)
//
result1 = result1 * getkp("c0_16O2"),
Diese greift ihrerseits auf die Interpreterfunktion “Gauss” zu:
x = arg1
sigma = arg2
mean = 0
result1 = exp(-0.5*(x-mean)^2/sigma^2)
Die weiteren Eingabedaten für das Modell im ASCII Code sind im Folgenden verkürzt
dargestellt:
SSITKA_h4-05s_35_05_345
Endzeit 1200.00
Genauigkeit 2.0000e-02