Past exam of the numerical solution of differential equations course of the University of Oxford / 2011 / A / 3

Let ${\delta }_{xx}{U}_{ij}=(U_{i+1,j}-2U_{ij}+U_{i-1,j})/(\Delta x{)}^2$, and define ${\delta }_{yy}$ analogously. The Crank-Nicolson diffusion scheme isUse $U_{ij}^0=u_0(x_i,y_j)$ and zero values on every boundary node for every time level. This is the theta method with $\theta =1/2$ applied to the spatial finite difference method for the heat equation.

Write $s={\mu }_x+{\mu }_y$. The update can be rearranged asWhen $s\le 1$, the right side is a convex combination of old values: its coefficients are nonnegative and sum to one. Assume all old values lie in $[U_{\min },U_{\max }]$. A new maximum $v>U_{\max }$ is interior because boundary values are zero and $U_{\max }\ge 0$. At its node the left side is at least $(1+s)v-sv=v$, while the right side is at most $U_{\max }$, a contradiction. A new minimum $v<U_{\min }$ similarly makes the left side at most $v$ and the right side at least $U_{\min }$, again a contradiction.

The initial values satisfy the bounds by construction; mathematical induction provesThe new-time matrix is a strictly diagonally dominant matrix with positive diagonal and nonpositive off-diagonal entries, so the linear update is uniquely solvable. This Crank-Nicolson maximum principle under a step restriction concerns preservation of value bounds, which is a stronger requirement than scalar A-stability.

For nonnegative initial data, $U_{\min }=0$. The proved Crank-Nicolson maximum principle under a step restriction immediately gives $\boxed{U_{ij}^m\ge 0}$ for all indices. For nonpositive initial data, $U_{\max }=0$ and the same bound gives $\boxed{U_{ij}^m\le 0}$. Equivalently, apply positivity to $-U$, since the finite difference method is linear and has homogeneous Dirichlet boundary conditions. Both conclusions use ${\mu }_x+{\mu }_y\le 1$.