Prove the augmented-boundary error estimate and recover fourth-order convergence
Establish a rigorous proof that the augmented boundary-variable error for the augmented high-order compact scheme satisfies the discrete estimate \(\|\mathbf{E}^q\|_2=O(h^{7/2})\), and determine whether a sharper treatment of the boundary contribution yields the optimal fourth-order error estimate suggested by the numerical results.
References
The convergence analysis yields a conditional $O(h{7/2})$ error estimate, provided that the augmented boundary variables satisfy an $O(h{7/2})$ discrete $L2$-error bound. Although a rigorous proof of this boundary estimate is not yet available, it is consistently supported by the numerical results, under this assumption, the remaining convergence analysis is rigorous. The resulting estimate does not recover the full fourth-order rate observed numerically. This loss of half an order appears to arise from a non-sharp treatment of the boundary contribution. The numerical results suggest that the theoretical bound is not optimal and that the actual convergence rate is fourth order.