Stable $C^1$-conforming finite element methods for a class of nonlinear fourth-order evolution equations

ملخص تقني: طرق العناصر المحدودة المستقرة والمتوافقة مع درجة الاستمرارية $C^1$ لفئة من المعادلات التطورية غير الخطية من الرتبة الرابعة

صياغة المسألة
يتناول البحث التحليل العددي لفئة عامة من معادلات التفاضل الجزئي التطورية غير الخطية من الرتبة الرابعة، والمعرفة على نطاق محدود $\Omega \subset \mathbb{R}^d$ ($d \in \{1, 2, 3\}$). تُعطى المعادلة الحاكمة بما يلي:
$$\frac{\partial u}{\partial t} - \beta_1 \Delta u + \beta_2 \Delta^2 u = \beta_3(1 - |u|^2)u - \beta_4 u \times \Delta u + \beta_5 \Delta(|u|^2 u) + \beta_6 (j \cdot \nabla)u$$
خاضعة لشروط نيومان الحدية ($\frac{\partial u}{\partial n} = 0, \frac{\partial (\Delta u)}{\partial n} = 0$) وشرط ابتدائي $u(0) = u_0$. حيث يمكن أن تكون $u$ ذات قيمة قياسية ($m=1$) أو ذات قيمة متجهة ($m=3$). وتتضمن المعادلة نماذج فيزيائية مختلفة اعتماداً على المعاملات $\beta_i$ والبعد $d$، بما في ذلك:

• المغناطيسية الدقيقة (Micromagnetics): معادلة لاندو-ليفيتز-باريياختر (LLB) (مع أو بدون عزم دوران سبيني) وتعديلات معادلة لاندو-ليفز-بلوخ.
• تكون الأنماط (Pattern Formation): معادلة سويفت-هوهنبرج.
• انفصال الطور والانتشار (Phase Separation and Diffusion): معادلات كاهن-هيليارد الحملية مع حدود مصدرية، ونماذج الانتشار غير الطبيعية ذات التدفق الثنائي، ونماذج انتشار الجمهرة المعممة.

يشير المؤلفون إلى أنه بينما توجد طرق عددية لنماذج محددة من هذه المعادلات (مثل معادلة سويفت-هوهنبرج القياسية أو معادلة كاهن-هيليارد القياسية)، إلا أن هناك نقصاً في الأدبيات العلمية لطريقة عددية موحدة وعامة للنظام الكامل (1.1) مع احتمال كون جميع المعاملات غير صفرية.

المنهجية
يقترح المؤلفون ويحللون طرق عناصر محدودة منفصلة زمنياً ومكانياً بالكامل تستخدم عناصر متوافقة مع درجة الاستمرارية $C^1$. وقد استلزم هذا الاختيار وجود رتبة رابعة من المشتقات المكانية ($\Delta^2 u$)، مما يتطلب أن يكون فضاء العناصر المحدودة فضاءً جزئياً من $H^2(\Omega)$.

• التقطيع المكاني: تستخدم الطريقة كثيرات حدود مجزأة من الدرجة $r$ (حيث $r=3$ للأبعاد $d=1,2$ و $r=9$ للبعد $d=3$ للحفاظ على استمرارية $C^1$). ومن العناصر المحددة المذكورة: عناصر هيرميت $C^1$-P3 (في البعد 1)، وعناصر ماكرو HCT (في البعد 2)، وعناصر Zenísek (في البعد 3).
• التقطيع الزمني: تم اقتراح مخططي خطية للخطوات الزمنية:
  1. طريقة أويلر شبه الضمنية (من الدرجة الأولى في الزمن).
  2. طريقة BDF2 شبه الضمنية (من الدرجة الثانية في الزمن).
     في كلا المخططين، يتم خطية الحدود غير الخطية عن طريق تقييمها عند الخطوة الزمنية السابقة، مما يؤدي إلى نظام من المعادلات الخطية في كل خطوة. وهذا يتجنب التكلفة الحسابية لحل الأنظمة غير الخطية عبر طريقة نيوتن.

