Optimal error estimates for a fully discrete, highly efficient decoupled scheme for the 2D/3D diffuse interface two-phase MHD flows

ملخص تقني مفصل لورقة بحثية بعنوان: "تقديرات الخطأ المثلى لمخطط منفصل تماماً، عالي الكفاءة، ومفكك لتدفقات الديناميكا المائية المغناطيسية (MHD) ثنائية الطور ذات الواجهة المنتشرة في بعدين/ثلاثة أبعاد"

1. بيان المشكلة

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

• معادلات كاهن-هيليارد (Cahn–Hilliard): تصف تطور حقل الطور (الواجهة) بين السائلين.
• معادلات نافييه-ستوكس (Navier–Stokes): تصف الهيدروديناميكا (السرعة والضغط).
• معادلات ماكسويل (Maxwell): تصف تطور المجال المغناطيسي.

التحديات الرئيسية:

• اللاخطية القوية: يتضمن النظام حدود اقتران معقدة (قوى لورنتز، قوى الشعيرية، الحمل الحراري) وجهوداً لاخطية (جهد جيجنغ-لاند ذو البئر المزدوج).
• الاقتران: المعادلات مقترنة بقوة، مما يجعل المخططات المتكاملة (fully implicit) مكلفة حاسوبياً، والمخططات المفككة (decoupled) صعبة الاستقرار.
• محدودية الدقة: غالباً ما كانت الدراسات السابقة على نماذج مماثلة تعطي تقديرات خطأ دون المستوى الأمثل (على سبيل المثال، أخطاء معيار $L^2$ للسرعة والمجال المغناطيسي كانت أقل بدرجة واحدة من خطأ الاستيفاء) بسبب "التلوث الاصطناعي" الناتج عن تقريب حقل الطور $\phi$. وغالباً ما تفشل مؤثرات الإسقاط القياسية في التعامل مع هذا التلوث بفعالية.

2. المنهجية

يقترح المؤلفون طريقة عناصر محدودة (FEM) منفصلة تماماً، ومقسمة تلاففياً (convex-splitting)، ومفككة.

التقطيع الزمني:

• يستخدم مخطط أويلر الخلفي شبه الضمني (semi-implicit backward Euler) للخطوات الزمنية.
• يعتمد نهج التقسيم التلافي (convex-splitting) لجهود الحدود غير الخطية لضمان استقرار الطاقة غير المشروط.

التقطيع المكاني:

• يستخدم عناصر تايلور-هود (Taylor–Hood) المستقرة وفق معيار (inf-sup) ($P_{r+1}$ للسرعة، $P_r$ للضغط)، وعناصر ميني (Mini) للحالات ذات الرتب المنخفضة.
• يقوم بتقطيع حقل الطور والجهد الكيميائي باستخدام عناصر متصلة.
• يقوم بتقطيع المجال المغناطيسي باستخدام عناصر مناسبة متوافقة مع الفضاء $H(\text{curl})$.

استراتيجية التفكيك (مخطط تصحيح الضغط):
لتجنب حل نظام مقترن ضخم في كل خطوة زمنية، يستخدم المؤلفون مخطط تصحيح الضغط:

1. الخطوة 1: حل حقل الطور ($\phi$)، والجهد الكيميائي ($\omega$)، والسرعة المتوسطة ($\tilde{u}$)، والمجال المغناطيسي ($B$) بشكل مفكك.
2. الخطوة 2: حل معادلة بواسون لتصحيح الضغط.
3. الخطوة 3: تحديث حقل السرعة باستخدام تصحيح الضغط لفرض شرط انعدام التباعد (divergence-free).

الأدوات الرياضية المبتكرة:
للتغلب على "التلوث الاصطناعي" الناتج عن تقريب حقل الطور، قدم المؤلفون تقنيات إسقاط مبتكرة:

• إسقاط ريتز شبه الإهليلجي (Elliptic Ritz Quasi-projection): للجهد الكيميائي $\omega$.
• إسقاط ستوكس شبه الإسقاطي (Stokes Quasi-projection): للسرعة $u$ والضغط $p$.
• إسقاط ماكسويل (Maxwell Projection): للمجال المغناطيسي $B$.
  تم تصميم هذه الإسقاطات خصيصاً للتعامل مع حدود الاقتران التي تتضمن $\nabla \phi$ و $u \times B$، مما يسمح بحدود خطأ مثالية.

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

1. تقديرات الخطأ المثلى: تستنتج الورقة تقديرات خطأ مثلى لمعايير $L^2$ و $H^1$ لحقل الطور، والسرعة، والضغط، والمجال المغناطيسي. وهذا يصحح النتائج السابقة حيث كانت أخطاء السرعة والمجال المغناطيسي دون المستوى الأمثل.
   • لعناصر تايلور-هود ($r \ge 1$)، تكون تقديرات الخطأ هي $O(\Delta t + h^{r+1})$ لحقل الطफा والـ $O(\Delta t + h^{r+2})$ (أو $h^{r+1}$ لـ $r=1$) للسرعة والمجال المغناطيسي في معيار $L^2$.
2. استقرار الطاقة غير المشروط: ثبت أن المخطط مستقر طاقياً بشكل غير مشروط، حيث يحقق قانون طاقة منفصل يحاكي تبدد الطاقة الفيزيائية المستمرة، بغض النظر عن حجم الخطوة الزمنية $\Delta t$.
3. كفاءة التفكيك: من خلال استخدام مخطط تصحيح الضغط، تقلل الطريقة التكلفة الحسابية بشكل كبير مقارنة بالمحللات المتكاملة (monolithic solvers) مع الحفاظ على دقة عالية.
4. التعامل مع التلوث الاصطناعي: يعد التطبيق الناجح لإسقاطات ريتز وستوكس للقضاء على تدهور الدقة الناتج عن تقريبات حقل الطور تقدماً نظرياً هاماً.

4. النتائج

النتائج النظرية:

• النظرية 2.1: تثبت معدلات التقارب المثلى.
  • حقل الطوت $\phi$: خطأ $L^2 \sim O(\Delta t + h^{r+1})$.
  • السرعة $u$ والمجال المغناطيسي $B$: خطأ $L^2 \sim O(\Delta t + \beta_h)$، حيث يحقق $\beta_h$ الرتبة المثلى ($h^{r+2}$ لـ $r \ge 2$).
• النظرية 2.2: تثبت استقرار الطاقة غير المشروط للمخطط المنفصل.

التحقق العددي:

• اختبارات التقارب: أكدت التجارب العددية في بعدين وثلاثة أبعاد باستخدام حلول دقيقة سلسة معدلات التقارب النظرية.
  • الحالة الأولى (عناصر ميني): لوحظ تقارب $O(h^2)$ للسرعة والمجال المغناطيسي في معيار $L^2$.
  • الحالة الثانية (عناصر تايلور-هود): لوحظ تقارب من رتبة أعلى مثالية متسق مع النظرية.
• المحاكاة الفيزيائية:
  • التحلل الإسْبيني (Spinodal Decomposition): أظهر عملية فصل الطور وتحقق من حفظ الكتلة وتبدد الطاقة.
  • عدم استقرار كيلفن-هلمهولتز (K–H Instability): تمت محاكاة اضطرابات جيبية أحادية النمط ومزدوجة النمط. نجحت النتائج في التقاط هياكل دوامات "عين القط" (cat's eye) والتفاف الواجهة، بالإضافة إلى تأثيرات التخميد لقوى لورنتز تحت تأثير نفاذية مغناطيسية ($\mu$) وموصلية ($\sigma$) متغيرة.

5. الأهمية

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