Error analysis of a divergence-preserving mixed finite element scheme for the incompressible Hall--magnetohydrodynamic equations

ملخص تقني: تحليل الخطأ لمخطط عناصر منتهية مختلطة حافظ للتباعد لنظام معادلات هول-MHD غير القابلة للانضغاط

بيان المشكلة
يعد نظام هول-MHD (Hall–MHD) غير القابل للانضغاط نموذجاً للتفاعل بين حركة المائع والمجالات المغناطيسية في الموائع الموصلة كهربائياً، مع دمج تأثير "هول" الذي يصبح جوهرياً عند مقاييس تقارب سمك طبقة غلاف الأيونات (ion skin depth). يُدخل هذا التأثير حداً شديد الصلابة وعالي اللاخطية ($\eta \text{curl}(\text{curl} B \times B)$) في معادلة الحث، مما يؤدي إلى سلوك تشتتي، وتدرجات حادة، وتذبذبات عالية التردد (موجات ويسلر/whistler waves). هذه السمات تفرض تحديات تحليلية وحسابية جسيمة؛ فوجود الحلول القوية عالمياً في الأبعاد الثلاثة لا يزال مسألة مفتوحة، كما تعاني الطرق العددية القياسية غالباً من عدم الاستقرار والاشتراط الصارم بأن يظل المجال المغناطيسي خالي التباعد ($\text{div } B = 0$). إن الانتهاكات لشرط خلو التباعد هذا يمكن أن تؤدي إلى ظهور أحادي أقطاب مغناطيسية غير فيزيائية وتسبب عدم استقرار عددي شديد. علاوة على ذلك، فإن النهج العددي الحالي، مثل تلك القائمة على طرق المتغير المساعد القياسي (SAV) أو التقنيات الطيفية، غالباً ما تشتت دوال الطاقة المعدلة بدلاً من طاقة النظام الفيزيائية الحقيقية، أو تفشل في فرض شرط خلو التباعد بدقة عند المستوى المنفصل، أو تفتقر إلى تحليل الخطأ الصارم والمرونة الهندسية.

المنهجية
لمعالجة هذه التحديات، يقترح المؤلفون طريقة عناصر منتهية مختلطة (FEM) خطية، كاملة التجزئة، تحافظ على البنية لنظام Hall–MHD مُنظم بـ "فويت" (Voigt-regularised). يقوم تنظيم "فويت" باستبدال المشتق الزمني في معادلة الزخم بنسخة مفلترة مكانياً $(I - \alpha_1 \Delta)\partial_t u$ ويحتفظ بحد عطالة الإلكترون $\alpha_2 \partial_t J$ في قانون أوم المعمم. يعمل هذا التعديل على تحسين جودة الحل (well-posedness) والاستقرار العددي مع الحفاظ على البنية الفيزيائية الأساسية وقوانين الطاقة للنظام الأصلي.

يتم بناء المخطط العددي باستخدام المكونات التالية:

1. الصياغة المختلطة: يتم إعادة صياغة المشكلة باستخدام المجال الكهربائي $E$ وكثافة التيار $J$ كمتغيرات أساسية، مما يؤدي إلى نظام مقترن من معادلات الزخم، والحث، وقانون أوم، وقانون أمبير.
2. التجزئة المكانية: يستخدم المخطط حساب العناصر المنتهية الخارجية (FEEC) لضمان استيفاء شرط خلو التباعد بدقة.
   • يتم تقريب السرعة $u$ والضغط $p$ باستخدام زوج العناصر المنتهية MINI (الأكثر استقراراً وفق شرط inf-sup من الدرجة الأدنى).
   • يتم تقريب المجال المغناطيسي $B$ وكثافة التيار $J$ باستخدام فضاءات Raviart–Thomas ($RT_0$) وNédélec ($X_0$) الدنيا، على التوالي.
   • تشكل هذه الفضاءات مجمع "دي رام" (de Rham complex) المنفصل، مما يضمن الحفاظ على $\text{div } B_h = 0$ بدقة عند المستوى المنفصل.
