Abstract
This paper investigates a solution technique for solving a two-dimensional Kuramoto-Sivashinsky equation discretized using a finite difference method. It consists of an order reduction method into a coupled system of second-order equations, and to formulate the fully discretized, implicit time-marched system as a Lyapunov-Sylvester matrix equation. Convergence and stability is examined using Lyapunov criterion and manipulating generalized Lyapunov-Sylvester operators. Some numerical implementations are provided at the end to validate the theoretical results.