Hallo ans Forum,
ich beschäftige mich gerade mit der Austreibung der Volatilen aus einem Ash-Volatile-Gemisch. Dabei gibt es in OpenFOAM ja die zwei Modelle constantRate- und singleKineticRate-Devolatilization (~/OpenFOAM-1.7.0/src/lagrangian/intermediate/submodels/ReactingMultiphase/DevolatilisationModel/).
Dazu 3 Sachen:
In der singleKineticRateDevolatlization.H ist definiert, dass E der prä-exponentielle Faktor ist, und A1 die Aktivierungsenergie.
In der singleKineticRateDevolatilization.C allerdings taucht die Formel (eine einfache Arrheniusgleichung) in folgender Form auf:
const scalar kappa = A1_*exp(-E_/(specie::RR*T));
(1) Sind dabei nicht A1 und E vertauscht?
Mein Testcase ist folgendes: einfache Rohrgeometrie, in welche ich zu einem bestimmten Zeitpunkt ein Partikel einspeise (Zusammensetzung: 0.7 ash und 0.3 Vol). Meine Ergebnisse bei constantRateDevolatilization sehen gut aus, der Antail an Asche (Y(ash)) erhöht sich wie erwartet nach einer gewissen Zeit. Bei der singleKineticRateDevolatilization allerdings "springt" der Anteil an Asche von einem Zeitschritt auf den nächsten sofort auf 1.
(2) Woran kann das liegen? Sind evtl. Meine Werte für Aktivierungsenergie und Prä-exponentiellen Faktor falsch?
(3) Wo bekommt man normalerweise die Werte für A0, A1 und E her?
MfG,
STefan
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP