مخطط عددي للدراسة الحاسوبية للانتشار ثنائي الأبعاد وأنظمة بورجر مع تقدير الاستقرار والخطأ
Numerical Scheme for the Computational Study of Two Dimensional Diffusion and Burgers’ Systems with Stability and Error Estimate

المجلة: Journal of Nonlinear Mathematical Physics، المجلد: 32، العدد: 1
DOI: https://doi.org/10.1007/s44198-025-00277-6
تاريخ النشر: 2025-04-22
المؤلف: Muhammad Bilal وآخرون
الموضوع الرئيسي: المعادلات التفاضلية والأساليب العددية

نظرة عامة

تقدم هذه الورقة نهجًا عدديًا لحل معادلات تفاضلية جزئية (PDEs) ثنائية الأبعاد مفردة ومترابطة باستخدام نسخة جديدة من موجات هار، وتحديدًا موجات هار من الدرجة الثالثة (S3HW)، بالتزامن مع طرق الفروق المحدودة. يتم هيكلة المنهجية في مرحلتين: تتضمن المرحلة الأولى التقدير العددي للمشتقات الزمنية من خلال الفروق المحدودة، مما يحول المشكلة إلى صيغة زمنية منفصلة. تركز المرحلة الثانية على تقريب المشتقات المكانية والحلول باستخدام S3HW، تليها تطبيق تقنية التوافق لاشتقاق نظام من المعادلات الجبرية الخطية. تؤدي حل هذه النظام إلى الحصول على معاملات الموجات غير المعروفة، والتي تُستخدم بعد ذلك لحساب الحلول العددية. تتضمن الورقة أيضًا تحليلًا للخطأ، والتقارب، والاستقرار، مما يؤدي إلى تقدير جديد للخطأ، وتحقق من الطريقة من خلال محاكاة عددية تُظهر فعاليتها مقارنة بالتقنيات الحالية، بما في ذلك موجات هار من الدرجة الثانية ودوال الأساس الشعاعية.

في الختام، أثبت تنفيذ S3HW وصيغة الفروق المحدودة نجاحه في حل PDEs ثنائية الأبعاد الخطية وغير الخطية عدديًا. تشير النتائج إلى أن الطريقة المقترحة تتفوق على التقنيات الحالية، لا سيما من حيث معايير الخطأ والاستقرار، بينما تعالج بفعالية القضايا المتعلقة بكثافة مصفوفة المعاملات. يُقترح العمل المستقبلي لتوسيع هذه الطريقة لتشمل مشاكل من درجات أعلى وأبعاد متعددة، بما في ذلك تلك التي تتضمن مشتقات كسرية، مع مراعاة تحسين التعقيد الحسابي من خلال اختيار حجم خطوة مناسب.

مقدمة

تسلط مقدمة هذه الورقة البحثية الضوء على أهمية المعادلات التفاضلية الجزئية (PDEs) في نمذجة ظواهر فيزيائية متنوعة، مثل توصيل الحرارة، وديناميات الموجات، والاضطراب. تركز بشكل خاص على معادلة الانتشار ومعادلة بورجر، اللتين تعتبران محوريتين في مجالات مثل ميكانيكا السوائل والعمليات الكيميائية. توضح الورقة الصيغ الرياضية لهذه المعادلات، مع التأكيد على أدوارها كمعايير لاختبار طرق عددية جديدة.

يستعرض المؤلفون التقنيات العددية الحالية لحل نماذج الانتشار ومعادلات بورجر ثنائية الأبعاد، بما في ذلك الفروق المحدودة، والعناصر المحدودة، والطرق المعتمدة على الموجات. بشكل ملحوظ، يناقشون مزايا استخدام موجات هار، لا سيما موجات هار من الدرجة الثالثة التي تم تقديمها حديثًا، والتي أظهرت معدلات تقارب محسنة مقارنة بنظيراتها من الدرجة الثانية. الهدف الرئيسي من هذه الدراسة هو توسيع تطبيق موجات هار من الدرجة الثالثة لمشاكل عالية الأبعاد، تحديدًا لمعادلة الانتشار ومعادلات بورجر. ستفصل الأقسام التالية من الورقة منهجية الموجات، وتحليل الخطأ، وتقديم النتائج العددية، مما يؤدي إلى استنتاج يلخص النتائج.

نقاش

يتناول قسم النقاش في الورقة تنفيذ طريقة موجة هار من الدرجة الثالثة (HW) لحل معادلات الانتشار ثنائية الأبعاد ومعادلات بورجر. تبدأ الطريقة بتقسيم الفترة \([A, B)\) إلى \(3Q\) فترات فرعية متساوية، وتعريف المعلمات للترجمة والتوسيع، وإنشاء دوال الموجات. يتم تعريف دالة التدرج \(H_1(x)\) والموجات المتماثلة وغير المتماثلة، مما يؤدي إلى صياغة تكاملات متكررة ضرورية للحل العددي. يتم اشتقاق التعبيرات التحليلية لهذه التكاملات، والتي تعتبر حاسمة لتقريب حل المعادلات المعنية.

