DOI: https://doi.org/10.3934/nhm.2025058
تاريخ النشر: 2025-01-01
المؤلف: Chuan Ran وآخرون
الموضوع الرئيسي: حلول المعادلات التفاضلية الكسرية
نظرة عامة
تقدم هذه الدراسة إطار عمل للاختلافات المحدودة/غاليركين غير المستمر المحلي (LDG) مصمم لحل معادلة انتشار-موجة كسرية من الرتبة الرابعة. يتم تقريب المشغل الكسر الزمني باستخدام طريقة الاختلافات المحدودة، مما يحقق دقة تقليم من \( O((\Delta t)^{3-\alpha}) \)، حيث \( \Delta t \) هو زيادة الوقت و \( \alpha \) هو الرتبة الكسرية. تستخدم التقطيع المكاني تقنية LDG، مما يؤدي إلى صياغة منفصلة بالكامل ضمنية. يثبت المؤلفون الاستقرار غير المشروط والتقارب للخوارزمية المقترحة من خلال الاستقراء الرياضي، ويحققون صحة حدود الخطأ النظرية من خلال سلسلة من التجارب الحسابية.
في الختام، تقدم البحث طريقة للاختلافات المحدودة/LDG من الرتبة العالية لنظام انتشار-موجة كسرية من الرتبة الرابعة، خاصة للمشتقات الكسرية مع \( 1 < \alpha < 2 \). يستنتج المؤلفون صيغة التقطيع ويثبتون بدقة الاستقرار غير المشروط والتقارب للنظام المنفصل بالكامل. تظهر الأمثلة العددية أداءً ممتازًا ودقة للطريقة. تهدف الأعمال المستقبلية إلى توسيع نظام LDG لعناصر مكانية من رتبة أعلى، مثل \( P_3 \)، وتطوير أنظمة LDG منفصلة بالكامل لنماذج ثنائية أو ثلاثية الأبعاد. تعتبر الصياغة المحلية لإطار عمل LDG مفيدة بشكل خاص للتطبيقات عالية الأبعاد، وقد يؤدي استكشاف طرق بدون شبكة، مثل نهج غاليركين الخالي من العناصر، إلى تعزيز الدقة وقابلية تطبيق الطريقة المقترحة.
مقدمة
تسلط مقدمة هذه الورقة البحثية الضوء على الاهتمام المتزايد في حساب التفاضل الكسرية، وخاصة تطبيقاته عبر مجالات مختلفة مثل المالية والهندسة وعلم الأحياء. تركز على معادلات انتشار كسرية زمنية (T-FDEs)، والتي هي مجموعة فرعية من المعادلات الخطية التكاملية التفاضلية. بينما توجد بعض الحلول التحليلية، لا تزال العديد من المعادلات التفاضلية الجزئية الكسرية (FPDEs) معقدة وغير قابلة للحل، مما يستدعي تطوير طرق عددية فعالة للتطبيقات العملية. تؤكد الورقة على أهمية الأساليب العددية، بما في ذلك الاختلافات المحدودة، والعناصر المحدودة، والأساليب الطيفية، من بين أمور أخرى.
تتناول الدراسة بشكل خاص معادلة انتشار-موجة كسرية من الرتبة الرابعة المعطاة بواسطة
\[
\frac{\partial^\alpha u(x, t)}{\partial t^\alpha} + \mu^2 \frac{\partial^4 u(x, t)}{\partial x^4} = f(x, t), \quad x \in (a, b), \, t \in (0, T], \, u(x, 0) = u_0(x), \, u_t(x, 0) = u_1(x),
\]
حيث \(1 < \alpha < 2\) تشير إلى الرتبة الكسرية في الزمن، ويتم استخدام المشتق الكسرية كابوتو. تشير المقدمة أيضًا إلى طرق عددية مختلفة تم اقتراحها سابقًا لحل معادلات مشابهة، بما في ذلك تلك التي قدمها هو وزانغ، وجعفري وآخرون، والتي استكشفت طرقًا مختلفة مثل مخططات كرانك-نيكولسون وطرق بدون شبكة. توضح الورقة هيكلها، مشيرة إلى أن الأقسام التالية ستتناول الأسس النظرية، وتطوير طريقة للاختلافات المحدودة/LDG، والتجارب العددية، والاستنتاجات.
النتائج
في هذا القسم، يحدد المؤلفون الرموز الأساسية والنتائج المساعدة المتعلقة بتحليلهم لنطاق حسابي ممثل بشبكة تعرف على أنها \( a = x_{1/2} < x_{3/2} < \ldots < x_{N+1/2} = b \). يتم الإشارة إلى خلايا الشبكة بواسطة \( I_j = [x_{j-1/2}, x_{j+1/2}] \) لـ \( j = 1, \ldots, N \)، مع أحجام خلايا متناسبة \( \Delta x_j = x_{j+1/2} - x_{j-1/2} \). يتم تعريف الفضاء متعدد الحدود من الدرجة القصوى \( k \) على أنه \( V^k_h = \{ v : v|_{I_j} \in P^k(I_j), x \in I_j, j = 1, \ldots, N \} \). يقدم المؤلفون نوعين من الإسقاطات، الإسقاط القياسي \( L^2 \) \( P \) والإسقاطات الخاصة \( P^+ \) و \( P^- \)، والتي تسهل تحليل الخطأ في النطاق المكاني أحادي البعد \( \Omega = [a, b] \). يتميز الإسقاطات بخصائص محددة تضمن التعامد بالنسبة لوظائف الاختبار متعددة الحدود. بشكل ملحوظ، يتم إعطاء خاصية التقريب لهذه الإسقاطات بواسطة \[ \| \omega_e \| + h \| \omega_e \|_\infty + h^{1/2} \| \omega_e \|_{\tau h} \leq C h^{k+1}, \] حيث تمثل \( \omega_e \) مصطلحات الخطأ المرتبطة بالإسقاطات، و \( C \) هو ثابت إيجابي عام. يختتم القسم بتعريفات للمنتج الداخلي القياسي والمعايير المرتبطة في \( L^2(D) \)، مع تبسيط الرموز عندما يكون \( D = \Omega \).
المناقشة
في هذا القسم، يوضح المؤلفون التقطيع العددي لمشتق زمني كسرية كما هو موضح في المعادلة (1.1)، مع التركيز على خصائص الاستقرار والتقارب لنظامهم المقترح. يستخدمون صيغة اختلاف خلفية ثلاثية النقاط لتقريب المشتق الزمني، مما يؤدي إلى خطأ تقليم من الرتبة $O((\Delta t)^2)$. يتم تقريب المشتق الكسرية كابوتو باستخدام سلسلة من النقاط الزمنية المنفصلة، مما ينتج عنه صياغة تتضمن معاملات $b_i$ التي تكون إيجابية وتتناقص إلى الصفر مع زيادة $n$. يستنتج المؤلفون نظامًا منفصلًا بالكامل بناءً على تمثيل نظام من الرتبة الأولى، مما يضمن اختيار التدفقات العددية للحفاظ على الاستقرار.
يتم إثبات استقرار النظام تحت ظروف حدود دورية أو مدعومة بشكل مضغوط، مما يظهر الاستقرار غير المشروط من خلال الاستقراء الرياضي. يثبت المؤلفون أن الحل العددي يبقى محدودًا من حيث الشروط الأولية والمشتق الزمني الأول. علاوة على ذلك، يقدمون تقديرات للخطأ تشير إلى معدلات التقارب المعتمدة على الرتبة الكسرية $\alpha$، مع حدود محددة لـ $1 < \alpha < 2$ وعندما تقترب $\alpha$ من 2. تؤكد الأمثلة العددية المقدمة النتائج النظرية، مما يظهر أن الطريقة تحقق معدلات التقارب المتوقعة وتحافظ على دقة عالية عبر مختلف الرتب الكسرية. يختتم المؤلفون باقتراح أعمال مستقبلية لتوسيع الطريقة لتشمل مشاكل ذات أبعاد أعلى وأطر عددية بديلة.
DOI: https://doi.org/10.3934/nhm.2025058
Publication Date: 2025-01-01
Author(s): Chuan Ran et al.
Primary Topic: Fractional Differential Equations Solutions
Overview
This study introduces a finite difference/local discontinuous Galerkin (LDG) framework designed to solve a fourth-order fractional diffusion-wave equation. The temporal fractional operator is approximated using a finite difference method, achieving a truncation accuracy of \( O((\Delta t)^{3-\alpha}) \), where \( \Delta t \) is the time increment and \( \alpha \) is the fractional order. The spatial discretization employs the LDG technique, resulting in a fully implicit discrete formulation. The authors establish the unconditional stability and convergence of the proposed algorithm through mathematical induction, and they validate the theoretical error bounds with a series of computational experiments.
In conclusion, the research presents a high-order finite difference/LDG method for a fourth-order fractional diffusion-wave system, specifically for fractional derivatives with \( 1 < \alpha < 2 \). The authors derive a discretization formula and rigorously prove the unconditional stability and convergence of the fully discrete scheme. Numerical examples demonstrate the method's excellent performance and accuracy. Future work aims to extend the LDG scheme to higher-order spatial elements, such as \( P_3 \), and to develop fully discrete LDG schemes for two- or three-dimensional models. The local formulation of the LDG framework is particularly advantageous for high-dimensional applications, and exploring meshless methods, like the element-free Galerkin approach, may further enhance the accuracy and applicability of the proposed method.
Introduction
The introduction of this research paper highlights the growing interest in fractional calculus, particularly its applications across various fields such as finance, engineering, and biology. It focuses on time-fractional diffusion equations (T-FDEs), which are a subset of linear integro-differential equations. While some analytical solutions exist, many fractional partial differential equations (FPDEs) remain complex and intractable, necessitating the development of efficient numerical methods for practical applications. The paper emphasizes the importance of numerical approaches, including finite difference, finite element, and spectral methods, among others.
The study specifically addresses the fourth-order fractional diffusion-wave equation given by
\[
\frac{\partial^\alpha u(x, t)}{\partial t^\alpha} + \mu^2 \frac{\partial^4 u(x, t)}{\partial x^4} = f(x, t), \quad x \in (a, b), \, t \in (0, T], \, u(x, 0) = u_0(x), \, u_t(x, 0) = u_1(x),
\]
where \(1 < \alpha < 2\) denotes the fractional order in time, and the Caputo fractional derivative is utilized. The introduction also references various numerical methods previously proposed for solving similar equations, including those by Hu and Zhang, Jafari et al., and others, which have explored different approaches such as Crank-Nicolson schemes and meshfree methods. The paper outlines its structure, indicating that subsequent sections will cover theoretical preliminaries, the development of a finite difference/LDG method, numerical experiments, and conclusions.
Results
In this section, the authors establish foundational notations and auxiliary results pertinent to their analysis of a computational domain represented by a mesh defined as \( a = x_{1/2} < x_{3/2} < \ldots < x_{N+1/2} = b \). The cells of the mesh are denoted by \( I_j = [x_{j-1/2}, x_{j+1/2}] \) for \( j = 1, \ldots, N \), with corresponding cell sizes \( \Delta x_j = x_{j+1/2} - x_{j-1/2} \). The piecewise-polynomial space of degree at most \( k \) is defined as \( V^k_h = \{ v : v|_{I_j} \in P^k(I_j), x \in I_j, j = 1, \ldots, N \} \). The authors introduce two types of projections, the standard \( L^2 \) projection \( P \) and special projections \( P^+ \) and \( P^- \), which facilitate error analysis in the one-dimensional spatial domain \( \Omega = [a, b] \). The projections are characterized by specific properties that ensure orthogonality with respect to polynomial test functions. Notably, the approximation property for these projections is given by \[ \| \omega_e \| + h \| \omega_e \|_\infty + h^{1/2} \| \omega_e \|_{\tau h} \leq C h^{k+1}, \] where \( \omega_e \) represents the error terms associated with the projections, and \( C \) is a generic positive constant. The section concludes with definitions of the scalar inner product and associated norms in \( L^2(D) \), with the notation simplified when \( D = \Omega \).
Discussion
In this section, the authors detail the numerical discretization of a time-fractional derivative as described in Eq (1.1), focusing on the stability and convergence properties of their proposed scheme. They employ a three-point backward difference formula to approximate the time derivative, leading to a truncation error of order $O((\Delta t)^2)$. The Caputo fractional derivative is approximated using a series of discrete time points, resulting in a formulation that incorporates coefficients $b_i$ which are positive and decrease to zero as $n$ increases. The authors derive a fully discrete scheme based on a first-order system representation, ensuring that the numerical fluxes are chosen to maintain stability.
The stability of the scheme is established under periodic or compactly supported boundary conditions, demonstrating unconditional stability through mathematical induction. The authors prove that the numerical solution remains bounded in terms of the initial conditions and the first time derivative. Furthermore, they provide error estimates that indicate convergence rates dependent on the fractional order $\alpha$, with specific bounds for $1 < \alpha < 2$ and as $\alpha$ approaches 2. The numerical examples presented validate the theoretical findings, showing that the method achieves the expected convergence rates and maintains high accuracy across various fractional orders. The authors conclude by suggesting future work to extend the method to higher-dimensional problems and alternative numerical frameworks.