A Discrete Model to Solve a Bifractional Dissipative Sine-Gordon Equation: Theoretical Analysis and Simulations

In this work, we consider a generalized form of the classical (2+1)-dimensional sine-Gordon system. The mathematical model considers a generalized reaction term, and the two-dimensional Laplacian includes the presence of space-fractional derivatives of the Riesz type with two different differentiation orders in general. The system is equipped with a conserved quantity that resembles the energy functional in the integer-order scenario. We propose a numerical model to approximate the solutions of the fractional sine-Gordon equation. A discretized form of the energy-like quantity is proposed, and we prove that it is conserved throughout the discrete time. Moreover, the analysis of consistency, stability, and convergence is rigorously carried out. The numerical model is implemented computationally, and some computer simulations are presented in this work. As a consequence of our simulations, we show that the discrete energy is approximately conserved throughout time, which coincides with the theoretical results.