An efficient two-grid fourth-order compact difference scheme with variable-step BDF2 method for the semilinear parabolic equation

Bingyin Zhang,Hongfei Fu
2023-10-14
Abstract:Due to the lack of corresponding analysis on appropriate mapping operator between two grids, high-order two-grid difference algorithms are rarely studied. In this paper, we firstly discuss the boundedness of a local bi-cubic Lagrange interpolation operator. And then, taking the semilinear parabolic equation as an example, we first construct a variable-step high-order nonlinear difference algorithm using compact difference technique in space and the second-order backward differentiation formula (BDF2) with variable temporal stepsize in time. With the help of discrete orthogonal convolution (DOC) kernels and a cut-off numerical technique, the unique solvability and corresponding error estimates of the high-order nonlinear difference scheme are established under assumptions that the temporal stepsize ratio satisfies rk < 4.8645 and the maximum temporal stepsize satisfies tau = o(h^1/2 ). Then, an efficient two-grid high-order difference algorithm is developed by combining a small-scale variable-step high-order nonlinear difference algorithm on the coarse grid and a large-scale variable-step high-order linearized difference algorithm on the fine grid, in which the constructed piecewise bi-cubic Lagrange interpolation mapping operator is adopted to project the coarse-grid solution to the fine grid. Under the same temporal stepsize ratio restriction rk < 4.8645 and a weaker maximum temporal stepsize condition tau = o(H^1.2 ), optimal fourth-order in space and second-order in time error estimates of the two-grid difference scheme is established if the coarse-fine grid stepsizes satisfy H = O(h^4/7). Finally, several numerical experiments are carried out to demonstrate the effectiveness and efficiency of the proposed scheme.
Numerical Analysis
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to develop an efficient two - grid fourth - order compact difference scheme, combined with the variable - step - size BDF2 method, for solving semi - linear parabolic equations. Specifically, the paper focuses on the following aspects: 1. **Lack of high - order two - grid difference algorithms**: - Due to the lack of analysis of appropriate mapping operators (such as the local bicubic Lagrange interpolation operator) between the coarse grid and the fine grid, there are few studies on high - order two - grid difference algorithms. 2. **Numerical stability and error estimation**: - The paper proposes and analyzes a new high - order two - grid difference algorithm and proves its numerical stability and error estimation. In particular, by introducing the discrete orthogonal convolution (DOC) kernel and truncation techniques, under the time - step - size ratio \( r_k := \frac{\tau_k}{\tau_{k - 1}} < 4.8645 \) and the maximum time - step - size condition \( \tau = o(h^{1/2}) \), the unique solvability and error estimation of the nonlinear compact difference scheme are established. 3. **Efficient computing methods**: - The paper combines a small - scale variable - step - size high - order nonlinear difference algorithm (on the coarse grid) and a large - scale variable - step - size high - order linearized difference algorithm (on the fine grid), and projects the coarse - grid solution onto the fine grid through the constructed piecewise bicubic Lagrange interpolation mapping operator, thereby developing an efficient two - grid high - order difference algorithm. 4. **Optimal error estimation**: - Under the weaker maximum time - step - size condition \( \tau = o(H^{1/2}) \) and the coarse - and - fine - grid - step - size relationship \( H = O(h^{4/7}) \), the fourth - order spatial and second - order temporal optimal error estimations of the two - grid difference scheme are established. ### Main contributions - A variable - step - size two - grid fourth - order compact difference method is proposed, which combines the compact difference scheme, the two - grid method, the variable - step - size BDF2 formula, and the piecewise bicubic Lagrange interpolation operator. - Under the condition that the nonlinear term satisfies the local continuity condition, by using the DOC kernel and truncation techniques, the unique solvability and error estimation of the nonlinear compact difference scheme are strictly proved. - Through the discussion of the boundedness of the proposed piecewise bicubic Lagrange interpolation operator, the optimal error estimation of the two - grid compact difference scheme is established. - The effectiveness and efficiency of the proposed method are verified by numerical experiments, and the computational efficiency and stability are compared with other schemes. These works provide an efficient and stable numerical method for solving semi - linear parabolic equations, especially suitable for cases requiring high precision and adaptive time - step - sizes.