Convergence of a CutFEM for fluid--structure interaction with a deforming interface

ملخص تقني: تقارب طريقة CutFEM للتفاعل بين المائع والمنشأ مع واجهة مشوهة

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

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

المنهجية
يقترح المؤلفون ويحللون طريقة CutFEM (طريقة العناصر المحدودة المقطوعة) شبه المنفصلة زمنياً لهذا النظام المتكامل تماماً من التفاعل بين المائع والمنشأ. تستخدم الطريقة المكونات الرئيسية التالية:

1. الشبكة غير المتوافقة: يتم تقسيم المائع على شبكة خلفية ثابتة لا تتوافق مع حدود المنشأ. يتم تقسيم المنشأ باستخدام شبكة متوافقة "إيزو-بارامتريّة" (isoparametric) على نطاقه المرجعي.
2. صيغة مضاعف لاغرانج: يتم فرض شروط الاقتران بين المائع والمنشأ بشكل ضعيف باستخدام مضاعف لاغرانج. يمثل هذا المضاعف قوة شد المائع المسحوبة إلى واجهة المنشأ المرجعية. يتم فرض الشرط الكينماتيكي ($u = \partial_t d \circ X^{-1}$) على الواجهة الحالية، بينما تُصاغ موازنة الشد عبر علاقة تكاملية تتضمن التشكيل المرجعي.
3. استقرار عقوبة الشبح (Ghost Penalty): لضمان الاستقرار والمتانة ضد العناصر المقطوعة الصغيرة بشكل عشوائي، تُضاف حدود عقوبة الشبح العالمية إلى صيغ السرعة والضغط. تعمل هذه الحدود على معاقبة القفزات في المشتقات العمودية عبر الأوجه الداخلية للشبكة الخلفية، مما يؤدي فعلياً إلى توسيع استقرار فضاء العناصر المحدودة إلى نطاق المائع الموسع.
4. المخطط شبه المنفصل زمنياً: المخطط منفصل مكانياً فقط، مما يعني أن الزمن مستمر. يتم تقريب سرعة المائع، والضغط، وإزاحة المنشأ، ومضاعف لاغرانج باستخدام فضاءات العناصر المحدودة من الدرجة $k$ (للسرعة والإزاحة) و $k-1$ (للضغط).

المساهمات الرئيسية والإطار النظري
تتمثل المساهمة المركزية لهذا البحث في كونه أول تحليل دقيق لخطأ "قبلي" (a priori) لطريقة عناصر محدودة غير متوافقة مطبقة على نظام FSI متكامل تماماً مع واجهة مشوهة حقاً. كانت التحليلات الدقيقة السابقة تقتصر إلى حد كبير على النماذج الخطية، أو تقريبات الإزاحة الصغيرة، أو طرق الشبكة المتوافقة.

يعتمد إثبات التقارب على عدة مكونات تحليلية مبتكرة:

