DOI: https://doi.org/10.4208/jcm.2510-m2025-0131
تاريخ النشر: 2026-03-15
المؤلف: Weiping Bu وآخرون
الموضوع الرئيسي: حلول المعادلات التفاضلية الكسرية
نظرة عامة
في هذه الدراسة، يبحث المؤلفون في معادلة تحت الانتشار تتميز بتأخير زمني ثابت، $\tau$. يقومون بتحليل انتظام الحل، كاشفين أن المشتق الزمني من الدرجة الأولى يظهر تفردًا عند $t = 0^+$، بينما يظهر المشتق الزمني من الدرجة الثانية تفردات عند كل من $t = 0^+$ و $t = \tau^+$. يتم تقسيم الحل إلى مكونات تفردية ومنتظمة، مما يسهل فهم سلوكه بشكل أوضح.
لمعالجة المشكلة عدديًا، يقوم المؤلفون بتطوير مخطط عنصر نهائي كامل الاستخدام يعتمد على طريقة غاليركين القياسية في الأبعاد المكانية وتقريب من نوع غرينوالد-ليتنكوف للتفريق الزمني. يتم إثبات استقرار هذا المخطط العددي، ويتم تقديم عدم المساواة الجديدة من نوع غروينوال لتقديم تقدير خطأ محلي في الزمن. تختتم الورقة باختبارات عددية تؤكد النتائج النظرية، مما يوضح أن الطريقة المقترحة تلتقط ديناميات معادلة الانتشار الزمنية الكسرية مع تأخير زمني ثابت بشكل فعال.
مقدمة
تسلط مقدمة الورقة الضوء على أهمية المعادلات التفاضلية، وخاصة المعادلات التفاضلية ذات التأخير، في نمذجة الظواهر المعقدة عبر مجالات مختلفة مثل الفيزياء، وعلم الأحياء، والطب. يؤكد المؤلفون على مزايا المعادلات التفاضلية الكسرية مقارنة بالنماذج ذات الترتيب الصحيح، خاصة في التقاط تأثيرات الذاكرة والخصائص الوراثية للأنظمة. يركز البحث على معادلة انتشار زمنية كسرية معينة مع تأخير زمني ثابت، تُعبر عنها كالتالي:
\[
\partial_t u(x, t) = \partial_t^{1-\alpha}(p \Delta u(x, t) + au(x, t)) + bu(x, t – \tau) + f(x, t),
\]
حيث يتم تعريف المعلمات والدوال ضمن مجالات محددة.
تستعرض الورقة التقدمات الحديثة في التحليل النظري والعددي لمعادلات الانتشار الكسرية مع التأخيرات، مشيرة إلى طرق ونتائج متنوعة من دراسات سابقة. ومن الجدير بالذكر أن المؤلفين يقدمون نتائج انتظام جديدة لمشكلتهم، مشيرين إلى أن المشتق الأول للحل يظهر انفجارًا فقط عند \(t = 0\)، بينما ينفجر المشتق الثاني عند كل من \(t = 0\) و \(t = \tau\). تشير هذه الانتظامية المحسنة إلى أن المشكلة قيد الدراسة لديها خصائص أفضل من تلك التي تم استكشافها في الأعمال السابقة. يقترح المؤلفون طريقة عددية جديدة تعتمد على تقريب غرينوالد-ليتنكوف ومخطط عنصر نهائي كامل، تهدف إلى تعزيز الاستقرار والتقارب، والتي سيتم مناقشتها في الأقسام التالية من الورقة.
مناقشة
في هذا القسم، يبحث المؤلفون في انتظام الحل للمشكلة المحددة بواسطة المعادلات (1.1)-(1.3) باستخدام طرق الفصل المتغير وتحويل لابلاس. يتم التعبير عن الحل من حيث دالة ميتاج-ليفلي، مما يؤدي إلى حل فريد في المجال \([0, L] \times (0, \tau]\). يؤسس المؤلفون ليلتين رئيسيتين توفران حدودًا على انتظام الحل، موضحين أنه تحت ظروف معينة على البيانات الأولية وعوامل القوة، يبقى الحل \(u_k^{\tau}(x, t)\) محدودًا في معايير مختلفة عبر الفترات المحددة.
تقدم الليمّا الأولى (ليمّا 2.1) علاقة تتضمن دالة ميتاج-ليفلي، بينما تؤسس الليمّا الثانية (ليمّا 2.2) شروطًا يتم بموجبها تحديد أن الحل ومشتقاته محدودة. على وجه التحديد، توضح أنه إذا كانت الشروط الأولية وعوامل القوة تلبي شروط نمو معينة، فإن الحل يظهر انتظامًا في كل من الزمن والمكان. تؤكد النتيجة الرئيسية (نظرية 2.1) أن الحل يحقق حدودًا على مشتقاته، مما يشير إلى أن الحل سلس ويتصرف بشكل جيد عبر الفترات المحددة. تمهد هذه التحليلات الطريق لاستكشاف المزيد من المخططات العددية واستقرارها في الأقسام التالية.
DOI: https://doi.org/10.4208/jcm.2510-m2025-0131
Publication Date: 2026-03-15
Author(s): Weiping Bu et al.
Primary Topic: Fractional Differential Equations Solutions
Overview
In this study, the authors investigate a subdiffusion equation characterized by a constant time delay, $\tau$. They analyze the regularity of the solution, revealing that the first-order time derivative exhibits a singularity at $t = 0^+$, while the second-order time derivative shows singularities at both $t = 0^+$ and $t = \tau^+$. The solution is decomposed into singular and regular components, facilitating a clearer understanding of its behavior.
To address the problem numerically, the authors develop a fully discrete finite element scheme utilizing the standard Galerkin method in spatial dimensions and a Grünwald-Letnikov type approximation for temporal discretization. The stability of this numerical scheme is established, and a new discrete Gronwall inequality is introduced to derive a local error estimate in time. The paper concludes with numerical tests that corroborate the theoretical findings, demonstrating that the proposed method effectively captures the dynamics of the time-fractional diffusion equation with constant time delay.
Introduction
The introduction of the paper highlights the significance of differential equations, particularly delay differential equations, in modeling complex phenomena across various fields such as physics, biology, and medicine. The authors emphasize the advantages of fractional differential equations over integer-order models, particularly in capturing memory effects and hereditary characteristics of systems. The focus of the study is on a specific time-fractional diffusion equation with constant time delay, expressed as:
\[
\partial_t u(x, t) = \partial_t^{1-\alpha}(p \Delta u(x, t) + au(x, t)) + bu(x, t – \tau) + f(x, t),
\]
where the parameters and functions are defined within specified domains.
The paper reviews recent advancements in the theoretical and numerical analysis of fractional diffusion equations with delays, citing various methods and results from previous studies. Notably, the authors present new regularity results for their problem, indicating that the first derivative of the solution exhibits a blow-up only at \(t = 0\), while the second derivative blows up at both \(t = 0\) and \(t = \tau\). This improved regularity suggests that the problem under consideration has better properties than those explored in earlier works. The authors propose a new numerical method based on the Grünwald-Letnikov approximation and a fully discrete finite element scheme, aiming to enhance stability and convergence, which will be discussed in subsequent sections of the paper.
Discussion
In this section, the authors investigate the regularity of the solution to the problem defined by equations (1.1)-(1.3) using the variable separation and Laplace transform methods. The solution is expressed in terms of the Mittag-Leffler function, leading to a unique solution in the domain \([0, L] \times (0, \tau]\). The authors establish two key lemmas that provide bounds on the regularity of the solution, showing that under certain conditions on the initial data and forcing terms, the solution \(u_k^{\tau}(x, t)\) remains bounded in various norms across the specified intervals.
The first lemma (Lemma 2.1) provides a relationship involving the Mittag-Leffler function, while the second lemma (Lemma 2.2) establishes conditions under which the solution and its derivatives are bounded. Specifically, it demonstrates that if the initial conditions and forcing terms satisfy certain growth conditions, then the solution exhibits regularity in both time and space. The main result (Theorem 2.1) asserts that the solution satisfies bounds on its derivatives, indicating that the solution is smooth and well-behaved across the defined intervals. This analysis lays the groundwork for further exploration of numerical schemes and their stability in subsequent sections.