Past exam of the mathematics course of the University of Cambridge / 2015 / iii / Paper 68 / 5

For energy stability for variable-coefficient reaction diffusion, let $h_x=1/(M+1)$, let $D_h$ be the Dirichlet discrete Laplacian, and set $V_h=\mathrm{diag}(a_m)$. Use the mesh-weighted inner product $(v,w{)}_h=h_x{\sum }_{m=1}^M{v}_m{w}_m$ with the corresponding discrete L2 norm. For zero boundary values, summation by parts givesFor a solution, or a difference of two solutions, the energy method givesThe Gronwall inequality therefore provesThe constant is independent of $M$, so this is stability of a numerical method on every fixed finite time interval. The bound allows physical growth if $a_{+}>0$; stability here does not mean uniform boundedness as $t\to \infty$.

There is also a uniform maximum-norm bound. At a positive spatial maximum the second difference is nonpositive. Applying this observation to $e^{-a_{+}t}v$ and to its negative gives $\Vert v(t){\Vert }_{\infty }\le e^{a_{+}t}\Vert v(0){\Vert }_{\infty }$. With a perturbing source $r(t)$, the Duhamel principle gives the corresponding initial-data bound plus ${\int }_0^t{e}^{a_{+}(t-s)}\Vert r(s)\Vert ds$ in either norm. Thus the estimate controls accumulated residuals as well as initial perturbations.

The original PDF has $u_m^{n+1}$ on the left of the update; the converted TeX omits its $+1$. Use the Forward Euler method with $k=\Delta t\ge 0$ and $\mu =k/h_x^2$. The explicit time stepping for bounded reaction diffusion update matrix isWhen $0\le \mu \le 1/2$, the heat update $H_h$ has nonnegative stencil weights $\mu ,1-2\mu ,\mu$ and row sums at most one. Hence its induced maximum operator norm is at most one. It is also a symmetric matrix with eigenvaluesall in $[-1,1]$, so its discrete L2 norm operator bound is at most one too. Let $A_{\ast }=\max (|a_0|,|a_{+}|)$. In either of these norms,Iteration gives the mesh-uniform stability estimateThe same estimate applies to differences of solutions. Adding a per-step defect $k{r}^n$ gives $\Vert e^n\Vert \le e^{A_{\ast }T}(\Vert e^0\Vert +k{\sum }_{j<n}\Vert r^j\Vert )$.

This proof only requires $\mu \le 1/2$ and bounded reaction coefficients. It does not assume the full update $E_h$ is nonnegative: a negative $a_m$ can make its middle stencil weight negative when $\mu =1/2$. Nor does it claim contractivity or an all-time bound for every coefficient. The required stability is a bound uniform in the discretization on fixed finite intervals.