ثم تُطبق الطريقة المقترحة على مشاكل عددية متنوعة، مما يظهر فعاليتها من خلال معايير الخطأ مثل \(L^\infty\)، \(L^2\)، الجذر التربيعي لمتوسط المربعات (RMS)، والخطأ النسبي (RE). تشير النتائج إلى أن دقة الطريقة تتحسن مع زيادة مستويات الدقة، متفوقة على التقنيات الحالية مثل HW من الدرجة الثانية ودوال الأساس الشعاعية. كما يتم تحليل استقرار الطريقة وتقاربها، مما يؤكد أن طريقة HW من الدرجة الثالثة تتقارب بمعدل \(O(1/Q^3)\) مع زيادة الدقة. تشير النتائج إلى أنه بينما تظهر الطريقة وعدًا للمعادلات التفاضلية الجزئية الخطية وغير الخطية، يمكن أن يركز العمل المستقبلي على توسيع تطبيقها لمشاكل من درجات أعلى وأبعاد متعددة، مع معالجة التعقيد الحسابي من خلال اختيار حجم خطوة استراتيجي.

Journal: Journal of Nonlinear Mathematical Physics, Volume: 32, Issue: 1
DOI: https://doi.org/10.1007/s44198-025-00277-6
Publication Date: 2025-04-22
Author(s): Muhammad Bilal et al.
Primary Topic: Differential Equations and Numerical Methods

Overview

This paper presents a numerical approach for solving two-dimensional single and coupled partial differential equations (PDEs) using a novel version of Haar wavelets, specifically the scale-3 Haar wavelets (S3HW), in conjunction with finite difference methods. The methodology is structured in two phases: the first phase involves the numerical estimation of temporal derivatives through finite differences, converting the problem into a time-discrete format. The second phase focuses on approximating spatial derivatives and solutions using S3HW, followed by the application of a collocation technique to derive a system of linear algebraic equations. The solution of this system yields the unknown wavelet coefficients, which are then used to compute numerical solutions. The paper also includes an analysis of error, convergence, and stability, leading to a new error estimate, and validates the method through numerical simulations that demonstrate its effectiveness compared to existing techniques, including scale-2 Haar wavelets and radial basis functions.

In conclusion, the implementation of the S3HW and finite difference formulation has proven successful for numerically solving two-dimensional linear and nonlinear PDEs. The results indicate that the proposed method outperforms existing techniques, particularly in terms of error norms and stability, while effectively addressing issues related to the sparsity of the coefficient matrix. Future work is suggested to extend this method to higher-order and multi-dimensional problems, including those involving fractional derivatives, with considerations for optimizing computational complexity through appropriate step size selection.

Introduction

The introduction of this research paper highlights the significance of partial differential equations (PDEs) in modeling various physical phenomena, such as heat conduction, wave dynamics, and turbulence. Specifically, it focuses on the diffusion equation and Burgers’ equation, which are pivotal in fields like fluid mechanics and chemical processes. The paper outlines the mathematical formulations of these equations, emphasizing their roles as benchmarks for testing new numerical methods.

The authors review existing numerical techniques for solving two-dimensional diffusion and Burgers’ models, including finite difference, finite element, and wavelet-based methods. Notably, they discuss the advantages of using Haar wavelets, particularly the newly introduced scale-3 Haar wavelets, which have shown improved convergence rates compared to their scale-2 counterparts. The primary aim of this study is to extend the application of scale-3 Haar wavelets to high-dimensional problems, specifically for the diffusion equation and Burgers’ equations. The subsequent sections of the paper will detail the wavelet methodology, error analysis, and present numerical results, culminating in a conclusion that summarizes the findings.

Discussion

The discussion section of the paper elaborates on the implementation of a scale-3 Haar wavelet (HW) method for solving two-dimensional diffusion equations and Burgers’ equations. The method begins by partitioning the interval \([A, B)\) into \(3Q\) equal subintervals, defining parameters for translation and dilation, and establishing wavelet functions. The scaling function \(H_1(x)\) and the symmetric and anti-symmetric wavelets are defined, leading to the formulation of repeated integrals necessary for the numerical solution. The analytical expressions for these integrals are derived, which are crucial for approximating the solution of the equations under consideration.

The proposed method is then applied to various numerical problems, demonstrating its effectiveness through error norms such as \(L^\infty\), \(L^2\), root mean square (RMS), and relative error (RE). The results indicate that the accuracy of the method improves with increasing resolution levels, outperforming existing techniques like scale-2 HW and radial basis functions. The stability and convergence of the method are also analyzed, confirming that the scale-3 HW method converges at a rate of \(O(1/Q^3)\) as the resolution increases. The findings suggest that while the method shows promise for linear and nonlinear PDEs, future work could focus on extending its application to higher-order and multi-dimensional problems, addressing computational complexity through strategic step size selection.