Time-dimensional reduction and a Carleman contraction principle for an inverse initial data problem for the incompressible Navier-Stokes equations with unknown body force

إليك ملخص تقني مفصل للورقة البحثية بعنوان "الاختزال في البعد الزمني ومبدأ التقلص لـ "كارليمان" لمسألة بيانات أولية عكسية لمعادلات نافيير-ستوكس غير القابلة للانضغاط مع قوة جسم مجهولة" من تأليف فوونج م. نغوين ولوك هـ. نغوين.

1. بيان المسألة

تتناول الورقة مسألة بيانات أولية عكسية (inverse initial data problem) صعبة لمعادلات نافيير-ستوكس غير القابلة للانضغاط في نطاق محدود $\Omega \subset \mathbb{R}^d$ عبر فاصل زمني $(0, T)$.

• المعادلات الحاكمة: يتكون النظام من توازن الزخم وشرط عدم الانضغاط:
  $$\begin{cases} u_t + (u \cdot \nabla)u = -\nabla p + \mu \Delta u + f & \text{في } \Omega \times (0, T), \\ \text{div } u = 0 & \text{في } \Omega \times (0, T), \\ u = 0 & \text{على } \partial\Omega \times (0, T), \\ u(x, 0) = u_0, \quad p(x, 0) = p_0 & \text{في } \Omega. \end{cases}$$
• المجاهيل: الهدف هو استعادة السرعة الأولية $u_0$ والضغط الأولي $p_0$.
• التحدي: قوة الجسم $f(x)$ هي مجهولة وغير ثابتة زمنياً. في معظم المسائل العكسية، تُفترض معاملات النموذج (مثل $f$) معلومة؛ أما هنا، فإن $f$ هي مجهول إضافي يزيد من تعقيد عملية إعادة البناء.
• البيانات المتاحة: المعلومات الوحيدة المتاحة هي بيانات الحدود الجانبية المقاسة على $\partial\Omega \times (0, T)$:
  • المشتقة العمودية للسرعة: $g = \partial_\nu u$.
  • الضغط ومشتقته العمودية: $h = [p, \partial_\nu p]^T$.

2. المنهجية

يقترح المؤلفون إطار عمل جديد يجمع بين الاختزال في البعد الزمني مع مبدأ تقلص "كارليمان" الموزون (Carleman-weighted contraction mapping). تسير المنهجية عبر ثلاث مراحل:

أ. التخلص من القوة المجهولة

بما أن قوة الجسم $f$ غير ثابتة زمنياً، يقوم المؤلفون باشتقاق معادلة الزخم بالنسبة للزمن ($t$). يؤدي هذا إلى إلغاء $f$ تماماً، مما ينتج عنه نظام جديد لـ مشتقة السرعة بالنسبة للزمن ($u_t$) ومشتقة الضغط بالنسبة للزمن ($p_t$):
$$u_{tt} + (u_t \cdot \nabla)u + (u \cdot \nabla)u_t = -\nabla p_t + \mu \Delta u_t.$$
ومن خلال أخذ تباعد (divergence) هذه المعادلة واستخدام $\text{div } u = 0$، يستنتجون نظاماً مقترناً يتضمن فقط $u_t$ و$p_t$، مما يفصل فعلياً قوة الجسم المجهولة عن عملية إعادة البناء.

ب. الاختزال في البعد الزمني باستخدام كثيرات حدود ليجاندر-الأسية

يقوم المؤلفون بتوسيع المجالات المجهولة $u(x,t)$ و$p(x,t)$ باستخدام أساس محدد:
$$u(x,t) = \sum_{n=0}^\infty u_n(x) \Psi_n(t), \quad p(x,t) = \sum_{n=0}^\infty p_n(x) \Psi_n(t),$$
حيث $\Psi_n(t) = e^t Q_n(t)$، و $Q_n$ هي كثيرات حدود ليجاندر المزاحة على الفترة $(0, T)$.

• القطع (Truncation): يتم قطع السلاسل اللانهائية عند مستوى محدد $N$.
• الإسقاط (Projection): يتم إسقاط المعادلات المشتقة زمنياً على دالات الأساس $\Psi_m(t)$ باستخدام حاصل ضرب داخلي موزون $\langle \cdot, \cdot \rangle_{e^{-2t}}$.
• النتيجة: يحول هذا المسألة العكسية المعتمدة على الزمن إلى نظام مقترن من معادلات تفاضلية جزئية إهليلجية (elliptic PDEs) محدودة الأبعاد لمعاملات المتجهات $U_N = [u_0, \dots, u_N]^T$ و $P_N = [p_0, \dots, p_N]^T$. كما يتم إسقاط بيانات الحدود لتوفير بيانات "كوشي" (Cauchy data) لهذه المعاملات.

