2025-11-19T11:46:14.009582

Mixing in viscoelastic fluids using elastic turbulence

van Buel, Stark
We investigate the influence of elastic turbulence on mixing of a scalar concentration field within a viscoelastic fluid in a two-dimensional Taylor-Couette geometry using numerical solutions of the Oldroyd-B model. The flow state is determined through the secondary-flow order parameter indicating the regime of elastic turbulence. When starting in the turbulent state and subsequently lowering the Weissenberg number, a weakly-chaotic flow occurs below $\mathrm{Wi_c}$. Advection in both the turbulent and weakly-chaotic flow states induces mixing, which we illustrate by the time evolution of the standard deviation of the solute concentration from the uniform distribution. In particular, in the elastic turbulent state mixing is strong and we quantify it by the mixing rate, the mixing time, and the mixing efficiency. All three quantities follow scaling laws. Importantly, we show that the order parameter is strongly correlated to the mixing rate and hence is also a good indication of mixing within the fluid.
academic

الخلط في السوائل اللزجة المرنة باستخدام الاضطراب المرن

المعلومات الأساسية

  • معرّف الورقة: 2409.06391
  • العنوان: الخلط في السوائل اللزجة المرنة باستخدام الاضطراب المرن
  • المؤلفون: Reinier van Buel, Holger Stark (جامعة برلين التقنية)
  • التصنيف: physics.flu-dyn
  • تاريخ النشر: سبتمبر 2024
  • رابط الورقة: https://arxiv.org/abs/2409.06391

الملخص

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

خلفية البحث والدافع

خلفية المشكلة

  1. تحديات الخلط عند أرقام Reynolds المنخفضة: في السوائل النيوتونية، تكون التدفقات طبقية عند ظروف Reynolds المنخفضة، ويكون الخلط بطيئاً ويسيطر عليه الانتشار بشكل أساسي. يشكل هذا قيداً مهماً في أجهزة الميكروفلويديكس.
  2. اكتشاف الاضطراب المرن: إضافة البوليمرات عالية الوزن الجزيئي إلى السائل يجعله لزجاً مرناً، مما يسمح بظهور ظاهرة الاضطراب المرن حتى عندما تكون القوى الخاملة مهملة.
  3. أهمية الخلط على المستوى الدقيق: في أجهزة المختبر على الرقاقة (lab-on-a-chip)، يعتبر الخلط الفعال للسوائل حاسماً، خاصة على مستوى الميكرومتر.

دافع البحث

  • استكشاف إمكانات الاضطراب المرن في تحسين كفاءة الخلط
  • فهم الطبيعة الفيزيائية لآليات الخلط في السوائل اللزجة المرنة
  • توفير إرشادات نظرية لتطبيقات الميكروفلويديكس

المساهمات الأساسية

  1. تحديد طبيعة الانتقال: تحديد الانتقال من التدفق الطبقي إلى الاضطراب المرن كانتقال حرج فرعي، واكتشاف وجود حالة تدفق فوضوية ضعيفة
  2. إنشاء طرق تحديد الخلط: إدخال ثلاثة مؤشرات كمية - معدل الخلط وزمن الخلط وكفاءة الخلط - التي تتبع جميعاً قوانين التحجيم
  3. اكتشاف الارتباط القوي بين معامل الترتيب والخلط: إثبات أن معامل ترتيب التدفق الثانوي يرتبط ارتباطاً قوياً بمعدل الخلط، ويمكن أن يكون مؤشراً جيداً لدرجة الخلط في السائل
  4. التحقق من آلية خلط Batchelor: ملاحظة الاضمحلال الأسي في حالة الاضطراب المرن، المقابل لمنطقة خلط Batchelor، وهو متسق مع النتائج التجريبية

شرح الطريقة

الإعداد الهندسي والمعادلات الحاكمة

يستخدم البحث هندسة Taylor-Couette ثنائية الأبعاد تحتوي على أسطوانتين متحدتي المركز:

  • نصف قطر الأسطوانة الخارجية: ro=10μmr_o = 10\mu m
  • نصف قطر الأسطوانة الداخلية: ri=2.5μmr_i = 2.5\mu m
  • تدور الأسطوانة الخارجية بسرعة زاوية Ω=2πs−1\Omega = 2\pi s^{-1} عكس اتجاه عقارب الساعة

