Numerical simulation of dilute polymeric fluids with memory effects in the turbulent flow regime

ملخص تقني: المحاكاة العددية للسوائل البوليمرية المخففة ذات تأثيرات الذاكرة في نظام التدفق المضطرب

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

1. الأبعاد العالية: تُعرف معادلة فوكر-بلانك (FP) على حاصل الضرب الديكارتي للمجال المكاني العياني ($\Omega \subset \mathbb{R}^d$) والمجال التكويني المجهري ($D \subset \mathbb{R}^d$)، مما ينتج عنه مسألة ذات أبعاد $2d$.
2. اللا-محلية في الزمن: يُدخل المشتق الكسري الزمني اعتماداً على التاريخ، مما يتطلب تخزين ومعالجة تاريخ الحل بأكمله، وهو أمر مكلف حاسوبياً.
3. دقة الاضطراب: يتطلب تحليل التدفق المضطرب دقة عالية في المكان والزمان، مما يضاعف التكلفة الحسابية للنظام عالي الأبعاد وغير المحلي.

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

المنهجية
يقترح المؤلفون مخططاً عددياً عيانياً بحتاً يختزل نظام (micro-macro) عالي الأبعاد إلى نظام من المعادلات التفاضلية الجزئية (PDEs) المعرفة فقط على المجال المكاني. تتكون المنهجية من المكونات الرئيسية التالية:

1. طريقة هيرميت الطيفية (Hermite Spectral Method): للتعامل مع الفضاء التكويني غير المحدود $D = \mathbb{R}^d$ المرتبط بنموذج النابض الهوكيني (Hookean spring)، يقوم المؤلفون بتجزئة معادلة (FP) باستخدام دوال هيرميت الطيفية. يحول هذا المسألة إلى نظام من المعادلات للأنماط الطيفية $\phi_{i_1, \dots, i_d}$.

   • القياس الأمثل (Optimal Scaling): المساهمة النظرية الحاسمة هي إثبات أن اختيار معامل القياس الهيرميتي $a = 1/\sqrt{2}$ يضمن أن موتر الإجهاد الإضافي العياني المستمد من التقريب الطيفي يطابق تماماً موتر الإجهاد الإضافي التحليلي للنظام الميكرو-ماكرو الأساسي. وهذا يلغي خيارات المعاملات التعسفية الموجودة في الأعمال السابقة.
   • القطع (Truncation): تم إثبات أن الأنماط الطيفية من الدرجة 0 و2 هي الوحيدة المطلوبة لحساب موتر الإجهاد الإضافي الناتج عن البوليمر، مما يقلص النظام إلى 7 معادلات فقط في الفضاء ثلاثي الأبعاد (4 في الفضاء ثنائي الأبعاد).

2. طريقة ضغط النواة (Kernel Compression Method): لمعالجة اللا-محلية للمشتق الكسري الزمني، يتم تقريب نواة تكامل ريمان-ليوفيل باستخدام مجموع أوزان من الدوال الأسية بناءً على تقريب (AAA algorithm). يحول هذا المعادلة التفاضلية الكسرية إلى نظام من المعادلات التفاضلية العادية (ODEs) للأنماط "الكسرية" المساعدة، مما يتجنب الحاجة لتخزين التاريخ الكامل.

3. التكامل الزمني والاقتران:

   • يتم حل معادلات نافيير-ستوكس باستخدام طريقة الإسقاط (Chorin-type) مع تكامل زمني من الدرجة الثانية.
   • يتم دمج معادلات (MFP) والأنماط الكسرية (ODEs) باستخدام صيغة الفرق الخلفي (BDF) من الدرجة الثانية.
   • يتم استقراء حدود الاقتران لتحقيق مخطط مقترن بالكامل. يثبت المؤلفون أن هذا النهج يعطي معدلات تقارب مستقلة عن رتبة المشتق الكسري الزمني $\alpha$.

