DOI: https://doi.org/10.3934/math.2026922
تاريخ النشر: 2026-01-01
المؤلف: Xinpei Wu وآخرون
الموضوع الرئيسي: حلول المعادلات التفاضلية الكسرية
نظرة عامة
في هذه الورقة، يستقصي المؤلفون طريقة عددية لمعادلة الانتشار الزمنية الكسرية العادية، مع التركيز على خصائصها الرياضية والعددية. تم تصميم الطريقة لتكون مستقرة بلا شروط وتلتزم بمبدأ الحد الأقصى، باستخدام دالة وزن للمشتق العادي تحافظ على قيمة تكامل واحدة عبر جميع الأوامر الكسرية والأوقات. تسهل هذه الخاصية مقارنة عادلة لسلوكيات الانتشار عبر أوامر كسرية متغيرة. يقوم المؤلفون بتنفيذ مخطط فرق محدود مع تمييز زمني ضمني، مما يؤدي إلى نظام جبر ثلاثي الأبعاد يتم حله بكفاءة باستخدام خوارزمية توماس.
تؤكد التحليل النظري أن المخطط العددي المقترح مستقر بلا شروط، ويحافظ على مبدأ الحد الأقصى، ويظهر تقاربًا صحيحًا، مع أوامر تقارب من \(O(\Delta t^{2-\alpha})\) في الزمن و \(O(h^2)\) في الفضاء. تؤكد الأمثلة العددية النتائج النظرية، مما يظهر استقرار الطريقة وتقاربها بينما توضح تأثيرات الأوامر الكسرية المختلفة على عمليات الانتشار. يعمل النموذج العادي كإطار قوي لتحليل مشاكل الانتشار الزمنية الكسرية، مع خطط للبحث المستقبلي لتوسيع النهج إلى نماذج ذات أبعاد أعلى، ومشاكل غير خطية، ومعادلات ذات أوامر متغيرة. بالإضافة إلى ذلك، يهدف المؤلفون إلى معالجة الانفرادات الضعيفة الأولية بالقرب من \(t = 0\) من خلال دمج شبكات زمنية متدرجة غير متساوية أو تقنيات تصحيح أولية لضمان تقارب مثالي من الخطوة الزمنية الأولى.
مقدمة
تناقش مقدمة الورقة أهمية معادلات الانتشار الكسرية في نمذجة عمليات الانتشار الشاذة التي تتميز بتأثيرات الذاكرة. ركزت التطورات الأخيرة على معادلات الانتشار الزمنية الكسرية ذات الأوامر المتغيرة، والتي تقدم تمثيلًا أكثر واقعية لتأثيرات الذاكرة غير المتجانسة مقارنة بالنماذج ذات الأوامر الثابتة. ومن الجدير بالذكر أن Zheng وآخرين تناولوا مشكلة الانفرادات غير الفيزيائية في الزمن الأولي، بينما استكشف Hong وآخرون المشاكل العكسية المتعلقة بالأوامر المتغيرة مكانيًا. على الرغم من مزايا نماذج الأوامر المتغيرة، لا تزال هناك تحديات في عزل تأثيرات الأمر الكسرى على ديناميات الانتشار بسبب اعتماد نواة الذاكرة على الأمر الكسرى.
لمعالجة هذه التحديات، اقترح Lee و Kim نموذج انتشار كسرية زمنية عادية ذات أوامر متغيرة، والذي يفصل المساهمة الإجمالية لنواة الذاكرة عن الأمر الكسرى، مما يسهل مقارنة أكثر اتساقًا لسلوكيات الانتشار. تهدف هذه الدراسة إلى تحليل معادلة الانتشار الزمنية الكسرية العادية المعرفة على مجال مكاني أحادي البعد، مقدمة نتائج نظرية حول الاستقرار بلا شروط، والحفاظ على مبدأ الحد الأقصى، وتقارب الطريقة العددية. كما تسلط الورقة الضوء على أهمية الحفاظ على الصلاحية الفيزيائية في الحلول العددية وتظهر دقة الطريقة من خلال اختبارات عددية متنوعة، مما يمهد الطريق لاستكشافات إضافية في الأقسام التالية.
نقاش
في هذا القسم، يقدم المؤلفون طريقة فرق محدود (FDM) لحل معادلة الانتشار الزمنية الكسرية العادية، مع التركيز على استقرارها بلا شروط، ومبدأ الحد الأقصى، وخصائص التقارب. تستخدم الطريقة تمييزًا زمنيًا ضمنيًا وتحل نظام مصفوفة ثلاثي الأبعاد باستخدام خوارزمية توماس. يقوم المؤلفون بإنشاء شبكة موحدة على المجال الحسابي $\Omega = (L, R)$، مع حجم شبكة $h = \frac{R – L}{N – 1}$، ويشيرون إلى الحل عند نقاط الشبكة المتقطعة كـ $u^n_i = u(x_i, t_n)$. يستنتجون صيغة تكاملية لتقريب المشتق الزمني الكسرى، مما يؤدي إلى خطأ اقتطاع $R \Delta t$ الذي يتم تحليله باستخدام توسع تايلور تحت فرضية وجود نعومة كافية للحل.
تؤكد تحليل الاستقرار أن المخطط الضمني مستقر بلا شروط، مع الحفاظ على مبدأ الحد الأقصى. يظهر المؤلفون أن معاملات النظام الثلاثي الأبعاد إيجابية وأن المخطط ينتج حلاً فريدًا في كل خطوة زمنية. يثبتون حد الخطأ العالمي، مما يظهر أن الحل العددي يتقارب بمعدل $O(\Delta t^{2-\alpha} + h^2)$، حيث $\alpha$ هو الأمر الكسرى. يختتم القسم باختبارات حسابية تؤكد أداء الطريقة، بما في ذلك اختبارات التقارب مع ظروف أولية سلسة وغير سلسة، ومقارنات مع مخطط صريح، مما يبرز مزايا النهج الضمني من حيث الكفاءة الحسابية والاستقرار.
DOI: https://doi.org/10.3934/math.2026922
Publication Date: 2026-01-01
Author(s): Xinpei Wu et al.
Primary Topic: Fractional Differential Equations Solutions
Overview
In this paper, the authors investigate a numerical method for a normalized time-fractional diffusion equation, emphasizing its mathematical and numerical properties. The method is designed to be unconditionally stable and adheres to the maximum principle, utilizing a weighting function for the normalized derivative that maintains an integral value of one across all fractional orders and times. This characteristic facilitates a fair comparison of diffusion behaviors across varying fractional orders. The authors implement a finite difference scheme with implicit time discretization, resulting in a tridiagonal algebraic system efficiently solved using the Thomas algorithm.
The theoretical analysis confirms that the proposed numerical scheme is unconditionally stable, preserves the maximum principle, and exhibits correct convergence, with convergence orders of \(O(\Delta t^{2-\alpha})\) in time and \(O(h^2)\) in space. Numerical examples corroborate the theoretical findings, demonstrating the method’s stability and convergence while illustrating the effects of different fractional orders on diffusion processes. The normalized model serves as a robust framework for analyzing time-fractional diffusion problems, with plans for future research to extend the approach to higher-dimensional models, nonlinear problems, and variable-order equations. Additionally, the authors aim to address initial weak singularities near \(t = 0\) by incorporating non-uniform graded time meshes or initial correction techniques to ensure optimal convergence from the first time step.
Introduction
The introduction of the paper discusses the significance of fractional diffusion equations in modeling anomalous diffusion processes characterized by memory effects. Recent advancements have focused on variable-order time-fractional diffusion equations, which offer a more realistic representation of heterogeneous memory effects compared to constant-order models. Notably, Zheng et al. addressed the issue of nonphysical initial-time singularities, while Hong et al. explored inverse problems related to spatially variable orders. Despite the advantages of variable-order models, challenges remain in isolating the effects of fractional order on diffusion dynamics due to the dependency of the memory kernel on the fractional order.
To address these challenges, Lee and Kim proposed a normalized variable-order time-fractional diffusion model, which decouples the total contribution of the memory kernel from the fractional order, thereby facilitating a more consistent comparison of diffusion behaviors. This study aims to analyze the normalized time-fractional diffusion equation defined on a one-dimensional spatial domain, presenting theoretical results on unconditional stability, maximum principle preservation, and convergence of the numerical method. The paper also highlights the importance of maintaining physical validity in numerical solutions and demonstrates the method’s accuracy through various numerical tests, ultimately setting the stage for further exploration in subsequent sections.
Discussion
In this section, the authors present a finite difference method (FDM) for solving the normalized time-fractional diffusion equation, emphasizing its unconditional stability, maximum principle, and convergence properties. The method employs an implicit time discretization and solves a tridiagonal matrix system using the Thomas algorithm. The authors construct a uniform grid over the computational domain $\Omega = (L, R)$, with mesh size $h = \frac{R – L}{N – 1}$, and denote the solution at discrete grid points as $u^n_i = u(x_i, t_n)$. They derive a quadrature formula to approximate the fractional time derivative, leading to a truncation error $R \Delta t$ that is analyzed using Taylor expansion under the assumption of sufficient smoothness of the solution.
The stability analysis confirms that the implicit scheme is unconditionally stable, with the maximum principle preserved. The authors demonstrate that the coefficients of the tridiagonal system are positive and that the scheme yields a unique solution at each time step. They establish a global error bound, showing that the numerical solution converges with a rate of $O(\Delta t^{2-\alpha} + h^2)$, where $\alpha$ is the fractional order. The section concludes with computational tests validating the method’s performance, including convergence tests with both smooth and non-smooth initial conditions, and comparisons with an explicit scheme, highlighting the advantages of the implicit approach in terms of computational efficiency and stability.