يتم وصف حركة السائل بمعادلة Navier-Stokes المعممة: ρ(∂u∂t+u⋅∇u)=−∇p+ηs∇2u+∇⋅τ\rho\left(\frac{\partial \mathbf{u}}{\partial t} + \mathbf{u} \cdot \nabla \mathbf{u}\right) = -\nabla p + \eta_s \nabla^2 \mathbf{u} + \nabla \cdot \boldsymbol{\tau}

نموذج Oldroyd-B التأسيسي

يتم اعتماد نموذج Oldroyd-B لموتر إجهاد البوليمر τ\boldsymbol{\tau}: τ+λτ∇=ηp[∇⊗u+(∇⊗u)T]\boldsymbol{\tau} + \lambda \stackrel{\nabla}{\boldsymbol{\tau}} = \eta_p[\nabla \otimes \mathbf{u} + (\nabla \otimes \mathbf{u})^T]

حيث τ∇\stackrel{\nabla}{\boldsymbol{\tau}} هي المشتقة العلوية المتحركة.

نمذجة الخلط

يتم التحكم في تطور مجال تركيز المذاب c(r,t)c(\mathbf{r}, t) بواسطة معادلة الانتقال الحراري والانتشار: ∂c∂t+(u⋅∇)c=D∇2c\frac{\partial c}{\partial t} + (\mathbf{u} \cdot \nabla)c = D\nabla^2 c

في هذا البحث، تم تعيين D=0D = 0، مع إهمال الانتشار الجزيئي.

المعاملات اللابعدية الرئيسية

  • رقم Reynolds: Re=ρΩro2/ηs=0.628×10−3Re = \rho\Omega r_o^2/\eta_s = 0.628 \times 10^{-3}
  • رقم Weissenberg: Wi=λΩWi = \lambda\Omega
  • نسبة اللزوجة: β=ηp/ηs=1.5\beta = \eta_p/\eta_s = 1.5

الإعداد التجريبي

الطريقة العددية

  • استخدام حل OpenFOAM® مفتوح المصدر بالحجم المحدود
  • استخدام حل RheoTool المتخصص للتعامل مع التدفقات اللزجة المرنة
  • الشبكة: شبكة ذات 100 خلية شعاعية و120 خلية زاوية
  • خطوة زمنية: δt=10−5s\delta t = 10^{-5}s

تعريف معامل الترتيب

يتم تعريف شدة التدفق الثانوي كمعامل ترتيب: Φ=σ=⟨[u(r,t)−u0(r)]2⟩/u0max⁡\Phi = \sigma = \sqrt{\langle[\mathbf{u}(\mathbf{r},t) - \mathbf{u}_0(\mathbf{r})]^2\rangle}/u_0^{\max}

مؤشرات تحديد الخلط

الانحراف المعياري المعياري: Δc=⟨[c(r,t)−c∞]2⟩⟨[c(r,0)−c∞]2⟩\Delta c = \sqrt{\frac{\langle[c(\mathbf{r},t) - c_\infty]^2\rangle}{\langle[c(\mathbf{r},0) - c_\infty]^2\rangle}}

حيث c∞c_\infty هو التركيز في حالة الخلط الكامل.

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

انتقالات حالة التدفق

  1. الانتقال الحرج الفرعي: اكتشاف أن الانتقال من التدفق الطبقي إلى الاضطراب المرن يظهر خصائص حرجة فرعية، مع وجود تأثير التخلف الزمني
  2. رقم Weissenberg الحرج:
    • الفرع الصاعد: Wic=9.95Wi_c = 9.95
    • الفرع الهابط: WicI=9.26Wi_c^I = 9.26
  3. حالة فوضوية ضعيفة: ملاحظة حالة تدفق فوضوية ضعيفة أسفل WicIWi_c^I

تحليل خصائص الخلط

حالة التدفق الفوضوي الضعيف (5.03<Wi≤8.175.03 < Wi \leq 8.17)

  • اضمحلال الانحراف المعياري وفقاً لقانون قوة بأس قريب من −1/6-1/6
  • يظهر الطيف الترددي قانون قوة بأس −5-5 أو أكبر

حالة الاضطراب المرن (Wi≥9.42Wi \geq 9.42)

  • اضمحلال قانون القوة في المرحلة الأولية (أس حوالي 1/51/5)، يتبعه اضمحلال أسي
  • شكل الاضمحلال الأسي: Δc∼exp⁡(−αmixt)\Delta c \sim \exp(-\alpha_{mix}t)
  • يظهر الطيف الترددي خصائص طيف Batchelor (تحجيم k−1k^{-1})

