DOI: https://doi.org/10.1038/s41598-026-53831-7
PMID: https://pubmed.ncbi.nlm.nih.gov/42174034
تاريخ النشر: 2026-05-22
المؤلف: Chaturya Karanam وآخرون
الموضوع الرئيسي: المعادلات التفاضلية والأساليب العددية
نظرة عامة
تتناول هذه الورقة البحثية مشاكل نقاط التحول المضطربة بشكل فردي (SPTPPs) مع شروط حدود من نوع روبن، والتي تتميز بطبقات حادة على الحدود والداخلية التي تعقد الطرق العددية القياسية. لمواجهة هذه التحديات، يقترح المؤلفون طريقة تجميعية باستخدام B-spline مكعب مثلثي (TCBS)، والتي تستفيد من سلاسة ومرونة دوال الأساس TCBS لالتقاط التغيرات السريعة المرتبطة بسلوك نقطة التحول بشكل فعال. الإطار النظري الذي تم تأسيسه في الدراسة يوضح أن الطريقة تحقق تقارب موحد بالنسبة للمعامل من الدرجة \( O(N^{-2} (\ln N)^3) \)، مما يشير إلى دقة من الدرجة الثانية مع عامل لوغاريتمي.
يؤكد المؤلفون نتائجهم النظرية من خلال تجارب عددية أجريت على شبكة شيشكين، التي تسمح بشبكة مصقولة في مناطق الطبقات، مما يعزز دقة النهج القائم على B-spline. تظهر النتائج أن الطريقة تحافظ على معدلات التقارب قريبة من اثنين عبر قيم مختلفة لمعامل الاضطراب، مما يؤكد تقديرات الخطأ المتوقعة. بالإضافة إلى ذلك، يتم تسليط الضوء على الكفاءة الحسابية لطريقة تجميع TCBS من خلال الأوقات المبلغ عنها لوحدة المعالجة المركزية، وتؤكد المقارنات مع النتائج الموجودة في الأدبيات موثوقية وفعالية النهج المقترح. بشكل عام، تقدم الدراسة مخططًا عدديًا قويًا لمعالجة SPTPPs، مما يمهد الطريق للبحوث المستقبلية في هذا المجال.
مقدمة
تناقش المقدمة المشاكل المضطربة بشكل فردي (SPPs)، التي تتميز بمعامل إيجابي صغير يؤثر على المشتق من أعلى درجة، مما يؤدي إلى سلوك متعدد المقاييس في الحلول. غالبًا ما تظهر هذه المشاكل طبقات حادة على الحدود أو الداخلية، مما يعقد الحلول العددية، خاصة مع الطرق الكلاسيكية التي تعتمد على الشبكات الموحدة. تنشأ مثل هذه التحديات في تطبيقات متنوعة، بما في ذلك ديناميات السوائل ونمذجة أشباه الموصلات. وبالتالي، تم تطوير تقنيات عددية متقدمة، مثل طرق المشغل الملائمة والشبكات المعدلة للطبقات، لالتقاط هذه الظواهر بدقة مع الحفاظ على الاستقرار والدقة الموحدة بالنسبة للمعامل.
تحدث مشاكل نقاط التحول (TPPs)، وهي فئة صعبة من SPPs، عندما يتلاشى معامل الحد المرافق أو يتغير الإشارة داخل المجال، مما يؤدي إلى تفاعلات معقدة بين الطبقات الحدودية والداخلية. تسلط المقدمة الضوء على التقدمات الأخيرة في الطرق العددية لمشاكل TPPs، بما في ذلك طرق الفرق المحدودة وطرق العناصر المحدودة، التي تهدف إلى تعزيز الدقة والكفاءة الحسابية. على الرغم من التقدم، لا يزال معالجة الأنظمة المترابطة من SPTPPs مع شروط حدود روبن غير مستكشفة بشكل كافٍ، مما يتطلب تطوير تقنيات عددية قوية تضمن الاستقرار والدقة عبر المجال بأكمله. تختتم هذه الفقرة بتحديد النظام المحدد من SPPs من الدرجة الثانية مع نقطة تحول داخلية التي ستتناولها الورقة، مع التأكيد على الحاجة إلى استراتيجيات عددية فعالة لحل الطبقات الحدودية المرتبطة.
النتائج
في هذا القسم، يؤسس المؤلفون مبدأ الحد الأقصى ونتائج الاستقرار المرتبطة لحل مشكلة محددة، مقسمين المجال إلى ثلاثة فترات فرعية: $\Omega_1 = [-1, -\tau]$, $\Omega_2 = [-\tau, \tau]$, و$\Omega_3 = [\tau, 1]$، حيث $0 < \tau \geq \frac{1}{2}$. يفترضون أن دالة المعامل $|c_r(s)|$ محصورة بشكل موحد بعيدًا عن الصفر، تحديدًا $|c_r(s)| \geq \gamma > 0$ لـ $\tau < |s| \leq 1$. تعتبر الشروط على معاملات الربط والحدود القطرية حاسمة لضمان التزام النظام بمبدأ الحد الأقصى. تؤكد النظرية 1.1 أنه إذا تم استيفاء شروط حدود معينة، فإن الحل $Z(s)$ يبقى غير سالب عبر المجال. لتحقيق نتائج نظرية، يتم إجراء تجارب عددية باستخدام MATLAB، مع التركيز على مجموعة من مشاكل الاختبار. يتم تقييم دقة الطريقة العددية المقترحة من خلال مبدأ الشبكة المزدوجة، مقارنة الحلول المحسوبة على الشبكات الخشنة والمصقولة. يتم تعريف الحد الأقصى للخطأ النقطي، وتقدير معدلات التقارب، مما يكشف أنه مع تحسين الشبكة من $N = 64$ إلى $N = 1024$، تنخفض الأخطاء النقطية القصوى بشكل كبير، مما يظهر استقرار وتقارب الطريقة. تشير النتائج إلى أن الطريقة تحقق تقريبًا من الدرجة الثانية، مع انخفاض الأخطاء باستمرار عبر مختلف معاملات الاضطراب، مما يؤكد قوة وسلوك موحد بالنسبة للمعامل للنهج العددي. تظهر المقارنات مع الطرق الموجودة تحسينًا في الدقة، مدعومة بشكل أكبر بالتمثيلات الرسومية للحلول العددية وتوزيعات الأخطاء.
المناقشة
في هذا القسم، يناقش المؤلفون التقدم في الطرق العددية لحل المشاكل المضطربة بشكل فردي (SPPs)، مع التركيز بشكل خاص على أنظمة المعادلات التفاضلية الجزئية الزمنية المضطربة بشكل فردي (SPTPPs) مع شروط حدود روبن. أثبتت الأعمال السابقة، بما في ذلك تلك التي قام بها جيثا وتاميلسيلفان، طرق الفرق المحدودة (FDMs) مع درجات متفاوتة من الدقة، لكن غالبًا ما تكافح هذه الطرق مع الطبقات الحدودية وتتطلب شبكات دقيقة لتحقيق نتائج مرضية. يقترح المؤلفون طريقة تجميعية باستخدام B-spline مكعب مثلثي (TCBS) على شبكة شيشكين، والتي تم تصميمها للتعامل بشكل فعال مع التدرجات الحادة والطبقات الحدودية، محققة تقارب موحد بالنسبة للمعامل من الدرجة \(O(N^{-2} (\ln N)^3)\). يمثل هذا تحسينًا كبيرًا مقارنة بالتقارب من الدرجة الأولى الذي يُرى عادةً في FDMs القياسية.
تستخدم طريقة TCBS دوال أساس مثلثية سلسة، والتي تناسب بشكل أفضل تقريب الحلول المتغيرة بسرعة مقارنةً بطرق B-spline متعددة الحدود الكلاسيكية. يوضح المؤلفون بناء شبكة شيشكين، التي تقسم المجال الحسابي إلى مناطق تسمح بدقة مصقولة في المناطق التي تحدث فيها التدرجات الحادة. تم صياغة المخطط العددي لضمان أن النظام المنفصل الناتج يكون نادرًا ومنظمًا بشكل جيد، مما يعزز الكفاءة الحسابية والاستقرار. تختتم الفقرة بمخطط هيكل الورقة، مشيرة إلى أن الأقسام التالية ستقدم تحليل تقارب صارم، ونتائج عددية، ومناقشات حول آثار الطريقة المقترحة.
DOI: https://doi.org/10.1038/s41598-026-53831-7
PMID: https://pubmed.ncbi.nlm.nih.gov/42174034
Publication Date: 2026-05-22
Author(s): Chaturya Karanam et al.
Primary Topic: Differential Equations and Numerical Methods
Overview
This research paper addresses singularly perturbed turning point problems (SPTPPs) with Robin-type boundary conditions, which are characterized by sharp boundary and interior layers that complicate standard numerical methods. To tackle these challenges, the authors propose a trigonometric cubic B-spline (TCBS) collocation method, which leverages the smoothness and flexibility of TCBS basis functions to effectively capture the rapid variations associated with turning point behavior. The theoretical framework established in the study demonstrates that the method achieves parameter-uniform convergence of order \( O(N^{-2} (\ln N)^3) \), indicating second-order accuracy with a logarithmic factor.
The authors validate their theoretical findings through numerical experiments conducted on a Shishkin mesh, which allows for refined mesh in layer regions, enhancing the accuracy of the spline-based approach. The results show that the method maintains convergence rates close to two across various perturbation parameter values, confirming the predicted error estimates. Additionally, the computational efficiency of the TCBS collocation method is highlighted through reported CPU times, and comparisons with existing results from the literature further affirm the reliability and effectiveness of the proposed approach. Overall, the study presents a robust numerical scheme for addressing SPTPPs, paving the way for future research in this area.
Introduction
The introduction discusses singularly perturbed problems (SPPs), characterized by a small positive parameter affecting the highest-order derivative, leading to multiscale behavior in solutions. These problems often exhibit sharp boundary or interior layers, complicating numerical resolution, especially with classical methods that rely on uniform meshes. Such challenges arise in various applications, including fluid dynamics and semiconductor modeling. Consequently, advanced numerical techniques, such as fitted operator methods and layer-adapted meshes, have been developed to accurately capture these phenomena while maintaining stability and parameter-uniform accuracy.
Turning point problems (TPPs), a challenging subclass of SPPs, occur when the convective term’s coefficient vanishes or changes sign within the domain, resulting in complex interactions between boundary and interior layers. The introduction highlights recent advancements in numerical methods for TPPs, including finite difference and finite element methods, which aim to enhance accuracy and computational efficiency. Despite progress, the treatment of coupled systems of SPTPPs with Robin boundary conditions remains underexplored, necessitating the development of robust numerical techniques that ensure stability and accuracy across the entire domain. The section concludes by outlining the specific system of second-order SPPs with an interior turning point that the paper will address, emphasizing the need for effective numerical strategies to resolve the associated boundary layers.
Results
In this section, the authors establish the maximum principle and associated stability results for the solution of a specified problem, partitioning the domain into three subintervals: $\Omega_1 = [-1, -\tau]$, $\Omega_2 = [-\tau, \tau]$, and $\Omega_3 = [\tau, 1]$, where $0 < \tau \geq \frac{1}{2}$. They assume that the coefficient function $|c_r(s)|$ is uniformly bounded away from zero, specifically $|c_r(s)| \geq \gamma > 0$ for $\tau < |s| \leq 1$. The conditions on the coupling coefficients and diagonal terms are crucial for ensuring the system adheres to the maximum principle. Theorem 1.1 asserts that if certain boundary conditions are satisfied, then the solution $Z(s)$ remains non-negative across the domain. To validate the theoretical findings, numerical experiments are conducted using MATLAB, focusing on a set of test problems. The proposed numerical method's accuracy is assessed through the double mesh principle, comparing solutions computed on coarse and refined meshes. The maximum pointwise error is defined, and convergence rates are estimated, revealing that as the mesh is refined from $N = 64$ to $N = 1024$, the maximum pointwise errors decrease significantly, demonstrating stability and convergence of the method. The results indicate that the method achieves nearly second-order convergence, with errors consistently decreasing across various perturbation parameters, confirming the robustness and parameter-uniform behavior of the numerical approach. Comparisons with existing methods show improved accuracy, further supported by graphical representations of the numerical solutions and error distributions.
Discussion
In this section, the authors discuss advancements in numerical methods for solving singularly perturbed problems (SPPs), particularly focusing on systems of singularly perturbed time-dependent partial differential equations (SPTPPs) with Robin boundary conditions. Previous works, including those by Geetha and Tamilselvan, have established finite difference methods (FDMs) with varying degrees of accuracy, but often these methods struggle with boundary layers and require fine meshes for satisfactory results. The authors propose a Trigonometric Cubic B-spline (TCBS) collocation method on a Shishkin mesh, which is designed to effectively handle sharp gradients and boundary layers, achieving a parameter-uniform convergence of order \(O(N^{-2} (\ln N)^3)\). This represents a significant improvement over the first-order convergence typically seen in standard FDMs.
The TCBS method utilizes smooth trigonometric basis functions, which are better suited for approximating rapidly varying solutions compared to classical polynomial spline methods. The authors detail the construction of the Shishkin mesh, which partitions the computational domain into regions that allow for refined resolution in areas where steep gradients occur. The numerical scheme is formulated to ensure that the resulting discrete system is sparse and well-structured, enhancing computational efficiency and stability. The section concludes with an outline of the paper’s structure, indicating that subsequent sections will provide a rigorous convergence analysis, numerical results, and discussions on the implications of the proposed method.