振荡反应是七年前上基础课的时候的一个实习,实验结果做出来了、也用Athena 模拟出来了,但那个时候还不会MATLAB,想不通Athena几行代码为什么就可以模拟这个反应。
最近因为考虑一个自催化过程,才知道数学上为什么会是这种结果。考虑一个反应的步骤是
x --> y
y --> z
z + x --> y
速度常数分别是k1, k2, k3。如果k2>>k1并且k3>>k1 那么这个反应就是个自催化反应。如果k3>>k2,用“稳态近似”处理z,得到的方程是一个关于y的二阶常微分方程:
y'' + k1y' + k1k2y = 0
它的解带有三角函数项,所以在合适的条件下会表现出振荡——也就是振荡反应。
顺便鄙视一下搞“全球变暖”气候模型预测的二流研究人员——对一个无比复杂的系统搞模拟和预测,不骗行吗?【补充一下:全球温度几十年尺度的升高和降低的一个可能原因就是“大洋涛动”——也是一种振荡。】
https://m.sciencenet.cn/blog-70036-515015.html
上一篇:
今年“全球变暖”研究四大科学事件下一篇:
全球气候协议根本不可能达成