Entropy-Audited Thermoviscoelastic Modelling of Polymer Relaxation Under Non-Isothermal Loading: A Reproducible Numerical Study

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

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

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

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

لم يكن الفرق بين النموذجين خطأً صغيراً أو طفيفاً؛ ففي نهاية الفترة التي تم فيها تثبيت المادة عند درجة الحرارة العالية، توقع النموذج المبسط إجهاداً أعلى بنسبة 68.0 بالمائة مما توقعه النموذج الأكثر دقة والمترابط كلياً. وهذا يعني أن النسخة المبسطة ظنت أن المادة تحتفظ بما يقرب من ضعف القوة التي كانت تحتفظ بها في الواقع. والأكثر أهمية من ذلك، أن هذا الخطأ لم يختفِ عندما بردت المادة مرة أخرى. فبعد عودة درجة الحرارة إلى نقطة البداية والسماح للمادة بالاستقرار لمدة 150 ثانية أخرى، ظل النموذج المبسط يتوقع مستوى إجهاد أعلى بنسبة 28.8 بالمائة. لقد تعرضت المادة لـ "خداع" من قبل الحرارة؛ فقد استرخت أكثر مما ينبغي خلال المرحلة الساخنة، ولم يستطع النموذج المبسط، الذي فشل في رؤية حدوث ذلك، تصحيح تنبؤه لاحقاً. لقد غير تاريخ الحرارة الحالة الداخلية للمادة بشكل دائم، والنموذج الذي تجاهل سرعة هذا التغيير لم يستطع التعافي.

بعيداً عن مجرد قياس الإجهاد، قدمت الدراسة طبقة جديدة من التحقق تسمى "تدقيق الإنتروبيا" (entropy audit). في الفيزياء، ينص القانون الثاني للديناميكا الحرارية على أن الطاقة في أي نظام يجب أن تتبدد دائماً أو تُفقد كحرارة أثناء العمليات غير العكسية مثل الاسترخاء؛ ولا يمكن أبداً خلقها من العدم. قام الباحثون ببناء محاكاتهم بحيث يمكنها حساب مقدار الطاقة التي تتبدد بدقة في كل لحظة. ووجدوا أن النموذج الدقيق أظهر باستمرار كمية موجبة من الطاقة المفقودة، وهو أمر صحيح فيزيائياً. أما النموذج المبسط، فرغم أنه قد يبدو جيداً في الظاهر، إلا أنه فشل في التقاط تدفق الطاقة الحقيقي. ومن خلال تتبع هذه الطاقة، استطاع الباحثون إثبات أن نموذجهم المعقد لم يكن مجرد تخمين، بل كان متسقاً رياضياً مع قوانين الفيزياء. وعلى مدار دورة الـ 800 ثانية بأكملها، أظهرت المحاكاة الدقيقة أن المادة تراكمت لديها 57.01 كيلوجول من الطاقة المتبددة لكل متر مكعب، وهو رقم دقيق يعمل كبصمة لسلوك المادة.

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

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

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