3. التجزئة الزمنية: يُستخدم مخطط خطوة زمنية شبه ضمني (semi-implicitly) متماثل مائل (skew-symmetric) للتعامل مع حدود الحمل واللورنتز بشكل شبه ضمني للحفاظ على الخطية والاستقرار الطاقي غير المشروط.
4. تحليل الخطأ: يضع المؤلفون معدلات تقارب مثلى من خلال تفكيك الخطأ إلى أخطاء التقريب (باستخدام إسقاطات ستوكس والإسقاطات المتعامدة) والأخطاء المنفصلة. ويستخدمون متباينات "جافني" (Gaffney) المنفصلة وتقديرات الطاقة للتعامل مع حدود الاقتران غير الخطية.

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

• مخطط حافظ للبنية: يقدم البحث طريقة عناصر منتهية خطية وكاملة التجزئة تفرض قيد خلو التباعد على المجال المغناطيسي بدقة، مما يتجنب الحاجة إلى خطوات الإسقاط أو طرق الجزاء (penalty methods) التي قد تضعف الدقة أو الاستقرار.
• استقرار طاقي غير مشروط: يحقق المخطط المقترح قانون تبدد طاقة منفصل غير مشروط، محاكياً قانون الطاقة الفيزيائي المستمر، حتى في وجود حد "هول" شديد الصلابة.
• تقديرات خطأ صارمة: يقدم المؤلفون تحليلاً كاملاً للخطأ لنظام "فويت" المنظم، موضحين معدلات تقارب مثلى من رتبة $O(h^s + \tau)$ تحت افتراضات انتظام مناسبة.
• التمديد إلى النظام غير المنظم: من الأهمية بمكان أن التحليل يمتد إلى نظام Hall–MHD غير المنظم (حيث $\alpha_1 = \alpha_2 = 0$). يثبت المؤلفون أنه يمكن أيضاً إثبات التقارب للنظام غير المنظم، بشرط أن تكون المقاومية $\sigma > 0$. وفي حالة عدم وجود مقاومية ($\sigma = 0$)، يظهر أن تنظيم "فويت" يلعب دوراً حاسماً في تقديرات الخطأ.
• التطبيق في أبعاد 2.5D و3D: تم التحقق من صحة الإطار في كل من بيئات 2.5D (الثبات الانتقالي) و3D الكاملة، مما يثبت قابلية تطبيق الطريقة على الهندسات المعقدة والظواهر الفيزيائية.

النتائج
تؤكد المحاكاة العددية في أبعاد 2.5D و3D النتائج النظرية:

• التقارب: يحقق المخطط معدلات التقارب المكاني والزمني المتوقعة للسرعة، والمجال المغناطيسي، وكثافة التيار. وبينما تنخفض معدلات التقارب للسرعة وكثافة التيار قليلاً في بعض الحالات (ويُعزى ذلك إلى محدودية انتظام الحل بسبب حد هول)، يحافظ المجال المغناطيسي على تقارب مثالي.
• الحفاظ على التباعد: تظل أقصى قيمة لخطأ التباعد ($\max | \text{div } B_h |$) عند مستوى سماحية الحل (solver tolerance) طوال عمليات المحاكاة، مما يؤكد الحفاظ الدقيق على شرط خلو التباعد.
• الاستقرار الطاقي: تضمحل الطاقة الكلية للنظام مونوتونياً (رتيباً) مع مرور الوقت، بما يتوافق مع قانون تبدد الطاقة النظري.
• الدقة الفيزيائية: نجحت محاكاة تدفق ABC (Arnold–Beltrami–Childress)، ودوامة Orszag–Tang، وإعادة الاتصال المدفوعة بهول في ورقة Harris (Harris sheet) في التقاط الظواهر الفيزيائية المعقدة، بما في ذلك إعادة الاتصال المغناطيسي، وتشكل الجزر المغناطيسية، وحبال التدفق، وأنماط التيار رباعية الأقطاب المميزة لديناميكيات Hall–MHD. يحل المخطط هذه السمات بدقة حتى على الشبكات الخشنة نسبياً.

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