🔬 condensed matter

Shear Banding in Simulations of Polymer Melts

تقارن هذه الورقة بين المحاكاة العددية للمصهورات البوليمرية المبتورة وبين نموذج نظري يربط بين ديناميكيات أنبوب "رولي-بولي" (Rolie-Poly) وتحرر القيد المنقول (Convected Constraint Release)، حيث تجد توافقاً شبه كمي يتنبأ بأن التدرج الطبقي للقص يحدث عندما يتجاوز عدد التشابك عتبة حرجة تعتمد على الصلابة، مع تسليط الضوء على قصور نماذج الأنبوب في السوائل ذات التدفق القوي والمنفصلة جزئياً.

المؤلفون الأصليون: Lucas L. Nelson, Gary S. Grest, Peter D. Olmsted

نُشر 2026-04-07

📖 5 دقيقة قراءة🧠 قراءة متعمّقة

CC BY 4.0

المؤلفون الأصليون: Lucas L. Nelson, Gary S. Grest, Peter D. Olmsted

البحث الأصلي مرخَّص بموجب CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ هذا شرح مولَّده بالذكاء الاصطناعي للبحث أدناه. لم يكتبه المؤلفون ولم يصادقوا عليه. وللتحقق من الدقة التقنية، يرجى الرجوع إلى البحث الأصلي. اقرأ إخلاء المسؤولية الكامل

الصورة الكبيرة: سباغيتي في خلاط

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

منطقة سريعة: حيث تندفع المعكرونة بجنون حول بعضها.
منطقة بطيئة: حيث تكاد المعكرونة لا تتحرك.

هذه الظاهرة تسمى "التنديد بالقص" (Shear Banding). الأمر يشبه قول السائل: "لا أستطيع التعامل مع هذه السرعة بشكل موحد، لذا سأنقسم إلى مسار سريع ومسار بطيء".

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

الشخصيات الرئيسية

لفهم هذه الدراسة، نحتاج للتعرف على "الممثلين" في هذه الدراما:

البوليمرات (السباغيتي): جزيئات طويلة تشبه السلاسل. في المصهور (مثل البلاستيك المنصهر)، تكون جميعها متشابكة مع بعضها البعض، مثل وعاء من السباغيتي تُرِك لبعض الوقت.
التشابكات (العُقد): النقاط التي تعلق فيها السلاسل ببعضها البعض. كلما زاد تشابكها، زادت صعوبة تحريكها.
القص (التحريك): القوة المطبقة لجعل السائل يتدفق.
"الأنبوب" (الممر): في الفيزياء، نتخيل أن كل خيط سباغيتي محبوس في ممر مصنوع من جيرانه. يمكنه التمايل للأمام والخلف (الانزلاق أو الـ reptation)، لكن لا يمكنه القفز بسهما بسهولة إلى الجوانب.
CCR (بطاقة "الخروج من السجن مجانًا"): يرمز هذا إلى "تحرير القيود الناجم عن النسخ" (Convected Constraint Release). تخيل أنه عندما تحرك السباغيتي بسرعة، يتم تدمير أو إعادة ترتيب "الممرات" (الأنابيب) لأن الخيوط تُسحب بقوة شديدة. هذا يحرر الخيوط لتتحرك بشكل أسرع. تدرس الورقة معاملًا محددًا يسمى $\beta$ ، والذي يتحكم في مدى سرعة تدمير هذه الممرات.

التجربة: مختبر رقمي

لم يستخدم الباحثون وعاءً حقيقيًا من السباغيتي، بل بنوا عالمًا افتراضيًا داخل سوبر كمبيوتر (باستخدام برنامج يسمى LAMMPS).

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

استخدموا طريقتين للتحريك:

"المعلم الصارم" (SLLOD): هذه الطريقة تجبر السائل على التحرك بسرعة موحدة تمامًا. إنها تشبه معلمًا يجبر كل طالب على المشي بنفس الوتيرة تمامًا. هذه الطريقة تمنع التنديد، لذا استخدموها لقياس كيفية تفاعل السلاسل مع الإجهاد.
"المجال الحر" (DPD): هذه الطريقة تسمح للسائل بالتحرك كما يشاء. إذا أراد السائل الانقسام إلى مسارات سريعة وبطيئة، فإنه يُسمح له بذلك. وهنا شاهدوا "التنديد" فعليًا.

