A splitting moving mesh method for reaction-diffusion equations of quenching type (Q2495754)

The authors study the following degenerate reaction-diffusion equation: \(\varphi(x,y)u_t=u_{xx}/a^2+u_{yy}/b^2+f(u)\) in \(Dx(0,T)\), \(T<\infty\), \(D=(0,1)^2\). Here \(\varphi(x,y)=(a^2x^2+b^2y^2)^{9/2}\) \((q\geq 0)\), \(f(u)\) is strictly increasing for \(0\leq u<1\), \(f(0)>0\), \(\lim_{u\to 1}=f(u)=\infty\). Initial and boundary conditions are taken as zero. Replacing the spatial derivatives by some finite differences on a nonuniform mesh, using a trapezoidal rule for an integral appearing in the formal solution to the obtained semi-discretized system and suitably approximating the matrix exponential they get a variable temporal step with an adaptive moving mesh in space Peaceman-Rachford splitting scheme combined with the bilinear interpolation. The temporal and spatial adaptation is implemented using an arc-length type of the monitor function on \(u_t\). Imitating the physical behaviour of the solution to the differential problem it is proved that the solution to the discrete problem increases monotonically in time before quenching when the temporal step is sufficiently small and further assumptions are satisfied. (A time \(0<t_0\leq T\) is called the quenching time if \(\sup\{|u_t(x,y,t)|:(x,y)\leq\overline D\}\to\infty\) as \(t\to t^-_0\).) Next a stability analysis for the discrete solution both away and near the quenching time is performed. Finally, a parallel implementation and numerical experiments for three interesting (also from physical point of view) problems are investigated.