DOI: https://doi.org/10.3934/era.2026109
تاريخ النشر: 2026-01-01
المؤلف: Man Luo وآخرون
الموضوع الرئيسي: حلول المعادلات التفاضلية الكسرية
نظرة عامة
تقدم هذه الورقة طريقة عددية لحل معادلات التطور الكسرية المميزة بالنوى الضعيفة التفرد. تعتمد الطريقة على استخدام طريقة كرانك-نيكولسون (CN) للتفريق الزمني، مع استخدام خطوات زمنية متدرجة، بينما يتم تحقيق التفريق المكاني من خلال طريقة الفرق المحدود الضمنية في اتجاهات متناوبة (ADI). يقوم المؤلفون بتحليل استقرار الطريقة وتقاربها، مشيرين إلى معدل التقارب \( O(k^2 + h_x^2 + h_y^2) \)، حيث يمثل \( k \) الحد الأقصى للخطوة الزمنية و \( h_x \)، \( h_y \) تشير إلى خطوات الشبكة المكانية الموحدة. تم تقديم ثلاثة أمثلة عددية لتوضيح دقة الطريقة وسلوك التقارب.
في الختام، تطور الدراسة بنجاح إطارًا عدديًا قائمًا على ADI لمعادلات التطور الكسرية ذات النوى الضعيفة التفرد، باستخدام شبكات زمنية متدرجة محددة على أنها \( t_n = (nk’)^\gamma \). تشير النتائج إلى أن الحد الأقصى لترتيب التقارب \( L^2 \) في الزمن يصل إلى 2. ستركز الأبحاث المستقبلية على توسيع هذه المنهجية لتشمل معادلات المشتقات الكسرية في الفضاء، بهدف بناء مخططات مكانية دقيقة من الدرجة العالية على شبكات غير متجانسة.
مقدمة
تركز مقدمة هذه الورقة البحثية على تحليل مشكلة قيمة حدودية أولية ثنائية الأبعاد ممثلة بالمعادلة
\[
u_t + \int_0^t \beta(ts) A u(x, y, s) ds = f(x, y, t), \quad (x, y) \in \Omega, \quad 0 < t \leq T,
\]
مع شروط حدودية \( u(x, y, t) = 0 \) على \( \partial \Omega \) وحالة أولية \( u(x, y, 0) = v(x, y) \) لـ \( (x, y) \in \Omega \). المجال المكاني \( \Omega = (0, 1)^2 \) محدود، و \( A \) هو مشغل خطي ذاتي التوافق، موجب التعريف. يُلاحظ أن النواة \( \beta(t) = t^{\alpha-1} \Gamma(\alpha) \) تتميز بموجبيتها وسلوكها التفردي، حيث \( \alpha \) هو مشغل تكامل كسري مع \( \alpha > 0 \). تسلط الورقة الضوء على الأهمية المتزايدة لمعادلات التفاضل الكسرية (FDEs) عبر مجالات متنوعة، بما في ذلك الاقتصاد، والأرصاد الجوية، وعلم الأحياء، والفيزياء، والكيمياء.
يقترح المؤلفون نهجًا لحل عددي لمعادلات التطور الكسرية ذات النوى الضعيفة التفرد، باستخدام طريقة الفرق المحدود الضمنية في اتجاهات متناوبة (ADI) للتفريق المكاني وطريقة كرانك-نيكولسون (CN) ذات الخطوات الزمنية المتغيرة للتفريق الزمني. يتناولون التحديات التي تطرحها السلوكيات التفردية للمشتقات الكسرية، والتي يمكن أن تؤدي إلى تدهور معدلات التقارب عند تطبيق الطرق العددية القياسية. تؤكد الورقة على استخدام الشبكات المتدرجة لإدارة التفردات الضعيفة بفعالية، مشيرة إلى أن طريقة CN على هذه الشبكات تحقق تقاربًا كاملًا. بالإضافة إلى ذلك، يقدم المؤلفون قاعدة تكامل متخصصة لتفريق مصطلح التطور الكسرية ويحددون استراتيجية لتقليل المشكلة ثنائية الأبعاد إلى سلسلة من المشكلات أحادية البعد، مما يخفف من التعقيد الحسابي. ستفصل الأقسام التالية من الورقة الرموز، وصياغة المخطط المنفصل، وتحليل الاستقرار والتقارب، والتجارب العددية التي تؤكد النتائج النظرية.
نقاش
في هذا القسم، يناقش المؤلفون تطوير مخطط عددي لمعادلات التطور الكسرية ذات النوى الضعيفة التفرد، باستخدام طريقة الاتجاه المتناوب الضمنية (ADI) مع دمجها مع مخطط كرانك-نيكولسون (CN) للتفريق الزمني. يتم تعريف المجال المكاني على أنه \(\Omega = (0, 1) \times (0, 1)\) مع شبكة موحدة، ويستخدم التفريق الزمني شبكات متدرجة من الشكل \(t_n = (nk’)^\gamma\)، حيث \(k’ = T/N\) و \(\gamma \geq 1\). يسمح هذا الاختيار بزيادة تجمع مستويات الزمن بالقرب من \(t = 0\) مع زيادة \(\gamma\)، مما يعالج بفعالية السلوك التفردي في الزمن الأولي.
يشتق المؤلفون المخطط المنفصل من خلال التكامل على فترات زمنية والتفريق المكاني باستخدام نهج ADI. يقدمون عدة ليمات تثبت العلاقات بين دوال الشبكة ومشتقاتها، والتي تعتبر حاسمة للتحليل النظري لاستقرار المخطط وتقاربه. يتم تأكيد استقرار الحلول العددية من خلال النظرية 4.1، التي تحدد حدًا لنور الحل العددي في الخطوة الزمنية النهائية من حيث الشروط الأولية ومصطلحات القوة. تشير تحليل التقارب، المضمن في النظرية 4.3، إلى أن الخطأ ينخفض مع اقتراب معلمات التفريق المكاني والزمني من الصفر، محققًا ترتيب تقارب أقصى من \(L^2\) يبلغ 2 في الزمن. يختتم المؤلفون بالإشارة إلى نيتهم توسيع هذا العمل ليشمل معادلات المشتقات الكسرية في الفضاء مع مخططات مكانية دقيقة من الدرجة العالية على شبكات غير متجانسة.
DOI: https://doi.org/10.3934/era.2026109
Publication Date: 2026-01-01
Author(s): Man Luo et al.
Primary Topic: Fractional Differential Equations Solutions
Overview
This paper presents a numerical method for solving fractional evolution equations characterized by weakly singular kernels. The approach employs the Crank-Nicolson (CN) method for temporal discretization, utilizing graded time steps, while spatial discretization is achieved through an alternating direction implicit (ADI) finite difference method. The authors analyze the stability and convergence of the method, reporting a convergence rate of \( O(k^2 + h_x^2 + h_y^2) \), where \( k \) represents the maximum time step and \( h_x \), \( h_y \) denote uniform spatial grid steps. Three numerical examples are provided to illustrate the method’s accuracy and convergence behavior.
In conclusion, the study successfully develops an ADI-based numerical framework for fractional evolution equations with weakly singular kernels, employing graded temporal meshes defined as \( t_n = (nk’)^\gamma \). The findings indicate that the maximal \( L^2 \) convergence order in time achieves 2. Future research will focus on extending this methodology to space-fractional derivative equations, aiming to construct high-order accurate spatial schemes on nonuniform grids.
Introduction
The introduction of this research paper focuses on the analysis of a two-dimensional initial-boundary value problem represented by the equation
\[
u_t + \int_0^t \beta(ts) A u(x, y, s) ds = f(x, y, t), \quad (x, y) \in \Omega, \quad 0 < t \leq T,
\]
with boundary conditions \( u(x, y, t) = 0 \) on \( \partial \Omega \) and initial condition \( u(x, y, 0) = v(x, y) \) for \( (x, y) \in \Omega \). The spatial domain \( \Omega = (0, 1)^2 \) is bounded, and \( A \) is a self-adjoint, positive definite linear operator. The kernel \( \beta(t) = t^{\alpha-1} \Gamma(\alpha) \) is noted for its positive definiteness and singular behavior, where \( \alpha \) is a fractional integration operator with \( \alpha > 0 \). The paper highlights the increasing relevance of fractional differential equations (FDEs) across various fields, including economics, meteorology, biology, physics, and chemistry.
The authors propose a numerical solution approach for fractional evolution equations with weakly singular kernels, employing an alternating direction implicit (ADI) finite difference method for spatial discretization and a variable time-step Crank-Nicolson (CN) method for temporal discretization. They address the challenges posed by the singular behavior of fractional derivatives, which can lead to deteriorating convergence rates when standard numerical methods are applied. The paper emphasizes the use of graded meshes to effectively manage weak singularities, asserting that the CN method on these meshes achieves full-order convergence. Additionally, the authors introduce a specialized quadrature rule for discretizing the fractional evolution term and outline a strategy to reduce the two-dimensional problem into a series of one-dimensional problems, thereby alleviating computational complexity. The subsequent sections of the paper will detail the notation, formulation of the discrete scheme, stability and convergence analysis, and numerical experiments that validate the theoretical results.
Discussion
In this section, the authors discuss the development of a numerical scheme for fractional evolution equations with weakly singular kernels, utilizing an Alternating Direction Implicit (ADI) method combined with a Crank-Nicolson (CN) scheme for temporal discretization. The spatial domain is defined as \(\Omega = (0, 1) \times (0, 1)\) with a uniform grid, and the temporal discretization employs graded meshes of the form \(t_n = (nk’)^\gamma\), where \(k’ = T/N\) and \(\gamma \geq 1\). This choice allows for increased clustering of time levels near \(t = 0\) as \(\gamma\) increases, effectively addressing singular behavior at the initial time.
The authors derive the discrete scheme through integration over time intervals and spatial discretization using the ADI approach. They present several lemmas that establish relationships between grid functions and their derivatives, which are crucial for the theoretical analysis of the scheme’s stability and convergence. The stability of the numerical solutions is confirmed through Theorem 4.1, which bounds the norm of the numerical solution at the final time step in terms of initial conditions and forcing terms. The convergence analysis, encapsulated in Theorem 4.3, indicates that the error decreases as the spatial and temporal discretization parameters approach zero, achieving a maximal L2 convergence order of 2 in time. The authors conclude by noting their intention to extend this work to space-fractional derivative equations with high-order accurate spatial schemes on nonuniform grids.