الاكتشاف الرئيسي: عتبة "التشابك"

قارن الباحثون نتائج الكمبيوتر الخاصة بهم بنموذج رياضي يسمى نموذج DO (نسبة إلى دولاتا وأولمستيد). يحاول هذا النموذج التنبؤ متى سينقسم السائل إلى حزم (Bands).

التنبؤ:
يقول النموذج: "إذا كانت السباغيتي الخاصة بك متشابكة بما يكفي (عدد كبير من العقد، يسمى $Z$ )، ولم يتم تدمير 'الممرات' بسرعة كبيرة (قيمة $\beta$ منخفضة)، فسيحدث التنديد".

النتيجة:
اتفقت محاكاة الكمبيوتر مع النموذج!

قصيرة/تشابك قليل: تدفق السائل بسلاسة. لا يوجد تنديد.
طويلة/تشابك كثير: انقسم السائل إلى مسارات سريعة وبطيئة. حدث التنديد.

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

المفاجأة: لماذا قد لا تندد البلاستيكيات الحقيقية؟

هذا هو الجزء الأكثر إثارة للاهتمام. يتنبأ النموذج بأن التنديد يجب أن يحدث في العديد من البلاستيكات الحقيقية. ولكن في الحياة الواقعية، نادرًا ما نرى ذلك في البلاستيك الصلب (المصاهير). لماذا؟

تقترح الورقة أن السبب هو الصلابة.

التشبيه: تخيل الفرق بين قطعة معكرونة لينة ومطبوخة، وبين عود سباغيتي صلب وجاف.
النتيجة: السلاسل الأكثر صلابة (مثل تلك الموجود، في العديد من البلاستيكات الحقيقية) تمتلك قيمة $\beta$ أعلى. وهذا يعني أن "ممراتها" تُدمر بسرعة كبيرة عند تحريكها.
النتيجة المترتبة: عندما تُدمر الممرات بسرعة، يمكن للسلاسل أن تسترخي وتتحرك بسهولة دون أن تعلق في "ازدحام مروري". هذا يجعل التدفق سلسًا ويمنع الانقسام إلى حزم.

يشير المؤلفون إلى أن البوليمرات الحقيقية (مثل البولي إيثيلين) صلبة بما يكفي لـ "الهروب" من حالة التنديد. ومع ذلك، فإن البوليمرات المرنة جدًا (مثل السيليكون، أو PDMS) قد تظهر بالفعل سلوك التنديد هذا إذا نظرنا إليه بدقة.

مشكلة "التاريخ": الأمر يعتمد على كيفية البدء

وجدت الورقة أيضًا أن التنديد يعتمد على التاريخ.

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

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

الملخص: ماذا يعني هذا بالنسبة لنا؟

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

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

إليك ملخص تقني مفصل للورقة البحثية بعنوان "التندب القصي في محاكاة مصهورات البوليمرات" (Shear Banding in Simulations of Polymer Melts) لـ نيلسون، غريست، وأولمستيد.

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

تتناول الورقة البحثية التناقض طويل الأمد بين التوقعات النظرية والملاحظات التجريبية فيما يتعلق بظاهرة التندب القصي (Shear Banding) في مصهورات البوليمرات المتشابكة.

السياق النظري: تتوقع نماذج الأنبوب الكلاسيكية (مثل نموذج Doi-Edwards) أنه تحت معدلات قص عالية، تصبح العلاقة بين إجهاد القص ( $\sigma_{xy}$ ) ومعدل القص ( $\dot{\gamma}$ ) غير رتيبة (non-monotonic). تؤدي هذه اللا-رتابة إلى عدم استقرار في التدفق يُعرف بالتندب القصي، حيث ينفصل السائل إلى مناطق ذات معدلات قص عالية ومنخفضة.
الصراع: بينما يُلاحظ التندب القصي بشكل متكرر في محاليل البوليمرات وفي محاكاة الديناميكا الجزيئية (MD) للمصهورات، فإن الأدلة التجريبية على التندب القصي في حالة الاستقرار لمصهورات البوليمرات شحيحة أو مثيرة للجدل. غالبًا ما تُظهر التجارب زيادة طفيفة في الإجهاد بدلاً من الهضبة (plateau)، وغالبًا ما تختلط الظواهر مع انزلاق الجدران وكسر الحواف.
الفجوة: هناك حاجة للتوفيق بين أدلة المحاكاة الواضحة على التندب وبين عدم الملاحظة التجريبية، وتحديداً من خلال فهم دور تحرير القيود الناجم عن الحمل (Convected Constraint Release - CCR) وكيف تؤثر صلابة البوليمر على الظروف الحرجة للتندب.

