2025-11-22T13:22:16.398648

A time-space B-spline integrator for the Burgers' equation

Dag, UÄurluoÄlu, Adar

The purpose of this paper is to propose a new algorithm for obtaining approximate solutions to the Burgers' equation (BE). Integration in time by a quadratic B-spline collocation method is shown. To the best of our knowledge, B-splines have not previously been used to integrate partial differential equations in both time and space. First, the BE is integrated using quadratic B-spline functions in time, and then the time-integrated BE is further solved in space via the cubic B-spline collocation method. The resulting recursive algebraic equation is used to obtain both shock wave and front propagation solutions of the BE, demonstrating the effectiveness of the space--time B-spline collocation method.

academic

محاكي B-spline في الزمان والمكان لمعادلة Burgers

المعلومات الأساسية

معرّف الورقة: 2510.09408
العنوان: A time-space B-spline integrator for the Burgers' equation
المؤلفون: Idris Dag, Serkan Uğurluoğlu, Nihat Adar (جامعة Eskisehir Osmangazi، تركيا)
التصنيف: math.NA cs.NA
تاريخ النشر: 10 أكتوبر 2025 (مسودة arXiv)
رابط الورقة: https://arxiv.org/abs/2510.09408

الملخص

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

الخلفية البحثية والدافع

خلفية المشكلة

أهمية معادلة Burgers: معادلة Burgers هي معادلة تفاضلية جزئية شبه خطية من النوع القطعي، وتُستخدم على نطاق واسع في العلوم الهندسية لنمذجة الاضطراب والديناميكا الغازية والتوصيل الحراري والمرونة وتدفق المرور والموجات الصوتية وغيرها. وهي توفر صيغة رياضية لوصف التفاعل بين الحمل والانتشار.
ضرورة الحل العددي: عادة ما لا توجد حلول تحليلية للمعادلات التفاضلية الجزئية غير الخطية، خاصة عندما تكون ثابت اللزوجة υ صغيراً، حيث تتقارب الحلول المتسلسلة ببطء وتفتقر إلى الفائدة العملية. تؤدي قيم υ الصغيرة إلى تكوين موجات صدمية، مما يجعل من الصعب على الطرق التحليلية والعددية توفير حلول معقولة.
قيود الطرق الموجودة: تتضمن محاكيات الزمن التقليدية تنسيق Crank-Nicolson وطرق Runge-Kutta والمحاكيات الهندسية والمحاكيات الأسية وغيرها، لكن هذه الطرق تواجه تحديات في الدقة والاستقرار عند التعامل مع مشاكل الموجات الصدمية والانتشار الأمامي لمعادلة Burgers.

الدافع البحثي

الدافع الأساسي لهذه الورقة هو تطوير طريقة عددية جديدة قادرة على:

استخدام B-spline للتكامل في كلا البعدين الزمني والمكاني
التعامل بفعالية مع مشاكل الموجات الصدمية في حالات قيم اللزوجة الصغيرة
توفير حلول عددية بدقة أعلى

المساهمات الأساسية

طريقة رائدة: أول اقتراح لاستخدام B-spline كمحاكي زمني ومكاني معاً لحل معادلة Burgers
استراتيجية B-spline المزدوجة: مزيج مبتكر من استخدام B-spline من الدرجة الثانية في الاتجاه الزمني و B-spline من الدرجة الثالثة في الاتجاه المكاني
محاكاة الموجات الصدمية والانتشار الأمامي: محاكاة ناجحة لحلول الموجات الصدمية وحلول الموجات المتنقلة لمعادلة Burgers
إطار الخوارزمية العودية: إنشاء نظام حل المعادلات الجبرية العودية القائم على تجميع B-spline

شرح الطريقة

تعريف المهمة

حل معادلة Burgers أحادية البعد المعتمدة على الزمن: $u_t + uu_x - \nu u_{xx} = 0, \quad x \in [a,b], \quad t \in (0,T]$

مع الشروط الابتدائية: $u(x,0) = f(x)$

الشروط الحدية: $u(a,t) = u_0, u_x(a,t) = u'_0, u(b,t) = u_m, u_x(b,t) = u'_m$

حيث $\nu$ هو ثابت اللزوجة.

معمارية النموذج

التقسيم الزمني (B-spline من الدرجة الثانية)

تقسيم الشبكة: تقسيم المجال الزمني $[c,d]$ بالتساوي إلى نقاط شبكة $c = t_0 < t_1 < \cdots < t_k = d$ ، حيث $t_j = t_0 + j\Delta t$
دوال أساس B-spline من الدرجة الثانية: تعريف $B_j(t) \in C^1[a,b]$ ، $j = -1,\ldots,n$ :

