إنشاء حلول انفجار من النوع I لمعادلة مكافئية شبه خطية من رتبة أعلى

الملخص

ندرس المعادلة المكافئية شبه الخطية من رتبة أعلى

tu=(Δ)mu+u|u|p1,

في الفضاء كله N، حيث p>1 وm1 عدد صحيح فردي. نعرض حلولا غير ذاتية التشابه من النوع I تنفجر لهذه المعادلة، ونحصل على وصف حاد لسلوكها التقاربي. وتعتمد طريقة الإنشاء على التحليل الطيفي للمؤثر الخطي الناتج عن الخطية، وهو غير ذاتي المرافقة، في إطار مناسب من المتغيرات المحجّمة. وبالنظر إلى الخواص الطيفية والقطاعية المعروفة للمؤثر الخطي الناتج عن الخطية التي حصل عليها Galaktionov [Galaktionov(2001)]، نعيد تناول التقنية التي طورها Merle-Zaag [Merle and Zaag(1997)] في الحالة الكلاسيكية m=1، والتي تتألف من خطوتين: اختزال المسألة إلى مسألة منتهية الأبعاد، ثم حل المسألة منتهية الأبعاد بحجة طوبولوجية كلاسيكية قائمة على نظرية الدرجة. ويقدم تحليلنا تبريرا صارما لنتيجة صورية في [Galaktionov(2001)].

الكلمات المفتاحية:
معادلة مكافئية من رتبة أعلى، حل انفجار، ملف الانفجار، الاستقرار
1991 Mathematics Subject Classification:
أولي: 35K50, 35B40; ثانوي: 35K55, 35K57.
May 13, 2026

Tej-Eddine Ghoul, Van Tien Nguyen and Hatem Zaag

New York University in Abu Dhabi, P.O. Box 129188, Abu Dhabi, United Arab Emirates.

Université Paris 13, Sorbonne Paris Cité, LAGA, CNRS (UMR 7539), F-93430, Villetaneuse, France.

1. المقدمة.

نهتم بالمعادلة المكافئية شبه الخطية

