2025-11-10T03:15:01.076525

Construction of Basis Functions for the Geometry Conforming Immersed Finite Element Method

Adjerid, Lin, Meghaichi
The Frenet apparatus is a new framework for constructing high order geometry-conforming immersed finite element functions for interface problems. In this report, we present a procedure for constructing the local IFE bases in some detail as well as a new approach for constructing orthonormal bases using the singular value decomposition of the local generalized Vandermonde matrix. A sample implementation in MATLAB is provided to showcase the simplicity and extensionability of the framework.
academic

بناء دوال الأساس لطريقة العناصر المحدودة المغمورة المطابقة للهندسة

المعلومات الأساسية

  • معرّف الورقة: 2510.12018
  • العنوان: Construction of Basis Functions for the Geometry Conforming Immersed Finite Element Method
  • المؤلفون: Slimane Adjeril (Virginia Tech)، Tao Lin (Virginia Tech)، Haroun Meghaichi (Ohio State University)
  • التصنيف: math.NA cs.NA
  • تاريخ النشر: 15 أكتوبر 2025
  • رابط الورقة: https://arxiv.org/abs/2510.12018

الملخص

تقترح هذه الورقة إطار عمل جديد قائم على جهاز Frenet لبناء دوال العناصر المحدودة المغمورة (IFE) عالية الرتبة المطابقة للهندسة لحل مسائل الواجهات. تقدم الورقة تفاصيل عملية بناء دوال IFE الأساسية المحلية، وطريقة جديدة لبناء أساس متعامد معياري باستخدام تحليل القيم الذاتية للمصفوفات Vandermonde المعممة المحلية. يتم توفير تطبيق مثال في MATLAB يوضح بساطة الإطار وقابليته للتوسع.

الخلفية البحثية والدافع

خلفية المشكلة

يتناول هذا البحث مسائل الواجهات الإهليلجية: −∇⋅(β∇u)=f, on Ω−∪Ω+-\nabla \cdot (\beta\nabla u) = f, \text{ on } \Omega^- \cup \Omega^+ حيث يتم تقسيم المجال الحسابي Ω⊂R2\Omega \subset \mathbb{R}^2 بواسطة الواجهة Γ\Gamma إلى نطاقين فرعيين Ω−\Omega^- و Ω+\Omega^+، ومعامل الانتشار β\beta هو دالة ثابتة متعددة التعريف.

دافع البحث

  1. قيود طرق IFE الموجودة: تستطيع طرق العناصر المحدودة المغمورة التقليدية فقط تقريب شروط القفز على الواجهة
  2. متطلبات الدقة العالية: الحاجة إلى بناء دوال IFE لمتعددات الحدود من رتب عشوائية
  3. متطلبات المطابقة الهندسية: الحاجة إلى دوال IFE تحقق شروط القفز على الواجهة بدقة
  4. الاستقرار الحسابي: الحاجة إلى نظام دوال أساسية بأرقام شرط جيدة

الأهمية

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

المساهمات الأساسية

  1. اقتراح إطار عمل بناء GC-IFE قائم على جهاز Frenet: استخدام الإطار المرجعي Frenet من الهندسة التفاضلية لتبسيط شروط القفز على الواجهة
  2. تطوير طريقتي بناء دوال أساسية:
    • طريقة البناء الأولي القائمة على صيغة خاصة
    • طريقة البناء المعممة القائمة على فكرة التوسيع
  3. اقتراح خوارزميتي إعادة تعامد:
    • الطريقة القائمة على تحليل SVD لمصفوفة الكتلة (الطريقة 1)
    • الطريقة القائمة على تحليل SVD لمصفوفة Vandermonde المعممة (الطريقة 2)
  4. توفير تطبيق MATLAB كامل: يوضح الجدوى العملية للطريقة وقابليتها للتوسع

شرح الطريقة

تحويل Frenet

الفكرة الأساسية هي استخدام تحويل Frenet لتحويل الإحداثيات الفيزيائية (x,y)(x,y) إلى إحداثيات محلية (η,ξ)(\eta,\xi): [ηξ]→PΓ[x(η,ξ)y(η,ξ)]=g(ξ)+ηn(ξ)\begin{bmatrix} \eta \\ \xi \end{bmatrix} \xrightarrow{P_\Gamma} \begin{bmatrix} x(\eta,\xi) \\ y(\eta,\xi) \end{bmatrix} = g(\xi) + \eta n(\xi)

حيث g(ξ)g(\xi) هي معاملة منحنى الواجهة، و n(ξ)n(\xi) هو متجه الوحدة الناظم.

تحويل شروط القفز على الواجهة

