A Hybrid DEC-SIE Framework for Potential-Based Electromagnetic Analysis of Heterogeneous Media

ملخص تقني: إطار عمل هجين (DEC-SIE) للتحليل الكهرومغناطيسي القائم على الجهد للوسط غير المتجانس

بيان المشكلة
يمثل تحليل المجالات الكهرومغناطيسية في البيئات المعقدة، متعددة المواد، والمفتوحة، تحديات حوسبية كبيرة. تتفوق طرق المعادلات التفاضلية (DE) التقليدية، مثل طريقة العناصر المحدودة (FEM)، في نمذجة المواد غير المتجانسة، لكنها تتطلب قطعاً اصطناعياً للنطاق (على سبيل المثال، عبر الطبقات المتطابقة تماماً PML أو شروط الحدود الماصة ABC) للتعامل مع المناطق المفتوحة، مما يعقد المسائل ذات النطاقات المفتوحة. وفي المقابل، تلبي طرق المعادلات التكاملية (IE) شروط الإشعاع طبيعياً في النطاقات غير المحدودة، لكنها تواجه صعوبة في التعامل مع الأوساط غير المتجانسة. علاوة على ذلك، غالباً ما تعاني حلول المعادلات التفاضلية القياسية من "انهيار التردد المنخفض"، حيث يصبح النظام سيئ الحالة (ill-conditioned) مع اقتراب التردد من الصفر بسبب عدم التوازن بين مقاييس الطاقة الكهربائية والمغناطيسية. وبينما توفر الصيغ القائمة على الجهد (باستخدام الجهد المتجه المغناطيسي $\mathbf{A}$ والجهد السلمي الكهربائي $\Phi$) حصانة ذاتية ضد انهيار التردد المنخفض، فإن ربط هذه الصيغ مع معادلات التكامل السطحي (SIE) للمناطق المفتوحة كان تاريخياً أمراً صعباً. تتضمن معادلات التكامل السطحي القياسية القائمة على الجهد للأوساط النفاذة هياكل عمليات معقدة تحتوي على ما يصل إلى 14 مؤثر تكاملياً مختلفاً، مما يخلق أعباءً عددية وتنفيذية جسيمة.

المنهجية
يقترح المؤلفون إطار عمل عددياً هجيناً يربط بين حساب التفاضل الخارجي المنفصل (DEC) للمنطقة الداخلية غير المتجانسة، ومعادلة التكامل السطحي (SIE) المبسطة للمنطقة الخارجية المتجانسة.

1. الصيغة القائمة على الجهد: تستخدم الطريقة صيغة "كامبس" (Lorenz-gauge) لمعادلات ماكسويل بدلالة الجهد المتجه المغناطيسي $\mathbf{A}$ والجهد السلمي الكهربائي $\Phi$. يتجنب هذا النهج بطبيعته انهيار التردد المنخفض ويتماشى بشكل طبيعي مع الارتباطات متعددة الفيزياء.
2. تجزئة النطاق:
   • الداخل (DEC): يتم تقسيم المنطقة المحدودة غير المتجانسة باستخدام DEC. وخلافاً لطريقة العناصر المحدودة التقليدية، يعد DEC إطار عمل خالٍ من الإحداثيات يحافظ على الثوابت الطوبولوجية ويسمح بالتمييز (discretization) المحافظ على البنية على شبكات غير منتظمة.
   • الخارج (SIE): يتم التعامل مع المنطقة الخارجية المتجانسة غير المحدودة باستخدام معادلات التكامل السطحي (SIEs)، والتي تلبي شرط الإشعاع بطبيعتها وتلغي الحاجة إلى القطع الاصطناعي للنطاق.
