On the Extended Simple Equations Method (SEsM) for Obtaining Numerous Exact Solutions to Fractional Partial Differential Equations: A Generalized Algorithm and Several Applications

In this article, we present the extended simple equations method (SEsM) for finding exact solutions to systems of fractional nonlinear partial differential equations (FNPDEs). The expansions made to the original SEsM algorithm are implemented in several directions: (1) In constructing analytical solutions: exact solutions to FNPDE systems are presented by simple or complex composite functions, including combinations of solutions to two or more different simple equations with distinct independent variables (corresponding to different wave velocities); (2) in selecting appropriate fractional derivatives and appropriate wave transformations: the choice of the type of fractional derivatives for each system of FNPDEs depends on the physical nature of the modeled real process. Based on this choice, the range of applicable wave transformations that are used to reduce FNPDEs to nonlinear ODEs has been expanded. It includes not only various forms of fractional traveling wave transformations but also standard traveling wave transformations. Based on these methodological enhancements, a generalized SEsM algorithm has been developed to derive exact solutions of systems of FNPDEs. This algorithm provides multiple options at each step, enabling the user to select the most appropriate variant depending on the expected wave dynamics in the modeled physical context. Two specific variants of the generalized SEsM algorithm have been applied to obtain exact solutions to two time-fractional shallow-water-like systems. For generating these exact solutions, it is assumed that each system variable in the studied models exhibits multi-wave behavior, which is expressed as a superposition of two waves propagating at different velocities. As a result, numerous novel multi-wave solutions are derived, involving combinations of hyperbolic-like, elliptic-like, and trigonometric-like functions. The obtained analytical solutions can provide valuable qualitative insights into complex wave dynamics in generalized spatio-temporal dynamical systems, with relevance to areas such as ocean current modeling, multiphase fluid dynamics and geophysical fluid modeling.