في النظام الإحداثي المحلي، تصبح شروط القفز على الواجهة:

  1. ⟦u^⟧Γ^KF=0\llbracket \hat{u} \rrbracket_{\hat{\Gamma}_{K_F}} = 0
  2. ⟦β^u^η⟧Γ^KF=0\llbracket \hat{\beta}\hat{u}_\eta \rrbracket_{\hat{\Gamma}_{K_F}} = 0
  3. ⟦β^∂j∂ηjL(u^)⟧Γ^KF=0,j=0,1,…,m−2\llbracket \hat{\beta}\frac{\partial^j}{\partial\eta^j}\mathcal{L}(\hat{u}) \rrbracket_{\hat{\Gamma}_{K_F}} = 0, j = 0,1,\ldots,m-2

طرق بناء دوال الأساس

الطريقة 1: البناء بصيغة خاصة

بناء دوال أساسية بالصيغة التالية: ϕ^i,j(η,ξ)=1β^(η,ξ)qj(ηηh)pi(ξ−ξmidξh)\hat{\phi}_{i,j}(\eta,\xi) = \frac{1}{\hat{\beta}(\eta,\xi)}q_j\left(\frac{\eta}{\eta_h}\right)p_i\left(\frac{\xi-\xi_{mid}}{\xi_h}\right)

بالنسبة للحالات 1≤j≤m,0≤i≤m1 \leq j \leq m, 0 \leq i \leq m، تحقق هذه الدوال شروط القفز تلقائياً.

الطريقة 2: البناء المعمم بالتوسيع

استخدام صيغة أكثر عمومية:

