Evolving finite elements for advection diffusion with an evolving interface

ملخص تقني: العناصر المحدودة المتطورة لنقل الانتشار مع واجهة متطورة

بيان المشكلة
تتناول هذه الورقة التقريب العددي للمعادلات المكافئة (parabolic equations) المعرفة على مجال ثابت $\Omega$ يحتوي على واجهة متحركة $\Gamma(t)$. تقسم الواجهة المجال إلى نطاقين فرعيين $\Omega_1(t)$ و $\Omega_2(t)$. تنمذج المعادلات الحاكمة ديناميكيات النقل والانتشار لـ حقول قياسية $u_1$ و $u_2$ داخل هذه النطاقات الفرعية، مع خضوعها لـ:

• ديناميكيات النقل-الانتشار-التفاعل داخل كل نطاق فرعي مع معاملات تعتمد على الزمن ($A_i, B_i, C_i$).
• شرط حدي من نوع ديريكليه (Dirichlet) على الحد الخارجي $\partial\Omega$.
• استمرارية الحل عبر الواجهة ($u_1 = u_2$).
• قفزة تدفق محددة عبر الواجهة ($A_1 \partial_{\nu} u_1 - A_2 \partial_{\nu} u_2 = g$).

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

المنهجية
يقترح المؤلفون نهج "لاغرانج-أويلر التعسفي" (Arbitrary Lagrangian-Eulerian - ALE) باستخدام فضاءات العناصر المحدودة المتطورة (ESFEM) بناءً على الإطار المجرد الذي وضعه المرجع [21].

1. الإعداد التحليلي الوظيفي: تمت صياغة المشكلة المستمرة ضمن سياق فضاءات هيلبرت المتطورة. يعرّف المؤلفون أزواجاً متوافقة من الفضاءات $(X(t), \phi_t)$ حيث $\phi_t$ هي خريطة تدفق تنقل الدوال من مجال مرجعي إلى المجال المتطور. يسمح هذا بتعريف مشتقات مادية ضعيفة واستخلاص صيغة ضعيفة جيدة التحديد (مبرهنة 2.16).
2. بناء الشبكة:
   • الشبكة الأولية: يتم تقسيم المجال في البداية باستخدام مضلعات بسيطة خطية. ومن الأهمية بمكان أن تقع رؤوس أوجه الواجهة تماماً على الواجهة الأولية $\Gamma(0)$ (الافتراض M5).
   • التحويل الأيزو-بارامتري (Isoparametric Mapping): لتحقيق دقة أعلى رتبة، يتم تحويل الشبكة الخطية إلى شبكة أيزو-بارامترية باستخدام التشاكل $\Psi_h$. تقوم هذه الخريطة بدفع أوجه الواجهة لتستقر تماماً على الواجهة المنحنية الدقيقة $\Gamma(0)$ مع الحفاظ على استمرارية الشبكة.
   • التطور الزمني: تتطور الشبكة عن طريق تحريك العقد اللاغرانجية وفقاً لمجال سرعة عالمي سلس ومعروف $w$. يتم بناء التدفق المنفصل $\Phi^h$ كاستيفاء حدودي (polynomial interpolation) للتدفق المستمر، مما يضمن بقاء الشبكة كتقسيم متوافق مقبول في جميع الأوقات.
3. مؤثر الرفع (Lift Mapping): تم إدخال مؤثر رفع $\Lambda_h$ لتعيين الدوال من المجال الحسابي المنفصل (الذي يقرب الواجهة المنحنية) إلى المجال المستمر الدقيق. هذا الأمر ضروري لمقارنة الحل المنفصل بالحل الدقيق في تحليل الخطأ.
4. المخطط المنفصل: تستخدم الطريقة عناصر أيزو-بارامترية من الرتبة $k$. يتم صياغة المسألة التباينية المنفصلة باستخدام المشتق المادي المنفصل والمعاملات المرفوعة. يتم حل المخطط باستخدام طريقة الخطوات الزمنية (تحديداً صيغ الفرق الخلفي BDF في التجارب العددية).

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

