Lipschitz Stability for an Inverse Problem of Biharmonic Wave Equations with Damping

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

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

تتناول الورقة مسألة عكسية متزامنة لـ معادلة موجة ثنائية هاملتون (biharmonic wave equation) مخمدة ذات كثافة متغيرة على نطاق محدود $\Omega \subset \mathbb{R}^n$. النظام الأمامي الحاكم هو:
$$\begin{cases} \rho(x)\partial^2_t u + \Delta^2 u + \gamma \partial_t u = 0, & (t, x) \in \mathbb{R}^+ \times \Omega, \\ u(0, x) = f(x), \quad \partial_t u(0, x) = g(x), & x \in \Omega, \\ u = 0, \quad \frac{\partial u}{\partial n} = 0, & (t, x) \in \mathbb{R}^+ \times \partial\Omega. \end{cases}$$
الهدف: استعادة مجهولين في آن واحد من قياسات حدودية:

1. معامل الكثافة المتغير $\rho(x)$.
2. الإزاحة الأولية $f(x)$.

البيانات: القياسات المتاحة هي بيانات كوشي الحدودية لـ "لابلاسيان" الحل خلال فترة زمنية محدودة $[0, T]$:
$$\Delta u|_{\partial\Omega} \quad \text{و} \quad \frac{\partial (\Delta u)}{\partial n}\bigg|_{\partial\Omega}.$$

الافتراضات:

• الكثافة $\rho(x)$ ثابتة جزئياً ومحدودة ($0 < \rho_{\min} \le \rho \le \rho_{\max}$).
• السرعة الأولية $g$ معلومة (أو متطابقة للحلين اللذين تتم مقارنتهما).
• يُفرض شرط عدم التدهور (non-degeneracy condition) على تسارع حل مرجعي ($\partial^2_t u_2$)، مما يضمن أن اللوح ليس ساكناً ($\|\partial^2_t u_2\|_{L^2} \ge \alpha > 0$).
• النطاق $\Omega$ متمركز حول نقطة (star-shaped) بالنسبة لنقطة $x_0$.
• زمن المراقبة $T$ يستوفي شرط التحكم الهندسي: $T > \frac{2 \text{diam}(\Omega)}{\sqrt{\rho_{\min}}}$.

2. المنهجية

يستخدم المؤلفون نهجاً تحليلياً وظيفياً صارماً مدمجاً مع تقنيات المضروبات (multiplier techniques). تسير المنهجية عبر ثلاث مراحل رئيسية:

أ. جودة التحديد (Well-Posedness) للمسألة الأمامية

لضمان أن المسألة العكسية محددة جيداً، يثبت المؤلفون أولاً جودة التحديد للنظام الأمامي.

• نظرية شبه المجموعة (Semigroup Theory): يعرّفون فضاء الطاقة $\mathcal{H} = (H^4(\Omega) \cap H^2_0(\Omega)) \times H^2_0(\Omega)$ ومعامل نظام $A$.
• مبرهنة هيل-يوسيدا (Hille-Yosida Theorem): يثبتون أن $A$ هو مبدد أقصى (maximally dissipative) ومعرّف بكثافة. وهذا يضمن أن $A$ يولد شبه مجموعة (semigroup) من نوع $C_0$ من التقلصات، مما يضمن وجود الحل، وفرادته، والاعتماد المستمر للحلول على البيانات الأولية.

ب. متباينة القابلية للملاحظة (Observability Inequality)

خطوة حاسمة هي اشتقاق متباينة قابلية الملاحظة التي تربط طاقة النظام الداخلية بالقياسات الحدودية.

• طريقة المضروب (Multiplier Method): باستخدام المجال المتجه $m(x) = x - x_0$ (بالاستفادة من الهندسة المتمركزة)، يقومون بضرب معادلة التفاضل الجزئي في $m \cdot \nabla u$ وإجراء التكامل بالتجزئة.
• النتيجة: يشتقون متباينة تظهر أن الطاقة الأولية $E(0)$ محكومة بالتكاملات الحدودية لـ $\Delta u$ و $\frac{\partial \Delta u}{\partial n}$.
• الاعتماد على التخميد: يعد تتبع معامل التخميد $\gamma$ إنجازاً تقنياً رئيسياً. تأخذ المتباينة الشكل التالي:
  $$E(0) \le C(1+\gamma) \int_0^T \int_{\partial\Omega} \left( \left|\frac{\partial \Delta u}{\partial n}\right|^2 + |\Delta u|^2 \right) dS dt.$$

ج. تحليل الاستقرار عبر نظام الفرق

لإثبات استقرار ليفشيتز، يحلل المؤلفون الفرق $u = u_1 - u_2$ بين حلين يتوافقان مع معلمات مختلفة $(\rho_1, f_1)$ و $(\rho_2, f_2)$.

