A fully decoupled and structure-preserving relaxation Crank--Nicolson finite element method for Gross--Pitaevskii--Poisson model

1. بيان المشكلة

تتناول الورقة البحثية المحاكاة العددية لنظام جروس-بيتايفسكي-بواسون (GPP) المقترن، والذي ينمذج ديناميكيات البلازما فائقة البرودة المكونة من نوعين متفاعلين من الجسيمات (مثل الكاتيونات والأنيونات). يُحكم هذا النظام بواسطة:

• معادلتين غير خطيتين لشرودنجر (معادلات غروس-بيتافسكي) تصفان الدوال الموجية $\psi_+$ و $\psi_-$.
• معادلة بواسون التي تصف الجهد الكهروستاتيكي المتسق ذاتياً $\phi$ الناتج عن كثافات الشحنة $|\psi_+|^2$ و $|\psi_-|^2$.

التحديات الرئيسية:

• اللاخطية القوية: يحتوي النظام على حدود تفاعل ذاتي تكعيبي ($|\psi|^2\psi$)، وحدود تفاعل متبادل ($|\psi_{\mp}|^2\psi_{\pm}$)، والاقتران بجهد بواسون.
• الاقتران: الدوال الموجية والجهد مقترنان بقوة، مما يجعل المخططات الصريحة أو الضمنية القياسية مكلفة حاسوبياً أو غير مستقرة.
• الحفاظ على البنية: النماذج الفيزيائية من هذا النوع تحفظ الكتلة (لكل نوع) والطاقة الإجمالية. غالباً ما تفشل الطرق العددية القياسية في حفظ هذه الثوابت على المستوى المنفصل، مما يؤدي إلى عدم استقرار طويل الأمد أو نتائج غير فيزيائية.
• الكفاءة: تتطلب الطرق الموجودة غالباً حل أنظمة جبرية غير خطية مقترنة ضخمة أو عدة أنظمة خطية في كل خطوة زمنية (على سبيل المثال، طرق التقسيم التي تتطلب ما يصل إلى 7 أنظمة خطية).

2. المنهجية

يقترح المؤلفون طريقة العناصر المحدودة (FEM) لكرانك-نيكلسون بالاسترخاء، المفككة تماماً والمحافظة على الب estrutura.

أ. إعادة الصياغة عبر المتغيرات المساعدة

للتعامل مع اللاخطية، أدخل المؤلفون متغيرين مساعدين:
$$Z_+ = |\psi_+|^2, \quad Z_- = |\psi_-|^2$$
يعيد هذا صياغة حدود التفاعل غير الخطية الأصلية (على سبيل المثال، $|\psi_+|^2\psi_+$) إلى أشكال خطية تتضمن المتغيرات المساعدة (على سبيل المثال، $Z_+\psi_+$). يحول هذا النظام إلى شكل مكافئ حيث يتم "إرخاء" اللاخطية.

ب. المخطط العددي

يجمع المخطط بين:

1. التقطيع الزمني: مخطط كرانك-نيكلسون بالاسترخاء. يستخدم تقريب نقطة المنتصف للمشتق الزمني، ويقوم بعمل متوسط للمتغيرات المساعدة والدوال الموجية.
2. التقطيع المكاني: طريقة العناصر المحدودة (FEM) القياسية باستخدام فضاءات متعددة الحدود الجزئية ($P_k$).

الميزات الرئيسية للمخطط:

• مفكك تماماً: في كل خطوة زمنية، ينقسم النظام إلى مسائل خطية مستقلة.
  • نظامان خطيان مستقلان لـ $\psi_+$ و $\psi_-$ (يتم حلهما بالتوازي).
  • معادلة بواسون واحدة لـ $\phi$.
  • تحديثان خطيان للمتغيرات المساعدة $Z_+$ و $Z_-$.
• لا يتطلب تكراراً: على عكس الطرق الضمنية التقليدية، لا يتطلب هذا المخطط حلولاً تكرارية (مثل نيوتن-رافسون) للتعامل مع اللاخطية لأن تقنية الاسترخاء تجعل النظام خطياً.
• الحفاظ على البنية: تم تصميم المخطط ليحفظ بدقة:
  • الكتلة المنفصلة: $M(\psi_{\pm}) = \int |\psi_{\pm}|^2 dx$ محفوظة لكل نوع.
  • الطاقة المنفصلة المعدلة: يتم حفظ دالة طاقة منفصلة محددة (تتضمن المتغيرات المساعدة) بدقة.

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

1. خوارزمية مبتكرة: تطوير طريقة عناصر محدودة لكرانك-نيكلسون قائمة على الاسترخاء خصيصاً لنظام GPP المقترن مع أنواع متعددة. وهي تقلل التعقيد الحسابي بشكل كبير مقارنة بطرق التقسيم أو الطرق التكرارية الموجودة.
2. الكفاءة الحسابية: يتطلب المخطط حل ثلاثة أنظمة خطية فقط في كل خطوة زمنية (اثنان للدوال الموجية، وواحد للجهد)، بينما قد تتطلب طرق التقسيم المماثلة ما يصل إلى سبعة أنظمة. كما يتجنب تكلفة التكرارات غير الخطية تماماً.
3. تحليل الخطأ الصارم:
   • إثبات الوجود والوحدانية (well-posedness and uniqueness) للمخطط المنفصل تماماً.
   • وضع تقديرات الخطأ المثلى:
     • دقة من الدرجة الثانية في الزمن ($O(\tau^2)$).
     • دقة مثلى من الدرجة $(k+1)$ في الفضاء ($O(h^{k+1})$) في معيار $L^2$ للعناصر $P_k$.
   • يستخدم الإثبات الاستقراء الرياضي، وتقسيم الخطأ القائم على الإسقاط، والتعامل الدقيق مع التفاعل بين الأنواع المتعددة واقتران بواسون.
4. قوانين الحفظ: إثبات صارم بأن المخطط المنفصل يحفظ الكتلة وثابت الطاقة المعدل، مما يضمن الاستقرار طويل الأمد والموثوقية الفيزيائية.

4. النتائج

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

• معدلات التقارب:
  • التقارب المكاني: أكدت التجارب باستخدام عناصر $P_1$ و $P_2$ معدلات التقارب النظرية $O(h^2)$ و $O(h^3)$ على التوالي.
  • التقارب الزمني: أكدت التجارب التقارب من الدرجة الثانية في الزمن ($O(\tau^2)$).
• خصائص الحفظ: أظهرت المحاكاة العددية أن أخطاء الكتلة والطاقة الإجمالية ظلت عند مستوى دقة الآلة (خطأ التقريب) عبر فترات زمنية طويلة، مما يؤكد الطبيعة الحافظة للبنية للمنهجية.
• الظواهر الفيزيائية:
  • نجحت في رصد أنماط موجات الكثافة وظواهر عدم الاستقرار (عدم الاستقرار التعديلي) المتأصلة في نظام GPP.
  • أظهرت القدرة على محاكاة الديناميكيات المعقدة، مثل انقسام السوليتونات وتشكيل حالات ثنائية القطب/رباعية الأقطاب تحت معاملات تفاعل متغيرة.

5. الأهمية

يوفر هذا العمل إطاراً عددياً قوياً، فعالاً، وسليماً من الناحية النظرية لمحاكاة البلازما فائقة البرودة ونماذج تكاثف بوز-أينشتاين متعددة المكونات.

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

باختاًصار، تقدم الورقة البحثية طفرة في الطرق العددية للأنظمة الكمومية المقترنة، حيث توازن بين الدقة العالية، الكفاءة الحسابية، والالتزام الصارم بقوانين الحفظ الفيزيائية.