Efficient generation of large-scale non-equilibrium distributions of particles

ملخص تقني: التوليد الفعال لتوزيعات الجسيمات غير المتوازنة واسعة النطاق

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

المنهجية: خوارزمية التورم والهجرة العشوائية (SRM)
تقترح الورقة البحثية خوارسة "التورم والهجرة العشوائية" (SRM)، وهي نهج هجين يجمع بين تقنية "التورم" المعدلة لـ ميلر-توركاتو وطريقة "الإضافة التسلسلية العشوائية" (RSA) المعدلة. تعمل الطريقة على صندوق وحدة يحتوي على شوائب أحادية التشتت، وتتبع الخطوات التكرارية التالية:

1. التهيئة: يتم توليد تكوين أولي غير متداخل عند كسر حجم منخفض باستخدام طريقة RSA الكلاسيكية.
2. التورم: يتم زيادة نصف قطر جميع الجسيمات في آن واحد بمقدار زيادة صغيرة ($c_w$).
3. الهجرة العشوائية: يتم إزاحة كل جسيم بمقدار متجه عشوائي صغير ($c_m$).
4. فحص الصلاحية المحلية: يتم تقييم المسافات مع الجيران باستخدام نظام تقسيم مكاني تكيفي يعتمد على الخلايا. إذا حدث تداخل للجسيم، يتم إعادته إلى موضعه السابق، وتتم محاولة إزاحة عشوائية جديدة. تستمر هذه الحلقة حتى يتم العثور على موضع صالح أو الوصول إلى الحد الأقصى من المحاولات ($N_l$).
5. خطوة الهز: إذا لم يتمكن جسيم من إيجاد موضع صالح بعد عدة محاولات، يتم إجراء استرخاء عالمي ("هز"). يتضمن ذلك إزاحة جميع الجسيمات عشوائياً واختبار التداخلات للتغلب على الاحتقان المحلي، خاصة في التكوينات المتكتلة.
6. الإنهاء: تتكرر الدورة حتى يتم الوصول إلى كسر الحجم المستهدف.

الكفاءة الحوسبية
المكون الحاسم للمنهجية هو تنفيذ بحث عن الجيران يعتمد على الخلايا. من خلال تقسيم الـ RVE إلى خلايا أكبر قليلاً من قطر الجسيم، تقوم الخوارمة بحصر فحوصات التداخل في الجسيمات الموجودة داخل الخلية الحالية وجيرانها المباشرين (كتلة $3 \times 3$ في البعدين). يقلل هذا من التعقيد الحوسبي لاكتشاف التصادم من التدرج التربيعي إلى التدرج القريب من الخطي ($O(N)$) بالنسبة لعدد الجسيمات. ومع تحسينات إضافية لأنماط الوصول إلى الجسيمات للحد من أخطاء ذاكرة الكاش (cache misses)، يتدرج زمن التوليد تقريباً وفق العلاقة $t \propto N^{1.15}$، مما يسمح بإجراء عمليات محاكاة تصل إلى $10^7$ جسيم في كل من البعدين والثلاثة أبعاد.

النتائج الرئيسية

• القابلية للتوسع: نجحت خوارزمية SRM المحسنة في توليد تكوينات تصل إلى $10^7$ جسيم. بالنسبة لكسور الحجم المنخفضة والمتوسطة، تحافظ الطريقة على كفاءة عالية.
• التحكم في عدم الاتزان: من خلال ضبط معدل التورم ($c_w$) ومعدل الهجرة ($c_m$)، يمكن للخوارزمية إنتاج طيف من البنى المجهرية يتراوح من التوزيعات الشبيهة بالاتزان إلى الترتيبات المتكتلة بشدة.
  • المقارنة مع الديناميكا الجزيئية (MD): على عكس الديناميكا الجزيئية الموجهة بالأحداث (MD)، والتي تتطور طبيعياً نحو الاتزان، تولد SRM تكوينات ذات توزيعات لكسر الحجم المحلي (LVF) أوسع بكثير وتوقيعات تكتل أقوى. تظل الحالات "المتكتلة" الناتجة عن MD أقرب إلى الاتزان مما تنتجه SRM.
  • حدود كسر الحجم: في البعد الثنائي، يمكن للطريقة تجاوز حد التعبئة العشوائية الوثيقة (RCP) الكلاسيكي (لتصل إلى ~87% للأقراص). في البعد الثلاثي، الحد العملي هو حوالي 60% للكرات؛ وما وراء ذلك، تصبح الهجرة العشوائية غير فعالة، ويُفضل استخدام الديناميكا الجزيئية الموجهة بالأحداث.
• التمدد إلى الشوائب غير الكروية: تم تطبيق الخوارزمية بنجاح على الصفائح الدائرية الرقيقة (spherodisks). ومن خلال تعديل المعاملات، تمكن المؤلفون من توليد:
  • شبكات "بيت الورق" (HoC): هياكل مترابطة للغاية تتشكل بمعدلات تورم عالية ومعدلات هجرة منخفضة.
  • نطاقات شبه نيماتية: مناطق مستقرة جزئياً ومنظمة، تتشكل من خلال إرخاء هياكل "بيت الورق" بمعدلات هجرة عالية.
• تأثير البنية-الخاصية: كشف تحليل التخلل (percolation analysis) لهياكل الصفائح أن تكوينات "بيت الورق" تظهر مسافة تخلل حرجة أقل بنحو عشرة أضعاف من الهياكل شبه النيماتية المسترخية. يوضح هذا أن الترتيب المجهري يؤثر بشكل كبير على خصائص النقل.

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

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