DOI: https://doi.org/10.1186/s13104-024-07005-1
PMID: https://pubmed.ncbi.nlm.nih.gov/39833977
تاريخ النشر: 2025-01-20
المؤلف: Zerihun Ibrahim Hassen وآخرون
الموضوع الرئيسي: المعادلات التفاضلية والأساليب العددية
نظرة عامة
في هذه الدراسة، يتناول المؤلفون مشكلة تأخير التوصيل-الانتشار البارابوليكية المعتمدة على الزمن والمضطربة بشكل فردي، والتي تخضع لظروف حدود ديريشليه، والتي تتميز بوجود طبقات حدودية تظهر تدرجات حادة أو تذبذبات. تكافح الطرق العددية التقليدية لتقديم حلول دقيقة في مثل هذه السيناريوهات. للتغلب على هذه التحديات، يقترح المؤلفون طريقة عددية تعتمد على الانحدار الموحد بدون تذبذبات. تستخدم هذه الطريقة طريقة أويلر الضمنية في الاتجاه الزمني وطريقة الانحدار الأسّي مع شبكة موحدة في الاتجاه المكاني، مع دمج عامل ملاءمة أسّي لإدارة معامل الاضطراب بشكل فعال.
يتم تحليل استقرار وتقديرات الخطأ الموحد للبرنامج العددي المقترح بدقة، مما يظهر أنه يحقق تقاربًا خطيًا موحدًا في المعيار الأقصى، بغض النظر عن شبكة المعطيات ومعامل الاضطراب. يتم تسليط الضوء على كفاءة الطريقة الحاسوبية، حيث إنها مناسبة للتنفيذ على أجهزة الكمبيوتر الشخصية مع الحد الأدنى من وقت وحدة المعالجة المركزية للنقاط الشبكية المطلوبة. ومع ذلك، تظل تقييمات الدوال الأسية مكثفة حسابيًا، خاصة للدقة العالية. يتم دعم النتائج النظرية من خلال أمثلة عددية، والتي توضح دقة الطريقة المتفوقة مقارنة بالأساليب الموجودة في الأدبيات.
مقدمة
في هذه المقدمة، يناقش المؤلفون المعادلات التفاضلية ذات التأخير المضطرب بشكل فردي (SPD-DEs)، والتي تتميز باعتمادها على الحالات والمشتقات السابقة، وهي مهمة في تطبيقات متنوعة عبر العلوم والهندسة، بما في ذلك أنظمة التحكم، والتفاعلات الكيميائية، وعلم الأوبئة. تركز الورقة بشكل خاص على مشاكل تأخير التوصيل-الانتشار البارابوليكية المضطربة بشكل فردي (SPDPCDPs)، التي تمثل أنظمة حيث يكون التفاعل بين التوصيل، والانتشار، وتأثيرات التأخير حاسمًا، مثل النقل التفاعلي في الوسائط المسامية. يبرز المؤلفون التحديات في حل SPDPCDPs بسبب وجود طبقات حدودية ناتجة عن معاملات الاضطراب، مما يؤدي إلى تباينات كبيرة في الحلول التي تكافح الطرق العددية الكلاسيكية للتعامل معها بشكل فعال.
لمعالجة هذه التحديات، يستعرض المؤلفون الطرق العددية الحالية، بما في ذلك تقنيات الشبكة الملائمة والمشغل الملائم، التي تهدف إلى الاستقرار والتقارب الموحد بغض النظر عن قيمة معامل الاضطراب. يلخصون مختلف الأساليب التي تم تطويرها سابقًا من قبل باحثين آخرين، مشيرين إلى نجاحاتهم في تحقيق التقارب الموحد للمعامل. ثم يقدم المؤلفون مخططهم العددي الجديد لحل SPDPCDPs مع تأخيرات كبيرة وظروف حدود ديريشليه، والذي يستخدم طريقة أويلر الضمنية للتفريق الزمني ومخطط انحدار أسّي للتفريق المكاني. يتميز هذا النهج الجديد باستقلاله عن الشبكات المصممة خصيصًا والمعرفة السابقة بخصائص طبقة الحدود، مما يعد بتحسين الدقة والاتساق في النتائج.
النتائج
في هذا القسم، يقدم المؤلفون أمثلة عددية ونتائج تتعلق بمشكلتين محددتين تتضمنان معادلات تفاضلية بارابولية مضطربة بشكل فردي تتميز بتأخيرات كبيرة. الهدف هو إظهار فعالية الطريقة المقترحة في معالجة هذه المعادلات المعقدة. تشير النتائج إلى أن الطريقة تلتقط بنجاح ديناميات الأنظمة قيد الدراسة، مما يوفر رؤى حول سلوك الحلول في وجود تأخيرات كبيرة. لا تقتصر هذه الاستكشافات على التحقق من قابلية تطبيق الطريقة فحسب، بل تبرز أيضًا إمكانياتها للبحث المستقبلي في المجالات ذات الصلة.
المناقشة
في هذا القسم، يناقش المؤلفون صياغة وتحليل طريقة عددية لحل مشاكل تأخير التوصيل-الانتشار البارابوليكية المضطربة بشكل فردي (SPDPCDPs) مع ظروف حدود ديريشليه. يتم تعريف المشكلة على مجالات مكانية وزمنية، مع تضمين المشغل التفاضلي لعنصر تأخير ومعاملات اضطراب. يثبت المؤلفون وجود وحيدة الحلول من خلال تطبيق مبادئ قصوى وليمات، مما يظهر أن الحلول تظهر طبقة حدودية مع اقتراب معامل الاضطراب $\epsilon$ من الصفر.
تستخدم الخطة العددية طريقة أويلر الضمنية لتفريق الزمن وتستخدم انحدارات أسّيّة لتفريق الفضاء. يستنتج المؤلفون عدة ليمات ونظريات لضمان استقرار وموثوقية وتقارب الطريقة، موضحين أن الحل يبقى محدودًا ومتقاربًا بشكل موحد بغض النظر عن معلمات الشبكة ومعامل الاضطراب. يتم التحقق من النتائج من خلال أمثلة عددية، والتي توضح فعالية الطريقة في التقاط سلوك طبقة الحدود وتحقيق حلول دقيقة. تشير المقارنات مع الطرق الموجودة إلى أن النهج المقترح يقدم دقة وكفاءة متفوقة، خاصة مع زيادة حجم الشبكة وانخفاض معامل الاضطراب.
DOI: https://doi.org/10.1186/s13104-024-07005-1
PMID: https://pubmed.ncbi.nlm.nih.gov/39833977
Publication Date: 2025-01-20
Author(s): Zerihun Ibrahim Hassen et al.
Primary Topic: Differential Equations and Numerical Methods
Overview
In this study, the authors address a singularly perturbed time-dependent delay parabolic convection-diffusion problem subject to Dirichlet boundary conditions, which is characterized by the presence of boundary layers exhibiting steep gradients or oscillations. Traditional numerical methods struggle to provide accurate solutions in such scenarios. To overcome this challenge, the authors propose an oscillation-free parameter uniform exponentially spline numerical method. This method employs the implicit Euler method in the temporal direction and an exponential spline method with a uniform mesh in the spatial direction, incorporating an exponential fitting factor to effectively manage the perturbation parameter.
The stability and uniform error estimates of the proposed numerical scheme are rigorously analyzed, demonstrating that it achieves uniform linear convergence in the maximum norm, independent of mesh and perturbation parameters. The method’s computational efficiency is highlighted, as it is suitable for implementation on personal computers with minimal CPU time for the required mesh points. However, the evaluation of exponential functions remains computationally intensive, particularly for high precision. The theoretical findings are substantiated through numerical examples, which illustrate the method’s superior accuracy compared to existing approaches in the literature.
Introduction
In this introduction, the authors discuss singularly perturbed delay differential equations (SPD-DEs), which are characterized by their dependence on past states and derivatives, and are significant in various applications across science and engineering, including control systems, chemical reactions, and epidemiology. The paper specifically focuses on singularly perturbed delay parabolic convection-diffusion problems (SPDPCDPs), which model systems where the interplay between convection, diffusion, and delay effects is critical, such as in reactive transport in porous media. The authors highlight the challenges in solving SPDPCDPs due to the presence of boundary layers caused by perturbation parameters, which lead to significant variations in solutions that classical numerical methods struggle to handle effectively.
To address these challenges, the authors review existing numerical methods, including fitted mesh and fitted operator techniques, which aim for stability and uniform convergence regardless of the perturbation parameter’s value. They summarize various approaches previously developed by other researchers, noting their successes in achieving parameter uniform convergence. The authors then introduce their own novel numerical scheme for solving SPDPCDPs with large delays and Dirichlet boundary conditions, which utilizes an implicit Euler method for temporal discretization and an exponential spline scheme for spatial discretization. This new approach is distinguished by its independence from specially designed meshes and prior knowledge of boundary layer characteristics, promising improved accuracy and consistency in results.
Results
In this section, the authors present numerical examples and results pertaining to two specific problems involving singularly perturbed parabolic differential equations characterized by significant delays. The aim is to demonstrate the effectiveness of the proposed method in addressing these complex equations. The findings indicate that the method successfully captures the dynamics of the systems under consideration, providing insights into the behavior of solutions in the presence of large delays. This exploration not only validates the method’s applicability but also highlights its potential for further research in related areas.
Discussion
In this section, the authors discuss the formulation and analysis of a numerical method for solving singularly perturbed delay parabolic convection-diffusion problems (SPDPCDPs) with Dirichlet boundary conditions. The problem is defined over spatial and temporal domains, with the differential operator incorporating a delay term and perturbation parameters. The authors establish the existence and uniqueness of solutions through the application of maximum principles and lemmas, demonstrating that the solutions exhibit a boundary layer as the perturbation parameter $\epsilon$ approaches zero.
The numerical scheme employs an implicit Euler method for time discretization and utilizes exponential splines for spatial discretization. The authors derive several lemmas and theorems to ensure the stability, consistency, and convergence of the method, showing that the solution remains bounded and uniformly convergent regardless of the mesh parameters and perturbation parameter. The results are validated through numerical examples, which illustrate the method’s effectiveness in capturing boundary layer behavior and achieving accurate solutions. Comparisons with existing methods indicate that the proposed approach offers superior accuracy and efficiency, particularly as the mesh size increases and the perturbation parameter decreases.