Finite Difference Scheme for One-Dimensional Coupled Parabolic System with Blow-up

This study aims to find an efficient technique to estimate the blow-up time (BUT) of one-dimensional semi-linear coupled parabolic systems. Firstly, a fully discrete finite difference formula is derived with a non-fixed time-stepping formula, based on the Crank-Nicolson method. In addition, the consistency, stability, and convergence of the proposed scheme are considered. Secondly, two numerical experiments are presented. For each experiment, we apply the proposed scheme to calculate the numerical blow-up time, error bounds, and the numerical order of convergence for blow-up times. The obtained results show that the proposed C.N scheme is consistent with the system considered. However, it is conditionally stable, and the Crank-Nicolson scheme converges in the stability region and achieves first-and second-order accuracy in temporal and spatial dimensions, respectively. Also, it helps to increase the order of numerical convergence. Furthermore, the numerical experiments demonstrate that the numerical blow-up simultaneously occurs at only the center point. Finally, the numerical blow-up time sequence is convergent. Moreover, convergence for the blow-up time agrees well with the theoretical order of convergence of the proposed scheme.