2. المنهجية

استخدم المؤلفون نهجاً متعدد الأوجه يجمع بين النمذجة المستمرة ومحاكاة الديناميكا الجزيئية واسعة النطاق.

أ. الإطار النظري (نموذج DO)

استخدموا نموذج دولاتا-أولمستيد (DO)، وهو امتداد متسق ديناميكياً حرارياً لنموذج Rolie-Poly (RP).
الآلية الرئيسية: يربط النموذج تطور موتر التشكيل للأنبوب ( $A$ ) مع العدد المختزل للتشابكات ( $\nu$ ). وهو يتضمن CCR، حيث يؤدي تمدد البوليمر الناتج عن التدفق إلى تسريع تحرير التشابكات.
المعلمات (Parameters): يحكم النموذج المعلمات التالية:
- $Z$ : عدد التشابك عند التوازن.
  
  $\beta$ : معلمة CCR (قوة التفكك الناتج عن التدفق).
- $\alpha$ : معلمة تباين الاسترخاء (تتحكم في فروق الإجهاد الطبيعي).
- $\lambda_{max}$ : نسبة التمدد القصوى.
التوقع: يتنبأ النموذج بوجود "مخطط طور" (phase diagram) حيث يحدث التندب القصي فقط إذا تجاوز $Z$ قيمة حرجة $Z_c$ ، والتي تعتمد على $\beta$ و $\alpha$ . وتحديداً، يُتوقع أن $\beta$ الأعلى (CCR أقوى) تعمل على تثبيط التندب من خلال القضاء على اللا-رتابة في منحنى الإجهاد-الانفعال.

ب. محاكاة الديناميكا الجزيئية

النظام: بوليمرات (Kremer-Grest - KG) مكونة من خرز ونابض شبه مرنة، تمت محاكاتها باستخدام برنامج LAMEDPS.
المتغيرات:
- أطوال السلاسل ( $N_{mon}$ ) من 200 إلى 2000 خرزة.
- صلابة الانحناء ( $k_\theta$ ) تم تغييرها (0، 1.5، 2.0) لضبط طول التشابك والصلابة.
- بروتوكولات القص: تم استخدام طريقتين:
  1. SLLOD: تفرض ملف سرعة موحد (تمنع التندب، وتُستخدم لقياس الخصائص التكوينية).
  2. Fix Deform (DPD Thermostat): تسمح لملف السرعة بالتطور طبيعياً، مما يسمح بحدوث التندب القصي.
التحليل:
- قياس التشابك: تم استخدام كود $Z1+$ لحساب عدد عُقد المسار الأولي ( $Z_k$ ) لتحديد نسبة التفكك $\nu = Z_k/Z_{k,eq}$ .
- توفيق المعلمات: تم تحديد معلمة CCR ( $\beta$ ) عن طريق مطابقة بيانات المحاكاة لـ $\nu$ مقابل رقم وايزنبرغ راوس ( $Wi_R$ ) مع توقعات نموذج DO.
- ملاحظة التندب: تم تحليل ملفات السرعة للكشف عن تشكل نطاقات (bands) سريعة وبطيئة متميزة.

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

تكميم $\beta$ : نجحت الدراسة في ربط معلمة CCR ( $\beta$ ) بصلابة البوليمر (المقاسة بنسبة طول كون (Kuhn length) إلى طول التعبئة $b_K/p$ ). ووجدوا أن $\beta$ تزداد مع الصلابة، متبعة قانون القوة $\beta \propto (b_K/p)^{1.6}$ .
التحقق من صحة نموذج DO: أكدت المحاكاة أن نموذج DO يتنبأ بشكل صحيح بوجود عدد تشابك حرج ( $Z_c$ ) يحدث فوقه التندب.
تحديد التناقضات: حدد المؤلفون مجالات محددة يبتعد فيها نموذج DO المبسط عن واقع المحاكاة، لا سيما فيما يتعلق باعتماد $Z_c$ على $\beta$ وانهيار النماذج الأنبوبية عند معدلات القص العالية.
الاعتماد على التاريخ: سلطت الدراسة الضوء على الطبيعة شبه المستقرة (metastable) للتندب القصي في المحاكاة، حيث أظهرت أن تاريخ القص (البدء من السكون مقابل خفض المعدل من قص عالٍ) يؤثر بشكل كبير على حدوث التندب.

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