علاقات التحجيم

  1. زمن الخلط: tmix∼(Wi−WicI)−0.52t_{mix} \sim (Wi - Wi_c^I)^{-0.52}
  2. تركيز المنطقة II: cII/c∞∼(Wi−WicI)0.28c_{II}/c_\infty \sim (Wi - Wi_c^I)^{0.28}
  3. علاقة التحجيم: cII∼tmix−1/2c_{II} \sim t_{mix}^{-1/2}

الارتباط بين معامل الترتيب ومعدل الخلط

تظهر النتائج التجريبية ارتباطاً قوياً بين معامل الترتيب Φd\Phi_d ومعدل الخلط αmix\alpha_{mix}، مما يشير إلى أن تعزيز تقلبات السرعة يحسن بشكل مباشر كفاءة الخلط في السائل.

الأعمال ذات الصلة

أبحاث الاضطراب المرن

  • اكتشاف Groisman & Steinberg الأول لظاهرة الاضطراب المرن في محاليل البوليمر
  • تحليل الاستقرار الخطي واللاخطي في هندسة Taylor-Couette
  • ملاحظة أنماط التدفق الحلزونية والشريطية في المحاكاة العددية ثلاثية الأبعاد

أبحاث تعزيز الخلط

  • خلط التدفق الممتد في آلة الدرفلة الرباعية
  • عدم الاستقرار المرن في هندسة الفتحة المتقاطعة
  • تعزيز نقل الحرارة بالحمل في القنوات الملتوية

الخلاصة والمناقشة

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

  1. طبيعة الانتقال: تأكيد الطبيعة الحرجة الفرعية لانتقال الاضطراب المرن في تدفق Taylor-Couette
  2. آلية الخلط: إثبات أن الاضطراب المرن يمكن أن يحقق خلطاً فعالاً عند أرقام Reynolds المنخفضة جداً
  3. منطقة Batchelor: ملاحظة خصائص خلط Batchelor في حالة الاضطراب المرن
  4. مؤشر التنبؤ: يمكن استخدام معامل الترتيب كمؤشر موثوق للتنبؤ بكفاءة الخلط

القيود

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

الاتجاهات المستقبلية

  1. تأثير الموجات المرنة: دراسة تأثير الموجات المرنة على كفاءة الخلط ومعدل الخلط وزمن الخلط
  2. التوسع ثلاثي الأبعاد: التحقق من قوانين التحجيم المكتشفة في الهندسة ثلاثية الأبعاد
  3. التطبيقات العملية: استكشاف إمكانات التطبيق العملي في أجهزة الميكروفلويديكس

التقييم المتعمق

المميزات

  1. البحث المنهجي: تحليل شامل لخصائص الخلط خلال عملية الانتقال من التدفق الطبقي إلى الاضطراب المرن
  2. التحديد الكمي: إدخال مؤشرات كمية متعددة لتحديد درجة الخلط والكفاءة
  3. التحقق النظري: تتوافق النتائج العددية مع نظرية خلط Batchelor والملاحظات التجريبية
  4. القيمة العملية: توفير إرشادات نظرية مهمة لتطبيقات الميكروفلويديكس

أوجه القصور

  1. مشاكل الدقة العددية: يشير الاعتماد على الشبكة لتأثير التخلف الزمني إلى قيود محتملة في الطريقة العددية
  2. فهم الآلية: شرح غير كافٍ للآلية الفيزيائية للاضمحلال وفقاً لقانون القوة في حالة التدفق الفوضوي الضعيف
  3. نطاق المعاملات: نطاق المعاملات المدروسة محدود نسبياً، مع نقص في استكشاف فضاء المعاملات الأوسع

التأثير

  1. المساهمة النظرية: توفير رؤى مهمة لفهم آليات الخلط في السوائل اللزجة المرنة
  2. آفاق التطبيق: ذات أهمية إرشادية كبيرة لتصميم أجهزة الميكروفلويديكس
  3. القيمة المنهجية: يمكن تعميم طرق التحديد المعتمدة على هندسات تدفق أخرى

السيناريوهات المعمول بها

  • تحسين الخلط في أجهزة الميكروفلويديكس
  • التحكم في الخلط في عمليات معالجة البوليمر
  • الخلط على المستوى الدقيق في التطبيقات الطبية الحيوية

المراجع

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