• جودة التحديد (Well-Posedness): توفر الورقة إعداداً تحليلياً وظيفياً دقيقاً يثبت وجود وتفرد الحل للمشكلة المستمرة (مبرهنة 2.16)، بما في ذلك نتائج انتظام إضافية تحت افتراضات محددة على البيانات.
• إطار عمل ESFEM الأيزو-بارامتري: يوسع المؤلفون إطار العناصر المحدودة المتطورة المجرد للتعامل مع مسائل الواجهة ثنائية الجانب مع قفزات التدفق. كما يفصلون بناء العناصر الأيزو-بارامترية التي تقرب الواجهة المنحنية بدقة، مما يتيح تقارباً من رتبة أعلى.
• حدود الخطأ المثلى: تثبت النتيجة النظرية الرئيسية (مبرهنة 5.1) حدود خطأ ذات رتبة مثلى. وتحديداً، يتم تحديد الخطأ في معيار $L^2$ بـ $O(h^{k+1})$، والخطأ في نصف معيار $H^1$ بـ $O(h^k)$، حيث $h$ هو حجم الشبكة و $k$ هي درجة كثير الحدود. وهذه الحدود تظل صالحة حتى ثابت معين وتطابق خطأ الاستيفاء لحل سلس معروف.
• التنفيذ والتحقق: تتضمن الورقة تنفيذاً قوياً باستخدام حزمة Firedrake و GMSH. وتتحقق التجارب العددية في بعدين وثلاثة أبعاد من معدلات التقارب النظرية لمختلف رتب كثيرات الحدود.

النتائج
تُظهر التجارب العددية فعالية الطريقة على مسألة نموذجية ذات حل دقيق معروف، وهو حل مستمر ولكنه غير قابل للاشتقاق عبر الواجهة.

• التقارب: بالنسبة للتقسيم المكاني من الرتبة $k$، تؤكد النتائج العددية تقارباً من الرتبة $k+1$ في معيار $L^2$ ومن الرتبة $k$ في نصف معيار $H^1$.
• المتانة: تظل النتائج صالحة لكل من النطاقات ثنائية وثلاثية الأبعاد، حتى عندما يتم إنشاء الشبكات بشكل مستقل (بدلاً من التكرار المنتظم)، مما يشير إلى استقرار الطريقة.
• توافر الكود: يوفر المؤلفون كود المحاكاة للسماح بإعادة الإنتاج والاستخدام من قبل المجتمع العلمي.

الأهمية والنطاق
تدعي الورقة أنها توفر إطاراً نظرياً وعملياً كاملاً لحل مسائل النقل-الانتشار مع الواجهات المتحركة باستخدام عناصر محدودة أيزو-بارامترية عالية الرتبة. تكمن الأهمية في سد الفجوة بين نظرية الفضاءات المتطورة المجردة والتطبيق العملي لمسائل الواجهة مع قفزات التدفق.

يحافظ المؤلفون على نطاق متواضع فيما يتعلق بمجال السرعة: حيث يفترضون أن مجال سرعة عالمي وسلس معطى مسبقاً. ويصرحون بوضوح أن العمل لا يتناول كيفية حساب مجال السرعة هذا (على سبيل المثال، عبر إعادة المعلمة لمنع تشابك الشبكة)، مشيرين إلى أن مثل هذه الطرق موجودة في الأدبيات [8, 18, 19]. علاوة على ذلك، بينما يمكن للطريقة نظرياً أن تنطبق على أنظمة المعادلات التفاضلية الجزئية المقترنة حيث يكون تطور المجال جزءاً من الحل، إلا أن التحليل الحالي لا يغطي مثل هذه المسائل المقترنة. يظل التركيز صارماً على تقريب مسألة الإرسال (transmission problem) المكافئة المعطاة لتطور مجال معلوم.