Formation of Implosion Singularities in 3D Compressible Navier-Stokes-Korteweg Equation

ملخص تقني: تشكل مفردات الانفجار الداخلي في معادلات نافييه-ستوكس-كورتويج القابلة للانضغاط ثلاثية الأبعاد

بيان المشكلة
تبحث هذه الورقة في تشكل المفردة في زمن محدد (finite-time singularity) لنظام نافييه-ستوكس-كورتويج (NSK) القابل للانضغاط ثلاثي الأبعاد، مع معاملات لزوجة وتوتر سطحي تعتمد على الكثافة. يُعرف النظام بـ:
$$\begin{cases} \rho_t + \text{div}(\rho u) = 0, \\ (\rho u)_t + \text{div}(\rho u \otimes u) + \nabla P = \text{div}(2\mu(\rho)D(u)) + \nabla(\lambda(\rho)\text{div} u) + \text{div} K, \end{cases}$$
حيث تحقق المعاملات علاقات قانون القوة المحددة:
$$\mu(\rho) = \nu\rho^\alpha, \quad \lambda(\rho) = 2\nu(\alpha - 1)\rho^\alpha, \quad \kappa(\rho) = \varepsilon^2\alpha^2\rho^{2\alpha-3}.$$
أثبتت أعمال سابقة بواسطة Gu–Huang–Meng–Zhou [19] و Huang–Lei–Zhou [21] وجود حلول قوية عالمية لأي بيانات أولية كبيرة (بعيداً عن الفراغ) عندما يقع الأس $\alpha$ ضمن نطاق مناسب محدد (على سبيل المثال، $\alpha \in (\frac{9\sqrt{3}-4\sqrt{2}}{9\sqrt{3}-2\sqrt{2}}, 1)$ في الحالة ثلاثية الأبعاد). السؤال المركزي الذي تعالجه هذه الورقة هو ما إذا كانت الحلول ستوجد عالمياً أم ستنفجر في زمن محدد عندما يقع $\alpha$ خارج هذا النطاق، وتحديداً بالنسبة للأسس الموجبة الصغيرة ($\alpha < 1/2$). في هذا النظام، تصبح اللزوجة الحجمية الفعالة سالبة ($\lambda(\rho) + \frac{2}{3}\mu(\rho) < 0$)، وهي سمة مرتبطة بالأوساط الصوتية غير المتوازنة والمادة النشطة، مما يتحدى نظريات الوجود العالمي القياسية.

المنهجية
يستخدم المؤلفون نهجاً اضطرابياً يعتمد على بناء بروفيلات انفجار داخلي (imploding profiles) منتظمة ذاتية التشابه مشتقة من معادلات أويلر القابلة للانضغاط. تمر المنهجية عبر عدة خطوات رئيسية:

1. إعادة صياغة السرعة الفعالة: للتعامل مع موتر كورتويج من الدرجة الثالثة، تمت إعادة صياغة النظام باستخدام سرعة فعالة $v = u + d\rho^{\alpha-2}\nabla\rho$ (حيث $d = \nu - \sqrt{\nu^2-\varepsilon^2}$). هذا التحويل يعيد صياغة النظام إلى بنية بارابولية-بارابولية لـ $(\rho, v)$، مما يكشف عن حد انتشار متدهور خفي $d\Delta(\rho^\alpha)$ في معادلة الكتلة.
2. القياس الذاتي المتشابه (Self-Similar Scaling): تم تحويل المشكلة إلى إحداثيات ذاتية التشابه $(\tau, y)$ حيث $\tau = -\log(T-t)/\Lambda$ و $y = x/(T-t)^{1/\Lambda}$. في هذه الإحداثيات، تظهر التأثيرات اللزجة والتوتر السطحي كاضطرابات تتلاشى أسياً ($e^{-f_1\tau}$) بشرط أن يحقق معامل القياس $\Lambda$ شرطاً محدداً ($f_1 > 0$).
3. تحليل استقرار الاضطرابات: يتمثل جوهر البرهان في تحليل الاستقرار غير الخطي للنظام حول بروفيل أويلر سلس، متماثل كروياً، وذاتي التشابه $(\bar{S}, \bar{V})$. يتم اشتقاق معادلات الاضطراب واستخدام حجة "بوتستراب" (bootstrap argument) للتحكم في الحل.
4. تقديرات الطاقة والتحلل الطيفي:
   • تقديرات الطاقة الموزونة: تم وضع تقديرات طاقة موزونة عالية الرتبة للتحكم في تلاشي الاضطرابات. استخدم المؤلفون دالة وزن محددة $\phi$ وهيكلية من المعاملات لإدارة التفاعل بين النقل والتبدد.
   • حجة البوتستراب: تم تصميم مخطط بوتستراب صارم لإغلاق التقديرات للمشتقات من الرتب الدنيا (تحديداً الحد الأدنى للكثافة) والمشتقات من الرتب العليا. الابتكار الحاسم هو استخدام كل من بنية النقل والتبدد لوضع حد أدنى للكثافة، للتغلب على الصعوبة المتمثلة في أن معادلة الكثافة تمتلك بنية تبدد بالإضافة إلى بنية النقل.
   • التحكم في النمط غير المستقر: يتم تحليل المؤثر الخطي إلى فضاءات مستقرة وغير مستقرة. يتم التحكم في المكونات غير المستقرة، والتي تقابل فضاءً محدود الأبعاد، من خلال اختيار معاملات البيانات الأولية المحددة. يضمن هذا الاختيار بقاء الحل على "متشعب" (manifold) ذي كوديم world-dimension محدد يؤدي إلى الانفجار.

