db2b03d2-0682-4488-87bf-e151e9dd4fae 29 3 5 decimal

&[DATE] &[TIME] - &[FILENAME]

http://twt.mpei.ac.ru/TTHB/1/Chem-Kin-Chapter-3.pdf

a 1 - 10001 - 10001 - 134 mat :

a transpose k1 X 1 el * k2 X 2 el * k3 X 3 el * 31 mat * k1 X 1 el * - k1 X 1 el * k2 X 2 el * - k2 X 2 el * k3 X 3 el * - k3 X 3 el * 41 mat

t X F k1 X 1 el * - k1 X 1 el * k2 X 2 el * - k2 X 2 el * k3 X 3 el * - k3 X 3 el * 41 mat :

k1 0.15 : k2 0.1 : k3 0.075 : 13 mat IC 400014 mat transpose : 21 sys

S IC 0100500 t X F rkfixed :

x S 1 col : S1 x S 2 col augment : S2 x S 3 col augment : S3 x S 4 col augment : S4 x S 5 col augment : 51 sys

S1 S2 S3 S4 4 1 sys

k0 0.5 : k1 0.1 : k2 0.2 : 13 mat k3 0.025 : k4 0.018 : k5 0.017 : 13 mat 21 sys

t X D k1 X 1 el 2^* k2 X 1 el * X 2 el * - k0 + k3 X 1 el 2^* k4 X 2 el 2^* - k5 X 1 el * X 2 el * - 21 mat :

S 22821 mat 02502500 t X D rkfixed :

x S 1 col : S1 x S 2 col augment : S2 x S 3 col augment : S3 S 2 col S 3 col augment : 41 sys

S1 S2 2 1 sys S3

ε1 31012^* : ε1 61014^* : 12 mat

IC 0.002 00 T0 14 mat transpose :

R 1.986 : T0 298 : γ 10 : 13 mat

t X F ε1 - (29000 R X 4 el * / - exp * X 1 el * ε1 (29000 R X 4 el * / - exp * X 1 el * ε1 (48000 R X 4 el * / - exp * X 2 el * - ε1 (48000 R X 4 el * / - exp * X 2 el * γ 41 mat :

N IC 0501000 t X F dn_GearsBDF :

x N 1 col : N1 x N 2 col augment : N2 x N 3 col augment : N3 x N 4 col augment : 41 sys

N1 N2 N3 3 1 sys

k 1011^: E 80000 : R 8.314 : T 293 : 14 mat

t Y F k E R Y 2 el * / - exp * 1 Y 1 el - (* Y 3 el 031 mat :

γ IC 0 T γ 31 mat :

γ sol γ IC 0801000 t Y F rkfixed :

S1 1 sol : S2 1.5 sol : S3 2 sol : x S1 1 col : 41 sys

sols x S1 2 col augment x S2 2 col augment x S3 2 col augment 31 sys :

other S1 3 col S1 2 col augment S2 3 col S2 2 col augment S3 3 col S3 2 col augment 31 sys :

sols other

k 105^: J 100 : ET 5000 : T0 295 : C 1 : 15 mat

X r k ET X 2 el / - exp * X 1 el * :

t X f X r - J X r * 21 mat :

t0 0 : t1 100 : N 1000 : 13 mat

S C T0 21 mat t0 t1 N t X f dn_AdamsMoulton :

t S 1 col : CA t S 2 col augment : T t S 3 col T0 - augment : Y T 1 col 1 J / T 2 col * augment : 41 sys

CA Y 2 1 sys

k 105^: J 100 : ET 5000 : TS 295 : C 1 : 15 mat

α 0.1 : V 1 : ρ 1 : Cp 4 : S 2 : 15 mat

X r k ET X 2 el / - exp * X 1 el * :

t X f X r - J X r * α S * TS X 2 el - (* ρ Cp * V * / + 21 mat :

t0 0 : t1 300 : N 1000 : 13 mat

S C TS 21 mat t0 t1 N t X f dn_AdamsMoulton :

t S 1 col : CA t S 2 col augment : T t S 3 col TS - augment : Y T 1 col 11 α / Cp + / (T 2 col * augment : Y2 t S 3 col TS - augment : 51 sys

CA Y 2 1 sys

1 α / Cp + 14