ج. إطار تقلص "كارليمان"

لحل النظام الإهليجي المقترن الناتج دون الحاجة إلى تخمين أولي جيد (وهي نقطة فشل شائعة لطرق التحسين القياسية)، يستخدم المؤلفون تكرار بيكارد الموزون بـ "كارليمان":

1. المؤثر المنتظم (Regularized Operator): يعرفون مؤثراً $T_{\lambda, \epsilon}$ يربط بين الزوج التجريبي $(V, Q)$ وبين مُقلل دالة التكلفة الموزونة (weighted least-squares functional). يتم "تجميد" اللاخطية باستخدام الزوج التجريبي.
2. المعيار الموزون: يتم إجراء التقليل في فضاء مزود بـ وزن "كارليمان" $e^{2\lambda \mu_\beta(x)}$، حيث $\mu_\beta(x) = |x-x_0|^{-\beta}$.
3. إثبات التقلص: باستخدام تقديرات "كارليمان" (متراجحات تحدد حل المعادلة التفاضلية الجزئية بناءً على مصدرها وبيانات حدودها)، يثبتون أنه بالنسبة لمعامل $\lambda$ كبير بما يكفي، يكون المؤثر $T_{\lambda, \epsilon}$ متقلصاً (contraction) على مجموعة مسموح بها.
4. التقارب العالمي: من خلال مبرهنة النقطة الثابتة لـ "باناخ"، يتقارب التكرار عالمياً نحو نقطة ثابتة وحيدة، وهي تمثل تقريباً لحل النظام الإهليجي المختزل.

3. المساهمات الرئيسية

• التعامل مع قوة الجسم المجهولة: تنجح الورقة في حل مسألة عكسية حيث تكون قوة الجسم مجهولة تماماً، وهو سيناريو يُعتبر غالباً مستعصياً في الصيغ القياسية.
• التقارب العالمي بدون تخمين أولي: على عكس مناهج التحسين التقليدية (مثل طريقة غاوس-نيوتن) التي تكون غير محدبة (non-convex) وحساسة للتخمينات الأولية، يضمن منهج "كارليمان-بيكارد" المقترح التقارب العالمي إلى حل النموذج المختزل.
• تقنية اختزال مبتكرة: إن الجمع بين الاشتقاق الزمني لإلغاء $f$ واختزال الأساس الأسسي-الليجاندر يحول مسألة عكسية معلمية (parabolic) معقدة إلى نظام إهليجي قابل للحل.
• التحليل الصارم: يقدم المؤلفون برهاناً نظرياً كاملاً لخاصية التقلص للمؤثر واتساق النقطة الثابتة مع النظام المختزل مع اقتراب معامل الانتظام $\epsilon \to 0$.

4. النتائج العددية

قام المؤلفون بتطبيق الخوارزمية 1 (منهج كارليمان-بيكارد) واختبارها على بيانات اصطناعية ناتجة عن معادلات نافيير-ستوكس الأمامية في نطاق مربع ثنائي الأبعاد ($\Omega = (-1, 1)^2$).

• الإعداد:
  • ثلاث حالات اختبار متميزة بقوى جسم مختلفة وحقول سرعة وضغط أولية مختلفة.
  • تم توليد بيانات الحدود وتلويثها بـ ضوضاء عشوائية بنسبة 10%.
  • المعاملات: مستوى القطع $N=35$، ومعاملات وزن "كارليمان" $\beta=20, \lambda=6$.
• الأداء:
  • الدقة: نجحت الطريقة في استعادة الميزات النوعية الرئيسية (المواقع، الإشارات، والأشكال) للسرعة والضغط الأوليين.
  • مقاييس الخطأ:
    • الاختبار 1: أخطاء $L^2$ النسبية $\approx 14.5\%$ (السرعة) و $23.9\%$ (الضغط).
    • الاختبار 2: أخطاء $L^2$ النسبية $\approx 17.6\%$ (السرعة) و $18.7\%$ (الضغط).
    • الاختبار 3 (الأكثر تعقيداً): ارتفعت الأخطاء إلى $\approx 36\%$ (السرعة) و $42.8\%$ (الضغط)، لكن الطريقة لا تزال تلتقط الهياكل الأساسية.
  • التقارب: أظهر تكرار "بيكارد" تقارباً سريعاً ومستقراً، مع تناقص مستمر في البواقي النسبية وتغيرات التكرار. ظلت الطريقة مستقرة حتى مع وجود بيانات مشوشة.

5. الأهمية والتوجهات المستقبلية

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

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