المساهمات والنتائج الرئيسية
تثبت الورقة وجود بيانات أولية سلسة، كثافتها منفصلة بانتظام عن الفراغ، تتطور إلى مفردات انفجار داخلي في زمن محدد لنظام NSK ثلاثي الأبعاد عندما يكون $\alpha$ صغيراً بما يكفي ($\alpha < 1/2$).

• النظرية الرئيسية (Theorem 1.1): للمعاملات التي تحقق إما:
  • (P1) $1 < \gamma < 1 + 2/\sqrt{3}$، $0 < \alpha < 1/2$، و $f_1 > 0$؛ أو
  • (P2) $\gamma \in (1, 1 + 2/\sqrt{3}) \setminus \mathcal{J}$ (حيث $\mathcal{J}$ هي مجموعة استثنائية معدودة) و $0 < \alpha < \alpha^*(\gamma)$؛
    توجد بيانات أولية $C^\infty$ $(\rho_0, u_0)$ بحيث يتطور الحل المقابل $(\rho, u)$ في الفترة $[0, T)$ ويطور مفردة انفجار داخلي عند $t=T$.
• خصائص المفردة: عند زمن الانفجار $T$:
  • تصبح الكثافة لانهائية عند الأصل: $\lim_{t \to T^-} \rho(t, 0) = +\infty$.
  • تكون السرعة الفعالة $v$ غير محدودة في كل جوار للأصل.
  • يتقارب الحل، بعد إعادة القياس، نحو بروفيل أويلر المتشابه ذاتياً المحدد مسبقاً.
• النتيجة (Corollary 1.1): يطور النظام الأصلي مفردة في زمن محدد، وتحديداً تتميز بانفجار الكثافة ($\lim_{t \to T^-} \rho(t, 0) = +\infty$). ومع ذلك، لا يمكن استنتاج انفجار السرعة الفيزيائية $u$ مباشرة من انفجار السرعة الفعالة، حيث قد يصبح التدرج $\nabla\rho$ أيضاً مفرداً بالقرب من الأصل، ولا يمكن استبعاد حدوث إلغاء محتمل في العلاقة $v = u + d\rho^{\alpha-2}\nabla\rho$.

الأهمية
تدعي هذه الورقة أنها تكمل نظرية الوجود العالمي الحالية لنظام NSK. بينما أظهرت النتائج السابقة وجود حلول عالمية لنطاق محدد من $\alpha$، تحدد هذه الورقة آلية انفجار داخل في زمن محدد لنطاق "$\alpha$ الصغير". يتميز هذا النطاق ببنية لزوجة حجمية فعالة سالبة، وهي ظاهرة ذات صلة بالفيزياء غير المتوازنة. تُظهر النتائج أنه حتى مع البيانات الأولية السلسة البعيدة عن الفراغ، فإن التفاعل بين اللزوجة المعتمدة على الكثافة، والتوتر السطحي، واللزوجة الحجمية السالبة يمكن أن يؤدي إلى مفردات انفجار داخلي حيث تنفجر الكثافة وحقول السرعة الفعالة نفسها، وليس فقط مشتقاتها (تشكل الصدمات). يعتمد البناء على استقرار بروفيلات أويلر ذاتية التشابه، مما يوسع آلية الانفجار الداخلي المعروفة لمعادلات أويلر لتشمل إطار NSK اللزج والمتوتر سطحياً.