DOI: https://doi.org/10.1186/s13661-025-02085-9
تاريخ النشر: 2025-07-01
المؤلف: H. M. Ahmed وآخرون
الموضوع الرئيسي: حلول المعادلات التفاضلية الكسرية
نظرة عامة
تقدم هذه الدراسة طريقة عددية جديدة لحل معادلات كاواهارا العامة غير الخطية ذات الفترات الزمنية (NTFGKE) تحت ظروف حدود أولية موحدة (IBCs) من خلال استخدام متعددة الحدود شرويدر المعدلة (MSPs). تستخدم الطريقة مصفوفات تشغيلية (OMs) لحساب كل من المشتقات الكسرية (FDs) والمشتقات العادية (ODs) لـ MSPs، المدمجة مع طريقة التجميع الطيفي (SCM) لتشكيل إطار حسابي شامل. يقدم المؤلفون تحليل التقارب والخطأ للطريقة المقترحة، موضحين دقتها من خلال ثلاث حالات اختبار عددية تقارن مع الحلول المعروفة. تشير النتائج إلى أن الطريقة تعزز بشكل كبير كل من الدقة والكفاءة الحسابية، مما يجعلها أداة قيمة لمعالجة المشكلات غير الخطية المعقدة.
في الختام، تؤكد الدراسة أن اختيار MSPs هو خيار نظري سليم وعملي مفيد لتحقيق حلول عددية دقيقة (NSols) لـ NTFGKE مع IBCs متجانسة. يضمن دمج SCM مع OMs لمشتقات MSPs تمثيل حلول موثوق. تحقق البحث فعالية الطريقة من خلال أمثلة عددية، موضحًا وجود علاقة قوية بين الحلول الدقيقة (ESol) والحلول التقريبية (ASol). تشير النتائج إلى أن هذه الطريقة لا تعزز فقط التقنيات العددية لمعادلات التفاضل الكسرية (FDEs) ولكنها تفتح أيضًا آفاقًا للبحث المستقبلي في تطبيقها في سيناريوهات أكثر تعقيدًا، مثل معادلات الانتشار-الموجة ذات الفترات الزمنية المتغيرة متعددة الحدود ونماذج غير خطية أخرى.
مقدمة
في مقدمة هذه الورقة البحثية، يناقش المؤلفون أهمية معادلات التفاضل الكسرية (FDEs) عبر مجالات علمية متنوعة، مؤكدين على ضرورة الطرق العددية للحصول على الحلول بسبب نقص الحلول التحليلية. يسلطون الضوء على منهجيات متنوعة مستخدمة في الأدبيات، بما في ذلك طرق الموجات، والمصفوفات التشغيلية، وتقنيات التنبؤ-التصحيح، وطرق التجميع، لمعالجة أنواع مختلفة من FDEs. ثم يتحول التركيز إلى معادلة كاواهارا (KE)، وهي معادلة زمنية كسرية حاسمة تمثل انتشار الموجات في الوسائط المشتتة، وتعميماتها، التي تم استكشافها من خلال مناهج عددية متنوعة.
يقدم المؤلفون إطارًا حسابيًا جديدًا يستخدم متعددة الحدود شرويدر المعدلة (MSPs) كدوال أساسية لاشتقاق حلول عددية دقيقة لمعادلة كاواهارا العامة غير الخطية ذات الفترات الزمنية (NTFGKE). يقترحون مصفوفات تشغيلية (OMs) مصممة لكل من المشتقات الكسرية والمشتقات العادية بمعنى ليوفيل-كابوتو، مما يعزز الاستقرار والدقة في التقريبات العددية. توضح الورقة مزايا طريقتهم مقارنة بالتقنيات العددية الحالية، خاصة في التعامل مع المشكلات غير الخطية المعقدة بكفاءة. كما يتم تقديم هيكل الورقة، موضحًا الأقسام التالية التي تغطي المفاهيم الأساسية، وتطوير OMs، والتقنيات العددية، وتقديرات الخطأ، والتحليلات المقارنة مع الطرق الحالية.
النتائج
في هذا القسم، يقدم المؤلفون تحليل التقارب وتقديرات الخطأ للطريقة المقترحة، مع التركيز على فضاء الدالة المرتبط $\mathcal{S}_\mathcal{Q} = \text{Span}\langle \mathcal{U}_r(z) \phi_s(t) | r, s = 0, 1, \ldots, Q \rangle$. يتم تعريف مصطلح الخطأ على أنه $\mathcal{E}_Q(z, t) = \|\mathcal{Y}_Q(z, t) – \mathcal{Y}(z, t)\|$، مع استخدام معيار L∞ لتقدير الخطأ، على وجه التحديد $\|\mathcal{E}_Q\|_\infty = \max_{(z, t) \in J} |\mathcal{Y}(z, t) – \mathcal{Y}_Q(z, t)|$، حيث $J = [0, 1] \times [0, 1]$. يثبت النظرية 6.1 حدًا أعلى للخطأ، المعطى بواسطة $\|\mathcal{E}_Q\|_\infty \leq C_4 Q \Gamma(Q + 1/2)$، حيث $C_4$ هو ثابت مشتق من المشتقات القصوى للدالة $\nu(z, t)$.
تشمل الإثبات تعريف متعددة الحدود $\nu_Q(z, t)$ التي تتداخل مع $\nu(z, t)$ عند نقاط محددة، مما يؤدي إلى تقدير للخطأ يتضمن مشتقات $\nu(z, t)$. تشير النتيجة النهائية إلى أن الخطأ بين الحل الفعلي $\mathcal{Y}(z, t)$ وأفضل تقريب له $\mathcal{Y}_Q(z, t)$ محصور بواسطة $\|\mathcal{Y} – \mathcal{Y}_Q\|_\infty \leq \frac{1}{4} Q \Gamma(Q + 1/2)$. هذا يثبت تقارب المخطط العددي مع زيادة $Q$، مما يؤكد فعالية الطريقة المقترحة للتوسع.
مناقشة
في هذا القسم، يناقش المؤلفون المفاهيم الأساسية لمشتقات ليوفيل-كابوتو الكسرية (FDs) ومتعددة الحدود شرويدر (SPs)، والتي تعتبر أساسية لطريقتهم المقترحة لحل معادلة كورتويغ-دي فريس العامة غير الخطية (NTFGKE). يتم تعريف المشتق الكسرى لليوفيل-كابوتو، مع تسليط الضوء على خطيته وخصائصه المحددة للدوال الثابتة والوظائف الشبيهة بالقوة. يقدم المؤلفون أيضًا SPs، موضحين صياغتها التكرارية وصيغتها العكسية، والتي تعتبر ضرورية لبناء المصفوفات التشغيلية (OMs) المستخدمة في عملية الحل العددي.
تتقدم المناقشة إلى اشتقاق المصفوفات التشغيلية للمشتقات الكسرية المرتبطة بالمتعددة الحدود \( \phi_j(t) \) و \( \mathcal{U}_i(z) \). يثبت المؤلفون النظريات والليمات الرئيسية التي تسهل تمثيل هذه المصفوفات، مع التأكيد على هيكلها المثلث السفلي، مما يبسط الحسابات ويعكس التفاعلات بين الحدود المتعددة. تظهر المصفوفات التشغيلية أنها نموذج فعال للسلوكيات الديناميكية في الأنظمة التي تحكمها معادلات التفاضل الكسرية، خاصة من حيث الذاكرة وعدم المحلية. يخلص المؤلفون إلى أن طريقتهم، المستندة إلى نهج التجميع باستخدام المصفوفات التشغيلية، تظهر دقة وكفاءة كبيرة في حل NTFGKE، كما يتضح من الأمثلة العددية والمقارنات مع الطرق الحالية.
DOI: https://doi.org/10.1186/s13661-025-02085-9
Publication Date: 2025-07-01
Author(s): H. M. Ahmed et al.
Primary Topic: Fractional Differential Equations Solutions
Overview
This study introduces a novel numerical method for solving nonlinear time-fractional generalized Kawahara equations (NTFGKE) under uniform initial boundary conditions (IBCs) by employing modified Schröder polynomials (MSPs). The method utilizes operational matrices (OMs) to compute both fractional derivatives (FDs) and ordinary derivatives (ODs) of MSPs, integrated with the spectral collocation method (SCM) to form a comprehensive computational framework. The authors provide a convergence and error analysis of the proposed approach, demonstrating its accuracy through three numerical test cases that benchmark against established solutions. Results indicate that the method significantly enhances both accuracy and computational efficiency, making it a valuable tool for addressing complex nonlinear problems.
In conclusion, the study asserts that the choice of MSPs is both theoretically sound and practically beneficial for achieving accurate numerical solutions (NSols) to NTFGKE with homogeneous IBCs. The integration of SCM with OMs for MSPs’ derivatives ensures reliable solution representation. The research validates the method’s effectiveness through numerical examples, showing a strong correlation between exact solutions (ESol) and approximate solutions (ASol). The findings suggest that this approach not only advances numerical techniques for fractional differential equations (FDEs) but also opens avenues for future research into its application in more complex scenarios, such as multiterm variable-order time-fractional diffusion-wave equations and other nonlinear models.
Introduction
In the introduction of this research paper, the authors discuss the significance of fractional differential equations (FDEs) across various scientific fields, emphasizing the necessity of numerical methods for obtaining solutions due to the lack of analytical solutions. They highlight diverse methodologies employed in the literature, including wavelet methods, operational matrices, predictor-corrector techniques, and collocation methods, to address different types of FDEs. The focus then shifts to the Kawahara equation (KE), a critical time-fractional equation that models wave propagation in dispersive media, and its generalizations, which have been explored through various numerical approaches.
The authors introduce a novel computational framework utilizing modified Schröder polynomials (MSPs) as basis functions to derive accurate numerical solutions for a specific nonlinear time-fractional generalized Kawahara equation (NTFGKE). They propose operational matrices (OMs) tailored for both fractional derivatives and ordinary derivatives in the Liouville-Caputo sense, which enhance stability and accuracy in numerical approximations. The paper outlines the advantages of their method over existing numerical techniques, particularly in handling complex nonlinear problems efficiently. The structure of the paper is also provided, detailing the subsequent sections that cover foundational concepts, the development of OMs, numerical techniques, error estimates, and comparative analyses with existing methods.
Results
In this section, the authors present a convergence analysis and error estimates for the proposed expansion method, focusing on the associated function space $\mathcal{S}_\mathcal{Q} = \text{Span}\langle \mathcal{U}_r(z) \phi_s(t) | r, s = 0, 1, \ldots, Q \rangle$. The error term is defined as $\mathcal{E}_Q(z, t) = \|\mathcal{Y}_Q(z, t) – \mathcal{Y}(z, t)\|$, with the L∞ norm used for error estimation, specifically $\|\mathcal{E}_Q\|_\infty = \max_{(z, t) \in J} |\mathcal{Y}(z, t) – \mathcal{Y}_Q(z, t)|$, where $J = [0, 1] \times [0, 1]$. Theorem 6.1 establishes an upper bound for the error, given by $\|\mathcal{E}_Q\|_\infty \leq C_4 Q \Gamma(Q + 1/2)$, where $C_4$ is a constant derived from the maximum derivatives of the function $\nu(z, t)$.
The proof involves defining the polynomial $\nu_Q(z, t)$ that interpolates $\nu(z, t)$ at specific points, leading to an error estimate that incorporates the derivatives of $\nu(z, t)$. The final result indicates that the error between the actual solution $\mathcal{Y}(z, t)$ and its best approximation $\mathcal{Y}_Q(z, t)$ is bounded by $\|\mathcal{Y} – \mathcal{Y}_Q\|_\infty \leq \frac{1}{4} Q \Gamma(Q + 1/2)$. This establishes the convergence of the numerical scheme as $Q$ increases, confirming the effectiveness of the proposed expansion method.
Discussion
In this section, the authors discuss the foundational concepts of Liouville-Caputo fractional derivatives (FDs) and Schröder polynomials (SPs), which are integral to their proposed method for solving the nonlinear time-fractional generalized Korteweg-de Vries equation (NTFGKE). The Liouville-Caputo fractional derivative is defined, highlighting its linearity and specific properties for constant and power-like functions. The authors also introduce SPs, detailing their recursive formulation and inversion formula, which are essential for constructing operational matrices (OMs) used in the numerical solution process.
The discussion progresses to the derivation of operational matrices for fractional derivatives associated with the polynomials \( \phi_j(t) \) and \( \mathcal{U}_i(z) \). The authors establish key theorems and lemmas that facilitate the representation of these matrices, emphasizing their lower triangular structure, which simplifies computations and captures the interactions among polynomial terms. The operational matrices are shown to effectively model dynamic behaviors in systems governed by fractional differential equations, particularly in terms of memory and nonlocality. The authors conclude that their method, based on the collocation approach using operational matrices, demonstrates significant accuracy and efficiency in solving the NTFGKE, as evidenced by numerical examples and comparisons with existing methods.