High-order exponential solver method for particle-in-cell simulations in cylindrical geometry

ملخص تقني: طريقة حل أسيّة عالية الرتبة لمحاكاة "الجسيم في خلية" (PIC) في الهندسة الأسطوانية

بيان المشكلة
تعتمد مخططات تسريع الجسيمات المدفوعة بالليزر، ولا سيما تسريع ليزر فجوة الاستيقاظ (LWFA)، بشكل كبير على محاكاة "الجسيم في خلية" (PIC) للفهم المجهري. وبينما تعد عمليات المحاكاة ثلاثية الأبعاد الكاملة دقيقة، إلا أنها مكلفة حاسوبياً. وللتخفيف من ذلك، تُستخدم محاكاة الهندسة الأسطوانية باستخدام تفكيك النمط المداري (AM)، مما يقلل التكاليف الحسابية والذاكرة بعدة مراتب. ومع ذلك، فإن الحلول الحالية، مثل طريقة "فورييه-بيسليل" الطيفية المستخدمة في كود FBPIC، تضحي بالدقة مقابل الكفاءة، خاصة فيما يتعلق بتمثيل الجسيمات بالقرب من المحور الأسطواني والحاجة إلى دوال أساس خاصة. علاوة على ذلك، غالباً ما تعاني طرق الفروق المحدودة القياسية في الإحداثيات الأسطوانية من عدم الاستقرار أو السلوك الطيفي غير الفيزيائي عند الترددات العالية. وهناك حاجة إلى معادل في الفضاء الحقيقي للطرق الطيفية يحافظ على دقة عالية دون الاعتماد على تحويلات أساس معقدة، مع معالجة العيوب العددية المحددة بالقرب من المحور ($\rho=0$).

المنهجية
يقدم المؤلفون حلاً زمنياً أسياً عالي الرتبة متكيفاً مع الهندسة الأسطوانية، بناءً على عملهم السابق في الإحداثيات الكارتيزية. تتضمن المنهجية عدة مكونات رئيسية:

1. حل ماكسويل الأسي في الفضاء الحقيقي:

   • يتم حل معادلات ماكسويل باستخدام نهج فروق محدودة في الفضاء الحقيقي مقترناً بمشغلات الخطوات الزمنية الأسية.
   • تُقرب المشتقات المكانية باستخدام فروق محدودة متدرجة عالية الرتبة (تتراوح من السادسة إلى الثانية والثلاثين).
   • تتم معالجة التطور الزمني عبر توسيعات تايلور الصريحة للمشغلات الأسية ($\exp(\Delta t \hat{H}_m)$)، حيث يمثل $\hat{H}_m$ مشغل ماكسويل لكل نمط مداري $m$.
   • تستخدم الطريقة تفكيك سلسلة فورييه حقيقية للمتغير المداري، مع إقران معاملات الجيب وجيب التمام لتكوين متجهات حقيقية، مما يتجنب التحويلات الطيفية المعقدة.

2. التجزئة الأسطوانية والشروط الحدية:

   • تُستخدم شبكة متدرجة (شبكة Yee) في الاتجاه الشعاعي لضمان الاستقرار وكبح إشعاع شيرينكوف العددي.
   • يُولى اهتمام خاص لسلوك أنماط المجال بالقرب من المحور ($\rho=0$). وقد اشتق المؤلفون شروطاً حدية متماثلة وغير متماثلة محددة لمعاملات نماذج الحقول القياسية والمتجهة بناءً على سلوكها التقاربي (مثل $mS$ للأنماط الشبيهة بالقياسية، و $mA$ للأنماط الشبيهة بالقطبية).
   • تُستخدم استكمال لاغرانج عالي الرتبة لتدرج وإلغاء تدرج المجالات بين أنواع الشبكات، مما يضمن أدنى فقد في الدقة.

3. روتين الجسيم في خلية (PIC):

   • إيداع التيار: تم تكييف خوارزمية "إسيركيبوف" المعدلة لمركبة $z$. بالنسبة للمكونات المستعرضة ($J_\rho, J_\theta$)، تم إدخال خطوة تصحيح التيار، والتي تتضمن حل معادلة بواسون مستوية لتصحيح أخطاء التباعد بالقرب من المحور، والتي يصعب تعريفها بسبب الاقتران الشعاعي-الزاوي.
   • حفظ الشحنة: قدم المؤلفون مخطط تربيع بدقة من الرتبة الرابعة لتحسين حفظ الشحنة بالقرب من $\rho=0$، مما يعالج المشكلات حيث تتجاوز التقلبات 10% في الطرق القياسية.
   • أخذ العينات الزاوية: تم تنفيذ مخطط يسمح بإزاحات زاوية عشوائية للجسيمات داخل الخلية لتقليل عيوب الشبكة المرتبة وتحسين أخذ عينات طور الفضاء الزاوي.

4. نموذج غلاف الليزر:

   • تم تطوير حل أسي لمعادلة غلاف جهد الليزر القياسي. يسمح هذا بنشر غلاف الليزر دون الحاجة لحل التردد الحامل، مما يقلل التكلفة الحسابية بشكل كبير للنبضات الطويلة.
   • يتضمن النموذج تقريب مركز التوجيه لـ "بوندرموتيف"، حيث تتفاعل الجسيمات مع الحقول متوسطة الدورة.

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

• حل أسي في الفضاء الحقيقي: يقدم البحث طريقة الفروق المحدودة في النطاق الزمني الأسي (FDETD) للهندسة الأسطوانية، والتي تعمل كمعادل في الفضاء الحقيقي للطرق الطيفية، مما يلغي الحاجة إلى تحويلات فورييه-بيسليل.
• دقة عالية الرتبة: من خلال استخدام فروق محدودة متدرجة عالية الرتبة (تصل إلى 32) مع الخطوات الزمنية الأسية، تحقق الطريقة دقة تشبه الدقة الطيفية في الفضاء الحقيقي.
• التعامل مع المحور: توفر الدراسة معالجة دقيقة للشروط الحدية بالقرب من المحور وتمثيل الجسيمات، بما في ذلك خوارزمية تصحيح تيار محددة للتعامل مع أخطاء التباعد عند $\rho=0$.
• تكامل غلاف الليزر: يتيح تطوير حل أسي لجهد غلاف الليزر محاكاة دقيقة وعالية الكفاءة لانتشار الليزر في البلازما منخفضة الكثافة، خاصة عندما يكون غلاف الليزر متماثلاً أسطوانياً.

النتائج والاختبارات المعيارية
تحقق المؤلفون من الطريقة من خلال عدة اختبارات معيارية:

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

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