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

الملخص

ندرس في هذه المذكرة نظام الحرارة شبه الخطي

tu=Δu+f(v),tv=μΔv+g(u),μ>0,

حيث لا تمتلك اللاخطية بنية تدرجية، وتتخذ الصورة الخاصة

f(v)=v|v|p1andg(u)=u|u|q1withp,q>1,

أو

f(v)=epvandg(u)=equwithp,q>0.

نبرهن وجود حلول انفجارية من النمط I لهذا النظام، ونقدم وصفا دقيقا لملامح انفجارها. وتعتمد الطريقة على إجراء من خطوتين: اختزال المسألة إلى مسألة منتهية البعد عبر تحليل طيفي، ثم حل المسألة المنتهية البعد بحجة طوبولوجية كلاسيكية مبنية على نظرية الدليل. وكنتيجة لتقنيتنا، تكون الحلول المنشأة مستقرة إزاء اضطراب صغير في المعطيات الابتدائية. ويمكن العثور على النتائج والحجج الرئيسية المعروضة في هذه المذكرة في بحثينا [Ghoul et al.(2018c)Ghoul, Nguyen, and Zaag, Ghoul et al.(2018d)Ghoul, Nguyen, and Zaag].

الكلمات المفتاحية:
حل انفجاري، ملمح الانفجار، الاستقرار، نظام قطع مكافئ شبه خطي
1991 Mathematics Subject Classification:
رئيسي: 35K50, 35B40; ثانوي: 35K55, 35K57.
H. Zaag is supported by the ANR project ANAÉ ref. ANR-13-BS01-0010-03.
——————–
May 13, 2026

Tej-Eddine Ghoul, Van Tien Nguyen و 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. المقدمة.

في [Ghoul et al.(2018c)Ghoul, Nguyen, and Zaag, Ghoul et al.(2018d)Ghoul, Nguyen, and Zaag]، ندرس النظام القطع المكافئ شبه الخطي

