DOI: https://doi.org/10.1016/j.cnsns.2025.108824
تاريخ النشر: 2025-04-07
المؤلف: Jun Ma وآخرون
الموضوع الرئيسي: حلول المعادلات التفاضلية الكسرية
نظرة عامة
في هذا البحث، تم تطوير مخطط خطي كامل التمييز لمعالجة معادلة شرودنجر الزمنية غير الخطية ذات البعدين المميزة بحلول ضعيفة التفرد. يستخدم المخطط طريقة L1 لمشتق كابوتو الكسري، وصيغة عكسية لتقريب الحد غير الخطي، ومخطط فرق من خمس نقاط للتفريق المكاني. يثبت المؤلفون بشكل صارم استقرار المخطط غير المشروط والتقارب في النقاط الزمنية، ولا سيما دون فرض أي قيود على نسبة الشبكة.
تشمل النتائج تقدير خطأ محلي لمخطط L1، مما يؤكد أن المخطط الكامل التمييز يحقق معدل تقارب قدره \( O(\tau_t^{\alpha-1} n + h^2) \). تؤكد هذه النتيجة فعالية الطريقة في التعامل مع الحلول غير السلسة مع الحفاظ على الاستقرار والتقارب، مما يوفر إطارًا قويًا لحل المعادلات التفاضلية الكسري المعقدة. تؤكد التجارب العددية النتائج النظرية، مما يبرز دقة المخطط وقابليته للتطبيق العملي.
مقدمة
تناقش مقدمة الورقة أهمية معادلة شرودنجر الزمنية غير الخطية (NTFS) في مجالات فيزيائية متنوعة، لا سيما في ميكانيكا الكم والديناميات الكسري. يتم تقديم المعادلة على النحو التالي
\[
iD^\alpha_t u + \Delta u + f(|u|^2)u = 0, \quad (x, y) \in \Omega, \quad 0 < t \leq T,
\]
مع شروط حدودية وبدائية محددة. هنا، \(D^\alpha_t\) يدل على مشتق كابوتو الكسري، وعندما تقترب \(\alpha\) من 1، تتقارب المعادلة إلى معادلة شرودنجر غير الخطية الكلاسيكية. تسلط المقدمة الضوء على التقدمات الأخيرة في الطرق العددية لحل NTFS، بما في ذلك المخططات المختلفة التي تم تطويرها لمعالجة التحديات التي تطرحها التفردات الأولية في الحلول.
يشير المؤلفون إلى أن الأبحاث السابقة افترضت في الغالب السلاسة في الحلول، وهو ما يكون غالبًا مقيدًا جدًا. يهدفون إلى معالجة هذه الفجوة من خلال إنشاء تقديرات خطأ محلية لـ NTFS دون فرض شروط صارمة على نسبة الشبكة، مما يستوعب الحلول غير السلسة. تتضمن الطريقة المقترحة مخطط فرق L1 خطي وتستهدف اشتقاق تقديرات خطأ مثالية تحت افتراضات انتظام أضعف للدالة \(f\). يتم توضيح هيكل الورقة، مما يشير إلى أن الأقسام اللاحقة ستفصل المخطط العددي، وتحليل التقارب، والاختبارات العددية للتحقق من النتائج النظرية.
مناقشة
في هذا القسم، يناقش المؤلفون تطوير وتحليل طريقة خطية كاملة التمييز لتقريب مشتق كابوتو الكسري لحل معادلة شرودنجر الزمنية غير الخطية ذات البعدين. تستخدم الطريقة شبكة محددة بواسطة حجم خطوة الزمن $\tau = T/N$ وحجم شبكة الفضاء $h = L/M$، مع تقريب العقد المعبر عنه بـ $U^n_{j,k}$. تقدم الورقة مخطط L1 لتقريب المشتق الكسري، وتؤسس عدة ليمات تظهر أحادية الاتجاه وحدود سلسلة المضاعفات $\theta_n$. كما يفترض المؤلفون شروط انتظام على الحل الدقيق، والتي تدعمها الأعمال السابقة حول المعادلات الزمنية الكسري.
يتم تفصيل تحليل الاستقرار والتقارب للمخطط المقترح بشكل صارم. يستخرج المؤلفون عدم المساواة المنفصلة لغروينوال ويستخدمونها لإظهار أن المخطط مستقر بشكل غير مشروط بالنسبة للقيمة الابتدائية. يثبتون وجود خطأ تقريبي ويقدمون تقدير خطأ محلي، مما يظهر أن معدل التقارب هو $O(\tau t^{\alpha-1} n + h^2)$، مستقل عن حجم شبكة الفضاء. تؤكد التجارب العددية هذه النتائج النظرية، مما يؤكد معدلات التقارب المتوقعة لقيم مختلفة من $\alpha$ ونسب شبكة مختلفة، مما يعزز قوة وفعالية الطريقة المقترحة لحل PDEs الزمنية غير الخطية الكسري.
DOI: https://doi.org/10.1016/j.cnsns.2025.108824
Publication Date: 2025-04-07
Author(s): Jun Ma et al.
Primary Topic: Fractional Differential Equations Solutions
Overview
In this research, a linearized fully discrete scheme is developed to address the two-dimensional nonlinear time fractional Schrödinger equation characterized by weakly singular solutions. The scheme employs the L1 method for the Caputo fractional derivative, a backward formula for approximating the nonlinear term, and a five-point difference scheme for spatial discretization. The authors rigorously demonstrate the unconditional stability and pointwise-in-time convergence of the proposed scheme, notably without imposing any restrictions on the grid ratio.
The findings include a local error estimate for the L1 scheme, confirming that the fully discrete scheme achieves a convergence rate of \( O(\tau_t^{\alpha-1} n + h^2) \). This result underscores the effectiveness of the method in handling non-smooth solutions while maintaining stability and convergence, thereby providing a robust framework for solving complex fractional differential equations. Numerical experiments corroborate the theoretical results, highlighting the scheme’s accuracy and practical applicability.
Introduction
The introduction of the paper discusses the significance of the nonlinear time fractional Schrödinger equation (NTFS) in various physical fields, particularly in quantum mechanics and fractional dynamics. The equation is presented as
\[
iD^\alpha_t u + \Delta u + f(|u|^2)u = 0, \quad (x, y) \in \Omega, \quad 0 < t \leq T,
\]
with specified boundary and initial conditions. Here, \(D^\alpha_t\) denotes the Caputo fractional derivative, and as \(\alpha\) approaches 1, the equation converges to the classical nonlinear Schrödinger equation. The introduction highlights recent advancements in numerical methods for solving NTFS, including various schemes that have been developed to address the challenges posed by initial singularities in solutions.
The authors note that previous research has primarily assumed smoothness in solutions, which is often too restrictive. They aim to address this gap by establishing local error estimates for NTFS without imposing stringent conditions on the grid ratio, thereby accommodating non-smooth solutions. The proposed approach involves a linearized L1 finite difference scheme and aims to derive optimal error estimates under weaker regularity assumptions for the function \(f\). The structure of the paper is outlined, indicating that subsequent sections will detail the numerical scheme, convergence analysis, and numerical tests to validate the theoretical findings.
Discussion
In this section, the authors discuss the development and analysis of a linearized fully discrete method for approximating the Caputo fractional derivative of a solution to a two-dimensional nonlinear time-fractional Schrödinger equation. The method employs a grid defined by time step size $\tau = T/N$ and spatial mesh size $h = L/M$, with nodal approximations denoted as $U^n_{j,k}$. The paper introduces the L1 scheme for approximating the fractional derivative, establishing several lemmas that demonstrate the monotonicity and boundedness of the sequence of multipliers $\theta_n$. The authors also assume regularity conditions on the exact solution, which are supported by previous works on time-fractional equations.
The stability and convergence analysis of the proposed scheme are rigorously detailed. The authors derive a discrete Grönwall inequality and utilize it to show that the scheme is unconditionally stable with respect to the initial value. They establish a truncation error and provide a local error estimate, demonstrating that the convergence rate is $O(\tau t^{\alpha-1} n + h^2)$, independent of the spatial mesh size. Numerical experiments validate these theoretical findings, confirming the expected convergence rates for various values of $\alpha$ and different grid ratios, thus reinforcing the robustness and efficacy of the proposed method for solving nonlinear time-fractional PDEs.