• التشكيل المستند إلى الاستيفاء (Interpolated Configuration): يقدم المؤلفون تشكيلاً مستنداً إلى الاستيفاء $X^*_H$ للحل الدقيق لمقارنة الواجهة الدقيقة $\Gamma[X]$، والواجهة المستوفاة $\Gamma[X^*_H]$، والواجهة العددية $\Gamma[X_H]$. يسمح هذا النظام بالقياس المنهجي للأخطاء الناجمة عن إزاحة الواجهة.
• التقديرات الموحدة: يثبت التحليل تكافؤات المعايير الموحدة، ومتراجحات الأثر (trace inequalities)، والتقديرات العكسية، وشروط inf–sup على النطاقات المقطوعة المتحركة. والأهم من ذلك، أن هذه التقديرات موحدة بالنسبة لموضع الواجهة داخل الشبكة الخلفية وبالنسبة للزمن، بشرما كان التشوه العددي قريباً بما يكفي من التشوه الدقيق.
• اقتران الأخطاء: إحدى الصعوبات التحليلية الكبرى التي تمت معالجتها هي الاعتماد المتبادل ثنائي الاتجاه بين خطأ المائع وخطأ الهندسة (إزاحة المنشأ). يؤدي الخطأ في إزاحة المنشأ إلى تغيير نطاق المائع والواجهة عددياً، وهو ما يؤثر بدوره على حل المائع. يقوم المؤلفون بإغلاق هذه الحلقة عبر اشتقاق تقديرات الاستقرار لعدم التطابق الكينماتيكي المسقط على الواجهة واستخدام حجة الاستمرار في الزمن.
• حجة الاستمرار: يفترض الإثبات وجود حدود "قبلية" (على وجه التحديد حدود $W^{1,\infty}$ و $L^\infty$) حتى زمن أقصى $t^*$. ومن خلال إثبات تقديرات الطاقة بثوابت مستقلة عن $t^*$، وإظهار أنه بالنسبة لأحجام شبكة صغيرة و $k \ge 3$، فإن معدلات التقارب المستمدة تعني ضمن الحدود المفروضة، يمدد المؤلفون صلاحية التقديرات لتشمل الفترة الزمنية الكاملة $[0, T]$.

النتائج الرئيسية
تحت افتراضات كفاية انتظام الحل الدقيق وأن أحجام الشبكة تحقق $c_1 H \le h \le c_2 H$ (حيث $h$ و $H$ هما حجم شبكة المائع والمنشأ على التوالي)، يثبت المؤلفون معدلات التقارب التالية للعناصر المحدودة من الدرجة $k \ge 3$:

• سرعة المائع ($u_h$) وإزاحة المنشأ ($d_H$): تقارب من الرتبة $k$ في معايير الطاقة الطبيعية.
• ضغط المائع ($p_h$): تقارب من الرتبة $k-1$ في معيار $L^2$ (المرتبط بمعيار الطاقة).
• مضاعف لاغرانج (شد الواجهة): تقارب من الرتبة $k-1$ في معيار $H^{-1/2}$.

وعلى وجه التحديد، يُعطى حد الخطأ كالآتي:
$$\max_{t \in [0,T]} \left( \|u_h - u\|_{L^2}^2 + \|d_H - d \circ \Phi_H\|_{H^1}^2 + \|\partial_t d_H - \partial_t d \circ \Phi_H\|_{L^2}^2 \right) + \int_0^T \|u_h - u\|_{H^1}^2 dt + H^2 \int_0^T (\|p_h - p\|_{L^2}^2 + \|\sigma_H - \sigma \circ \Phi_H\|_{H^{-1/2}}^2) dt \le C(h^{2k} + H^{2k})$$

الأهمية
يزعم هذا البحث أن هذا العمل يمثل أول تحليل تقارب دقيق لطريقة عناصر محدودة غير متوافقة لنظام FSI متكامل تماماً مع واجهة مشوهة. وبينما حللت أعمال سابقة طرقاً غير متوافقة لنطاقات متحركة محددة مسبقاً أو مسائل مقترنة مبسطة، فإن اللاخطية وحلقة التغذية الراجعة بين حل المائع والهندسة المتحركة في نظام FSI متكامل تماماً لم يتم التعامل معها بصرامة في السياق غير المتوافق. توفر النتائج أساساً نظرياً لاستخدام CutFEM في تطبيقات FSI المعقدة التي تنطوي على تشوهات كبيرة، مما يثبت دقة وطريقة واستقرار المنهج دون الاعتماد على تقريبات الإزاحة الصغيرة أو النطاق الثابت. وتجدر الإشارة إلى أن القيد $k \ge 3$ هو قيد تقني يتعلق بتقديرات عكسية محددة مستخدمة في الإثبات، ويشير المؤلفون إلى أن النتيجة قد تمتد لتشمل $k=2$ مع تحسينات مستقبلية.