A stabilized dual-SAV parametric finite element framework for constrained planar geometric flows with mesh regularization

ملخص تقني: إطار عمل العناصر المحدودة البارامترية ثنائي الـ SAV المستقر للتدفقات الهندسية المستوية المقيدة

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

المنهجية
يقترح البحث إطار عمل العناصر المحدودة البارامية (PFEM) ثنائي المتغير المساعد القياسي (SAV) المستقر للمنحنيات المستوية المغلقة. الابتكار الجوهري هو فصل الطاقة الهندسية الفيزيائية عن طاقة تنظيم الشبكة الاصطناعية إلى متغيرين مساعدين قياسيين منفصلين ($r_g$ و $r_m$).

1. بنية الطاقة المزدوجة:

   • الطاقة الهندسية الفيزيائية ($E_{geom}$): تقود التطور العمودي للمنحنى. يتم إدخال متغير مساعد $r_g$ لخطية القوى الهندسية غير الخطية.
   • طاقة تنظيم الشبكة ($E_{mesh}$): طاقة اصطناعية (عادةً ما تكون طاقة ديريكليه موزونة) تُستخدم فقط لإعادة توزيع نقاط الشبكة مماسياً. يتم إدخال متغير مساعد منفصل $r_m$. ومن الأهمية بمكان أن يتم إسقاط القوة الناجمة عن $E_{mesh}$ على الاتجاه المماسي فقط؛ وهذا يضمن أن إعادة توزيع الشبكة تعمل على تحسين المعلمة (parametrization) دون تغيير التطور العمودي الهندسي المقرر.

2. التقطيع شبه الضمني ذو المقياس المجمد (Frozen-Metric):
   يستخدم المخطط تقطيعًا زمنيًا شبه ضمني حيث يتم تقييم جميع المصطلحات المعتمدة على الهندسة (المقاييس، القوى، القيود) على المنحنى المعروف عند الخطوة الزمنية $n$ ("المقياس المجمد").

   • الاستقرار من الدرجة صفر: تعمل أشكال البيلينير (bilinear forms) للاستقرار (اللايبلاسي أو البيهالمورفيك) على السرعة الجديدة نفسها بدلاً من فروق السرعة. هذا الخيار ضروري للحفاظ على الإلغاء الجبري المطلوب لتقديرات تبدد الطاقة وتجنب نقل تاريخ السرعة.
   • الاستجابة المكانية الخطية: بمجرد تحديد المتغيرات المساعدة القياسية ومعاملات لاغرانج، تصبح المسائل المكانية للسرعات العمودية والمماسية خطية.

3. الاختزال الكتلي الجبري للقيود العالمية:
   يتعامل المنهج مع $K$ من القيود الهندسية العالمية (مثل المساحة، الطول) عبر اختزال كتلي جبري.

   • يتم حل مسائل الاستجابة المكانية الخطية واسعة النطاق أولاً لتوليد "متجهات الاستجابة".
   • يتضمن النظام غير الخطي المتبقي المتغير المساعد الهندسي ($r_g$) ومعاملات لاغرانج $K$ فقط.
   • بناءً على ذلك، يكون بُعد النظام غير الخطي المختزل هو $K+1$، وهو مستقل عن عدد رؤوس الشبكة $N$. بالنسبة للقيود المتزامنة للمساحة والطول ($K=2$)، يؤدي هذا إلى نظام غير خطي ثلاثي الأبعاد.

المساهمات الرئيسية
يحدد البحث خمس مساهمات محددة:

1. صياغة Dual-SAV: تقديم متغيرات مساعدة مستقلة للطاقات الفيزيائية وطاقة الشبكة، مع ربط قوة الشبكة حصريًا بالحركة المماسية.
2. مخطط المقياس المجمد شبه الضمني: بناء تقطيع زمني ينتج مسائل استجابة مكانية خطية مع تلبية تقديرات تبدد الطاقة المعدلة.
3. التقطيع الكامل (PFEM): صياغة تحافظ على بنية الإلغاء الجبري عند المستوى المقطع، مما يضمن تبدد الطاقة المعدلة.
4. التعامل مع القيود المستقل عن الأبعاد: اشتقاق اختزال كتلي جبري يفصل حلول الاستجابة الخطية الكبيرة عن فرض القيود غير الخطية، مما يقلل حل القيود غير الخطي إلى البعد $K+1$.
5. التحقق العددي: اختبار شامل لتدفقات الرتبة الثانية والرابعة (تقصير المنحنى، تقصير المنحنى مع الحفاظ على المساحة، انتشار المنحنى، وتدفقات نوع هيلفريدش - Helfrich).

النتائج
تظهر التجارب العددية ما يلي:

• الدقة: يُظهر تدفق تقصير المنحنى (الدائرة المتقلصة) دقة زمنية من الدرجة الأولى تقريبًا، وهو ما يتوافق مع بنية "الباكورد-أويلر" شبه الضمنية.
• جودة الشبكة: في تدفق تقصير المنحنى مع الحفاظ على المساحة، يمنع إعادة التوزيع المماسي المقترح التكتل الشديد للرؤوس الملاحظ في حالات المرجع اللاغرانجية، مما يحافظ على نسبة شبكة منتظمة تقريبًا ($Q_{mesh} \approx 1$) مع فرض قيود المساحة بصرامة.
• المقاييس غير الموضعية: يدمج الإطار بنجاح مقياس $H^{-1}$ غير الموضعي لتدفقات انتشار المنحنى دون الحاجة إلى قيود مساحة صريحة عبر معاملات لاغرانج (حيث يتم الحفاظ على المساحة طبيعيًا من خلال خاصية المتوسط الصفري لسرعة $H^{-1}$).
• القيود عالية الرتبة: بالنسبة لتدفقات نوع هيلفريدش مع قيود المساحة والطول المتزامنة ($K=2$)، يتم حل النظام غير الخطي المختزل (البعد 3) بكفاءة. يتم فرض القيود بدقة تصل إلى دقة الآلة (أخطاء $\sim 10^{-15}$)، وتتبدد طاقة SAV الهندسية المعدلة بشكل أحادي الاتجاه.
• القابلية للتوسع: يظهر تحليل التكلفة الحسابية أنه مع زيادة عدد الرؤوس $N$، يظل عدد تكرارات نيوتن للنظام غير الخطي المختزل ثابتًا (حوالي 2.2)، بينما تتناسب التكلفة خطيًا مع عدد حلول الاستجابة الخطية.

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

• الاستقرار: التعامل القوي مع الصلابة عالية الرتبة عبر الاستقرار من الدرجة صفر والخطوات الزمنية شبه الضمنية.
• الكفاءة: تجنب حلول سادل-بوينت (saddle-point) غير الخطية الكبيرة؛ حيث تكون التكلفة غير الخطية مستقلة عن دقة الشبكة.
• الحفاظ على البنية: الحفاظ على تبدد الطاقة المعدلة وفرض القيود الهندسية العالمية بدقة.

يشير المؤلفون إلى أن النطاق الحالي يقتصر على المنحنيات المستوية المغلقة. وقد تم تحديد توسيع الإطار ليشمل الأسطح المتطورة في $\mathbb{R}^3$ كعمل مستقبلي يتطلب آليات إضافية لتنظيم شبكة السطح وتحليل الاستقرار. كما يقر البحث بأن العلاقة بين طاقات SAV المعدلة والطاقات الهندسية الأصلية تستحق مزيدًا من الاستقصاء، خاصة بالنسبة لخطوات الزمن الكبيرة.