الاتفاق على وجود التندب: تُظهر المحاكاة تندباً قصياً مستقراً واضحاً للسلاسل شديدة التشابك ( $Z > Z_c$ ). وتتفق منطقة التندب في فضاء المعلمات $(Z, \beta)$ شبه كمياً مع توقعات نموذج DO.
التناقض بين $\beta$ و $Z_c$ :
- توقع النموذج: من المفترض أن $\beta$ الأعلى (CCR أقوى) يجب أن تزيد من عدد التشابك الحرج $Z_c$ ، مما يجعل التندب أصعب.
- ملاحظة المحاكاة: أظهرت المحاكاة اتجاهاً طفيفاً حيث ينخفض $Z_c$ مع زيادة $\beta$ (السلاسل الأكثر صلابة ذات الـ $\beta$ الأعلى تندبت عند $Z$ أقل). وهذا يتناقض مع آلية النموذج التي تفترض أن CCR يثبط التندب.
الاعتماد على معدل القص: لوحظ التندب في نافذة متوسطة من معدلات القص. عند معدلات القص العالية جداً، أصبح النظام موحداً مرة أخرى (رب-ما بسبب التفكك الكامل/انهيار الأنبوب)، وعند المعدلات المنخفضة جداً، كان أيضاً موحداً.
نسبة الإجهاد الطبيعي ( $-N_2/N_1$ ): تشير المحاكاة والبيانات التجريبية لمصهورات البوليمر إلى أن نسبة الإجهاد الطبيعي أعلى في المصهورات منها في المحاليل، مما يعني وجود $\alpha$ (تباين استرخاء) أعلى (غالباً $\alpha > 0.5$ ). ومع ذلك، فإن قياس حد النظام الخطي صعب، مما يترك $\alpha$ غير مؤكدة نوعاً ما.
تأثيرات الحجم المحدود: لاحظ المؤلفون أن أحجام صندوق المحاكاة المحدودة يمكن أن تحجب التندب، خاصة عندما يصبح عرض نطاق القص العالي مقارباً للحجم الجزيئي (بضع مونومرات).

5. الأهمية والآثار المترتبة

تفسير التندب "المفقود" تجريبياً: تشير النتائج إلى أن البوليمرات الفيزيائية ذات الهياكل المجهرية المحددة (مثل البولي إيثيلين الموحد الذرات مع زوايا ثنائية السطح) قد تمتلك قيم $\beta$ عالية جداً. إذا كان هذا الاتجاه صحيحاً بأن $\beta$ العالية ترفع $Z_c$ بشكل كبير، فإن العديد من مصهورات البما بوليمر التجريبية قد يكون لديها ببساطة $Z < Z_c$ في الظروف العادية، مما يفسر ندرة ملاحظة التندب في التجارب مقارنة بالمحاكاة.
محدودية النماذج: تسلط الدراسة الضوء على أنه بينما يلتقط نموذج DO الفيزياء النوعية، فإنه يفشل كمياً عند معدلات القص العالية حيث ينهار مفهوم "الأنبوب" بسبب التفكك الكبير. إن افتراض نمط استرخاء واحد والشكل المحدد لـ CCR قد يكون غير كافٍ للتدفقات غير الخطية القوية.
الاتجاهات المستقبلية: يشير العمل إلى أنه لملاحظة التندب القصي في المصهورات الفيزيائية، يجب البحث عن بوليمرات مرنة ذات صلابة منخفضة (low $\beta$ ) وكثافة تشابك عالية ( $Z$ ). كما تؤكد الدراسة على الحاجة إلى قياسات أكثر دقة لنسبة الإجهاد الطبيعي في النظام الخطي لتحديد معلمة $\alpha$ في النماذج النظرية.

باختصار، توفر هذه الورقة جسراً صارماً بين المحاكاة المجهرية وعلم الريولوجيا المستمر، مؤكدة أن التندب القصي هو ظاهرة حقيقية في مصهورات البوليمرات تحت ظروف معينة، بينما تكشف في الوقت نفسه عن حدود النماذج الأنبوبية الحالية في وصف التفاعل المعقد بين الصلابة، والتشابك، وتاريخ التدفق.