4. التجزئة المكانية: يتم تجزئة نظام المعادلات التفاضلية الجزئية (PDE) الناتج باستخدام عناصر منتهية مستمرة من الدرجة الثانية على شبكات غير منتظمة.

النتائج الرئيسية
قام المؤلفون بتنفيذ المخطط باستخدام مكتبة MFEM وأجروا تجارب عددية مكثفة في بعدين وثلاثة أبعاد:

• التقارب: يحقق المخطط رتب تقارب مثالية. بالنسبة لمعادلة (FP) الكسرية الزمنية، يتلاشى الخطأ بمعدل $O(\Delta t^{1+\alpha})$. بالنسبة لنظام (NSFP) المقترن بالكامل، يظهر موتر الإجهاد الإضافي تقارباً من الدرجة الأولى، بما يتوافق مع الاعتماد على المشتق الكسري الزمني.
• التدفق حول أسطوانة: تظهر محاكاة التدفق حول عائق أن رتبة المشتق الكسري الزمني ($\alpha$) تؤثر بشكل كبير على استقرار التدفق. بينما تظل التدفقات ذات الرتبة الصحيحة ($\alpha=1$) لاميّة (laminar) في ظل الظروف المختبرة، فإن الحالات الكسرية ($\alpha=0.5, 0.8$) تنتقل إلى حالة الاضطراب. تؤدي قيم $\alpha$ الأقل إلى تطور أسرع للاضطراب وقوة بوليمرية أضعف.
• تدفق القناة ذو الجدار الخشن: في هندسة تحاكي خط أنابيب متآكل، تحدد الدراسة ثلاثة أنظمة تدفق بناءً على انتشار مركز الكتلة ($\epsilon$) ومساهمة لزوجة البوليمر: لامي مستقر، اضطراب مستدام، ونظام وسيط حيث يتولد الاضطراب أولاً ثم يخمد لاحقاً.
  • تأثيرات الذاكرة: تكشف المحاكاة أن تأثيرات الذاكرة (المشتقات الكسرية) تضعف قدرة البوليمرات على تقليل السحب. وتحديداً، تكون شدة القوة الناتجة عن البوليمر أقل في الحالات الكسرية مقارنة بالحالات ذات الرتبة الصحيحة، مما يؤدي إلى كفاءة أقل في كبح الاضطراب.
• تخفيف إجهاد الأنابيب (3D): تُظهر محاكاة ثلاثية الأبعاد لتدفق عبر أنبوب بإنعطافات قائمة المدى الذي تصل إليه الطريقة (باستخدام ما يصل إلى 256 نواة). تؤكد النتائج أنه بينما يصبح تدفق المذيب النقي مضطرباً بعد الانعطافات، فإن السائل البوليمري المخفف ذو تأثيرات الذاكرة يحافظ على حالة لاميّة، وإن كان تأثير تقليل السحب يتأثر بمعامل الذاكرة.

الأهمية والادعاءات
تزعم الورقة أنها تقدم أول دراسة عددية صارمة للعوامل المخفضة للسحب ذات تأثيرات الذاكرة في النظام المضطرب. وتكمن أهميتها الأساسية في:

1. الكفاءة: من خلال اشتقاق نموذج عياني بحت يكافئ نظام (micro-macro) الكامل (تحت معامل القياس الأمثل)، تتجاوز الطريقة لعنة الأبعاد المتأصلة في المحاكاة العددية المباشرة لمعادلة (FP).
2. الدقة النظرية: يضمن إثبات معامل القياس الأمثل $a = 1/\sqrt{2}$ أن النموذج العياني ليس تقريباً من حيث التكوين، بل هو تمثيل دقيق للاقتران الإجهادي لنموذج هوك.
3. الرؤية الفيزيائية: تتحدى النتائج العددية الافتراض القائل بأن تأثيرات الذاكرة تعزز دائماً تقليل السحب؛ بل تظهر أن تأثيرات الذاكرة يمكن أن تضعف القوة الناتجة عن البوليمر وتقلل من كفاءة تقليل السحب في التدفقات المضطربة.

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