• تحديد الحد المصدر: تصبح معادلة الفرق نظاماً غير متجانس حيث يكون الحد المصدر هو $-\rho \partial^2_t u_2$.
• تقديرات الطاقة: يشتقون تقديرات الطاقة لحل الفرق، مستفيدين من عدم تدهور $\partial^2_t u_2$ لحصر فرق الكثافة $\|\rho_1 - \rho_2\|_{L^\infty}$ من خلال تكامل الحد المصدر.
• الاقتران (Coupling): من خلال دمج تقديرات الطاقة مع متباينة القابلية للملاحظة، يربطون فروق المعلمات المجهولة مباشرة بالبيانات الحدودية.

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

1. الاستعادة المتزامنة لمتعدد المعلمات: بخلاف العديد من الأعمال الموجودة التي تركز على استعادة معلمة واحدة، توفر هذه الورقة إطاراً لـ الاستعادة المتزامنة لكل من الكثافة المتغيرة مكانياً $\rho$ والإزاحة الأولية $f$.
2. استقرار ليفشيتز صريح: تثبت الورقة تقديرات استقرار ليفشيتز (وهي أقوى من الاستقرار اللوغاريتمي) لكلا المعلمتين. وهذا يعني أن الأخطاء الصغيرة في القياسات الحدودية تؤدي إلى أخطاء متناسبة وصغيرة في المعلمات المستعادة.
3. الاعتماد الصريح على التخميد: تم تحديد ثوابت الاستقرار صراحةً بدلالة معامل التخميد $\gamma$ عبر العامل $(1+\gamma)^{1/2}$. وهذا يوفر رؤية حول كيفية تأثير التخميد على ضبط المسألة العكسية.
4. تحليل المشغل من الرتب العليا: تتعامل الورقة بصرامة مع تعقيدات المشغل ثنائي هاملتون ($\Delta^2$)، وهو أصعب بكثير في التحليل من مؤثر لابلاس القياسي بسبب متطلبات الانتظام الأعلى والشروط الحدودية.

4. النتائج الرئيسية

تثبت الورقة مبرهنتين رئيستين:

• المبرهنة 1.1 (استقرار الكثافة):
  تحت شرط عدم التدهور على التسارع، يتم حصر معيار $L^\infty$ لفرق الكثافة بواسطة البيانات الحدودية:
  $$\|\rho_1 - \rho_2\|_{L^\infty(\Omega)} \le C (1+\gamma)^{1/2} \left( \int_0^T \int_{\partial\Omega} \left( \left|\frac{\partial \Delta u}{\partial n}\right|^2 + |\Delta u|^2 \right) dS dt \right)^{1/2}.$$

• المبرهنة 1.2 (استقرار الإزاحة الأولية):
  بافتراض أن السرعات الأولية متطابقة ($g_1 = g_2$) وأن زمن المراقبة $T$ كبير بما يكفي، يتم حصر معيار $H^4$ لفرق الإزاحة الأولية:
  $$\|f_1 - f_2\|_{H^4(\Omega)} \le C (1+\gamma)^{1/2} \left( \int_0^T \int_{\partial\Omega} \left( \left|\frac{\partial \Delta u}{\partial n}\right|^2 + |\Delta u|^2 \right) dS dt \right)^{1/2}.$$

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

• التقدم النظري: تسد هذه الورقة فجوة حرجة في الأدبيات المتعلقة بـ معادلات الموجة ذات الكثافة المتغيرة والمخمدة وذات الرتبة الرابعة. فهي تتجاوز النماذج المبسطة لمعالجة المسألة العكسية المقترنة بضمانات استقرار صارمة.
• الاختبار غير التدميري (NDT): توفر النتائج أساساً نظرياً للكشف عن العيوب الداخلية (التي تُنمذج كتغيرات في الكثافة) وتقييم حالات الإجهاد الأولية في الألواح المرنة الرقيقة (مثل الهياكل في مجال الطيران أو الهندسة المدنية) باستخدام الاهتزازات الحدودية فقط.
• القلب الديناميكي (Dynamic Inversion): يشير الاعتماد الصريح على $\gamma$ إلى أن وجود التخميد، الذي يُنظر إليه غالباً كعامل معقد، يمكن التحكم فيه رياضياً، مما يقدم آفاقاً جديدة لاستراتيجيات القلب في الأوساط المبددة للطاقة.
• المتانة (Robustness): تضمن طبيعة "ليفشيتز" لتقديرات الاستقرار أن خوارزميات الاستعادة المستمدة من هذه النظرية ستكون متينة ضد ضجيج القياس، وهو شرط مسبق للتطبيق العددي العملي.