Fractional generalization of the fermi–Pasta–Ulam–Tsingou media and theoretical analysis of an explicit variational scheme

• Published in: Communications in nonlinear science & numerical simulation Vol. 88; p. 105158
• Main Author: Macías-Díaz, J.E.
• Format: Journal Article
• Language: English
• Published: Amsterdam Elsevier B.V 01.09.2020; Elsevier Science Ltd
• Subjects: Approximation; Boundary conditions; Computer simulation; Conservative fractional wave equation; Continuous fermi–Pasta–Ulam–Tsingou media; Damping; Derivatives; Discrete element method; Energy conservation; Energy dissipation; Energy-preserving method; Finite difference method; Fractional centered differences; Mathematical models; Nonlinear dynamics; Numerical efficiency analysis; Partial differential equations; Polynomials; Riesz space-fractional equations; 65Qxx; Energy-preserving method; Conservative fractional wave equation; Riesz space-fractional equations; Numerical efficiency analysis; Fractional centered differences; Continuous fermi–Pasta–Ulam–Tsingou media; 65Mxx
• ISSN: 1007-5704; 1878-7274
• DOI: 10.1016/j.cnsns.2019.105158

•A new fractional generalization of the FPUT is proposed.•We prove that the model has an associated variational structure.•A variational numerical scheme is proposed to solve the model.•The numerical properties of the scheme are thoroughly established.•A Matlab implementation of the scheme is provided in the appendix. In this work, we consider a partial differential equation that extends the well-known Fermi–Pasta–Ulam–Tsingou chains from nonlinear dynamics. The continuous model under consideration includes the presence of both a damping term and a polynomial function in terms of Riesz space-fractional derivatives. Initial and boundary conditions on a closed and bounded interval are considered in this work. The mathematical model has a fractional Hamiltonian which is conserved when the damping coefficient is equal to zero, and dissipated otherwise. Motivated by these facts, we propose a finite-difference method to approximate the solutions of the continuous model. The method is an explicit scheme which is based on the use of fractional centered differences to approximate the fractional derivatives of the model. A discretized form of the Hamiltonian is also proposed in this work, and we prove analytically that the method is capable of conserving or dissipating the discrete energy under the same conditions that guarantee the conservation or dissipation of energy of the continuous model. We show that solutions of the discrete model exist and are unique under suitable regularity conditions on the reaction function. We establish rigorously the properties of consistency, stability and convergence of the method. To that end, novel technical results are mathematically proved. Computer simulations that assess the capability of the method to preserve the energy are provided for illustration purposes.