{tu=𝐀mu+u|u|p1,u(0)=u0𝐀m(Δ)m, (1.1)

حيث u(t):N مع N1، ويرمز Δ إلى مؤثر لابلاس القياسي في N، أما الأسّان p وm فثابتان،

p>1andm,m1odd.

إن المعادلة المكافئية شبه الخطية من رتبة أعلى (1.1) هي تعميم طبيعي لمعادلة الحرارة شبه الخطية الكلاسيكية (m=1). وتظهر في كثير من التطبيقات الفيزيائية مثل نظرية الأغشية الرقيقة، والتشحيم، والحمل-الانفجار، وانتقال الطور، أو التطبيقات في الميكانيكا الإنشائية (انظر كتاب Petetier-Troy [Peletier and Troy(2001)] والمراجع الواردة فيه).

بواسطة النتائج القياسية يمكن حل مسألة كوشي المحلية للمعادلة (1.1) في L1L بفضل التمثيل التكاملي

u(t)=𝒦m(t)u0+0t𝒦m(ts)u(s)|u(s)|p1𝑑s, (1.2)

حيث إن 𝒦m(t) هو الحل الأساسي للمعادلة المكافئية الخطية t𝒦m=𝐀m𝒦m، والمعرّف عن طريق تحويل فورييه العكسي

𝒦m(x,t)=1(2π)NNe|ξ|2mtı(ξx)𝑑ξ,𝒦m(x,0)=δ(x).

من Fujita [Fujita(1966)] (m=1) وGalaktionov-Pohozaev [Galaktionov and Pohozaev(2002)] (m>1)، نعلم أن

pF1+2mN,

هو أس فوجيتا الحرج للمسألة بالمعنى الآتي. إذا كان p>pF، فإنه من أجل أي معطى ابتدائي صغير بما يكفي u0L1(N)L(N)، تقبل مسألة كوشي (1.1) حلا عالميا يحقق u(t)0 عندما t+ بانتظام في N. وإذا كان p(1,pF] وكانت المعطيات الابتدائية u00,Nu0𝑑x0 من أجل m>1 أو u00 من أجل m=1، فإن الحل الموافق للمسألة (1.1) ينفجر في زمن منته ما T>0، أي إن

limtTu(t)L(N)=+.

هنا يسمى T زمن الانفجار، وتسمى النقطة aN نقطة انفجار إذا وفقط إذا وجدت متتالية (an,tn)(a,T) بحيث |u(an,tn)|+ عندما n+. ويسمى حل (1.1) انفجارا من النوع I إذا حقق

c(Tt)1p1u(t)LC(Tt)1p1, (1.3)

وإلا فهو انفجار من النوع II. إضافة إلى ذلك، نسمي حل انفجار ذاتي التشابه إذا كان من الشكل

u(x,t)=(Tt)1p1Φ(y),y=x(Tt)12m, (1.4)

حيث إن Φ ليس ثابتا على نحو مطابق. ومن الواضح أن حل الانفجار ذاتي التشابه هو من النوع I.

عندما يكون m=1، تختزل المسألة (1.1) إلى معادلة الحرارة شبه الخطية الكلاسيكية

tu=Δu+u|u|p1, (1.5)

التي دُرست على نطاق واسع في العقود الأربعة الأخيرة، ولا يمكن لأي مراجعة أن تكون مستوعبة. وبالنظر إلى اهتمامنا بإنشاء حلول ذات سلوك انفجار موصوف، نذكر فقط الأعمال السابقة في هذا الاتجاه. فقد أعطى Bricmont-Kupiainen [Bricmont and Kupiainen(1994)] أول نتيجة بنائية، إذ أثبت وجود حل انفجار من النوع I للمعادلة (1.5) وفق الدينامية التقاربية

supxN|(Tt)1p1u(x,t)κ(1+(p1)4p|x|2(Tt)|log(Tt)|)1p1|0astT, (1.6)

لبعض ثابت كوني موجب κ=κ(p). نلاحظ أن مؤلفي [Bricmont and Kupiainen(1994)] عرضوا أيضا حلولا تنفجر في زمن منته وتحقق سلوكيات تقاربية أخرى يُتوقع أن تكون غير مستقرة. ونلاحظ أيضا أن Bressan [Bressan(1990), Bressan(1992)] أنجز بناء مماثلا في حالة لاخطية أسية. لاحقا، اقترح Merle-Zaag [Merle and Zaag(1997)] تعديلا لحجة [Bricmont and Kupiainen(1994)] وحصل على استقرار الحل المنشأ الذي يحقق (1.6) تحت اضطرابات صغيرة للمعطيات الابتدائية. وقد لوحظ استقرار السلوك التقاربي (1.6) عدديا بواسطة Berger-Kohn [Berger and Kohn(1988)] (انظر أيضا Nguyen [Nguyen(2017)] لتحليلات عددية أخرى). وبوجه خاص، أثبت Herrero-Velázquez [Herrero and Velázquez(1992)] أن دينامية الانفجار (1.6) نموذجية في الحالة أحادية البعد، وأعلنا الأمر نفسه للحالة ذات الأبعاد الأعلى (لكنهما لم ينشراه قط).

تعتمد طريقة [Bricmont and Kupiainen(1994)] و[Merle and Zaag(1997)] على فهم الخاصية الطيفية للمؤثر الخطي الناتج عن الخطية حول ملف متوقع في إطار متغيرات التشابه. وبصورة تقريبية، يملك المؤثر الخطي عددا منتهيا من القيم الذاتية الموجبة، وقيمة ذاتية صفرية، وطيفا سالبا؛ ثم يسيرون في خطوتين:

  • اختزال مسألة لا نهائية الأبعاد إلى مسألة منتهية الأبعاد، بمعنى أن ضبط الخطأ يختزل إلى ضبط المركبات الموافقة للقيم الذاتية الموجبة.

  • حل المسألة منتهية الأبعاد بفضل حجة طوبولوجية كلاسيكية قائمة على نظرية الدرجة.

وقد مُدّد هذا الإجراء العام ذو الخطوتين إلى أوضاع متعددة، مثل حالة معادلة Ginzgburg-Landau المركبة بواسطة Masmoudi-Zaag [Masmoudi and Zaag(2008)] وNouaili-Zaag [Nouaili and Zaag(2018)] (انظر أيضا Zaag [Zaag(1998)] لعمل سابق)؛ ومعادلة الحرارة شبه الخطية المركبة التي لا تملك بنية تغيرية بواسطة Duong [Duong()] وNouaili-Zaag [Nouaili and Zaag(2015)]؛ ومعادلات الحرارة شبه الخطية غير الثابتة بالتحجيم بواسطة Ebde-Zaag [Ebde and Zaag(2011)] وNguyen-Zaag [Nguyen and Zaag(2016)] وDuong-Nguyen-Zaag [Duong et al.(2018)Duong, Nguyen, and Zaag]. ونذكر أيضا عمل Tayachi-Zaag [Tayachi and Zaag(2015), Tayachi and Zaag(2016)] وGhoul-Nguyen-Zaag [Ghoul et al.(2017)Ghoul, Nguyen, and Zaag] الذي يتناول معادلة حرارة لاخطية ذات مصدر مزدوج يعتمد على الحل وعلى تدرجه في إطار حرج ما. في [Ghoul et al.(2018a)Ghoul, Nguyen, and Zaag, Ghoul et al.(2018b)Ghoul, Nguyen, and Zaag]، وفقنا في تكييف الطريقة لإنشاء حل انفجار مستقر لنظام مكافئي شبه خطي غير تغيري.

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

مبرهنة 1.1 (حلول انفجار من النوع I للمعادلة (1.1) ذات سلوك موصوف).

ليكن m1 عددا صحيحا فرديا. توجد معطيات ابتدائية u0L(N) بحيث إن الحل الموافق u(t) للمعادلة (1.1) يحقق
(i) ينفجر u(t) في زمن منته T=T(u0) عند الأصل فقط.
(ii) (السلوك التقاربي)

supxN|(Tt)1p1u(x,t)Φ(x[(Tt)|log(Tt)|]12m)|C|log(Tt)|12m, (1.7)

حيث

Φ(ξ)=κ(1+Bm,p|ξ|2m)1p1, (1.8)

مع

κ=(p1)1p1,Bm,p=(1)m+12(p1)(2m)!p((4m)!2[(2m)!]2). (1.9)

(iii) (ملف الانفجار النهائي) يوجد u𝒞(N\{0},) بحيث إن u(x,t)u(x) عندما tT بانتظام على المجموعات المدمجة من N\{0}، حيث

u(x)κ(Bm,p|x|2m2m|log|x||)1p1as|x|0. (1.10)
ملاحظة 1.2.

نعتقد أن ملف انفجار كهذا (1.8) موجود لكل m. ونلاحظ أن الثابت Bm,p<0 عندما يكون m زوجيا (انظر (1.9))، ومن ثم فإن الملف Φ ينفجر على الواجهة المنتهية |ξ|=ξ=(Bm,p)12m. وهذا يعني أن الحالة التي يكون فيها m زوجيا ستؤدي إلى حلول انفجار من النوع II بالمعنى الوارد في (1.3). وعلى الرغم من أن الأفكار الرئيسة للتبرير الكامل لمثل هذا السلوك الانفجاري عندما يكون m زوجيا تبقى هي نفسها، فإن البرهان سيكون دقيقا للغاية وسنعالجه في عمل مستقل.

ملاحظة 1.3.

حل الانفجار الموصوف في المبرهنة 1.1 ليس ذاتي التشابه بالمعنى الوارد في (1.4). نلاحظ أنه، خلافا لحلول الانفجار في معادلة الحرارة شبه الخطية الكلاسيكية من الرتبة الثانية (1.5)، خمّن Budd-Galaktionov-Williams [Budd et al.(2004)Budd, Galaktionov, and Williams]، من خلال حسابات عددية وتقاربية، وجود ما لا يقل عن 2m2 حلول انفجار ذاتية التشابه غير تافهة للمعادلة (1.1)، وأن الملفات ذات القيمة العظمى الوحيدة تقابل حلولا ذاتية التشابه مستقرة (نموذجية) للانفجار.

ملاحظة 1.4.

يتضمن برهان المبرهنة 1.1 (في البعد N=1 تبسيطا) وصفا مفصلا لمجموعة المعطيات الابتدائية التي تؤدي إلى الدينامية التقاربية (1.7). وبوجه خاص، تكون معطياتنا الابتدائية تقريبا من الشكل (انظر الصيغة (3.13) أدناه)

u0(x)=es0[Φ(xes0/2m)+As02χ(2xes0/2m,s0)k=12m1dkψk(xes0/2m)],

حيث إن A وs0 ثابتان كبيران مثبتان، وΦ هو الملف المعرف بواسطة (1.8)، وχ دالة قطع ملساء ما، و(d0,,d2m1)2m معاملات حرة، وψk,0k2m1 هي الدوال الذاتية للمؤثر الخطي (انظر القضية 2.1 لتعريف دقيق) الموافقة للقيمة الذاتية الموجبة λk=1k2m. ومن خلال حجة طوبولوجية، نثبت أنه توجد صيغة ملائمة للمعاملات (d0,,d2m1) بحيث إن حل المعادلة (1.1) ذي المعطى الابتدائي u0 يحقق خلاصة المبرهنة 1.1. وبمعنى ما، فإن الحل الذي أنشأناه مستقر بترافق 2(m1) بالمعنى الآتي. إن المركبات 2m من الحل الخطي الموافقة لـ λk=1k2m لها نمو أسي eλks. غير أن النمطين الأولين الموافقين لـ λ0 وλ1 يمكن إزالتهما باستعمال ثبات المسألة تحت ترجمتي الزمن والفضاء. ومن ثم، بتثبيت الاتجاهات 2(m1) ψ2,,ψ2m1 واضطراب المركبات الباقية (في L)، ما زلنا نحصل على الدينامية التقاربية نفسها (1.7) للحل المضطرب. ويتطلب برهان الاستقرار بترافق 2(m1) بعض انتظام Lipschitz لمجموعة المعطيات الابتدائية المعتبرة، وسيعالج ذلك على حدة في عمل آخر.

ملاحظة 1.5.

وفقا لإنشائنا، تقع الدينامية التقاربية (1.7) على المتشعب المركزي المولد بالدوال الذاتية الموافقة للقيمة الذاتية الصفرية λ2m=0. ويمكن توسيع تحليلنا لإنشاء حل للمعادلة (1.1) ينفجر في زمن منته ويمتلك دينامية تقاربية مختلفة عن (1.7). وتمتلك مثل هذه الحلول على وجه الخصوص ديناميات تقاربية واقعة على المتشعب المستقر المولد بالدوال الذاتية الموافقة للقيمة الذاتية السالبة λk=1k2m<0 مع k2m+1. وكما سيُشرح في الملاحظة 1.4، فإن المعطيات الابتدائية الموافقة التي تؤدي إلى مثل هذه الحلول ستتضمن k معاملات مع k2m+1 (اعتبر N=1)، بحيث تُسند حجة طوبولوجية لضبط المركبات k الأولى الموافقة للمتجهات الذاتية ψj من أجل 0jk1. وعلى الرغم من أن طريقة الإنشاء متشابهة في جميع الحالات، فإننا نقرر أن نتناول فقط حالة المتشعب المركزي، لأن برهانها هو الأدق بمعنى أنه يتطلب تحليلا أكثر تنقيحا في منطقة الانفجار، مؤديا إلى تصحيح لوغاريتمي ما لمتغير الانفجار كما يظهر في (1.7).

استراتيجية البرهان.

لنوجز الأفكار الرئيسة في برهان المبرهنة 1.1.

- متغيرات التشابه والمسألة الخطية. وفقا لثبات المسألة (1.1) تحت التحجيم، ندخل متغيرات التشابه

u(x,t)=(Tt)1p1w(y,s),y=x(Tt)12m,s=ln(Tt), (1.11)

مما يؤدي إلى المعادلة الجديدة

sw =𝐀mw12myw1p1w+w|w|p1
mwpp1w+w|w|p1, (1.12)

حيث يعطى المؤثر الخطي m بواسطة

m=𝐀m12my+Id. (1.13)

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

supyN|w(y,s)Φ(ys12m)|0ass+. (1.14)

ومن الطبيعي بعد ذلك خطية المعادلة (1.12) حول الملف المتوقع Φ بإدخال

q=wΦ.

مما يؤدي إلى معادلة من الشكل

sq=(m+V(y,s))q+B(q)+R(y,s), (1.15)

حيث يُبنى B(q) ليكون تربيعيا، ويقيس R الخطأ المتولد عن Φ ويكون محدودا بانتظام بـ 𝒪(s1)، وV هو الجهد المعرف كما يأتي

V(y,s)=p(Φp1κp1).

هدفنا هو أن ننشئ للمعادلة (1.15) حلا q معرفا لكل (y,s)N×[s0,+) بحيث

supyN|q(y,s)|0ass+.

- خواص المؤثر الخطي. بالنظر إلى المعادلة (1.15)، نرى أن الحد التربيعي اللاخطي وحد الخطأ صغيران ويمكن إهمالهما مقارنة بالحد الخطي. وبصورة تقريبية، سيؤدي الجزء الخطي دورا مهما في دينامية الحل. ومن الضروري تحديد الطيف والدوال الذاتية الموافقة لكل من m ومرافقه m. ووفقا لـ [Galaktionov(2001)]، يتألف طيف المؤثر الخطي m من قيم ذاتية حقيقية بسيطة فقط،

spec(m)={λβ=1|β|2m,β=(β1,,βn)n,|β|=β1++βn},

والدالة الذاتية الموافقة ψβ مع |β|=n هي كثيرة حدود من الرتبة n (انظر القضية 2.1 أدناه). علاوة على ذلك، تشكل عائلة الدوال الذاتية {ψβ}βn جزءا كاملا في Lρ2(N) حيث ρ دالة وزن متناقصة أسيا.

اعتمادا على السلوك التقاربي للجهد V، نلاحظ أن

  • داخل منطقة الانفجار، |y|Ks12m لبعض K كبير، يعد تأثير V اضطرابا لـ m.

  • خارج منطقة الانفجار، |y|Ks12m، يتصرف الجزء الخطي الكامل m+V مثل mpp1، وهو يملك طيفا سالبا صرفا. لذلك يكون ضبط الحل في هذه المنطقة بسيطا.

- تفكيك الحل والاختزال إلى مسألة منتهية الأبعاد. وفقا لطيف m، نفكك

q(y,s)=|β|=02mqβ(s)ψβ(y)+q(y,s),

حيث إن q هو إسقاط q على الفضاء الجزئي لـ m الموافق للقيم الذاتية السالبة تماما. وبما أن طيف الجزء الخطي من المعادلة التي يحققها q سالب، فإنه قابل للضبط إلى الصفر. ونود أن ننبه إلى أننا لا نستعمل تمثيل Feymann-Kac11 1 في [Bricmont and Kupiainen(1994)] و[Merle and Zaag(1997)]، تُعرّف النواة 𝒦 لشبه زمرة الحرارة المرتبطة بالمؤثر الخطي 1+V من خلال تمثيل Feymann-Kac 𝒦(s,σ,y,x)=e(sσ)1(y,x)𝑑μy,xsσ(ω)e0sσV(ω(τ),σ+τ)𝑑τ, حيث يعطى et1(y,x) بصيغة Mehler وdμy,xt هو قياس المذبذب على المسار المستمر: ω:[0,t]N مع ω(0)=x وω(t)=y. كما في الحالة m=1 المعالجة في [Merle and Zaag(1997)]، وذلك بسبب تعقيد تطبيقه في حالات الرتب الأعلى m2. ولتفادي مثل هذه الصيغة، نفكك كذلك

q(y,s)=|β|=2m+1Mqβ(s)ψβ+qM,(y,s),

لبعض M كبير بما يكفي (نموذجيا M4VLy,s). ويعطي الإسقاط المباشر

qβ=(1|β|2m)qβ+𝒪(s|β|+22m),|β|=2m+1,,M,

ومنها نحصل على التقدير الخشن |θβ(s)|=𝒪(s|β|+22m) من أجل |β|=2m+1,,M. أما بالنسبة إلى الجزء اللانهائي qM,، فنستكشف خواص شبه الزمرة esm ومتباينة Gronwall قياسية لإغلاق التقدير لهذا الجزء.

إن ضبط النمط الصفري q2m دقيق، إذ إن للجهد حجما حرجا، بمعنى ما، في تحليلنا. وبوجه خاص، نحتاج إلى تنقيح حذر للسلوك التقاربي لـ V لاشتقاق المعادلة التفاضلية العادية الحادة

qβ=2sqβ+𝒪(1s3)for|β|=2m,

التي تُظهر طيفا سالبا بعد تغيير المتغير τ=lns، ومن ثم التقدير الخشن |qβ(s)|=𝒪(logss2) من أجل |β|=2m. وهنا تكون القيمة الدقيقة لـ Bm,p المعطاة في (1.9) حاسمة في كثير من الهويات الجبرية لاشتقاق هذه المعادلة التفاضلية العادية الحادة. وفي هذه المرحلة، نختزل المسألة اللانهائية الأبعاد إلى مسألة منتهية الأبعاد بمعنى أنه يبقى ضبط عدد منته من الأنماط الموجبة qβ من أجل |β|2m1. ويتم ذلك عبر حجة طوبولوجية كلاسيكية قائمة على نظرية الدرجة.

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

2. مقاربة صورية عبر التحليل الطيفي.

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

2.1. الخواص الطيفية للمؤثر الخطي.

في هذا القسم الفرعي، نستذكر من [Galaktionov(2001)] الخواص الطيفية الأساسية للمؤثر الخطي m ومرافقه m. إن الحالة m=1 معروفة جيدا، إذ يمكننا إعادة كتابة

1=1ρ1.(ρ1)+Idwithρ1(y)=(4π)N2e|y|24,

وهو مؤثر ذاتي المرافق في فضاء Hilbert الموزون L2(e|y|2/4dy) ذي المجال 𝒟()=H2(e|y|2/4dy). وله طيف حقيقي متقطع، وتشتق الدوال الذاتية الموافقة من كثيرات حدود Hermite.

بالنسبة إلى m2، فإن المؤثر m ليس متناظرا ولا يقبل امتدادا ذاتي المرافق. لذلك نرمز بـ m إلى المرافق الصوري لـ m كما يأتي

mf=𝐀mf+12m.(yf)+f, (2.1)

بمعنى أن

<mf,g>=<f,mg>for(f,g)𝒟(m)×𝒟(m). (2.2)

من [Èĭ delman(1969)] و[Fedorjuk(1977)]، نعلم أن المعادلة الإهليلجية الآتية

𝐀mF+12m.(yF)=0inNwithNF(y)𝑑y=1. (2.3)

لها حل شعاعي وحيد يعطى بالصيغة الصريحة

F(y)=(2π)N20es2msN2JN22(s|y|)𝑑s, (2.4)

حيث إن Jν هي دالة Bessel. وبوجه خاص، يحقق F التقدير

|F(y)|<Ded|y|νinN,whereν=2m2m1, (2.5)

لبعض الثوابت الموجبة D وd المعتمدة على m وN.
ندخل دوال الوزن

ρ=ρ1,ρ(y)=ea|y|ν, (2.6)

حيث إن a=a(m,N)(0,d] ثابت صغير، وd وν معرفان في (2.5).

قضية 2.1 (الخواص الطيفية لـ m وm، [Galaktionov(2001)]).

ليكن m، لدينا
(i) إن m:Hρ2m(N)Lρ2(N) مؤثر خطي محدود طيفه

spec(m)={λβ=1|β|2m,β=(β1,,βN)N,|β|=0,1,2,}. (2.7)

ومجموعة الدوال الذاتية {ψβ}|β|0 كاملة في Lρ2(N)، حيث يملك ψβ التفكيك القابل للفصل

ψβ(y)=ψβ1(y1)ψβN(yN)withψk(ξ)=1k!j=0k2m(1)jmj!ξ2jmξk. (2.8)

(ii) إن m:Hρ2m(N)Lρ2(N) مؤثر خطي محدود طيفه

spec(m)={λβ=1|β|2m,β=(β1,,βN)N,|β|=0,1,2,}. (2.9)

ومجموعة الدوال الذاتية {ψβ}|β|0 كاملة في Lρ2(N)، حيث تعطى ψβ كما يأتي

ψβ(y)=(1)|β|β!yβF(y), (2.10)

حيث β!=β1!βN!، وyβ=β1y1βNyN والدالة F معرفة بواسطة (2.4).

(iii) (التعامد)

ψβ,ψγ=δβ,γ, (2.11)

حيث δβ,γ هو دلتا Kronecker.

ملاحظة 2.2.

نلاحظ من التعامد (2.11) ومن تعريف ψk أنه لكل كثيرة حدود Pn(y) من الدرجة n<|γ|، لدينا <Pn,ψγ>=0.

من أجل N=1، لدينا

ψk(y) =akykfor  0k<2m,ak=1k!,
ψ2mj(y) =i=0jc2miy2miandy2mj=i=0jc2miψ2mi(y)forj.

نختم هذا القسم الفرعي باستذكار خواص أساسية لشبه الزمرة esm من أجل s>0.

لمّة 2.3 (خواص شبه الزمرة esm).

تعطى نواة شبه الزمرة esm بواسطة

esm(y,x)=es[(1es)]N2mF(yes2mx[2m(1es)]12m),s>0, (2.12)

حيث يعرف F كما في (2.4). ويعرف فعل esm بواسطة

esmg(y)=Nesm(y,x)g(x)𝑑x. (2.13)

ولدينا أيضا التقديرات الآتية:
(i) esmgLCesgL لكل gL(N).
(ii) esmdiv(g)LCes(1es)12mgL لكل gL(N).
(iii) إذا كان |f(x)|η(1+|x|M+1) لكل xN، فعندئذ

|esmΠM,f(y)|CηeMs2m(1+|y|M+1),yN, (2.14)

حيث يعرف ΠM, كما في (3.10).

Proof.

يمكن التحقق من الصيغة (2.12) بحساب مباشر بفضل المعادلة (2.3). أما التقديرات (i)-(iii) فهي مباشرة من التعريفين (2.12) و(2.13). ∎

2.2. ملف انفجار تقريبي.

في هذا القسم الفرعي نستذكر المقاربة الصورية في [Galaktionov(2009)] (انظر أيضا [Budd et al.(2004)Budd, Galaktionov, and Williams]) لاستخلاص ملف انفجار ملائم لتحليلنا لاحقا. وقد استعملت هذه المقاربة في مسائل عدة تتضمن لابلاسيان من الرتبة الثانية، انظر مثلا [Bricmont and Kupiainen(1994)] و[Tayachi and Zaag(2016)] و[Ghoul et al.(2017)Ghoul, Nguyen, and Zaag, Ghoul et al.(2018a)Ghoul, Nguyen, and Zaag, Ghoul et al.(2018b)Ghoul, Nguyen, and Zaag]. وتعتمد الحجة على أساس الخواص الطيفية المعروفة للمؤثر المعاد تحجيمه m ومرافقه m المعطاة في القسم الفرعي السابق. وللتبسيط، نعتبر الحالة أحادية البعد والحلول الموجبة المتناظرة.

لندخل

w¯=wκ,

حيث إن κ=(p1)1p1 هو توازن ثابت للمعادلة (1.12). وهذا يعطي المعادلة المضطربة الآتية

sw¯=mw¯+p2κw¯2+R(w¯), (2.15)

حيث |R(w¯)|C|w¯|3 من أجل w¯1.
من القضية 2.1، نعلم أن ψn مع n2m+1 تقابل قيما ذاتية سالبة لـ m. لذلك يمكن أن نعتبر

w¯(y,s)=i=0mw¯2i(s)ψ2i(y), (2.16)

حيث i=0m|w¯2i(s)|0 عندما s+. وبإدخال هذه الفرضية في المعادلة (2.15)، وأخذ الجداء السلمي مع ψ2j، واستعمال القضية 2.1، نجد أنه من أجل 0jm،

w¯2j(s)=(1jm)w2j(s)+p2κ(j=0mw¯2i(s)ψ2i(y))2,ψ2j+𝒪(i=0m|w¯2i(s)|3). (2.17)

وبافتراض أن w¯2m هو الحد المهيمن، أي

for  0jm1,|w¯2j(s)||w¯2m(s)|ass+, (2.18)

فإن النظام (2.17) يختزل إلى

w¯2j(s) =(1jm)w2j(s)+𝒪(|w¯2m(s)|2)for0jm1,
w¯2m(s) =p2κμ¯2mw¯2m2+o(|w¯2m(s)|2),whereμ¯2m=ψ2m2,ψ2m.

وحل هذا النظام يعطي

w¯2m(s)=2κpμ¯2m1s+𝒪(ln2ss),i=0m1|w¯2i(s)|=𝒪(1s2)ass+,

وهو متوافق مع الفرضية (2.18).
ومن ثم، من (2.16) وتعريف ψ2m، نستنتج السلوك التقاربي الآتي

w(y,s)=κ2κpμ¯2ms(y2m+(1)m(2m)!(2m)!)+𝒪(lnss2), (2.19)

حيث يحدث التقارب في Lρ2() وكذلك بانتظام على المجموعات المدمجة بفعل الانتظام المكافئي القياسي.

يوفر التوسع (2.19) متغيرا ملائما للانفجار، هو z=ys12m، يتحكم في السلوك في المنطقة المتوسطة. وبوجه خاص، نحاول صوريا البحث عن حل w للمعادلة (1.12) من الشكل

w(y,s)=Φ(z)+(1)m+12κ(2m)!pμ¯2ms+𝒪(1s1+ϵ), (2.20)

لبعض ϵ>0، مع شرط الحد Φ(0)=κ.
وبإدخال هذه الفرضية في المعادلة (1.12) ومقارنة الحدود الرئيسة، نصل إلى

z2mΦ1p1Φ+Φp=0,Φ(0)=κ. (2.21)

وحل هذه المعادلة التفاضلية العادية يعطي

Φ(z)=κ(1+Bm,pz2m)1p1for some Bm,p.

وبمطابقة التوسعين (2.19) و(2.20)، نجد أن

Bm,p=2(p1)pμ¯2m(2m)!,μ¯2m=ψ2m2,ψ2m.

لنحسب μ¯2m فيما يأتي. باستعمال (2.8)، نكتب

ψ2m(y) =1(2m)!(y2m+(1)m(2m)!),
ψ4m(y) =1(4m)!(y4m+(1)m(4m)!(2m)!y2m+(4m)!),
ψ2m2(y) =1(2m)![(4m)!ψ4m+(1)m+1(2m)!((4m)!(2m)!2(2m)!)ψ2m+c0(m)ψ0].

وباستخدام علاقة التعامد (2.11)، والتعريف (2.10) لـ ψ2m، وحقيقة أن F(y)𝑑y=1، نحسب بالتكامل بالأجزاء

μ¯2m =(1)m+1(2m)!((4m)!(2m)!2(2m)!)ψ2m(y)ψ2m(y)𝑑y
=(1)m+1(2m)!((4m)![(2m)!]22)(y2m+c2m,0)y2mF(y)dy
=(1)m+1(2m)(2m)!((4m)![(2m)!]22)F(y)𝑑y
=(1)m+1(2m)!((4m)![(2m)!]22). (2.22)

وخلاصة ذلك أننا اشتققنا المرشح الآتي لملف الانفجار في متغيرات التشابه:

w(y,s)φ(y,s):=κ(1+Bm,py2ms)1p1+Am,ps, (2.23)

حيث

Bm,p=2(p1)pμ¯2m(2m)!,Am,p=(1)m+12κ(2m)!pμ¯2m. (2.24)

وبما أننا نريد الملف φ محدودا، فإن هذا يتطلب Bm,p>0 أو μ¯2m>0، ولا يحدث ذلك إلا عندما يكون m فرديا.

3. برهان المبرهنة 1.1 دون تفاصيل تقنية.

في هذا القسم نقدم جميع حجج برهان المبرهنة 1.1. ونعالج فقط الحالة أحادية البعد N=1 تبسيطا، إذ إن التحليل للحالات ذات الأبعاد الأعلى N2 هو نفسه تماما، باستثناء بعض الحسابات المعقدة لإسقاط (1.15) على الفضاءات الذاتية لـ m. ويكتمل برهان المبرهنة 1.1 في ثلاثة أجزاء:

  • في الجزء الأول نصوغ المسألة بخطية المعادلة المعاد تحجيمها (1.12) حول الملف التقريبي φ المعطى بواسطة (2.23). كما ندخل مجموعة متقلصة يُحبس فيها الحل المنشأ للمعادلة الخطية.

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

  • يخصص الجزء الأخير لبرهان البندين (i) و(iii) من المبرهنة 1.1. وعلى الرغم من أن حجة البرهان تكاد تكون هي نفسها في الحالة الكلاسيكية m=1، فإننا نرغب في إيجاز الأفكار الرئيسة تسهيلا على القارئ.

3.1. صياغة المسألة.

وفقا للتحليل الصوري المعطى في القسم 2.2، ندخل

q(y,s)=w(y,s)φ(y,s), (3.1)

ونكتب من (1.12) المعادلة التي يقودها q،

sq=(m+V(y,s))q+B(q)+R(y,s), (3.2)

حيث إن m هو المؤثر الخطي المعرف بواسطة (1.13) و

V(y,s) =p(φp1κp1), (3.3)
B(q) =(q+φ)|q+φ|p1φppφp1q, (3.4)
R(y,s) =sφ+𝐀mφ12myφφp1+φp. (3.5)

لتكن χ0𝒞0(+) دالة قطع مع supp(χ0)[0,2] وχ01 على [0,1]. ندخل

χ(y,s)=χ0(|y|Ks12m)withK1fixed, (3.6)

ونعرّف

qe(y,s)=q(y,s)(1χ(y,s)). (3.7)

نفكك

q(y,s)=k=0Mqk(s)ψk(y)+qM,(y,s), (3.8)

حيث

  • qk=Πk(q) هو إسقاط q على النمط الذاتي الموافق للقيمة الذاتية λk=1k2m، ويعرف كما يأتي

    qk(s)=<q,ψk>with ψk being introduced in (2.10). (3.9)
  • يسمى qM,=ΠM,(q) الجزء اللانهائي الأبعاد من q، حيث ΠM,(q) هو إسقاط q على الفضاء الجزئي الذاتي الذي يكون فيه طيف m أصغر من 1M2m. نلاحظ أن لدينا التعامد

    <qM,,ψk>=0forkM. (3.10)
  • إن M ثابت كبير عادة

    M=4mNwithN,NVLy,s. (3.11)

    وهو ما يسمح لنا بتطبيق متباينة Gronwall قياسية بنجاح لضبط الجزء اللانهائي الأبعاد qM,.

نهدف إلى أن ننشئ للمعادلة (1.15) حلا عالميا في الزمن q بحيث

q(s)L(N)0ass+. (3.12)

ووفقا للتفكيك (3.8)، يكفي أن نثبت وجود حل q ينتمي إلى المجموعة الآتية.

تعريف 3.1 (مجموعة متقلصة لحبس الحلول).

لكل A>0، ولكل s>0، نرمز بـ 𝒱A(s) إلى مجموعة جميع الدوال q(y,s) في L() التي تحقق

|qk(s)|As2for  0k2m1,|q2m(s)|A2logss2,
|qk(s)|Aksk+12mfor  2m+1kM,qe(s)L()AM+2s12m,
y,|qM,(y,s)|AM+1sM+22m(|y|M+1+1),

حيث إن qk وqM, وqe معرفة كما في التفكيك (3.8) و(3.7).

ملاحظة 3.2.

من التعريف، نرى أنه إذا كان q(s)𝒱A(s)، فإن الحد الآتي صحيح

|q(y,s)| CAM+1(k=0M+1(|y|k+1)sk+12m)𝟏{|y|2Ks12m}+qeL(N)
CAM+2s12mfor ally.

ومن ثم يختزل هدفنا (3.12) إلى أن ننشئ للمعادلة (1.15) حلا q(s) ينتمي إلى المجموعة المتقلصة 𝒱A(s) لكل s[s0,+).

3.2. وجود حلول محبوسة في 𝒱A.

في هذه الخطوة نهدف إلى إثبات أنه توجد فعلا معطيات ابتدائية ϕ(y)=q(y,s0) بحيث إن الحل الموافق q(y,s) لـ (1.15) ينتمي إلى المجموعة المتقلصة 𝒱A(s). وبإعطاء A1 وs0e، نعتبر معطيات ابتدائية من الشكل

ϕA,s0,𝐝(y)=As02χ(2y,s0)k=02m1dkψk(y), (3.13)

حيث إن 𝐝=(d1,,d2m1)2m معاملات حقيقية ستحدد لاحقا، وχ معرف في (3.6). وبوجه خاص، تنتمي المعطيات الابتدائية (3.13) إلى 𝒱A(s0) كما تبين القضية الآتية.

قضية 3.3 (خواص المعطيات الابتدائية (3.13)).

لكل A1، يوجد s0=s0(A)1 ومتوازي مستطيلات 𝒟s0[A,A]2m بحيث تتحقق الخواص الآتية لكل (d0,,d2m1)𝒟s0:

  • (i)

    تنتمي المعطيات الابتدائية ϕA,s0,𝐝 المعرفة في (3.13) إلى 𝒱A(s0) مع متباينات صارمة، باستثناء الأنماط الموجبة ϕk مع 0k2m1، حيث ϕk=Πk(ϕA,s0,𝐝).

  • (ii)

    التطبيق Γ:𝒟s02m، المعرف كما في Γ(d0,,d2m1)=(ϕ0,,ϕ2m1)، خطي وواحد لواحد من 𝒟s0 إلى [As02,As02]2m، ويصور 𝒟s0 في ([As02,As02]2m). علاوة على ذلك، فإن درجة Γ على الحد تختلف عن الصفر.

Proof.

يتبع البرهان مباشرة من تعريف الإسقاط Πk(ϕA,s0,𝐝)، وهو مباشر، لذلك نحذفه هنا. ∎

انطلاقا من معطيات ابتدائية ϕA,s0,𝐝 تنتمي إلى 𝒱A(s0)، ندعي أنه يمكننا ضبط المعاملات (d0,,d2m1)𝒟s0 بدقة بحيث يبقى الحل الموافق q(s) للمعادلة (1.15) في 𝒱A(s) لكل ss0. وبصورة أدق، ندعي ما يأتي.

قضية 3.4 (وجود حلول محبوسة في 𝒱A(s)).

توجد A1 وs0=s0(A)1 و(d0,,d2m1)𝒟s0 بحيث إذا كان q(s) هو حل المعادلة (1.15) ذي المعطيات الابتدائية (3.13)، فعندئذ q(s)𝒱A(s) لكل ss0(A).

Proof.

من مسألة كوشي المحلية لـ (1.1) في L()، نرى أنه لكل معطيات ابتدائية ϕA,s0,𝐝𝒱A(s0)، تملك المعادلة (1.15) حلا وحيدا q(s)𝒱A(s) لكل s[s0,s) مع s=s(𝐝). فإذا كان s=+ لبعض 𝐝=(d0,,d2m1)𝒟s0، فقد انتهينا. وإلا فإننا نسير بالمناقضة ونفترض أن s(𝐝)<+ لكل 𝐝𝒟s0. بالاستمرارية وتعريف s، نلاحظ أن q(s)𝒱A(s). وندعي الآتي.

قضية 3.5 (اختزال منته الأبعاد).

توجد A1 وs0=s0(A)1 و(d0,,d2m1)𝒟s0 بحيث تتحقق الخواص الآتية: إذا كان q(s) هو حل المعادلة (1.15) ذي المعطيات الابتدائية (3.13) وكان q(s)𝒱A(s) لكل s[s0,s1] وq(s1)𝒱A(s1)، فعندئذ

  • (i)

    (الاختزال إلى مسألة منتهية الأبعاد) (q0,,q2m1)(s1)([As12,As12]2m).

  • (ii)

    (التقاطعية) يوجد μ0>0 بحيث q(s1+μ)𝒱A(s1+μ) لكل μ(0,μ0).

لنؤجل برهان القضية 3.5 إلى القسم التالي، ونواصل حجتنا. من (i) في القضية 3.5، نحصل على (q0,,q2m1)(s)([As2,As2]2m)، ويكون التطبيق الآتي معرفا تعريفا جيدا

Θ:𝒟s0 ([1,1]2m)
(d0,,d2m1) s2A(q0,,q2m1)(s).

ومن التقاطعية المعطاة في البند (ii) من القضية 3.5، فإن (q0,,q2m1) يقطع حده فعلا عند s=s، مما يؤدي إلى استمرارية s وΘ. وباستعمال التقاطعية مرة أخرى، نرى أنه إذا كان (d0,,d2m1)𝒟s0، فإن q(s) يغادر 𝒱A(s) عند s=s0، ومن ثم s=s0 وΘ|𝒟s0=Γ، وهو التطبيق المعرف في البند (ii) من القضية 3.3. وباستعمال ذلك البند، نرى أن درجة Θ لا تساوي الصفر. وبما أن Θ مستمر، فهذا تناقض. وهذا ينهي برهان القضية 3.4، على افتراض صحة القضية 3.5. ∎

3.3. خاتمة برهان المبرهنة 1.1.

من القضية 3.4، يوجد حل q للمعادلة (1.15) بحيث q(s)𝒱A(s) لكل ss0. ومن الملاحظة 3.2، نستنتج أن q(s)L()C(A)s12m لكل ss0. ثم تتبع خلاصة البند (ii) من (1.11) و(3.1). أما البند (i) من المبرهنة 1.1 فهو نتيجة مباشرة للبندين (ii) و(iii). وبما أن برهان البند (iii) مماثل للحالة الكلاسيكية m=1، فإننا نكتفي بإيجاز الأفكار الرئيسة تسهيلا على القارئ. ويتبع وجود ملف الانفجار النهائي u𝒞({0}) من تقنية Merle [Merle(1992)]. وهنا نركز على وصف دقيق لملف الانفجار النهائي u في جوار المفردة. وللقيام بذلك، نتبع تقنية Herrero-Velázquez [Herrero and Velázquez(1993)] (انظر أيضا Bebernes-Bricher [Bebernes and Bricher(1992)] وZaag [Zaag(1998)] لمقاربة مشابهة) بإدخال الدالة المساعدة

h(ξ,τ;x0)=(Tt0(x0))1p1u(x,t), (3.14)

حيث

ξ=xx0[Tt0(x0)]12m,τ=tt0(x0)Tt0(x0),

وt0(x0) يحدد وحيدا بواسطة

|x0|=K[(Tt0(x0))|log(Tt0(x0))|]12m,K1fixed. (3.15)

نلاحظ أن h(ξ,τ;x0) هو أيضا حل لـ (1.1) بسبب ثبات (1.1) تحت التمددات. ومن (1.7)، لدينا

sup|ξ||log(Tt0(x0))|14m|h(ξ,0,x0)Φ(K)|C(Tt0(x0))12m0as|x0|0.

ليكن h^K(τ) حل (1.1) ذي المعطى الابتدائي الثابت Φ(K)، المعرف كما يأتي

h^K(τ)=κ(1τ+Bm,pK2m)1p1,τ[0,1).

وبالاستمرارية بالنسبة إلى المعطيات الابتدائية للمعادلة (1.1)، يمكن إثبات أن

sup|ξ||log(Tt0(x0))|14m,0τ<1|h(ξ,τ;x0)h^K(τ)|ϵ(x0)0,as|x0|0.

وبتمرير τ1 إلى النهاية نحصل على

u(x0)=(Tt0(x0))1p1limτ1h(0,τ;x0)(Tt0(x0))1p1h^K(1).

ومن التعريف (3.15)، نحسب

|log(Tt0(x0))|2m|log|x0||,Tt0(x0)|x0|2m2mK2m|log|x0||as|x0|0,

ومن ذلك نحصل على السلوك التقاربي

u(x0)κ(Bm,p|x0|2m2m|log|x0||)1p1as|x0|0.

وهذا ينهي برهان المبرهنة 1.1، على افتراض صحة القضية 3.5.

4. الاختزال إلى مسألة منتهية الأبعاد.

هذا القسم هو الجزء المركزي في تحليلنا، حيث نقدم جميع تفاصيل برهان القضية 3.5، وبذلك نكمل برهان المبرهنة 1.1. والفكرة الجوهرية هي إسقاط المعادلة (1.15) على مركبات مختلفة وفقا للتفكيك (3.8). وبوجه خاص، ندعي أن القضية 3.5 نتيجة مباشرة لما يأتي.

قضية 4.1 (ديناميات المعادلة (1.15)).

لكل A1، يوجد s0=s0(A)1 بحيث إذا كان q(s)𝒱A(s) لكل s[τ,τ1] مع τ1τs0، فإن التقديرات الآتية تتحقق لكل s[τ,τ1]:

  • (i)

    (ضبط الجزء منتهي الأبعاد)

    |qk(s)(1k2m)qk|Cs2for0k2m1. (4.1)
    |q2m(s)+2sq2m|CA3s3, (4.2)
    |qj(s)|Ce(j2m1)(sτ)|qj(τ)|+CAj1sj+12mfor2m+1jM. (4.3)
  • (ii)

    (ضبط الجزأين اللانهائي الأبعاد والخارجي)

    qM,(y,s)1+|y|M+1L()CeM(sτ)2mqM,(y,τ)1+|y|M+1L()+CAMsM+22. (4.4)
    qe(s)L()Ce(sτ)2(p1)qe(τ)L()+CAM+1s12m(1+sτ). (4.5)
ملاحظة 4.2.

نلاحظ أن المعادلة التفاضلية العادية الحادة (4.2) تأتي من الاختيار الدقيق للثابت Bm,p الظاهر في الملف Φ المعرف في (1.8). أما الاختيارات الأخرى فتعطي خطأ من الحجم 𝒪(1s2)، وهو كبير جدا لإغلاق التقدير لـ q2m.

لنؤجل برهان القضية 4.1 ونشرع في برهان القضية 3.5.

برهان القضية 3.5 بافتراض القضية 4.1.

بما أن حجة البرهان مشابهة لما أُنجز في [Merle and Zaag(1997)]، فإننا نوجز البرهان فقط تسهيلا على القارئ. نذكر من الفرض أن

q(s)𝒱A(s)for alls[s0,s1]andq(s1)𝒱A(s1). (4.6)

ومن ثم يثبت الجزء (i) من القضية 3.5 إذا أثبتنا أن جميع الحدود المعطاة في التعريف 3.1 يمكن تحسينها، باستثناء أول 2m مركبة qk مع 0k2m1، وذلك بمعنى أنه لكل s[s0,s1]:

|q2m(s)|<A2logss2,|qk(s)|<Aksk+12mfor2m+1kM, (4.7)
qM,(y,s)1+|y|M+1L<AM+1sM+22m,qe(s)L<AM+2s12m. (4.8)

نعالج أولا q2m. نجادل بالمناقضة مفترضين وجود s[s0,s1] بحيث لكل s[s0,s)،

|q2m(s)|<A2logss2,|q2m(s)|=A2logss2,

(لاحظ أن وجود s مضمون بالبند (i) من القضية 3.3). وباعتبار q2m(s)>0 واستعمال خاصية الأصغرية (الحالة q2m(s)<0 مشابهة)، لدينا من جهة

q2m(s)A2dds(logss2)=A2s32A2logss3,

وهذا صحيح بفضل (4.6). ومن جهة أخرى، نحصل من المعادلة التفاضلية العادية (4.2) على

q2m(s)2A2logss3+CAs3,

وينتج التناقض إذا كان A2C+1.
أما بالنسبة إلى qj مع 2m+1jM وqM, وqe، فنجادل كما يأتي: لتكن λ=logA ولنفرض أن s0λ بحيث لكل τs0 وs[τ,τ+λ] لدينا

τsτ+λτ+s02τ,hence,1τ1s.

نميز بين حالتين:
- من أجل ss0λ، نستعمل التقديرات (4.3) و(4.4) و(4.5) مع τ=s0 مع القضية 3.3 لنجد أن

|qj(s)|C(1+Aj1)sj+12m<Ajsj+12mfor2m+1jM,
qM,(y,s)1+|y|M+1LC(1+AM)sM+22m<AM+1sM+22m,
qe(s)LC(1+AM+1logA)s12m<AM+2s12m,

من أجل A كبير بما يكفي. ومن ثم، فإن التقديرين (4.7) و(4.8) يتحققان من أجل ss0λ.
- من أجل ss0>λ، نستعمل التقديرات (4.3) و(4.4) و(4.5) مع τ=sλ>s0 مع (4.6) لنكتب (تذكر أن sλ/2)

|qj(s)|e(j2m1)λAj(s/2)j+12m+CAj1sj+12m<Ajsj+12mfor  2m+1jM,
qM,(y,s)1+|y|M+1LCeMλ4mAM+1(s/2)M+22m+CAMsM+22<AM+1sM+22m,
qe(s)LCeλ2(p1)AM+2(s/2)12m+CAM+1(1+logA)s12m<AM+2s12m.

لذلك فإن التقديرين (4.7) و(4.8) يتحققان لكل s[s0,s1]، ومن ثم تتبع خلاصة الجزء (i) من القضية 3.5.

أما الجزء (ii) من القضية 3.5 فهو نتيجة مباشرة لديناميات qk المعطاة في (4.1) و(4.6). فعلا، من الجزء (i) من القضية 4.1، نعلم أن qk(s1)=ϵAs12 لبعض k{0,,2m1} وϵ=±1. وباستعمال التقدير (4.1) نحصل على

ϵqk(s1)(1k2m)ϵqk(s1)Cs12((1k2m)AC)1s12.

إذن، من أجل 0kk2m وA كبيرين بما يكفي، لدينا ϵqk(s1)>0. ومن ثم يكون qk خارجا عرضيا على منحنى الحد sϵAs2 عند s=s1. وهذا ينهي برهان القضية 3.5، على افتراض صحة القضية 4.1. ∎

نقدم الآن برهان القضية 4.1 لإكمال برهان القضية 3.5. نقسم البرهان إلى قسمين فرعيين وفقا لجزأي القضية 4.1.

4.1. ضبط الجزء منتهي الأبعاد.

نثبت في هذا الجزء البند (i) من القضية 4.1. إن إسقاطا مباشرا للمعادلة (1.15) على الدالة الذاتية ψk من أجل 0kM يعطي

qk(s)(1k2m)qk(s)=Πk(Vq+B(q)+R). (4.9)

تقدير Πk(Vq).

ندعي ما يأتي.

لمّة 4.3 (توسع V).

إن الجهد V المعرف بواسطة (3.3) يحقق التقدير

|V(y,s)|C(1+|y|2m)s,y,s1, (4.10)

ويقبل التوسع المنتظم الآتي لكل n،

V(y,s)=j=1n1sjVj(y)+𝒪(1+|y|2m(n+1)sn+1),|y|s12m,s>1, (4.11)

حيث إن Vj كثيرات حدود زوجية من الدرجة 2mj. وبصورة أدق، لدينا

V1(y)=2μ¯2mψ2mwithμ¯2m=<ψ2m2,ψ2m>. (4.12)
Proof.

التقدير (4.10) بديهي. ويتبع التقدير (4.11) من توسع Taylor في المتغير z=|y|2ms. أما الصيغة (4.12) فتأتي من التعريفين (2.23) و(2.8) لـ φ وψ2m. ∎

من اللمّة 4.3، نحصل على التقدير الآتي لـ Πk(Vq).

لمّة 4.4 (تقدير Πk(Vq)).

تحت فرضية القضية 4.1، لدينا

|Πk(Vq)|Cs2for0k2m1,
|Π2m(Vq)+2sq2m|CAs3,
|Πk(Vq)|CAk2sk+12mfor2m+1kM.
Proof.

بالتعريف 3.1، لدينا

|q(y,s)|CAM+1s1+1m(1+|y|M+1)for ally. (4.13)

وباستخدام هذا، والتقدير (4.10)، وملاحظة من التعريفين (2.10) و(2.5) أن ψk متناقصة أسيا، نحصل من أجل 0k2m1 على

|Πk(Vq)|CAM+1s2+1m(1+|y|M+1+2m)ψk(y)𝑑yCs2.

من أجل 2m+1kM، نكتب من التفكيك (3.8

Πk(Vq)=i=0Mqi(s)V(y,s)ψi(y)ψk(y)𝑑y+V(y,s)qM,(y,s)ψk(y)𝑑y.

وباستخدام (4.10) والحد المعطى لـ qM, في التعريف 3.1 نحصل على

|V(y,s)qM,(y,s)ψk(y)𝑑y|CAM+1s1+M+12m(1+|y|M+1+2m)ψk(y)𝑑yCAM+1s1+M+12m.

ومن التعامد (2.11)، نلاحظ أن <Pn,ψk>=0 لكل كثيرة حدود Pn من الدرجة nk1. ثم نستعمل (4.11) و(4.13) لتقدير

i=0Mqi(s)V(y,s)ψi(y)ψk(y)𝑑y =i=0Mj=1nqi(s)sj|y|s12mVj(y)ψi(y)ψk(y)𝑑y
+𝒪(AM+1sn+1+1m|y|s12m(1+|y|2m(n+1)+M+1)|ψk(y)|𝑑y)
+𝒪(AM+1s2+1m|y|s12m(1+|y|2m+M+1)|ψk(y)|𝑑y).

والحد الأخير محدود بـ 𝒪(ecs) بسبب التناقص الأسي لـ ψk. وبأخذ n بحيث n+1+1mk+22m، يكون الحد الثاني محدودا بـ 𝒪(AM+1sk+22m). وباستعمال حقيقة أن Vjψi كثيرة حدود من الدرجة 2mj+i، نرى أن <Vjψi,ψk>=0 من أجل 2mj+ik1. وبدمج هذا مع حدود qi المعطاة في التعريف 3.1 لـ 𝒱A، نحصل على التقدير الخشن

|i=0Mj=1nqi(s)sj|y|s12mVj(y)ψi(y)ψk(y)𝑑y|i=0Mj=1,2jm+iknCAisj+i+12m
CAk2sk+12m+i=k1MCAMs1+i+12mCAk2sk+12m.

أما بالنسبة إلى k=2m، فنحتاج إلى استعمال التعريف الدقيق (4.12) لـ V1 والعمل على نحو مماثل لما فعلناه من أجل k2m+1. فعلا، من التفكيك (3.8) والتوسع (4.11) مع n=2، لدينا

Π2m(Vq)q2ms|y|s12mV1(y)ψ2mψ2m𝑑y=i=0,i2mMqi(s)s|y|s12mV1ψiψ2m+𝒪(AM+1s3+1m).

وباستخدام التعريف (4.12) لـ V1، والحد المعطى لـ qi في التعريف 3.1، والحقيقة الآتية من التعريف (2.8) لـ ψβ ومن التعامد (2.11) وهي أن

ψ2mψiψ2m𝑑y=(c0ψ2m+i+c1ψi+c2ψi2m)ψ2m=0for  2m+1i4m1,

نصل إلى

|Π2m(Vq)+2sq2m|i=02m1CAs3+4mMCAMs1+i+12m+CAM+1s3+1mCAs3.

وهذا ينهي برهان اللمّة 4.4. ∎

ننتقل الآن إلى تقدير المساهمة الرئيسة للحد اللاخطي B(q) تحت الإسقاط Πk. نبدأ بالتوسع الآتي.

لمّة 4.5 (توسع B(q)).

لكل |q|<1 وss01 و|y|s12m، تقبل الدالة B(q) المعرفة بواسطة (3.4) التوسع المنتظم

|B(q)j=2M+1qjl=0M1sl(Bj,l(ys12m)+B~j,l(y,s))|C|q|M+2+CsM+1, (4.14)

حيث إن Bj,l(z) كثيرات حدود زوجية من الدرجة l و

|B~j,l(y,s)|C(1+|y|M+1)sM+12m.

وفضلا عن ذلك، لدينا التقدير لكل y وs1،

|B(q)|C|q|p¯withp¯=min{2,p}. (4.15)
Proof.

نلاحظ أولا من التعريف (2.23) لـ φ أن هناك ثابتين موجبين c0,C0 بحيث c0φ(y,s)C0 من أجل |y|s12m. لذلك يعطي توسع Taylor لـ B(q,φ) بدلالة q ما يأتي

|B(q)j=2M+1Bj(φ)qj|C|q|M+1,

حيث تقبل Bj(φ) التوسع الآتي بدلالة 1s،

|Bj(φ)l=0M1slBj,l(Φ)|CsM+1.

والآن، يعطي توسع Taylor لـ Bj,l(Φ) بدلالة المتغير z=ys12m النتيجة المطلوبة. وهذا ينهي برهان اللمّة 4.5. ∎

ومع اللمّة 4.5، نقدر المساهمة الرئيسة لـ B(q) تحت الإسقاط Πk. ندعي الآتي.

لمّة 4.6 (تقدير لـ Πk(B(q))).

تحت فرضية القضية 4.1، لدينا

|Πk(B(q))|Cs2for  0k2m1,|Π2m(B(q))|Cs3,
|Πk(B(q))|CAksk+22mfor  2m+1kM.
Proof.

من (4.14)، والتناقص الأسي لـ ψk و(4.13)، لدينا

Πk(B(q))=j=2M+1l=0M1sl|y|s12mqjBj,lψk𝑑y+𝒪(AM+1s2(1+1/m)+M+12m)+𝒪(ecs).

ومن التفكيك (3.8)، نكتب

qj=(i=0Mqiψi+qM,)j(q<+qM,)jforj2.

ومن التعريف 3.1 لـ 𝒱A والحدودية المنتظمة لـ q، نحصل على الحد الآتي لـ j2،

|qjq<j|C(|qM,|j+|q<|j1|qM,|)CAj(M+1)sM+4+2m2m(1+|y|j(M+1)).

لذلك،

Πk(B(q))=j=2M+1l=0M1sl|y|s12mq<jBj,lψk𝑑y+𝒪(Aj(M+1)sM+4+2m2m).

نعالج فقط الحالة k=2m، وهي الأدق. أما الحالة k2m فيمكن معالجتها على نحو مماثل، فنحذفها. من (4.13)، لدينا

j=3Ml=0M1sl|y|s12m|q<jBj,lψ2m|𝑑y+l=1M1sl|y|s12m|q<2Bj,lψ2m|𝑑yCAj(M+1)s3+1mCs3.

ويبقى ضبط الحد I:=|y|s12mq<2ψ2m𝑑y. وبحسب تعريف q<، نوسع

q<2 =(i=02mqiψi)2+2(i=02mqiψi)(i=2m+1Mqiψi)+(i=2m+1Mqiψi)2
=I1+I2+I3.

ومن التعريف 3.1، لدينا |I1|CA4log2s4(1+|y|4m) و|I2|CAM+4logss3+1/m(1+|y|M+2m). ومن ثم،

|y|s12m|(I1+I2)ψ2mdy|CA4log2s4+CAM+4logss3+1/m1s3.

لاحظ من علاقة التعامد (2.11) أنه إذا كان i+lJ2m من أجل J، فلدينا

<ψiψl,ψ2m> =(i=0i2mci,iy2imyi)(l=0l2mcl,ly2lmyl),ψ2m
=j=0i2m+l2mc^i+l,jyi+l2jm,ψ2m=j=0i2m+l2mc~i+l,jψi+l2jm,ψ2m=0.

لذلك نقدر

|y|s12mI3ψ2m𝑑y =i,l=2m+1Mci,lqi(s)ql(s)ψiψlψ2m𝑑y+𝒪(ecs)
=i,l=2m+1,i+l=J2mMci,lqi(s)ql(s)ψiψlψ2m𝑑y+𝒪(ecs).

من أجل J. وبما أن i+l4m+2، فإنه ينتج أن J3. ومن التعريف 3.1، لدينا الحد |qiql|Ai+lsi+l+22mAi+ls3+1m1s3. لذلك ||y|s12mI3ψ2m𝑑y|1s3. وهذا ينهي برهان اللمّة 4.6. ∎

ننتقل الآن إلى حد الخطأ المتولد R. نبدأ بالتوسع الآتي.

لمّة 4.7 (توسع R).

تحقق الدالة R المعرفة بواسطة (3.5) ما يأتي

R(s)LCs,

وتقبل التوسع الآتي: لكل n و|y|s12m وs1:

|R(y,s)j=1n11sj+1Rj(y)|C(1+|y|2mn)sn+1, (4.16)

حيث إن Rj كثيرات حدود من الدرجة 2mj من الشكل Rj(y)=i=0jdiy2mi. علاوة على ذلك، فإن معامل الدرجة 2m في R1 منعدم مطابقيا، ومن ثم <R1,ψ2m>=0.

Proof.

ليكن z=ys12m، ولاحظ أن Φ(z) يحقق المعادلة (2.21). لذلك نعيد كتابة (3.5) كما يأتي:

R(y,s)=z2msΦ+Am,ps2+1s𝐀mΦAm,p(p1)s+Q(Φ+Am,ps)Q(Φ), (4.17)

حيث Q(h)=hp. وبما أن c0Φ(z)C0 لكل |z|1 مع c0,C0 ثوابت موجبة ما، فإن لدينا التوسع الآتي لـ Q،

|Q(Φ+Am,ps)Q(Φ)j=1J1sjQj(Φ)|CsJ+1forJ.

ثم نوسع Qj وجميع الحدود المتبقية في (4.17) في متسلسلات قوى لـ Z=z2m لنحصل على النتيجة المطلوبة. نلاحظ أن معامل 1s في توسع R (بعد حساب أولي) يعطى بـ

Am,p(2m)!κBm,pp1=0,

حيث استعملنا (2.24). علاوة على ذلك، R1(y)=C1y2m+C2، حيث C2=C2(p,m) و

C1=κBm,pp1((4m)!(2m)!pBm,p2(p1)pAm,pκ1)=0.

ومرة أخرى، تكون القيم الدقيقة لـ Bm,p وAm,p المعطاة في (2.24) حاسمة لاستنتاج أن C1 منعدم مطابقيا. لذلك تعطي علاقة التعامد (2.11) أن <R1,ψ2m>=0. وهذا ينهي برهان اللمّة 4.7. ∎

وكنتيجة مباشرة للّمة 4.7، لدينا ما يأتي.

لمّة 4.8 (تقدير Πk(R)).

تحت فرضية القضية 4.1، لدينا

Πk(R)=𝒪(sM)for(kmod 2m)0,
|Π0(R)|Cs2,|Π2m(R)|Cs3,|Π2mi(R)|Csi+1fori.
Proof.

يتبع البرهان مباشرة من التوسع (4.16) ومن حقيقة أن

<y2mj,ψk>=i=0jai<ψ2mi,ψk>=0for(kmod 2m)0.

ومن أجل k مع (kmod 2m)0، نستعمل التوسع (4.16) مع n=M1 (يمكننا أن نستبدل M بأي عدد صحيح موجب L1) ونكتب

|Πk(R)| =|R(y,s)ψk(y)𝑑y|=||y|s12mR(y,s)ψk(y)𝑑y+|y|s12mR(y,s)ψk(y)𝑑y|
j=1M2|y|s12m|Rj(y)|sj+1|ψk(y)|𝑑y+CsM(1+|y|2m(M1))|ψk(y)|𝑑y
+Cs|y|s12m|ψk(y)|𝑑yCsM+CecsCsM.

ومن أجل (kmod 2m)0، أي k=2mi لبعض i، نستعمل (4.16) مع n=i+1 للحصول على الخلاصة. وهذا ينهي برهان اللمّة 4.8. ∎

إن جمع كل التقديرات المعطاة في اللمّات 4.4 و4.6 و4.8 والمعادلة (4.9) يعطي خلاصة الجزء (i) من القضية 4.1.

4.2. ضبط الجزأين اللانهائي الأبعاد والخارجي.

نثبت في هذا الجزء البند (ii) من القضية. نعالج أولا الجزء اللانهائي الأبعاد qM,، ثم الجزء الخارجي qe.

ضبط qM,:

بتطبيق ΠM, على المعادلة (1.15) واستعمال حقيقة أن ΠM,ψn=0 لكل nM (انظر (3.10))، نحصل على

sqM,=mqM,+ΠM,(Vq+B(q)+R). (4.18)

وبالتعريف، لدينا ΠM,(yk)=0 لكل kM. وباستعمال الحد المنتظم R(s)LCs والتوسع (4.16) مع n=M2m (انظر (3.11) لتعريف M)، وملاحظة أن ΠM,(Rj)=0 من أجل jM2m، نحصل على

|ΠM,(R)| C(1+|y|M+2m)sM2m+2𝟏|y|s12m+|R(y,s)|𝟏|y|s12mC(1+|y|M+1)sM+12m+1. (4.19)

أما تقدير ΠM,(Vq)، فندعي ما يأتي.

لمّة 4.9 (تقدير ΠM,(Vq)).

تحت فرضية القضية 4.1، لدينا

ΠM,(Vq)1+|y|M+1L()(V(s)L()+Cs)qM,1+|y|M+1L()+CAMsM+22m. (4.20)
Proof.

بالتعريف، نكتب

ΠM,(Vq)=VqM,ΠM,<(VqM,)+ΠM,(VqM,<),

حيث qM,<=ΠM,<(q)=(IdΠM,)(q). وباستعمال التقدير (4.10)، نستنتج

VqM,1+|y|M+1L+ΠM,<(VqM,)1+|y|M+1L(V(s)L()+1s)qM,1+|y|M+1L.

أما ضبط ΠM,(VqM,<)=iMΠM,(qiVψi)، فنحاجج كما يأتي:
- إذا كان Mi=0mod 2m، نأخذ n=Mi2m في التوسع (4.11) ونكتب

ΠM,(Vqiψi)=j=1nqisjΠM,(Vjψi)+ΠM,(V~nqiψi),

حيث |V~n(y,s)|C(1+|y|2m(n+1))sn+1. وباستعمال حقيقة أن ΠM,(yn)=0 من أجل nM، نحصل على j=1nΠM,(Vjqi)=0. ومن ثم،

|ΠM,(Vqiψi)|C|qi|sMi2m+1CAMsM+22m.

- إذا كان Mi=lmod 2m لبعض l(1,2m)، نأخذ n=Mil2m في التوسع (4.11) لنستنتج أن

|ΠM,(Vqiψi)|C|qi|sMil2m+1CAMsM+22m.

وهذا ينهي برهان اللمّة 4.9. ∎

وبخصوص ضبط ΠM,(B(q))، ندعي الآتي.

لمّة 4.10 (تقدير ΠM,(B(q))).

تحت فرضية القضية 4.1، لدينا

ΠM,(B(q))1+|y|M+1L()CA(M+2)2sM+1+p¯2mwithp¯=min{2,p}. (4.21)
Proof.

لتكن BM معرفة بواسطة

BM(y,s)=ΠM,<[j=2M+1qjl=0M1slBj,l(ys12m)], (4.22)

حيث إن ΠM,<=IdΠM, وBj,l معرفة في اللمّة 4.5. ثم ندعي الآتي: لكل y وs1،

|B(q)BM|CA(M+2)2sM+1+p¯2m(1+|y|M+1). (4.23)

يتبع التقدير (4.21) ببساطة من (4.23) لأن ΠM,(BM)=0. لنثبت (4.23). في المنطقة |y|s12m، نستعمل (4.15) والحدود المعطاة في التعريف 3.1 لتقدير

|B(q)|CAp¯(M+2)sp¯2mCAp¯(M+2)sM+1+p¯2m(1+|y|M+1)for|y|s12m.

وفي المنطقة |y|s12m، نكتب من اللمّة 4.5 التوسع

|B(q)ΠM,<(j=2M+1qjl=0M1slBj,l(ys12m))|
|ΠM,(j=2M+1qjl=0M1slBj,l(ys12m))|+|j=2M+1qjl=0M1slB~j,l(y,s)|+|q|M+2+CsM+1.

وتبين حجة مماثلة لما في برهان اللمّة 4.6 أن معامل الدرجة kM+1 في كثيرة الحدود j=0M+1qjl=0M1slBj,l(ys12m) محدود بـ Aksk+22m، ومن ثم

|ΠM,(j=2M+1qjl=0M1slBj,l(ys12m))|CA2(M+1)sM+32m(1+|y|M+1),|y|s12m.

ومن (4.13) والملاحظة 3.2، نحصل على

|j=2M+1qjl=0M1slB~j,l(y,s)|CA(M+2)(M+1)s1m+M+12m(1+|y|M+1)CA(M+2)2sM+1+p¯2m(1+|y|M+1),

و

|q|M+2C(AM+2s12m)M+1AM+2s1+1m(1+|y|m+1)CA(M+2)2sM+1+p¯2m(1+|y|M+1).

وهذا يكمل برهان (4.23) وكذلك برهان اللمّة 4.10. ∎

بإدخال (4.19) و(4.20) و(4.21) في المعادلة (4.18) نحصل على

sqM,=mqM,+GM,,

حيث يحقق GM,، من أجل s كبير بما يكفي،

GM,(y,s)1+|y|M+1L()V(s)L()qM,(y,s)1+|y|M+1L()+CAMsM+22m.

وباستخدام تمثيل شبه الزمرة بواسطة m، نكتب لكل s[τ,τ1]،

qM,(s)=e(sτ)mqM,(τ)+τse(ss)mGM,(s)𝑑s,

حيث يعرف esm في اللمّة 2.3. وبأخذ λ(s)=qM,(s)1+|y|M+1L() واستعمال (iii) من اللمّة (2.3)، لدينا

λ(s) eM2m(sτ)λ(τ)+τseM2m(ss)GM,(s)1+|y|M+1L()𝑑s
eM2m(sτ)λ(τ)+τseM2m(ss)(V(s)L()λ(s)+CAMsM+22m)𝑑s.

وبما أننا ثبتنا M كبيرا بحيث VLy,sM4m (انظر (3.11))، فإن حجة Gronwall قياسية مطبقة على الدالة eMs2mλ(s) تعطي

eMs2mλ(s)eM(sτ)4meMτ2mλ(τ)+CeMs2mAMsM+22,

ومنها نستنتج برهان (4.4).

ضبط qe.

نكتب من (1.15) المعادلة التي يحققها qe=(1χ(y,s))q،

sqe =(mpp1)qe+(1χ)(Q+R)+E1+E2,

حيث يعرف R بواسطة (3.5

Q =|q+φ|p1(q+φ)φp,E1=q(sχ+𝐀mχ+12myχ),
E2 =𝐀m(χq)+q𝐀mχχ𝐀mq=j=12m1cjj(q2mjχ).

وباستخدام تمثيل شبه الزمرة بواسطة m والبنود (i)(ii) من اللمّة 2.3، نكتب لكل s[τ,τ1]،

qe(s)L esτp1qe(τ)L+τsessp1(1χ)(Q(s)+R(s))+E1(s)L𝑑s
+Cj=12m1τsessp1(1e(ss))j2mq(s)2mjχ(s)L𝑑s.

ومن أجل K كبير بما يكفي، لدينا

(1χ(y,s)Q(s))LCφ(s)Lp1qe(s)L12(p1)qe(s)L.

نذكر من اللمّة 4.7 أن R(s)LCs. أما بالنسبة إلى E1، فنستعمل تعريف χ المعطى في (3.6) والحدود المعطاة في التعريف 3.1 للحصول على التقدير

E1(s)LCq(s)L(Ks12m|y|2Ks12m)CAM+1s12m.

وبما أن 2mjχ(s)LCs2mj2m، وباستعمال حقيقة أن 2mjχ ذات حامل مدمج، نقدر

j=12m1q(s)2mjχ(s)LCAM+1s1m.

ويعطي جمع التقديرات أعلاه

qe(s)L esτp1qe(τ)L
+τsessp1(12(p1)qe(s)L+CAM+1s12m+CAM+1[s1e(ss)]1m)𝑑s.

وبتطبيق متباينة Gronwall القياسية على الدالة esp1qe(s)L نحصل على

esp1qe(s)Lesτ2(p1)eτp1λ(τ)+CAM+1s12m(sτ+1),

ومنها يتبع التقدير (4.5). وهذا يكمل برهان القضية 4.1.

References

  • [Bebernes and Bricher(1992)] J. Bebernes and S. Bricher. Final time blowup profiles for semilinear parabolic equations via center manifold theory. SIAM J. Math. Anal., 23(4):852–869, 1992. ISSN 0036-1410. 10.1137/0523045. URL http://dx.doi.org/10.1137/0523045.
  • [Berger and Kohn(1988)] M. Berger and R. V. Kohn. A rescaling algorithm for the numerical calculation of blowing-up solutions. Comm. Pure Appl. Math., 41(6):841–863, 1988. ISSN 0010-3640. 10.1002/cpa.3160410606. URL http://dx.doi.org/10.1002/cpa.3160410606.
  • [Bressan(1990)] A. Bressan. On the asymptotic shape of blow-up. Indiana Univ. Math. J., 39(4):947–960, 1990. ISSN 0022-2518. 10.1512/iumj.1990.39.39045. URL http://dx.doi.org/10.1512/iumj.1990.39.39045.
  • [Bressan(1992)] A. Bressan. Stable blow-up patterns. J. Differential Equations, 98(1):57–75, 1992. ISSN 0022-0396. 10.1016/0022-0396(92)90104-U. URL http://dx.doi.org/10.1016/0022-0396(92)90104-U.
  • [Bricmont and Kupiainen(1994)] J. Bricmont and A. Kupiainen. Universality in blow-up for nonlinear heat equations. Nonlinearity, 7(2):539–575, 1994. ISSN 0951-7715. URL http://stacks.iop.org/0951-7715/7/539.
  • [Budd et al.(2004)Budd, Galaktionov, and Williams] C. J. Budd, V. A. Galaktionov, and J. F. Williams. Self-similar blow-up in higher-order semilinear parabolic equations. SIAM J. Appl. Math., 64(5):1775–1809, 2004. ISSN 0036-1399. 10.1137/S003613990241552X. URL http://dx.doi.org/10.1137/S003613990241552X.
  • [Duong()] G. K. Duong. Profile for the imaginary part of a blowup solution for a complex-valued seminar heat equation. arXiv:1712.07183. URL https://arxiv.org/abs/1712.07183.
  • [Duong et al.(2018)Duong, Nguyen, and Zaag] G. K. Duong, V. T. Nguyen, and H. Zaag. Construction of a stable blowup solution with a prescribed behavior for a non-scaling invariant semilinear heat equation. Tunisian Journal of Mathematics (to appear), 2018. URL https://arxiv.org/abs/1704.08580.
  • [Ebde and Zaag(2011)] M. A. Ebde and H. Zaag. Construction and stability of a blow up solution for a nonlinear heat equation with a gradient term. SeMA J., (55):5–21, 2011. ISSN 1575-9822.
  • [Èĭ delman(1969)] S. D. Èĭ delman. Parabolic systems. Translated from the Russian by Scripta Technica, London. North-Holland Publishing Co., Amsterdam-London; Wolters-Noordhoff Publishing, Groningen, 1969.
  • [Fedorjuk(1977)] M. V. Fedorjuk. Singularities of the kernels of Fourier integral operators, and the asymptotic behavior of the solution of a mixed problem. Uspehi Mat. Nauk, 32(6(198)):67–115, 287, 1977. ISSN 0042-1316.
  • [Fujita(1966)] H. Fujita. On the blowing up of solutions of the Cauchy problem for ut=Δu+u1+α. J. Fac. Sci. Univ. Tokyo Sect. I, 13:109–124 (1966), 1966. ISSN 0040-8980.
  • [Galaktionov(2009)] V. A. Galaktionov. Five types of blow-up in a semilinear fourth-order reaction-diffusion equation: an analytic-numerical approach. Nonlinearity, 22(7):1695–1741, 2009. ISSN 0951-7715. 10.1088/0951-7715/22/7/012. URL http://dx.doi.org/10.1088/0951-7715/22/7/012.
  • [Galaktionov and Pohozaev(2002)] V. A. Galaktionov and S. I. Pohozaev. Existence and blow-up for higher-order semilinear parabolic equations: majorizing order-preserving operators. Indiana Univ. Math. J., 51(6):1321–1338, 2002. ISSN 0022-2518. URL https://doi.org/10.1512/iumj.2002.51.2131.
  • [Galaktionov(2001)] Victor A. Galaktionov. On a spectrum of blow-up patterns for a higher-order semilinear parabolic equation. R. Soc. Lond. Proc. Ser. A Math. Phys. Eng. Sci., 457(2011):1623–1643, 2001. ISSN 1364-5021. 10.1098/rspa.2000.0733. URL http://dx.doi.org/10.1098/rspa.2000.0733.
  • [Ghoul et al.(2017)Ghoul, Nguyen, and Zaag] T. Ghoul, V. T. Nguyen, and H. Zaag. Blowup solutions for a nonlinear heat equation involving a critical power nonlinear gradient term. J. Differential Equations, 263(8):4517 – 4564, 2017. ISSN 0022-0396. http://dx.doi.org/10.1016/j.jde.2017.05.023. URL http://www.sciencedirect.com/science/article/pii/S0022039617302838.
  • [Ghoul et al.(2018a)Ghoul, Nguyen, and Zaag] T. Ghoul, V. T. Nguyen, and H. Zaag. Construction and stability of blowup solutions for a non-variational parabolic system. Ann. Inst. H. Poincaré Anal. Non Linéaire, to appear, 2018a. URL https://arxiv.org/abs/1610.09883.
  • [Ghoul et al.(2018b)Ghoul, Nguyen, and Zaag] T. Ghoul, V. T. Nguyen, and H. Zaag. Blowup solutions for a reaction-diffusion system with exponential nonlinearities. J. Differential Equations, to appear, 2018b. URL https://arxiv.org/abs/1707.08447.
  • [Herrero and Velázquez(1992)] M. A. Herrero and J. J. L. Velázquez. Generic behaviour of one-dimensional blow up patterns. Ann. Scuola Norm. Sup. Pisa Cl. Sci. (4), 19(3):381–450, 1992. ISSN 0391-173X. URL http://www.numdam.org/item?id=ASNSP_1992_4_19_3_381_0.
  • [Herrero and Velázquez(1993)] M. A. Herrero and J. J. L. Velázquez. Blow-up behaviour of one-dimensional semilinear parabolic equations. Ann. Inst. H. Poincaré Anal. Non Linéaire, 10(2):131–189, 1993. ISSN 0294-1449.
  • [Masmoudi and Zaag(2008)] N. Masmoudi and H. Zaag. Blow-up profile for the complex Ginzburg-Landau equation. J. Funct. Anal., 255(7):1613–1666, 2008. ISSN 0022-1236. 10.1016/j.jfa.2008.03.008. URL http://dx.doi.org/10.1016/j.jfa.2008.03.008.
  • [Merle(1992)] F. Merle. Solution of a nonlinear heat equation with arbitrarily given blow-up points. Comm. Pure Appl. Math., 45(3):263–300, 1992. ISSN 0010-3640. 10.1002/cpa.3160450303. URL http://dx.doi.org/10.1002/cpa.3160450303.
  • [Merle and Zaag(1997)] F. Merle and H. Zaag. Stability of the blow-up profile for equations of the type ut=Δu+|u|p1u. Duke Math. J., 86(1):143–195, 1997. ISSN 0012-7094. 10.1215/S0012-7094-97-08605-1. URL http://dx.doi.org/10.1215/S0012-7094-97-08605-1.
  • [Nguyen(2017)] V. T. Nguyen. Numerical analysis of the rescaling method for parabolic problems with blow-up in finite time. Phys. D, 339:49–65, 2017. ISSN 0167-2789. URL https://doi.org/10.1016/j.physd.2016.09.002.
  • [Nguyen and Zaag(2016)] V. T. Nguyen and H. Zaag. Construction of a stable blow-up solution for a class of strongly perturbed semilinear heat equations. Ann. Scuola Norm. Sup. Pisa Cl. Sci., 16(4):1275–1314, 2016. 10.2422/2036-2145.201412_001. URL http://dx.doi.org/10.2422/2036-2145.201412_001.
  • [Nouaili and Zaag(2015)] N. Nouaili and H. Zaag. Profile for a simultaneously blowing up solution to a complex valued semilinear heat equation. Comm. Partial Differential Equations, 40(7):1197–1217, 2015. ISSN 0360-5302. 10.1080/03605302.2015.1018997. URL http://dx.doi.org/10.1080/03605302.2015.1018997.
  • [Nouaili and Zaag(2018)] N. Nouaili and H. Zaag. Construction of a blow-up solution for the Complex Ginzburg-Landau equation in some critical case. Arch. Ration. Mech. Anal., to appear, 2018. URL https://arxiv.org/abs/1703.00081.
  • [Peletier and Troy(2001)] L. A. Peletier and W. C. Troy. Spatial patterns, volume 45 of Progress in Nonlinear Differential Equations and their Applications. Birkhäuser Boston, Inc., Boston, MA, 2001. ISBN 0-8176-4110-6. URL https://doi.org/10.1007/978-1-4612-0135-9. Higher order models in physics and mechanics.
  • [Tayachi and Zaag(2015)] S. Tayachi and H. Zaag. Existence and stability of a blow-up solution with a new prescribed behavior for a heat equation with a critical nonlinear gradient term. Actes du Colloque EDP-Normandie, Le Havre, arXiv:1610.01289, 2015.
  • [Tayachi and Zaag(2016)] S. Tayachi and H. Zaag. Existence of a stable blow-up profile for the nonlinear heat equation with a critical power nonlinear gradient term. arXiv:1506.08306, 2016.
  • [Zaag(1998)] H. Zaag. Blow-up results for vector-valued nonlinear heat equations with no gradient structure. Ann. Inst. H. Poincaré Anal. Non Linéaire, 15(5):581–622, 1998. ISSN 0294-1449. 10.1016/S0294-1449(98)80002-4. URL http://dx.doi.org/10.1016/S0294-1449(98)80002-4.