3. الابتكار الرئيسي: إعادة صياغة المكون السلمي:
   • تبسيط الـ SIE: يعيد المؤلفون صياغة معادلات SIE من خلال التعبير عن المسألة بدلالة المكونات الكارتيزية للجهد المتجه ($A_x, A_y, A_z$) ومشتقاتها العمودية. يقلل هذا من عدد مؤثرات التكامل السطحي المطلوبة من أربعة عشر (في صيغ الجهد المتجه التقليدية) إلى اثنين فقط: جهدي الطبقة المفردة ($S$) والطبقة المزدوجة ($D$) الكلاسيكيين.
   • تكييف الـ DEC: لضمان التوافق مع الـ SIE المبسط، تم تكييف إطار عمل DEC. بدلاً من معاملة الجهد المتجه $\mathbf{A}$ كـ "1-form" (مرتبط بحواف الشبكة)، يتم التعامل مع كل مكون كارتيزي ($A_x, A_y, A_z$) كـ "0-form" منفصل (مرتبط بعقد الشبكة). يحاكي هذا التعامل مع الحقول السلمية الراسخة في DEC ويتماشى طبيعياً مع تمثيل المكونات السلمية للـ SIE عند الواجهة.
4. آلية الربط: يربط النظام الهجين معادلات DEC الداخلية بمعادلات SIE الحدودية من خلال فرض استمرارية الجهود ومشتقاتها العمودية. يتضمن النظام الخطي الناتج مصفوفات متفرقة وجيدة الحالة (well-conditioned).

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

• تقليل المؤثرات: المساهمة الأساسية هي تقليل مجموعة مؤثرات التكامل السطحي من 14 إلى 2. ويتحقق ذلك من خلال قصر التحليل على المواد غير المغناطيسية ذات النفاذية المتجانسة مكانياً واستخدام تمثيل المكونات السلمية للجهد المتجه.
• الهجنة المحافظة على البنية: نجح العمل في توحيد DEC وSIE عبر تكييف تمثيل الـ $\mathbf{A}$ في DEC لتصبح "0-forms". يتجنب هذا التعقيد في مطابقة كميات الـ "1-form" المرتبطة بالحواف مع مؤثرات التكامل السطحي، مما يبسط الربط عند الواجهة مع الحفاظ على البنية الهندسية والطوبولوجية لإطار عمل DEC.
• الاستقرار عند التردد المنخفض: من خلال استخدام الصيغة القائمة على الجهد، تحافظ الطريقة على الاستقرار والدقة عبر نطاق واسع من الترددات، بما في ذلك نظام التردد المنخفض حيث تفشل الطرق التقليدية غالباً.

النتائج العددية
تم التحقق من صحة الطريقة المقترحة من خلال ثلاث تجارب عددية:

1. الكرة العازلة (Dielectric Sphere): أظهرت الطريقة تقارباً مع تحسين الشبكة، محققة أخطاء $L_2$ نسبية أقل من 1% للشبكات الدقيقة. كما أظهرت متانة ضد تباينات السماحية العالية (تصل إلى $\epsilon_r = 45$) وحافظت على حالة عددية جيدة (أرقام الحالة $\sim 10^3$).
2. الكرة متعددة الطبقات: لوحظ تقارب ودقة مماثلان لكرات ذات طبقات عازلة متعددة.
3. مشتت ثلاثي الأبعاد معقد: تم تحليل هيكل يتكون من أربعة أسطوانات عازلة بيضاوية على بلاطة عازلة. تم التحقق من الدقة باستخدام نظرية الانقراض (extinction theorem)، والتي أكدت أن الحقول التي تمت إعادة بنائها من البيانات الحدودية ألغت الحقل الساقط داخل المشتت بشكل صحيح، مع أخطاء انقراض نسبية في حدود $10^{-4}$ إلى $10^{-5}$.
4. الاستقرار واسع النطاق: أظهرت الط طريقة استقراراً مع اقتراب التردد من الصفر ($k_0 a \to 0$)، حيث ظلت أرقام الحالة محدودة، مما يؤكد الحصانة ضد انهيار التردد المنخفض.

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