A Scalable OpenLB LAMMPS Framework for Fully Resolved Simulations of Hindered Settling of Arbitrary Non-Spherical Particles
تقدم هذه الورقة إطار عمل OpenLB-LAMMPS قابل للتوسع يربط بين طرق "لاتيس بولتزمان" (Lattice Boltzmann) و"العناصر المنفصلة" (Discrete Element) لمحاكاة الترسيب المعاق المكتمل الدقة لما يصل إلى 100,000 جسيم غير كروي عشوائي، مما يكشف أن الجسيمات المكعبة تظهر ديناميكيات غلاف تنسيق واختلافات في ارتباطات السرعة تعتمد على حجم النظام مقارنة بالجسيمات الكروية.
المؤلفون الأصليون:Varghese Babu, Adrian Kummerländer, Mathias J. Krause, Santosh Ansumali, Dilip Gersappe
تخيل دلوًا من الرمل يسكب في الماء. لا تسقط الحبيبات بنفس السرعة التي كانت ستسقط بها لو كانت بمفردها؛ فهي تبطئ بعضها البعض، وتتصادم وتدفع السائل في رقصة معقدة ومزدحمة. هذه الظاهرة، المعروفة باسم "الترسب المعاق" (hindered settling)، هي عملية أساسية في الطبيعة والصناعة، تحكم كل شيء بدءًا من كيفية بناء الرواسب لدلتا الأنهار وصولاً إلى كيفية خلط الأدوية في المصانع. لعقود من الزمن، فهم العلماء هذا السلوك جيدًا عندما تكون الجسيمات كرات مثالية، مثل الكرات الزجاجية الصغيرة. ومع ذلك، فإن العالم الحقيقي نادرًا ما يتكون من كرات مثالية؛ فمعظم الجسيمات غير منتظمة، أو مسننة، أو مربعة الشكل، وتخلق أشكالها الغريبة طرقًا فريدة للتلامس والانزلاق فوق بعضها البعض. لقد كان التنبؤ بكيفية ترسب هذه الأشكال الغريبة تحديًا حوسبيًا مستعصيًا، لأن رياضيات تصادماتها والطريقة التي تسحب بها السائل حولها صعبة الحل في آن واحد.
ولمعالجة ذلك، بنى فريق من الباحثين إطارًا رقميًا جديدًا يربط بين أداتين برمجيتين مفتوحتي المصدر لمحاكاة هذه التفاعلات المعقدة. فقد دمجوا برنامج حل السوائل المسمى OpenLB، الذي ينمذج كيفية تحرك السوائل، مع برنامج محاكاة الجسيمات المسمى LAMMPS، الذي يتتبع كيفية اصطدام الأجسام الصلبة واحتكاكها ببعضها البعض. ويكمن الابتكار الرئيسي في كيفية تمثيل الأجسام الصلبة؛ فبدلاً من محاولة نمذجة شكل معقد مثل المكعب باستخدام تعريف رياضي واحد صلب، قاموا بتفكيك الجسم إلى مجموعة من الكرات المتداخلة، تمامًا مثل رجل ثلج مكون من ثلاث كرات. تسمح هذه "الكتلة" للكمبيوتر بحساب كيفية تلامس الجسم مع جيرانه بسهولة. وفي الوقت نفسه، يرى برنامج حل السوائل الجسم ليس كمجموعة من الكرات، بل ككتلة صلبة ذات شكل محدد، باستخدام تقنية تسمى "التحويل إلى فوكسل" (voxelization) لرسم سطح الجسم على شبكة السائل. يسمح هذا الفصل للفريق بمحاكاة ميكانيكا تلامس الجسيمات وتدفق السائل بدقة عالية، حتى عند التعامل مع أعداد هائلة من الأجسام.
اختبر الباحثون نظامهم من خلال محاكاة ترسب مكعبات في سائل، ومقارنة النتائج مع السلوكيات المعروفة للكرات. أجروا عمليات محاكاة تضم ما يصل إلى 103,823 مكعبًا، وهو مقيء يدفع حدود القدرة الحوسبية الحالية. وكشفت نتائجهم عن اختلاف واضح بين سلوك المكعبات والكرات أثناء الترسب. فبينما تميل الكرات إلى التكتل معًا من خلال التكتل الناتج عن سيطرة التلامس، شكلت المكعبات أغلفة تنسيق واضحة دون تلامس متوازي الأوجه. يشير هذا إلى أن الحواف والزوايا الحادة للمكعبات تخلق نوعًا معينًا من التباعد يمنعها من التراص بنفس كثافة الكرات. علاوة على ذلك، اكتشف الفريق أن تقلبات سرعة الجسيمات لم تكن عشوائية أو محلية؛ بل إن حركة جسيم واحد أثرت على حركة الآخرين عبر النظام بأكمله. وحتى في أكبر عمليات المحاكاة، امتد الارتباط بين سرعات الجسيمات ليشمل الحاوية بأكملة، مما يشير إلى أن الحركة الجماعية للحشد هي ظاهرة موحدة وليست مجرد سلسلة من الأحداث المنعزلة.
من خلال التحقق من صحة نموذجهم مقابل التجارب على الكرات، ثم تطبيق ذلك على المكعبات، أثبت المؤلفون أن شكل الجسيم يلعب دورًا حاسمًا في ديناميكيات التعلق. وأظهرت عمليات المحاكاة أن سرعة ترسب المكعبات تختلف عن الكرات، وأن البنية الداخلية لسحابة الترسب مختلفة جوهريًا. ويؤكد هذا العمل أن تجاهل شكل الجسيمات يمكن أن يؤدي إلى تنبؤات غير دقيقة في مجالات تتراوح من الهندسة الجيوتقنية إلى علوم البيئة. يوفر إطار عمل الفريق أداة قوية لدراسة هذه التدفقات واسعة النطاق، مما يوفر وسيلة لرؤية النظام الخفي داخل فوضى الجسيمات المترسبة. ورغم أن عمليات المحاكاة الحالية لم تتضمن الطبقة الرقيقة من السائل التي تعمل كمزلق للفجوة بين الجسيمات عندما تكون قريبة جدًا من بعضها، إلا أن الطريقة نجحت في التقاط الديناميكيات الأوسع للنظام. وتشير النتائج إلى أنه بالنسبة للجسيمات غير الكروية، فإن طريقة ترسبها تحكمها علاقة معقدة بين الشكل، والتلامس، وتفاعلات السويد طويلة المدى التي لا تستطيع النماذج البسيطة للكرات استيعابها.
ملخص تقني: إطار عمل OpenLB–LAMMPS قابل للتوسع للمحاكاة كاملة الدقة للترسيب المعاق لجسيمات عشوائية غير كروية
بيان المشكلة يُعد الترسيب المعاق (hindered settling) ظاهرة حيث تؤدي التفاعلات بين الجسيمات والجسيمات، وبين الجسيمات والسوائل، إلى إبطاء السرعة الهابطة للمعلقات ذات التركيز العالي، وهي ظاهرة أقل فهمًا للجسيمات غير الكروية مقارنة بالجسيمات الكروية. وبينما تصف علاقة ريتشاردسون-زاكي (Richardson-Zaki) المعلقات الكروية بفعالية، تواجه النماذج العددية للأشكال غير الكروية تحديات حوسبية كبيرة. تنشأ هذه الصعوبات أساسًا من ضرورة حل ميكانيكا التلامس المعقدة (التعاشق، والديناميكيات المعتمدة على التوجه) والتفاعلات الهيدروديناميكية في آن واحد. غالبًا ما تعاني أطر عمل ديناميكا السوائل الحسابية وطريقة العناصر المنفصلة (CFD–DEM) الحالية من صعوبة في محاكة الأنظمة واسعة النطاق (O(105) جسيم) ذات الأشكال الصلبة العشوائية بسبب القيود في اكتشاف التلامس ودقة التفاعل بين السائل والمنشأ (FSI).
المنهجية يقدم المؤلفون إطار عمل (FSI) قابل للتوسع وكامل الدقة يربط بين حزمتين من البرمجيات مفتوحة المصدر: OpenLB (محلل طريقة لاتيس بولتزمان - LBM) و LAMMPS (محلل طريقة العناصر المنفصلة - DEM).
محلل السائل (OpenLB): يتم نمذجة الطور السائل باستخدام طريقة لاتيس بولتزمان (LBM) مع نهج متجانس يعتمد على معادلات برينكمان-نافييه-ستوكس (BNSE). يتم تمثيل الأجسام الصلبة العشوائية على شبكة أوغلر (Eulerian lattice) ثابتة عبر مجال مسامية يعتمد على الزمن، d(x,t). ويتم فرض شرط عدم الانزلاق من خلال حد جزاء الحجم (سحب برينكمان). ويتم حساب القوى والعزوم الهيدروديناميكية عبر تبادل الزخم عبر الواجهات بين السائل والصلب.
محلل الجسيمات (LAMMPS): يتم نمذجة الأجسام الصلبة باستخدام طريقة العناصر المنفصلة (DEM). وللتعامل مع الأشكال العشوائية، يستخدم المؤلفون تمثيل "الكتل" (clump representation)، حيث يتم تفكيك الهندسات المعقدة (مثل المكعبات) إلى كرات متداخلة. يُستخدم هذا التمثيل حصريًا لاكتشاف التلامس وحساب القوة داخل LAMMPS.
استراتيجية الربط: يعمل إطار العمل وفق مخطط خطوة زمنية متسلسلة. يقوم محلل السائل بحساب القوى الهيدروديناميكية بناءً على مواضع وتوجهات الجسيمات الحالية، والتي تُطبق بعد ذلك كقوى خارجية في محاكاة DEM. وعلى العكس من ذلك، يتم رسم خرائط حالات الجسيمات المحدثة (المواضع، السرعات، السرعات الزاوية) مرة أخرى إلى شبكة السائل كجال مسامية.
والأهم من ذلك، أن محلل السائل يحل الهندسة السطحية الدقيقة للجسيمات عبر عملية "الفوكسلة" (voxelization)، بشكل مستقل عن تقريب "الكتل" متعدد الكرات المستخدم للتلامسات في LAMMPS. وهذا يميز هذا العمل عن الطرق السابقة حيث كان السائل "يرى" أيضًا تقريب الكتل، مما يتجنب أي أخطاء هندسية إضافية في التفاعل الهيدروديناميكي.
إعداد المحاكاة: تركز الدراسة على الترسيب المعاق في نطاق ثلاثي التكرار (triply-periodic domain). ولمنع التسارع اللانهائي الناتج عن الجاذبية، يتم تطبيق قوة تعويضية على عقد السائل للحفاظ على متوسط سرعة سائل يساوي صفرًا.
المساهمات الرئيسية
إطار عمل قابل للتوسع: يثبت المؤلفون القدرة على محاكاة أنظمة تحتوي على ما يصل إلى O(105) من الأجسام الصلبة العشوائية (تحديدًا المكعبات)، وهو ما يزيد بمقدار مرتبة عشرية عن العروض السابقة مفتوحة المصدر (مثل openHFDIB-DEM).
التمثيل الهجين: يسمح دمج الفوكسلة السطحية الدقيقة لحل السائل مع كتل متعددة الكرات لميكانيكا التلامس بمحاكة فعالة للهندسات المعقدة، مثل المكعبات، دون التضحية بالدقة الهيدروديناميكية.
التحقق من الصحة: تمت عملية التحقق من إطار العمل مقابل المعايير التجريبية والتحليلية، بما في ذلك ترسيب الكرة المنفردة (Ten Cate et al.)، وتصادمات الكرة والجدار، والسحب على مصفوفات دورية من الكرات، وسرعة ترسيب الجسيمات المكعبة (Kwon et al.).
النتائج
ديناميكيات الترسيب: تكشف عمليات محاكاة الترسيب المعاق للمكعبات والكرات متكافئة الحجم عن اختلافات واضحة. فبينما تميل الكرات إلى تشكيل عناقيد يهيمن عليها التلامس، تشكل المكعبات أغلفة تنسيق بارزة بدون تلامس موازٍ للوجوه، مما يشير إلى أن مورفولوجيا الجسيمات تغير البنية المجهرية للمعلق بشكل كبير.
دوال التوزيع الشعاعي (RDF): تظهر دالة التوزيع الشعاعي للمكعبات أغلفة أولى وثانية متميزة، بينما تظهر الكرات ميلًا للمس بعضها البعض مباشرة، مع زيادة عدد التلامسات مع زيادة كسر التعبئة.
ارتباطات السرعة: تدرس الدراسة الارتباط المكاني لتقلبات السرعة. وتتمثل النتيجة الرئيسية في أن طول ارتباط تقلبات السرعة يتناسب مع حجم النظام، حتى بالنسبة لأكبر الأنظمة المدروسة (∼105 جسيم). وخلافًا لبعض الملاحظات التجريبية حيث تتشبع التقلبات عند طول مقياس معين، تظهر المحاكاة ارتباطات تمتد عبر النطاق بأكمله.
الاعتماد على حجم النظام: بينما تظل متوسط السرعة النهائية وترتيب الجسيمات المحلي مستقرين عبر أحجام الأنظمة المختلفة (من ∼3,000 إلى ∼100,000 مكعب)، فإن الطبيعة طويلة المدى لارتباطات السرعة تستمر، مما يشير إلى أن الشروط الحدودية الدورية أو آلية قوة التوازن المحددة قد تؤثر على الارتباطات المرصودة.
الأهمية والادعاءات يزعم البحث توفير إطار عمل قوي ومفتوح المصدر لمحاكاة التدفقات الجيوتقنية والجسيمية واسعة النطاق التي تتضمن جسيمات غير كروية عشوائية. ومن خلال الربط الناجح بين OpenLB و LAMMPS، يردم المؤلفون الفجوة بين ميكانيكا التلامس الدقيقة والهيدروديناميكا كاملة الدقة للجسيمات غير الكروية. تسلط النتائج الضوء على الدور الحاسم لمورفولوجيا الجسيمات في ديناميكيات المعلقات، وتظهر تحديدًا أن المكعبات تظهر سلوكيات تكتل مختلفة مقارنة بالكرات، بما في ذلك ارتباطات تقلبات السرعة. يشير المؤلفون بتواضع إلى أنه بينما يعد إطار العمل قابلًا للتوسع، إلا أن هناك حاجة لعمل مستقبلي لمعالجة ارتباطات السرعة طويلة المدى المرصودة (والتي قد تكون ناتجة عن آثار الحدود) ولإدراج قوى التشحيم (lubrication forces) للتفاعلات الأكثر قربًا بين الجسيمات. يضع هذا العمل أساسًا لدراسة الأنظمة الجسيمية المعقدة بمقاييس كانت غير قابلة للوصول سابقًا لطرق CFD-DEM كاملة الدقة.