\hat{\lambda}_j^-(\eta,\xi) = \sum_{i=1}^{(m+1)^2} C_{i,j}^- R_i(\eta,\xi), & \eta < 0 \\ \hat{\lambda}_j^+(\eta,\xi) = \sum_{i=1}^{(m+1)^2} C_{i,j}^+ R_i(\eta,\xi), & \eta > 0 \end{cases}$$ يتم تحديد مصفوفة المعاملات $C^{\pm}$ من خلال حل نظام معادلات خطية. ### إعادة التعامد #### الطريقة 1: تحليل SVD لمصفوفة الكتلة بفرض $M_q = V_1\Lambda V_1^T$، نختار $Q_1 = V_1\Lambda^{-1/2}$ بحيث تصبح مصفوفة الكتلة المعاد بناؤها مصفوفة الوحدة. #### الطريقة 2: تحليل SVD لمصفوفة Vandermonde إجراء تحليل SVD على $\tilde{V} = \sqrt{W}V$: $\tilde{V} = U_2\Sigma V_2^T$، واختيار $Q_2 = V_2\Sigma^{-1}$. يشير التحليل النظري إلى أن الطريقة 2 أكثر استقراراً في التطبيق العددي، لأنها تتجنب تحليل SVD لمصفوفة برقم شرط مرفوع للقوة الثانية. ## إعداد التجارب ### مسائل الاختبار النظر في المجال $\Omega = [-1,1]^2$، مع واجهة $x^2 + y^2 = r_0^2$ (حيث $r_0 = 1/\sqrt{3}$)، دالة الاختبار: $$u(x,y) = \begin{cases} \frac{1}{\beta^+}\cos(2\pi r^2), & r > r_0 \\ \frac{1}{\beta^-}\cos(2\pi r^2) + \cos(2\pi r_0^2)\left(\frac{1}{\beta^+} - \frac{1}{\beta^-}\right), & r < r_0 \end{cases}$$ ### مؤشرات التقييم - رتبة التقارب لخطأ الإسقاط $L^2$ - رقم الشرط لمصفوفة الكتلة - الاستقرار العددي لبناء دوال الأساس ## نتائج التجارب ### التحقق من التقارب بالنسبة لرتب متعددات حدود مختلفة $m = 1,2,3,4$، يظهر خطأ الإسقاط $L^2$ رتبة تقارب مثلى: | $m$ | $N = 16$ | $N = 32$ | $N = 64$ | $N = 128$ | $N = 256$ | رتبة التقارب | |-----|----------|----------|----------|-----------|-----------|--------| | 1 | 8.14E-02 | 2.08E-02 | 5.23E-03 | 1.31E-03 | 3.28E-04 | ≈2.0 | | 2 | 9.29E-03 | 1.19E-03 | 1.50E-04 | 1.88E-05 | 2.36E-06 | ≈3.0 | | 3 | 8.80E-04 | 5.74E-05 | 3.62E-06 | 2.27E-07 | 1.42E-08 | ≈4.0 | | 4 | 7.55E-05 | 2.39E-06 | 7.49E-08 | 2.34E-09 | 7.33E-11 | ≈5.0 | ### تأثير تحسين رقم الشرط تظهر المقارنة بين أرقام الشرط لمصفوفة الكتلة قبل وبعد إعادة البناء: - يزداد رقم الشرط للبناء الأولي بشكل أسي مع رتبة متعددة الحدود - تستطيع الطريقة 2 تقليل رقم الشرط بشكل مستقر إلى 1.0 - تظهر الطريقة 1 عدم استقرار عددي في الحالات عالية الرتبة ### تأثير المعالجات المسبقة بالنسبة لمصفوفة المعاملات $A$ و $\tilde{A}$، يستطيع معالج Jacobi المسبق والمعالج المسبق بالتطبيع الصفي تحسين رقم الشرط بشكل كبير، خاصة في حالات الرتبة العالية وأحجام الشبكة الصغيرة. ## الأعمال ذات الصلة ### تطور طريقة IFE - طرق IFE التقليدية: تستطيع فقط تقريب شروط الواجهة - طرق IFE الموسعة: طرق عالية الرتبة قائمة على توسيع Cauchy - طريقة GC-IFE في هذه الورقة: طريقة مطابقة هندسية تحقق شروط الواجهة بدقة ### تقنيات التعامد - بناء متعددات حدود متعامدة في طرق العناصر الافتراضية - دوال أساسية متعامدة في طرق Galerkin المتقطعة - تقديم فكرة التعامد في إطار عمل IFE في هذه الورقة ## الخلاصة والمناقشة ### الاستنتاجات الرئيسية 1. **المساهمة النظرية**: إنشاء إطار عمل نظري GC-IFE قائم على جهاز Frenet 2. **الابتكار الخوارزمي**: اقتراح طريقتي بناء واستراتيجيتي تعامد 3. **التحقق العددي**: إثبات التقارب الأمثل والاستقرار العددي للطريقة 4. **القيمة العملية**: توفير تطبيق MATLAB كامل ### القيود 1. **القيود الهندسية**: تتطلب منحنى الواجهة أن يكون ناعماً بدرجة كافية ($C^3$) 2. **التعقيد الحسابي**: تكلفة تجميع المصفوفة أعلى في الحالات عالية الرتبة 3. **قابلية التوسع**: تنطبق حالياً فقط على مسائل الواجهات الإهليلجية ثنائية الأبعاد ### الاتجاهات المستقبلية 1. التوسيع إلى مسائل ثلاثية الأبعاد والهندسات المعقدة 2. التطبيق على مسائل الواجهات المعتمدة على الزمن 3. التوسيع إلى ميكانيكا المرونة وأنظمة Stokes 4. تطوير خوارزميات تطبيق أكثر كفاءة ## التقييم المتعمق ### المميزات 1. **الصرامة النظرية**: يوفر جهاز Frenet من الهندسة التفاضلية أساساً نظرياً متيناً 2. **ابتكار الطريقة**: أول تحقيق لطريقة IFE عالية الرتبة تحقق شروط القفز على الواجهة بدقة 3. **الاستقرار العددي**: تحسن إعادة التعامد بشكل كبير الاستقرار الحسابي 4. **التطبيق الكامل**: توفير وصف خوارزمي مفصل وكود MATLAB 5. **التحقق الكافي**: تؤيد التحليلات النظرية والتجارب العددية بعضها البعض ### أوجه القصور 1. **نطاق التطبيق محدود**: يقتصر حالياً على مسائل الواجهات الإهليلجية ثنائية الأبعاد 2. **تحليل التعقيد الحسابي**: لم يكن تحليل التعقيد الحسابي في الحالات عالية الرتبة عميقاً بما يكفي 3. **متطلبات هندسية صارمة**: متطلبات عالية لنعومة منحنى الواجهة 4. **حالات التطبيق العملي**: نقص التحقق من التطبيق على مسائل هندسية معقدة ### التأثير 1. **القيمة الأكاديمية**: توفير إطار نظري جديد وطريق تقني لطرق IFE 2. **القيمة العملية**: توفير أداة فعالة لحل مسائل الواجهات عالية الدقة 3. **قابلية التكرار**: يسهل التطبيق الكامل للكود على الباحثين الآخرين استخدامه وتوسيعه 4. **إمكانية التطور**: يضع أساساً متيناً للأبحاث اللاحقة ### السيناريوهات القابلة للتطبيق 1. **التدفقات متعددة الأطوار**: محاكاة عددية لواجهات السوائل المختلفة 2. **علوم المواد**: حسابات دقيقة عالية لمسائل واجهات المواد المركبة 3. **التطبيقات الطبية الحيوية**: مسائل نقل الحرارة والكتلة على واجهات الأنسجة 4. **الجيوفيزياء**: مسائل انتشار الموجات على واجهات الطبقات الجيولوجية ## المراجع تستشهد الورقة بـ 22 مرجعاً ذا صلة، تغطي أعمالاً مهمة في الهندسة التفاضلية والتحليل العددي وطرق العناصر المحدودة وغيرها، مما يوفر أساساً نظرياً متيناً لهذا البحث.