(t_{j+2}-t)^2 - 3(t_{j+1}-t)^2 + 3(t_j-t)^2, & t_{j-1} \leq t \leq t_j \\ (t_{j+2}-t)^2 - 3(t_{j+1}-t)^2, & t_j \leq t \leq t_{j+1} \\ (t_{j+2}-t)^2, & t_{j+1} \leq t \leq t_{j+2} \\ 0, & \text{غير ذلك} \end{cases}$$ 3. **تمثيل الحل التقريبي**: $$U(x,t) = \sum_{j=-1}^n \delta(x,t_j)B_j(t)$$ #### التقسيم المكاني (B-spline من الدرجة الثالثة) 1. **شبكة المكان**: تقسيم المجال المكاني $[a,b]$ بالتساوي، $x_i = x_0 + ih$، $i = 0,\ldots,m$ 2. **دوال أساس B-spline من الدرجة الثالثة**: $B_i(x) \in C^2[a,b]$، $i = -1,\ldots,m$ 3. **التوسع البارامتري**: $$\delta(x,t^j) = \sum_{i=-1}^{m+1} \sigma_i^j B_i(x)$$ ### نقاط الابتكار التقني 1. **استراتيجية التكامل الطبقية**: تكامل أولاً في الاتجاه الزمني باستخدام B-spline من الدرجة الثانية، ثم حل في الاتجاه المكاني باستخدام B-spline من الدرجة الثالثة 2. **طريقة نقاط التجميع**: فرض تحقيق المعادلة التفاضلية في نقاط الشبكة 3. **المعالجة الخطية**: خطية الحدود غير الخطية من خلال اختيار قيمة الخطوة الزمنية السابقة 4. **التحسين التكراري**: استخدام عملية تكرارية لتحسين الدقة في كل خطوة زمنية ## إعداد التجارب ### مشاكل الاختبار #### المشكلة 1: حل الموجة الصدمية الحل التحليلي: $$u(x,t) = \frac{x/t}{1 + \sqrt{t/t_0}\exp(x^2/(4\nu t))}, \quad t \geq 1, \quad 0 \leq x \leq 1$$ حيث $t_0 = \exp(1/(8\nu))$ #### المشكلة 2: حل الموجة المتنقلة الأمامية الحل التحليلي: $$u(x,t) = \frac{\alpha + \mu + (\mu - \alpha)\exp(\eta)}{1 + \exp(\eta)}$$ حيث $\eta = \frac{\alpha(x - \mu t - \gamma)}{\nu}$، $\alpha = 0.4$، $\mu = 0.6$، $\gamma = 0.125$ ### مؤشرات التقييم - **معيار $L_\infty$**: $|U - U_N|_\infty = \max|U_i^j - (U_N)_i^j|$ - **معيار $L_2$**: $|U - U_N|_2 = \sqrt{\sum_i |U_i^j - (U_N)_i^j|^2}$ - **القيمة القصوى**: القيمة العظمى للموجة الصدمية ### تفاصيل التنفيذ - خطوات الزمان والمكان: $h = \Delta t = 0.01, 0.005, 0.001$ - معاملات اللزوجة: $\nu = 0.01, 0.005$ - نطاق الزمن: $t = 3.25$ (مشكلة الموجة الصدمية)، $t = 1.2$ (مشكلة الموجة الأمامية) ## نتائج التجارب ### النتائج الرئيسية #### نتائج مشكلة الموجة الصدمية | $\nu=0.01$ | $h=\Delta t$ | $L_\infty$ | $L_2$ | القيمة القصوى | |-------------|--------------|------------|-------|------------| | | 0.01 | 0.0498 | 0.0211| 0.2233 | | | 0.005 | 0.0508 | 0.0207| 0.2225 | | | 0.001 | 0.0516 | 0.0206| 0.2203 | | $\nu=0.005$ | $L_\infty$ | $L_2$ | القيمة القصوى | |--------------|------------|-------|------------| | | 0.0582 | 0.0166| 0.2503 | | | 0.0571 | 0.0159| 0.2480 | | | 0.0560 | 0.0155| 0.2467 | #### نتائج مشكلة الموجة الأمامية | $\nu=0.01$ | $h=\Delta t$ | $L_\infty$ | $L_2$ | |-------------|--------------|------------|-------| | | 0.01 | 0.4875 | 0.1138| | | 0.005 | 0.2964 | 0.0638| | | 0.001 | 0.1514 | 0.0354| | $\nu=0.005$ | $L_\infty$ | $L_2$ | |--------------|------------|-------| | | 0.8841 | 0.2532| | | 0.6906 | 0.1329| | | 0.2427 | 0.0367| ### الاكتشافات التجريبية 1. **العلاقة بين الدقة والشبكة**: استخدام خطوات زمنية ومكانية أصغر يحقق حلولاً أكثر سلاسة ودقة أفضل 2. **تأثير اللزوجة**: ثوابت اللزوجة الأصغر تؤدي إلى أشكال موجية أكثر حدة وأخطاء أعلى 3. **توزيع الخطأ**: يحدث الخطأ الأقصى عادة بالقرب من قمة الموجة الصدمية أو منطقة الموجة الأمامية 4. **الاستقرار**: تحافظ الطريقة على استقرار عددي جيد في جميع إعدادات المعاملات المختلفة ## الأعمال ذات الصلة تتضمن الاتجاهات البحثية الرئيسية في هذا المجال: 1. **الطرق العددية التقليدية**: الفروقات المحدودة والعناصر المحدودة والطرق الطيفية وغيرها 2. **طرق B-spline**: تُستخدم بشكل أساسي للتقسيم المكاني، نادراً ما تُستخدم للتكامل الزمني 3. **طرق التجميع**: طرق Galerkin وطرق التجميع والطرق الخالية من الشبكات وغيرها 4. **المحاليل المتخصصة**: طرق عددية خاصة لمعادلة Burgers يكمن الابتكار في هذه الورقة في تطبيق B-spline للتكامل الزمني والمكاني معاً للمرة الأولى، مما يملأ فجوة في هذا المجال. ## الخلاصة والمناقشة ### الاستنتاجات الرئيسية 1. طريقة التجميع الزمكانية B-spline المقترحة قادرة على حل معادلة Burgers بفعالية 2. تُظهر الطريقة دقة واستقراراً جيداً في مشاكل الموجات الصدمية والانتشار الأمامي 3. يمكن التحكم في دقة الحل والسلاسة من خلال تعديل خطوات الزمان والمكان ### القيود 1. **التعقيد الحسابي**: يتطلب حل نظام معادلات خطية في كل خطوة زمنية 2. **معالجة عدم الخطية**: يتطلب عملية خطية وتكرارية 3. **حساسية المعاملات**: حساسة نسبياً لمعاملات اللزوجة والشبكة 4. **نطاق التطبيق**: يركز بشكل أساسي على معادلة Burgers أحادية البعد ### الاتجاهات المستقبلية يقترح المؤلفون أن تركز الأعمال المستقبلية على: - استخدام B-spline من درجات أعلى لتحسين الدقة - التوسع إلى معادلات Burgers متعددة الأبعاد - تحسين كفاءة الخوارزمية ## التقييم المتعمق ### المميزات 1. **ابتكار الطريقة**: أول اقتراح لطريقة التكامل الزمكانية B-spline، يتمتع بالأصالة 2. **الاكتمال النظري**: الاشتقاق الرياضي والتنفيذ العددي كلاهما شامل نسبياً 3. **كفاية التجارب**: اختبار مشاكل نموذجية من نوعين: الموجات الصدمية والانتشار الأمامي 4. **تحليل الدقة**: توفير تحليل خطأ مفصل ودراسة التقارب ### أوجه القصور 1. **نقص التحليل النظري**: افتقار إلى إثبات نظري صارم للتقارب والاستقرار 2. **تجارب المقارنة محدودة**: عدم إجراء مقارنة كافية مع الطرق المتقدمة الأخرى 3. **غياب تحليل الكفاءة الحسابية**: عدم توفير تحليل وقت الحساب واستخدام الذاكرة 4. **نقص التوجيه في اختيار المعاملات**: افتقار إلى التوجيه النظري لاختيار درجة B-spline ### التأثير 1. **المساهمة الأكاديمية**: توفير أفكار جديدة لحل المعادلات التفاضلية الجزئية 2. **القيمة العملية**: الطريقة بسيطة وسهلة التنفيذ، مناسبة للتطبيقات الهندسية 3. **القابلية للتوسع**: يمكن تعميم الإطار على معادلات تفاضلية جزئية غير خطية أخرى ### السيناريوهات المناسبة تناسب هذه الطريقة بشكل خاص: - الحسابات العلمية التي تتطلب حل معادلة Burgers بدقة عالية - محاكاة الظواهر الفيزيائية للموجات الصدمية والانتشار الأمامي - التطبيقات التي تتطلب دقة حسابية عالية لكن متطلبات الكفاءة معتدلة ## المراجع تستشهد الورقة بـ 16 مرجعاً ذا صلة، تغطي الطرق العددية لمعادلة Burgers ونظرية B-spline وطرق التجميع والمجالات ذات الصلة الأخرى، مما يوفر أساساً نظرياً متيناً لهذا البحث.