[2696] | 1 | % User-defined Rate Law functions |
---|
| 2 | % Note: insert this file at the end of Update_RCONST |
---|
| 3 | |
---|
| 4 | %--- Arrhenius |
---|
| 5 | function [rate] = ARR( A0,B0,C0 ) |
---|
| 6 | global TEMP CFACTOR |
---|
| 7 | rate = (A0) * exp(-(B0)/TEMP) * (TEMP/300.0)^(C0) ; |
---|
| 8 | return % ARR |
---|
| 9 | |
---|
| 10 | %--- Simplified Arrhenius, with two arguments |
---|
| 11 | %--- Note: The argument B0 has a changed sign when compared to ARR |
---|
| 12 | function [rate] = ARR2( A0,B0 ) |
---|
| 13 | global TEMP CFACTOR |
---|
| 14 | rate = (A0) * exp( (B0)/TEMP ) ; |
---|
| 15 | return % ARR2 |
---|
| 16 | |
---|
| 17 | function [rate] = EP2(A0,C0,A2,C2,A3,C3) |
---|
| 18 | global TEMP CFACTOR |
---|
| 19 | K0 = (A0) * exp(-C0/TEMP); |
---|
| 20 | K2 = (A2) * exp(-C2/TEMP); |
---|
| 21 | K3 = (A3) * exp(-C3/TEMP); |
---|
| 22 | K3 = K3*CFACTOR*1.0e+6; |
---|
| 23 | rate = K0 + K3/(1.0+K3/K2) ; |
---|
| 24 | return % EP2 |
---|
| 25 | |
---|
| 26 | function [rate] = EP3(A1,C1,A2,C2) |
---|
| 27 | global TEMP CFACTOR |
---|
| 28 | K1 = (A1) * exp(-(C1)/TEMP); |
---|
| 29 | K2 = (A2) * exp(-(C2)/TEMP); |
---|
| 30 | rate = K1 + K2*(1.0e+6*CFACTOR); |
---|
| 31 | return % EP3 |
---|
| 32 | |
---|
| 33 | function [rate] = FALL ( A0,B0,C0,A1,B1,C1,CF) |
---|
| 34 | global TEMP CFACTOR |
---|
| 35 | K0 = A0 * exp(-B0/TEMP)* (TEMP/300.0)^(C0); |
---|
| 36 | K1 = A1 * exp(-B1/TEMP)* (TEMP/300.0)^(C1); |
---|
| 37 | K0 = K0*CFACTOR*1.0e+6; |
---|
| 38 | K1 = K0/K1; |
---|
| 39 | rate = (K0/(1.0+K1))*(CF)^(1.0/(1.0+(log(K1))^2)); |
---|
| 40 | return % FALL |
---|
| 41 | |
---|