离散型化学动力学标定模型的还原能力检验
标签:
动力学标定最小二乘法 |
分类: 语言程序 |
和之前检验高斯数值模型的还原精度一样,检验离散型动力学模型的还原精度依然根据已知的动力学模型生成不同的升温速率曲线,然后用模型标定,看结果是否能够还原这些动力学参数。为此,我们选了一个极端的情况,活化能分布比较散(表1),然后,转化率曲线也比较奇怪(图1)。图1是标定结果与原始模型的对比,表2是标定出来的活化能分布及指前因子。可见选择的数值模型是能够还原真实的化学动力学参数。
| 表1 测试动力学参数 | |
| A=1e14 | |
| Ea(kj) | Frequency |
| 160 | 35 |
| 221 | 45 |
| 250 | 15 |
| 260 | 5 |
表2 标定的动力学参数
| Ei(KJ) | Xi(%) | Ai(/s) |
| 141 | 0 | 1.20E+14 |
| 145 | 0 | 1.20E+14 |
| 149 | 0 | 1.20E+14 |
| 153 | 0 | 1.20E+14 |
| 157 | 2.173 | 1.20E+14 |
| 161 | 32.795 | 1.20E+14 |
| 165 | 0 | 1.20E+14 |
| 169 | 0 | 1.20E+14 |
| 173 | 0 | 1.20E+14 |
| 177 | 0 | 1.20E+14 |
| 181 | 0 | 1.20E+14 |
| 185 | 0 | 1.20E+14 |
| 189 | 0 | 1.20E+14 |
| 193 | 0 | 1.20E+14 |
| 197 | 0 | 1.20E+14 |
| 201 | 0 | 1.20E+14 |
| 205 | 0 | 1.20E+14 |
| 209 | 0 | 1.20E+14 |
| 213 | 0 | 1.20E+14 |
| 217 | 0 | 1.20E+14 |
| 221 | 32.828 | 1.20E+14 |
| 225 | 12.477 | 1.20E+14 |
| 229 | 0 | 1.20E+14 |
| 233 | 0 | 1.20E+14 |
| 237 | 0 | 1.20E+14 |
| 241 | 0 | 1.20E+14 |
| 245 | 0 | 1.20E+14 |
| 249 | 3.94 | 1.20E+14 |
| 253 | 12.234 | 1.20E+14 |
| 257 | 0.378 | 1.20E+14 |
| 261 | 1.599 | 1.20E+14 |
| 265 | 1.576 | 1.20E+14 |
| 269 | 0 | 1.20E+14 |
| 273 | 0 | 1.20E+14 |
| 277 | 0 | 1.20E+14 |
| 281 | 0 | 1.20E+14 |
| 285 | 0 | 1.20E+14 |
| 289 | 0 | 1.20E+14 |
| 293 | 0 | 1.20E+14 |
| 297 | 0 | 1.20E+14 |
| 301 | 0 | 1.20E+14 |
| 305 | 0 | 1.20E+14 |

加载中…