Grünwald-Letnikov Fractional Derivative Applied to First-Order Ordinary Differential Equations

First-order ordinary differential equations (ODEs) are widely used in various fields of science and engineering to model natural phenomena. This study proposes an extension of these equations using fractional derivatives, specifically the Grünwald-Letnikov definition, to explore their impact on the...

Descripción completa

Detalles Bibliográficos
Autor: Willian Scotton, Jaque
Tipo de recurso: artículo
Estado:Versión publicada
Fecha de publicación:2024
País:Brasil
Institución:Universidade Estadual de Londrina (UEL)
Repositorio:Revista Semina: Ciências Exatas e Tecnológicas (Online)
Idioma:inglés
OAI Identifier:oai:ojs2.ojs.uel.br:article/51533
Acceso en línea:https://ojs.uel.br/revistas/uel/index.php/semexatas/article/view/51533
Access Level:acceso abierto
Palabra clave:differential equations
fractional calculus
Grünwald-Letnikov
equações diferenciais
cálculo fracionário
Descripción
Sumario:First-order ordinary differential equations (ODEs) are widely used in various fields of science and engineering to model natural phenomena. This study proposes an extension of these equations using fractional derivatives, specifically the Grünwald-Letnikov definition, to explore their impact on the behavior of solution curves. The fractional ODE considered is discretized using the finite difference method and solved numerically for different values of the derivative order (α). Tests were conducted to verify mesh independence and the quality of the computational implementation of the method, through which the accuracy and the absence of implementation errors were confirmed. The behavior of the solution curves for different values of α was analyzed, revealing a sharp decrease near the initial point (t = 0) and an almost linear growth at higher values of t, within the considered domain. Additionally, in solving a specific initial value problem with a known analytical solution, it was discovered that the accuracy of the numerical solutions for higher values of α was more dependent on the mesh refinement than the solutions for lower values.