我想使用 FEM 方法（例如deal.ii ）求解以下非线性反应扩散方程组（Schnakenberg Turing ）：

$$\partial_{t} u = \Delta u + \gamma\left(a-u+u²v\right)$$ $$\partial_{t} v = d\Delta v + \gamma\left(b-u²v\right)$$

其中 d, $\gamma$ , a, b 是常数。

可能我需要先应用一个时间积分方案（例如 Crank-Nicholson，隐式 Runge-Kutta），然后是 FE 空间离散化。

问题：我应该如何选择和执行这个非线性系统的时间积分？非常相似的例子或特定文献的任何提示？