المساهمات الرئيسية والنتائج النظرية
يقدم البحث تحليلاً دقيقاً للخطأ وإثباتاً للاستقرار للمخططات المقترحة:

1. الاستقرار:

   • ثبت أن المخططات مستقرة غير مشروطة في معيار $H^2$ للأبعاد $d=1$ و $d=2$. وهذا يعني أن الاستقرار يتحقق دون قيود على حجم الخطوة الزمنية $k$ بالنسبة لحجم الشبكة $h$، ودون افتراض وجود تقسيمات شبه منتظمة.
   • بالنسبة للبعد $d=3$، فإن المخططات مستقرة مشروطة. حيث يتحقق الاستقرار إما إذا كانت البيانات الأولية صغيرة بما يكفي أو إذا كانت الخطوة الزمنية تحقق قيداً معيناً بالنسبة للطاقة الأولية. ويقر المؤلفون بأن العناصر المتوافقة مع $C^1$ في البعد 3 (التي تتطلب كثيرات حدود من الدرجة 9) مكلفة حسابياً، مما يشير إلى أن الطرق غير المتوافقة قد تكون أكثر عملية للبعد $d=3$، وإن كان التحليل النظري قد أُدرج لأغراض الاكتمال.

2. تقديرات الخطأ والتقارب:

   • يحدد المؤلفون معدلات تقارب مثالية في معايير $L^\infty(0, T; H^\beta)$ حيث $\beta = 0, 1, 2$.
   • بالنسبة للمخطط شبه المنفصل مكانياً، يكون الخطأ مقيداً بـ $O(h^{r+1-\beta})$.
   • بالنسبة للمخططات المنفصلة تماماً:
     • تحقق طريقة أويلر شبه الضمنية خطأً قدره $O(h^{r+1-\beta} + k)$.
     • تحقق طريقة BDF2 شبه الضمنية خطأً قدره $O(h^{r+1-\beta} + k^2)$.
   • يعتمد التحليل على تطوير إسقاط إهليلجي متخصص متكيف مع المسألة غير الخطية وتقدير دقيق للحدود غير الخطية باستخدام تضمينات سوبوليف وتدابير العكس.

3. التجارب العددية:

   • تم التحقق من معدلات التقارب النظرية باستخدام الحزمة مفتوحة المصدر FreeFEM مع عناصر HCT في البعد 2.
   • أُجريت التجارب على معادلة لاندو-ليفيتز-باريياختر (متجهة) ومعادلة كاهن-هيليارد مع حد مصدري (قياسية).
   • أكدت النتائج رتب التقارب المتوقعة: الرتبة الرابعة في $L^2$، والرتبة الثالثة في $H^1$، والرتبة الثانية في $H^2$ للتقطيع المكاني، وهو ما يتوافق مع استخدام عناصر الدرجة الثالثة ($r=3$).

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

• الخطية: من خلال جعل الحدود غير الخطية خطية، تتجنب الطريقة تعقيد وعدم استقرار حل الأنظمة غير الخطية في كل خطوة زمنية، مما يجعلها مفيدة لتطبيقات مثل محاكاة المغناطيسية الدقيقة.
• الاستقرار: يضمن الاستقرار غير المشروط في $H^2$ للأبعاد $d=1,2$ أداءً قوياً دون قيود صارمة على الخطوة الزمنية.
• الدقة: يسمح استخدام العناصر المتوافقة مع $C^1$ بدقة أعلى في التقاط تدرجات الحل مقارنة بالطرق غير المتوافقة، رغم التكلفة الحسابية الإضافية لاستمرارية $C^1$.

يحافظ المؤلفون على نبرة متواضعة فيما يتعلق بالبعد 3، حيث ذكروا صراحة أنه بينما التحليل النظري كامل، فإن التطبيق العملي للعناصر المتوافقة مع $C^1$ في البعد 3 يعتبر حالياً غير ممكن حسابياً بسبب درجة كثيرات الحدود العالية المطلوبة. ويعمل هذا العمل كإطار نظري تأسيسي لحل المعادلات التطورية المعقدة من الرتبة الرابعة الناشئة في الفيزياء، وعلوم المواد، والبيولوجيا.