Iterative approximate factorization for difference operators of high-order bicompact schemes for multidimensional nonhomogeneous hyperbolic systems

An iterative method for solving equations of multidimensional bicompact schemes based on an approximate factorization of their difference operators is proposed for the first time. Its algorithm is described as applied to a system of two-dimensional nonhomogeneous quasilinear hyperbolic equations. The convergence of the iterative method is proved in the case of the two-dimensional homogeneous linear advection equation. The performance of the method is demonstrated on two numerical examples. It is shown that the method preserves a high (greater than the second) order of accuracy in time and performs 3–4 times faster than Newton’s method. Moreover, the method can be efficiently parallelized.