{tu=Δu+f(v),tv=μΔv+g(u),(x,t)+×N, (1.1)

مع N1، وμ>0 و(u(0),v(0))=(u0,v0)، حيث (u,v)(t):xN2 ولا تمتلك اللاخطية بنية تدرجية، وتتخذ الصورة الخاصة

f(v)=v|v|p1andg(u)=u|u|q1withp,q>1, (1.2)

أو

f(v)=epvandg(u)=equwithp,q>0. (1.3)

يمثل النظام (1.1) نموذجا بسيطا لنظام تفاعل-انتشار يصف انتشار الحرارة في مزيج قابل للاحتراق ذي مكونين، ولذلك كان موضوعا لدراسة مكثفة خلال العقدين الماضيين (انظر [Souplet(2005)]، و[Zheng et al.(2002)Zheng, Zhao, and Chen] والمراجع الواردة فيهما). نهتم هنا أساسا بإثبات وجود حلول انفجارية في زمن منته واستقرارها، مع تحقق سلوك تقاربي محدد سلفا. ونقصد بالانفجار في زمن منته أن T=T(u0,v0)، زمن الوجود الأعظمي للحل الكلاسيكي (u,v) للمسألة (1.1)، منته، وأن الحل ينفجر في الزمن المنتهي T بالمعنى الآتي:

limtT(u(t)L(N)+v(t)L(N))=+.

فضلا عن ذلك، يسمى حل الانفجار المنتهي (u,v) للنظام (1.1) حلا من النمط I إذا وجدت ثابتة موجبة C بحيث

u(t)L(N)Cu¯(t),v(t)L(N)Cv¯(t), (1.4)

حيث إن (u¯,v¯)(t) هو حل الانفجار الموجب الوحيد للنظام التفاضلي العادي المرتبط بـ (1.1)، أي إن

u¯(t)=Γ(Tt)α,v¯(t)=γ(Tt)β for (1.2),
u¯(t)=ln[(p(Tt))1q],v¯(t)=ln[(q(Tt))1p] for (1.3),

حيث يحدد (Γ,γ) بواسطة

γp=αΓ,Γq=γβ,α=p+1pq1,β=q+1pq1. (1.5)

وخلاف ذلك، يكون حل الانفجار من النمط II.

أما بالنسبة إلى النظام (1.1)-(1.2) مع μ=1، فقد اشتق فريدمان-غيغا [Friedman and Giga(1987)] وإسكوبيدو-هيريرو [Escobedo and Herrero(1991a)] وجود حلول انفجارية في زمن منته (انظر أيضا [Escobedo and Herrero(1991b)]، و[Escobedo and Herrero(1993)]، إلخ). ومن عمل أندروتشي-هيريرو-فيلاسكيز [Andreucci et al.(1997)Andreucci, Herrero, and Velázquez]، نعلم أن التقدير (1.4) صحيح إذا كان

pq>1,q(pN2)+<N+2orp(qN2)+<N+2.

انظر أيضا كاريستي-ميتيدييري [Caristi and Mitidieri(1994)]، ودنغ [Deng(1996)]، وفيلا-سوبليه [Fila and Souplet(2001)] لمزيد من النتائج المتعلقة بالتقدير (1.4). وبمعرفة أن الحل يبدي انفجارا من النمط I، تمكن مؤلفو [Andreucci et al.(1997)Andreucci, Herrero, and Velázquez] من الحصول على معلومات أدق عن السلوك التقاربي للحل قرب المفردة. ثم حسّن زاغ [Zaag(2001)] نتائجهم لاحقا. وعندما μ1، عرفت نتائج أقل بكثير، باستثناء محمودي-سوبليه-تاياشي [Mahmoudi et al.(2015)Mahmoudi, Souplet, and Tayachi] الذي يثبت نتيجة انفجار عند نقطة وحيدة تحسن النتيجة المحصل عليها في [Friedman and Giga(1987)]. أما بالنسبة إلى النظام (1.1) المقترن باللاخطية (1.3)، فالنتيجة الوحيدة المعروفة ترجع إلى سوبليه-تاياشي [Souplet and Tayachi(2016)]، إذ كيّفا التقنية المطورة في [Mahmoudi et al.(2015)Mahmoudi, Souplet, and Tayachi] للحصول على نتيجة الانفجار عند نقطة وحيدة لصنف من الحلول المتناقصة شعاعيا. وعلى حد علمنا، لا توجد نتائج تتعلق بالسلوك التقاربي، حتى في حالة تساوي الانتشار، أي μ=1. ونذكّر أيضا بأن دراسة النظام القطع المكافئ غير متساوي الانتشار (1.1) (قد يساوي μ أو لا يساوي 1) هي، بوجه خاص، أشد تعقيدا بكثير، سواء من حيث سلوك الحلول أو على المستوى التقني.

في هذه المذكرة نبرهن وجود حلول انفجارية من النمط I للنظام (1.1)، ونقدم أول وصف كامل لسلوكها التقاربي. وبصورة أدق، نثبت في [Ghoul et al.(2018c)Ghoul, Nguyen, and Zaag] المبرهنة الآتية.

مبرهنة 1.1 (حلول انفجارية من النمط I لـ (1.1)-(1.2) وسلوكها التقاربي، [Ghoul et al.(2018c)Ghoul, Nguyen, and Zaag]).

لتكن aN وT>0. توجد معطيات ابتدائية (u0,v0)L(N)×L(N) يكون من أجلها للنظام (1.1)-(1.2) الحل الوحيد (u,v) المعرف على N×[0,T) بحيث

  • (i)

    ينفجر الحل (u,v) في زمن منته T عند النقطة الوحيدة a.

  • (ii)

    (الملمح التقاربي) يتحقق لكل t[0,T) أن

    (Tt)αu(x,t)Φ(z)L(N)+(Tt)βv(x,t)Ψ(z)L(N)C(Tt), (1.6)

    حيث z=xa(Tt)|log(Tt)| والملمحان Φ0 وΨ0 يعطيان صراحة كما يأتي:

    zN,Φ(z)=Γ(1+b|z|2)α,Ψ(z)=γ(1+b|z|2)β, (1.7)

    مع كون Γ,γ,α,β معرفا في (1.5)، ومع

    b=b(μ,p,q)=(pq1)(2pq+p+q)4pq(p+1)(q+1)(μ+1)>0. (1.8)
  • (iii)

    (ملمح الانفجار النهائي) لكل xa، لدينا (u(x,t),v(x,t))(u(x),v(x))[𝒞2(N{0})]2 مع

    u(x)Γ(b|xa|22|log|xa||)p+1pq1,v(x)γ(b|xa|22|log|xa||)q+1pq1 (1.9)

    عندما |xa|0.

ملاحظة 1.2.

الملمح التقاربي المعرف في (1.7) مع μ=1 هو أحد الملامح ضمن نتيجة التصنيف المثبتة في [Andreucci et al.(1997)Andreucci, Herrero, and Velázquez] (انظر أيضا [Zaag(2001)]). وهذا يعني أننا نستطيع إنشاء حلول انفجارية من النمط I لـ (1.1)-(1.2) تحقق الملامح التقاربية الأخرى الموصوفة كما في [Andreucci et al.(1997)Andreucci, Herrero, and Velázquez]. غير أن تلك الإنشاءات ستكون أبسط من الحالة التي ندرسها (1.6)، إذ تتضمن تصحيحا لوغاريتميا لمتغير الانفجار.

ملاحظة 1.3.

التقدير (1.9) حاد بالمقارنة مع النتيجة المثبتة في [Mahmoudi et al.(2015)Mahmoudi, Souplet, and Tayachi] (انظر المبرهنة 1.3)، حيث لم يتمكن المؤلفون إلا من الحصول على تقديرات نقطية سفلية من دون التصحيح اللوغاريتمي.

أما بالنسبة إلى النظام (1.1) المقترن باللاخطية (1.3)، فإننا ندرسه في الفضاء التآلفي الخاص α من أجل ثابتة موجبة ما α،

α={(u,v)(ϕ¯,ψ¯)+L(N)×L(N)whereqϕ¯(x)=pψ¯(x)=ln(1+α|x|2)},

ونثبت في [Ghoul et al.(2018d)Ghoul, Nguyen, and Zaag] النتيجة الآتية:

مبرهنة 1.4 (حلول انفجارية من النمط I لـ (1.1)-(1.3) وسلوكها التقاربي، [Ghoul et al.(2018d)Ghoul, Nguyen, and Zaag]).

لتكن aN وT>0. توجد معطيات ابتدائية (u0,v0)α يكون من أجلها للنظام (1.1)-(1.3) الحل الوحيد (u,v) المعرف على N×[0,T) بحيث

  • (i)

    تنفجر الدالة (equ,epv) في زمن منته T عند النقطة الوحيدة a.

  • (ii)

    (الملمح التقاربي) يتحقق لكل t[0,T) أن

    (Tt)equ(x,t)Φ(z)L(N)+(Tt)epv(x,t)Ψ(z)L(N)C(Tt), (1.10)

    حيث z=xa(Tt)|log(Tt)| والملمحان Φ وΨ يعطيان صراحة كما يأتي:

    zN,pΦ(z)=qΨ(z)=11+b|z|2withb=12(μ+1). (1.11)
  • (iii)

    (ملمح الانفجار النهائي) لكل xa، لدينا (u(x,t),v(x,t))(u(x),v(x))[𝒞2(N{0})]2 مع

    u(x)1qln(2bp|log|xa|||xa|2),v(x)1pln(2bq|log|xa|||xa|2)as|xa|0. (1.12)
ملاحظة 1.5.

نستطيع إنشاء حلول انفجارية من النمط I لـ (1.1)-(1.3) تحقق ملامح انفجار مختلفة لا تتضمن تصحيحا لوغاريتميا لمتغير الانفجار الموصوف كما في (1.10). وهذا يعني أنه يمكننا الحصول على نتيجة تصنيف مماثلة للحلول الانفجارية من النمط I لـ (1.1)-(1.3) بتكييف تقنية [Andreucci et al.(1997)Andreucci, Herrero, and Velázquez] مع بعض الصعوبات التقنية الإضافية.

يعتمد برهان المبرهنة 1.1 والمبرهنة 1.4 على إجراء من خطوتين:

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

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

لقد طبق هذا الإجراء ذو الخطوتين بنجاح على معادلات تطور لاخطية متنوعة لإنشاء حلول انفجارية من النمطين I وII. وكان ذلك في حالة معادلة الحرارة شبه الخطية المعالجة في [Bricmont and Kupiainen(1994)]، و[Merle and Zaag(1997)]، و[Nguyen and Zaag(2017)] (انظر أيضا [Nguyen and Zaag(2016)]، و[Duong et al.(2018)Duong, Nguyen, and Zaag] لحالة الاضطرابات اللوغاريتمية، و[Bressan(1990)]، و[Bressan(1992)] و[Ghoul et al.(2017)Ghoul, Nguyen, and Zaag] للمصدر الأسي، و[Nouaili and Zaag(2015)] للحالة ذات القيم العقدية)، ومعادلة Ginzburg-Landau في [Masmoudi and Zaag(2008)]، و[Nouaili and Zaag(2018)] (انظر أيضا [Zaag(1998)] لعمل سابق). كما ظهر في معادلة Schrödinger اللاخطية في حالتي حرج الكتلة [Merle and Raphael(2003), Merle and Raphael(2004), Merle and Raphael(2005a), Merle and Raphael(2005b)] وفوق حرج الكتلة [Merle et al.(2015)Merle, Raphaël, and Rodnianski]؛ ومعادلة الموجة ذات الطاقة الحرجة [Duyckaerts et al.(2013)Duyckaerts, Kenig, and Merle]، و[Hillairet and Raphaël(2012)] وفوق الحرجة [C.(2016b)]؛ ومعادلة gKdV ذات الكتلة الحرجة [Martel et al.(2014)Martel, Merle, and Raphaël, Martel et al.(2015a)Martel, Merle, and Raphaël, Martel et al.(2015b)Martel, Merle, and Raphaël]؛ ونموذج Keller-Segel ثنائي البعد [Raphaël and Schweyer(2014b)]؛ والمعادلات الهندسية ذات الطاقة الحرجة وفوق الحرجة: خرائط الموجة [Raphaël and Rodnianski(2012)] و[Ghoul et al.(2018b)Ghoul, Ibrahim, and Nguyen]، وخرائط Schrödinger [Merle et al.(2013)Merle, Raphaël, and Rodnianski] وتدفق الحرارة التوافقي [Raphaël and Schweyer(2013), Raphaël and Schweyer(2014a)] و[Ghoul et al.(2018a)Ghoul, Ibrahim, and Nguyen]؛ ومعادلة الحرارة شبه الخطية في حالتي الطاقة الحرجة [Schweyer(2012)] وفوق الحرجة [C.(2016a)].

وكنتيجة لتقنيتنا، حصلنا على نتيجة الاستقرار الآتية.

مبرهنة 1.6.

الحلول المنشأة الموصوفة في المبرهنة 1.1 والمبرهنة 1.4 مستقرة بالنسبة إلى المعطيات الابتدائية.

ملاحظة 1.7.

يمكن فهم الفكرة الكامنة وراء نتيجة الاستقرار صوريا من ثبات المسألة تحت انتقال المكان-الزمان وتغيير القياس على النحو الآتي: للمؤثر الخطي حول الملمح المتوقع قيمتان ذاتيتان موجبتان λ0=1,λ1=12، وقيمة ذاتية صفرية λ2=0، ثم طيف سالب متقطع لا نهائي. ومن تحليل استقرار مسائل الانفجار، فإن المركبة الموافقة لـ λ0=1 لها نمو أسي es، ويمكن إلغاؤه بتغيير زمن الانفجار؛ وبالمثل بالنسبة إلى النمط λ1=12 عبر إزاحة نقطة الانفجار؛ أما النمط الحيادي λ2=0 فله عادة نمو كثير حدودي، ويمكن إلغاؤه باستعمال ثبات المسألة تحت تغيير القياس. وبما أن الأنماط المتبقية للمؤثر الخطي، الموافقة للطيف السالب، تضمحل أسيا، فإن المرء يستنتج استقرار الحلول المنشأة الموصوفة في المبرهنة 1.1 والمبرهنة 1.4. ومن نتيجة الاستقرار نتوقع أن ملامح الانفجار (1.7) و(1.11) عامة، أي إن ملامح الانفجار الأخرى غير مستقرة. وفي رأينا، هذه مسألة مفتوحة صعبة أعطى هيريرو-فيلاسكيز [Herrero and Velázquez(1992)] جوابا جزئيا خاصا عليها في حالة معادلة الحرارة شبه الخطية أحادية البعد.

2. حساب صوري لملمح الانفجار.

نستذكر بإيجاز في هذا القسم المقاربات الصورية في [Ghoul et al.(2018c)Ghoul, Nguyen, and Zaag, Ghoul et al.(2018d)Ghoul, Nguyen, and Zaag] لإنشاء ملمح انفجار تقريبي مناسب للنظام (1.1). ويمكن العثور على مقاربات مشابهة في [Tayachi and Zaag(2015), Tayachi and Zaag(2016)]، و[Ghoul et al.(2017)Ghoul, Nguyen, and Zaag]، و[Nouaili and Zaag(2018)] والمراجع الواردة فيها. وتعتمد الطريقة على التوسيعات التقاربية المتطابقة التي ترتكز أساسا على الخواص الطيفية للمؤثر الخطي حول ملمح متوقع.

متغيرات التشابه:

نجري تغيير المتغيرات المعروف

Φ(y,s)=(Tt)αu(x,t),Ψ(y,s)=(Tt)βv(x,t) for(1.2), (2.1)
Φ(y,s)=(Tt)equ(x,t),Ψ(y,s)=(Tt)epv(x,t) for(1.3), (2.2)

حيث أدخلت α,β في (1.5) و

y=xTt,s=log(Tt).

وبهذه الطريقة، تحقق (Φ,Ψ) النظام الجديد

{sΦ=1ΦαΦ+|Ψ|p1Ψ,sΨ=μΨβΨ+|Φ|p1Φ, for (1.2), (2.5)
{sΦ=1ΦΦ|Φ|2Φ+qΦΨ,sΨ=μΨΨμ|Ψ|2Ψ+pΦΨ, for (1.3), (2.8)

حيث

ηf=ηΔfy2f=ηρη(ρηf)withη{1,μ}, (2.9)

هو المؤثر ذاتي المرافقة بالنسبة إلى فضاء هيلبرت Lρη2(N,) المجهز بالجداء الداخلي

<f,g>Lρη2=Nf(y)g(y)ρη(y)dywithρη(y)=1(4π)N/2e|y|24η.

المسألة الخطية والخواص الطيفية للمؤثر الخطي المرتبط بها:

نلاحظ أن الحلول الثابتة غير الصفرية للنظامين (2.5) و(2.8) هي (Γ,γ) و(1/p,1/q) على الترتيب. وهذا يقترح التخطية

(Φ¯,Ψ¯)=(ΦΓ,Ψγ)for (1.2)and(Φ¯,Ψ¯)=(Φ1/p,Ψ1/q)for (1.3), (2.10)

حيث تحل (Φ¯,Ψ¯) النظام

i=1,2,s(Φ¯Ψ¯)=(+Mci)(Φ¯Ψ¯)+(Qi,1Qi,2), (2.11)

حيث ترمز i=1 إلى اللاخطية كثيرة الحدود (1.2)، وترمز i=2 إلى الحالة الأسية (1.3)؛ كما تبنى Qi,1 وQi,2 لتكونا تربيعيتين، ويعرف المؤثر الخطي والمصفوفات Mci كما يأتي:

=(100μ),Mc1=(αpγp1qΓq1β),Mc2=(0q/pp/q0). (2.12)

تعطي اللمّة الآتية الخواص الطيفية لـ +Mci.

لمّة 2.1 (قطرنة +Mci).

لكل n، توجد كثيرتا حدود fn,gn,f~n وg~n من الدرجة n بحيث

(+Mci)(fngn)=(1n2)(fngn),(+Mci)(f~ng~n)=λi,n(f~ng~n), (2.13)

حيث

λ1,n=(n2+(p+1)(q+1)pq1),λ2,n=(1+n2).
Proof.

انظر اللمّة 3.2 في [Ghoul et al.(2018c)Ghoul, Nguyen, and Zaag] للحالة كثيرة الحدود (1.2)، واللمّة 2.2 في [Ghoul et al.(2018d)Ghoul, Nguyen, and Zaag] للحالة الأسية (1.3). وندعو القارئ إلى الاطلاع على الصيغ الدقيقة للدوال الذاتية، وكذلك على تعريف مناسب للإسقاط وفق هذه الأنماط الذاتية في هذين البحثين. ∎

التوسع الداخلي:

من اللمّة 2.1، نعلم أن (fngn)n3 و(f~ng~n)n0 توافقان قيما ذاتية سالبة لـ +Mci؛ لذلك يمكننا، تحت فرضية التناظر الشعاعي للحل، اعتبار التوسع الصوري الآتي:

(Φ¯Ψ¯)(y,s)=a0(s)(f0g0)(y)+a2(s)(f2g2), (2.14)

حيث |a0(s)|+|a2(s)|0 عندما s+. وبإدخال هذه الفرضية الشكلية في (2.11) ثم الإسقاط على (fkgk),k=0,2 نحصل على النظام التفاضلي العادي

{a0=a0+𝒪(|a0|2+|a2|2),a2=ca22+𝒪(|a2|3+|a0a2|+|a0|3), (2.15)

حيث

c=2pq+p+q4pq(p+1)(q+1)(μ+1)for (1.2)andc=2pq(μ+1)for (1.3).

وبافتراض أن |a0(s)|=o(|a2(s)|) عندما s+، نحصل على

a2(s)=1cs+𝒪(logss2)and|a0(s)|=𝒪(1s2)ass+.

ومن (2.14)، و(2.10) وتعريف الدالة الذاتية (f2g2)، ننتهي إلى السلوك التقاربي

{Φ(y,s)=Γ[1p+1c|y|2s2p(1μ)cs]+𝒪(logss2),Ψ(y,s)=γ[1q+1c|y|2s2q(μ1)cs]+𝒪(logss2), for (2.5), (2.18)
{Φ(y,s)=1p[1pqc|y|2s+2μpqcs]+𝒪(logss2),Ψ(y,s)=1q[1pqc|y2|s+2pqcs]+𝒪(logss2), for (2.8), (2.21)

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

التوسع الخارجي:

توفر التوسعات التقاربية أعلاه متغير انفجار ملائما

z=ys=x(Tt)|log(Tt)|.

ثم نحاول البحث عن حل تقريبي لـ (2.5) (على الترتيب (2.8)) من الشكل

(ΦΨ)(y,s)=(Φ0Ψ0)(z)+1s(Φ1Φ1)(z)+, (2.22)

وبإدخال هذه الفرضية الشكلية في (2.5) (على الترتيب (2.8)) نحصل على النظام من الرتبة الرئيسية

z2Φ0αΦ0+Φ0p=0,z2Ψ0βΨ0+Ψ0q=0, for (2.5), (2.23)
z2Φ0Φ0+qΦ0Ψ0=0,z2Ψ0Ψ0+pΦ0Ψ0=0, for (2.8), (2.24)

خاضعا للشرط الابتدائي

(Φ0,Ψ0)(0)=(Γ,γ)for (2.18)and(Φ0,Ψ0)(0)=(1/p,1/q)for (2.21).

وتعطى حلول هذا النظام صراحة بواسطة

Φ0(z)=Γ(1+b|z|2)α,Ψ0(z)=γ(1+b|z|2)β for (2.23), (2.25)
Φ0(z)=1p(1+b|z|2),Ψ0(z)=1q(1+b|z|2) for (2.24), (2.26)

حيث b>0 ثابت تكامل. وبمطابقة التوسعات التقاربية (2.25) مع (2.18)، و(2.26) مع (2.21)، نحصل بدقة على قيمة الثابت b كما ورد في المبرهنتين 1.1 و1.4 على الترتيب.

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

{Φ(y,s)φ(y,s)=Φ0(ys)2Γp(1μ)cs,Ψ(y,s)ψ(y,s)=Ψ0(ys)2γq(μ1)cs, for (2.5), (2.29)
{Φ(y,s)=φ(y,s)=Φ0(ys)+2μqcs,Ψ(y,s)=ψ(y,s)=Ψ0(ys)+2pcs, for (2.8). (2.32)

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

نعرض الحجج الرئيسية كلها لبرهان الوجود من دون تفاصيل تقنية، ونحيل القارئ المهتم بشأنها إلى بحثينا [Ghoul et al.(2018c)Ghoul, Nguyen, and Zaag, Ghoul et al.(2018d)Ghoul, Nguyen, and Zaag]. نعالج أولا الحالة كثيرة الحدود (1.2)، أي برهان المبرهنة 1.1، ثم الحالة الأسية (1.3)، أي برهان المبرهنة 1.4، وهي أدق بسبب وجود الحدين |Φ|2Φ و|Ψ|2Ψ في صياغة متغيرات التشابه (انظر (2.8)).

3.1. الحالة كثيرة الحدود (1.2).

يكرس هذا القسم الفرعي لبرهان الجزء (ii) من المبرهنة 1.1. والجزآن (i) و(iii) نتيجتان للجزء (ii). ويمكن للقارئ أن يجد جميع تفاصيل البرهان في [Ghoul et al.(2018c)Ghoul, Nguyen, and Zaag].

صياغة المسألة:

بالنظر إلى متغيرات التشابه (2.1)، نرى أن إنشاء حلول انفجارية لـ (1.1) المقترنة بـ (1.2) وتحقق السلوك التقاربي (1.6) يكافئ، من أجل (2.5)، إنشاء حل عام في الزمن (Φ,Ψ) بحيث

supyN(|Φ(y,s)Φ(y/s)|+|Ψ(y,s)Ψ(y/s)|)0ass+, (3.1)

حيث Φ وΨ هما الملمحان المعرفان في المبرهنة 1.1. وانطلاقا من الحساب الصوري لملمح انفجار تقريبي المعروض في القسم السابق، نقوم بتخطية (2.5) حول (φ,ψ) المعرف في (2.29) بدلا من (Φ,Ψ)، أي إننا ندخل

(ΛΥ)=(ΦΨ)(φψ), (3.2)

مما يؤدي إلى النظام الخطي

s(ΛΥ)=(+Mc1+V(y,s))(ΛΥ)+(F1(Υ,y,s)F2(Λ,y,s))+(R1(y,s)R2(y,s)), (3.3)

حيث تعرف وMc1 في (2.12

V(y,s)=(0p(ψp1γp1)q(φq1Γq1)0)(0V1V20), (3.4)
(F1(Υ,y,s)F2(Λ,y,s))=(|Υ+ψ|p1(Υ+ψ)ψppψp1Υ|Λ+φ|q1(Λ+φ)φqqφq1Λ), (3.5)

و

(R1(y,s)R2(y,s))=(sφ+Δφ12yφ(p+1pq1)φ+ψpsψ+μΔψ12yψ(q+1pq1)ψ+φq). (3.6)

ويصبح هدفنا إنشاء، للنظام (3.3)، حل عام في الزمن (Λ,Υ) يحقق

supyN(|Λ(y,s)|+|Υ(y,s)|)0ass+. (3.7)

وبما أن الحل (Λ,Υ) يؤول إلى الصفر عندما s+، وأن الحد اللاخطي (F1,F2) بني ليكون تربيعيا وأن حد الخطأ (R1,R2) من الحجم s1، نرى أن ديناميات (3.3) تتأثر بقوة بالجزء الخطي +Mc1+V. وهنا يتصرف الكمون V على نحو مختلف كما يأتي:
- المنطقة الخارجية، أي |y|s: لكل ϵ>0، يوجد Kϵ>0 وsϵ>0 بحيث

sup|y|Kϵs,ssϵ|V(y,s)|ϵ.

من اللمّة 2.1، نرى أن المؤثر الخطي +Mc1+V يتصرف كمؤثر ذي طيف سالب تماما في المنطقة الخارجية، مما يجعل التحليل في هذه المنطقة أبسط.
- المنطقة الداخلية، أي |y|s: يعد الكمون V اضطرابا للجزء الخطي +Mc1.

وبما أن سلوك V في المنطقتين الداخلية والخارجية مختلف، فهذا يقترح النظر في ديناميات (3.3) من أجل |y|s و|y|s على حدة. ولهذه الغاية، ندخل دالة القطع

χ(y,s)=χ0(|y|Ks),χ0𝒞0(+,[0,1]),χ0(r)={1forr[0,1],0forr2, (3.8)

حيث K ثابت موجب سيختار كبيرا بما يكفي. ثم نعرف

(ΛeΥe)=(1χ(y,s))(ΛΥ) (3.9)

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

(ΛΥ)(y,s)=nM[θn(s)(fngn)+θ~n(f~ng~n)]+(ΛΥ)(y,s), (3.10)

حيث θn=Πn(ΛΥ) وθ~n=Π~n(ΛΥ) مع كون Πn وΠ~n الإسقاطين على النمطين (fngn) و(f~ng~n) على الترتيب، ويسمى (ΛΥ)=Π,M(ΛΥ) الجزء اللا نهائي البعد، مع كون Π,M المسقط على الفضاء الذاتي الجزئي الموافق لطيف الأدنى من 1M2. ونلاحظ أن التفكيك 3.10 وحيد.

تحضير المعطيات الابتدائية وتعريف المجموعة المتقلصة:

بالنظر إلى A>1 وs0e، نعتبر المعطيات الابتدائية للنظام (3.3) من الشكل

(ΛΥ)A,s0,d0,d1(y)=As02[d0(f0g0)+d1(f1g1)]χ(y,s0), (3.11)

حيث d0 وd1N معلمتان للمسألة. وهدفنا أن نبين أنه من أجل ثابت كبير مثبت A، ثم يثبت s0=s0(A) كبيرا أيضا، توجد (d0,d1)1+N بحيث إن النظام (3.3) ذا المعطيات الابتدائية عند s=s0 المعطاة بـ (3.11) له الحل الوحيد (Λ,Υ) الذي يحقق (3.7). وبصورة أدق، سنبرهن أن الحل (Λ,Υ) ينتمي إلى المجموعة المتقلصة الآتية:

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

لكل A1 وse، نعرف 𝒱A(s) على أنها مجموعة كل (Λ,Υ)L(N)×L(N) بحيث

|θ0(s)|As2,|θ1(s)|As2,|θ2(s)|A4logss2,
|θj(s)|Ajsj+12,|θ~j(s)|Ajsj+12for  3jM,|θ~i(s)|A2s2fori=0,1,2,
Λ(y,s)1+|y|M+1L(N)AM+1sM+22,Υ(y,s)1+|y|M+1L(N)AM+1sM+22,
Λe(s)L(N)AM+2s,Υe(s)L(N)AM+2s,

حيث تعرف Λe,Υe بواسطة (3.9)، وتعرف Λ,Υ، وθn، وθ~n كما في التفكيك (3.10).

ملاحظة 3.2.

يمكننا التحقق من أنه إذا كان (ΛΥ)𝒱A(s) من أجل se، فعندئذ

Λ(s)L()+Υ(s)L()CAM+2s, (3.12)

من أجل ثابت موجب ما C، ومن ثم يثبت التقدير (3.7).

فيما يأتي نتأكد من أن المعطيات الابتدائية (3.11) تنتمي إلى 𝒱A(s0).

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

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

  • (i)

    تنتمي المعطيات الابتدائية (3.11) إلى 𝒱A(s0) مع متراجحات صارمة، باستثناء تقديري θ0(s0) وθ1(s0).

  • (ii)

    التطبيق Θ:𝒟s01+N، المعرف بـ Θ(d0,d1)=(θ0(s0),θ1(s0))، خطي وواحد لواحد من 𝒟s0 إلى [As02,As02]1+N، ويرسل 𝒟s0 إلى ([As02,As02]1+N). إضافة إلى ذلك، فإن درجة Θ على الحد مختلفة عن الصفر.

Proof.

انظر القضية 3.3 في [Ghoul et al.(2018c)Ghoul, Nguyen, and Zaag]. ∎

وجود حل لـ (3.3) محصور في 𝒱A(s):

من الملاحظة 3.2، نهدف إلى إثبات الآتي.

قضية 3.4 (وجود حل لـ (3.3) محصور في 𝒱A(s)).

يوجد A1 بحيث إنه لكل AA1، يوجد s0,1(A) بحيث إنه لكل s0s0,1، يوجد (d0,d1) بحيث إذا كان (ΛΥ) هو حل (3.3) بالمعطيات الابتدائية عند s0 المعطاة بـ (3.11)، فعندئذ (Λ(s)Υ(s))𝒱A(s) لكل ss0.

Proof.

من أجل ثابت مثبت A1 وs0(A)1، نلاحظ من مسألة كوشي المحلية للنظام (1.1)-(1.2) في L(N)×L(N) أن لكل معطيات ابتدائية (3.11)، يملك النظام (3.3) حلا وحيدا يبقى في 𝒱A(s) حتى زمن أعظمي ما s=s(d0,d1). فإذا كان s(d0,d1)=+ من أجل (d0,d1)𝒟s0 ما، اكتمل البرهان. وإلا فإننا نبرهن بالخلف ونفترض أن s(d0,d1)<+ لأي (d0,d1)𝒟s0. وبالاستمرارية وتعريف s، نلاحظ أن الحل عند الزمن s يقع على حد 𝒱A(s). ومن ثم، تكون واحدة على الأقل من المتراجحات في تعريف 𝒱A(s) مساواة. وفي القضية الآتية، نبين أن هذا لا يمكن أن يحدث إلا للمركبتين θ0(s) وθ1(s).

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

افترض أن (ΛΥ) حل لـ (3.3) بالمعطيات الابتدائية عند s=s0 المعطاة بـ (3.11) مع (d0,d1)𝒟s0، وأن (Λ(s)Υ(s))𝒱A(s) لكل s[s0,s1] من أجل بعض s1s0 و(Λ(s1)Υ(s1))𝒱A(s1). عندئذ

  • (i)

    (θ0(s1),θ1(s1))([As12,As12])1+N.

  • (ii)

    يوجد ν0>0 بحيث

    ν(0,ν0),(Λ(s1+ν)Υ(s1+ν))𝒱A(s1+ν).
Proof.

برهان القضية 3.5 نتيجة مباشرة لديناميات النظام (3.3). والفكرة هي إسقاط النظام (3.3) على المركبات المختلفة للتفكيكين (3.10) و(3.9). ولجميع تفاصيل البرهان، انظر القسم 5.2 في [Ghoul et al.(2018c)Ghoul, Nguyen, and Zaag]. ∎

من الجزء (i) من القضية 3.5، نرى أن

(θ0(s),θ1(s))([As2,As2])1+N.

ومن ثم، يمكننا تعريف التدفق المعاد قياسه Θ عند s=s كما يأتي:

Θ:𝒟s0 ([1,1]1+N)
(d0,d1) s2A(θ0,θ1)d0,d1(s),

وهو مستمر بفضل الجزء (ii) من القضية 3.5. ومن جهة أخرى، من القضية 3.3، لدينا متراجحات صارمة للمركبات الأخرى من أجل (d0,d1)𝒟s0. وبتطبيق الجزء (ii) من القضية 3.5، نرى أن (Λ(s)Υ(s)) لا بد أن يغادر 𝒱A(s) عند s=s0، ومن ثم s(d0,d1)=s0. وباستحضار الجزء (ii) من القضية 3.3، فإن درجة Θ على الحد مختلفة عن الصفر. ثم ينتج تناقض من نظرية الدليل. وهذا يثبت أنه لا بد من وجود (d0,d1)𝒟s0 بحيث لكل ss0، (Λ(s)Υ(s))𝒱A(s). وبهذا يكتمل برهان القضية 3.4 وكذلك الجزء (ii) من المبرهنة 1.1. ∎

تكافؤ ملمح الانفجار النهائي:

نعرض الحجة الرئيسية لبرهان الجزء (iii) من المبرهنة 1.1. ولكل x00 مع |x0|1، ندخل لكل (ξ,τ)×[t0(x0)Tt0(x0),1) الدوال المساعدة

g(x0,ξ,τ)=(Tt0(x0))αu(x,t),h(x0,ξ,τ)=(Tt0(x0))βv(x,t),

حيث

x=x0+ξTt0(x0),t=t0(x0)+τ(Tt0(x0)), (3.13)

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

|x0|=K(Tt0(x0))|log(Tt0(x0))|for a fixed constant K1. (3.14)

ومن ثبات النظام (1.1)-(1.2) تحت تغيير القياس، تحقق (g(x0,ξ,τ),h(x0,ξ,τ)) أيضا (1.1)-(1.2). ومن (3.13)، و(3.13) والسلوك التقاربي (1.6)، لدينا

sup|ξ|2|log(Tt0(x0))|1/4|g(x0,ξ,0)Φ(K)|C|log(Tt0(x0))|1/40,

و

sup|ξ|2|log(Tt0(x0))|1/4|h(x0,ξ,0)Ψ(K)|C|log(Tt0(x0))|1/40,

عندما |x0|0. ومن الاستمرارية بالنسبة إلى المعطيات الابتدائية للنظام (1.1)-(1.2)، المقترنة بتموضع مكاني في الكرة B(0,|ξ|<|log(Tt0(x0))|1/4)، يمكننا أن نبين أن

sup|ξ|2|log(Tt0(x0))|1/4,0τ<1|g(x0,ξ,0)g^K(τ)|ϵ(x0)0,

و

sup|ξ|2|log(Tt0(x0))|1/4,0τ<1|h(x0,ξ,0)h^K(τ)|ϵ(x0)0,

عندما x00، حيث

g^K(τ)=Γ(1τ+bK2)α,h^K(τ)=γ(1τ+bK2)β,

هو حل النظام (1.1)-(1.2) بالمعطيات الابتدائية الثابتة (Φ(K),Ψ(K)).

وبأخذ τ1 واستعمال (3.13) نحصل على

u(x0) =limtTu(x,t)=(Tt0(x0))αlimτ1g(x0,0,τ)(Tt0(x0))αg^K(1),
v(x0) =limtTv(x,t)=(Tt0(x0))βlimτ1h(x0,0,τ)(Tt0(x0))βh^K(1),

عندما |x0|0. وباستعمال العلاقة (3.14)، نحصل على

|log(Tt0(x0))|2log|x0|,Tt0(x0)|x0|22K2|log|x0||as|x0|0,

ومن ثم

u(x0)Γ(b|x0|22|log|x0||)α,v(x0)γ(b|x0|22|log|x0||)β,

عندما |x0|0. وهذا يختتم برهان الجزء (iii) من المبرهنة 1.1. ونلاحظ أن الجزء (iii) يعطي مباشرة الانفجار عند نقطة وحيدة، وهو خلاصة الجزء (i). وبهذا يكتمل برهان المبرهنة 1.1. أما برهان المبرهنة 1.4، فنحيل إلى [Ghoul et al.(2018c)Ghoul, Nguyen, and Zaag].

3.2. الحالة الأسية (1.3).

في هذا القسم سنوجز تلك التنويعات على الحجج السابقة المطلوبة لبرهان المبرهنة 1.4. ويمكن العثور على جميع تفاصيل البرهان في [Ghoul et al.(2018d)Ghoul, Nguyen, and Zaag]. والفرق الرئيسي بين الحالتين هو وجود حدي التدرج اللاخطيين |Φ|2Φ و|Ψ|2Ψ في (2.8) بعد إجراء تغيير المتغيرات (2.2). وبالنظر إلى الملمح التقريبي (2.32)، فإن التحكم في هذين الحدين دقيق، ولا سيما عندما يؤول الحل إلى الصفر في المنطقة المتوسطة. ولمعالجتهما، ندخل تحكما دقيقا جدا في الحل داخل مجموعة متقلصة ذات 3 طبقات، معرفة كما يأتي: من أجل K0>0، وϵ0>0 وt[0,T)، نضع

𝒟1(t) ={x||x|K0|ln(Tt)|(Tt)}
{x||y|K0s}{x||z|K0},
𝒟2(t) ={x|K04|ln(Tt)|(Tt)|x|ϵ0}
{x|K04s|y|ϵ0es2}{x|K04|z|ϵ0ses2},
𝒟3(t) ={x||x|ϵ04}{x||y|ϵ04es2}{x||z|ϵ04ses2}.

- في منطقة الانفجار 𝒟1، نخطّي (2.21) حول الملمح التقريبي (2.32)، أي إن (Λ,Υ)=(Φ,Ψ)(φ,ψ) يحل النظام

s(ΛΥ)=(+Mc2+V(y,s))(ΛΥ)+(qp)ΛΥ+(R1R2)+(G1G2), (3.15)

حيث تعرف وMc2 بواسطة (2.12

V(y,s)=(qψ1q(ϕ1/p)p(ψ1/q)pϕ1)=(V1V2V3V4), (3.16)
(G1G2)=(|(Λ+ϕ)|2(Λ+ϕ)1+|ϕ|2ϕ1μ|(Υ+ψ)|2(Υ+ψ)1+μ|ψ|2ψ1), (3.17)

و

(R1R2)=(sϕ+Δϕ12yϕϕ+qϕψ|ϕ|2ϕ1sψ+μΔψ12yψψ+pϕψμ|ψ|2ψ1). (3.18)

والتحليل مشابه للحالة كثيرة الحدود وفق التفكيك (3.10) والتعريف (3.9).

- في المنطقة المتوسطة 𝒟2، نتحكم في (u,v) بإدخال الدوال المساعدة الآتية (u~,v~) المعرفة من أجل x0،

{u~(x,ξ,τ)=1qlnσ(x)+u(x+ξσ(x),t(x)+τσ(x)),v~(x,ξ,τ)=1plnσ(x)+v(x+ξσ(x),t(x)+τσ(x)), (3.19)

حيث يعرف t(x) وحيدا من أجل |x| صغير بما يكفي بواسطة

|x|=K04σ(x)|lnσ(x)|withσ(x)=Tt(x). (3.20)

وبثبات المسألة تحت تغيير القياس، نرى أن (u~,v~) تحقق أيضا النظام (1.1)-(1.3). ونثبت أن (u~,v~) تتصرف من أجل

|ξ|α0|lnσ(x)|andτ[t0t(x)σ(x),1)

من أجل بعض t0<T وα0>0، مثل حل النظام التفاضلي العادي

τu^=epv^,τv^=equ^, (3.21)

الخاضع للمعطيات الابتدائية

u^(0)=1qln[p(1+K02/162(μ+1))],v^(0)=1pln[q(1+K02/162(μ+1))].

ويعطى الحل الصريح بواسطة

u^(τ)=1qln[p(1τ+K02/162(μ+1))],v^(τ)=1pln[q(1τ+K02/162(μ+1))]. (3.22)

والتحليل في 𝒟2 يعطي مباشرة خلاصة الجزء (iii) من المبرهنة 1.4.

- في 𝒟3، نتحكم مباشرة في (u,v) باستعمال سلامة الطرح المحلية في الزمن لمسألة كوشي للنظام (1.1).

تعريف المجموعة المتقلصة الآتي لحصر الحل هو الفرق الجوهري بالمقارنة مع برهان الوجود في الحالة كثيرة الحدود.

تعريف 3.6 (تعريف مجموعة متقلصة).

لكل t0<T، وK0>0، وϵ0>0، وα0>0، وA>0، وδ0>0، وη0>0، وC0>0، ولكل t[t0,T)، نعرف 𝒮(t0,K0,ϵ0,α0,A,δ0,η0,C0,t) بأنها مجموعة كل الدوال (u,v) بحيث

  • (i)

    (التحكم في 𝒟1) (Λ(s)Υ(s))𝒱A(s)، حيث أدخلت 𝒱A(s) في التعريف 3.1.

  • (ii)

    (التحكم في 𝒟2) لكل |x|[K04|ln(Tt)|(Tt),ϵ0]، وτ=τ(x,t)=tt(x)σ(x) و|ξ|α0lnσ(x)،

    |u~(x,ξ,τ)u^(τ)|δ0,|ξu~(x,ξ,τ)|C0|lnσ(x)|,
    |v~(x,ξ,τ)v^(τ)|δ0,|ξv~(x,ξ,τ)|C0|lnσ(x)|,

    حيث تعرف u~,v~، وu^، وv^، وt(x) وσ(x) في (3.19)، و(3.22) و(3.20) على الترتيب.

  • (iii)

    (التحكم في 𝒟3) لكل |x|ϵ04،

    |xiu(x,t)xiu(x,t0)|η0and|xiv(x,t)xiv(x,t0)|η0fori=0,1.
ملاحظة 3.7.

بالمقارنة مع التعريف 3.1، تملك المجموعة المتقلصة 𝒮 تقديرات إضافية في المجالين 𝒟2 و𝒟3. وهذه التقديرات لازمة بصورة حاسمة لتحقيق التحكم في حد التدرج اللاخطي (G1G2) الظاهر في (3.15).

بعد تعريف المجموعة المتقلصة 𝒮 لحصر الحل، نحتاج إلى معطيات ابتدائية مناسبة لـ (3.15) بحيث ينتمي الحل الموافق تدريجيا إلى 𝒮(t) لكل t[t0,T). ولهذه الغاية، نعتبر الدوال الآتية المعتمدة على (N+1) معلمات (d0,d1)1+N:

(qupv)d0,d1(x,t0) =(u^(x)v^(x))(1χ1(x,t0))+{(11)s0+ln[(ϕψ)(y0,s0)]}χ1(x,t0)
+ln{(d0(f0(y0)g0(y0))+d1.(f1(y0)g1(y0)))A2s02χ(16y0,s0)}χ1(x,t0), (3.23)

حيث تعرف s0=ln(Tt0)، وy0=xes02، وϕ وψ بواسطة (2.32)، و(f0g0) و(f1g1) هي الدوال الذاتية المدخلة في اللمّة 2.1، وχ هي دالة القطع المعرفة بواسطة (3.8

χ1(x,t0)=χ0(|x||ln(Tt0)|Tt0)=χ0(y0s0),

و(u^,v^)𝒞(N{0})×𝒞(N{0}) تعرف بواسطة

u^(x)={ln(4(μ+1)|ln|x||p|x|2)for|x|C(a),ln(1+a|x|2)for|x|1,
v^(x)={ln(4(μ+1)|ln|x||q|x|2)for|x|C(a),ln(1+a|x|2)for|x|1.

بعد الحصول على تعريف مناسب للمعطيات الابتدائية والمجموعة المتقلصة، تبقى الخطوة الأخيرة هي أن نبين وجود (d0,d1)1+N بحيث إن للنظام (1.1)-(1.3) ذي المعطيات الابتدائية (3.23) حلا وحيدا (u,v)𝒮(t) لكل t[t0,T). والحجة الرئيسية في هذه الخطوة هي نفسها تماما كما في الحالة كثيرة الحدود، أي برهان القضية 3.4. ونحيل القارئ المهتم بالاختزال إلى مسألة منتهية البعد إلى القسم 4 في [Ghoul et al.(2018d)Ghoul, Nguyen, and Zaag] لجميع التفاصيل. وهذا يختتم برهان المبرهنة 1.4.

References