Co-moving volumes and the Reynolds transport theorem for two-phase flows

إليك ملخص تقني مفصل للورقة البحثية بعنوان "الأحجام المتحركة معاً ونظرية رينولدز للنقل لتدفقات ثنائية الطور" (Co-moving Volumes and the Reynolds Transport Theorem for Two-Phase Flows) من تأليف ديتر بوثي وماثياس كونه.

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

تتناول الورقة البحثية فجوة جوهرية في النمذجة الرياضية لـ التدفقات ثنائية الطور مع تغير الطور والانزلاق البيني ضمن إطار الواجهة الحادة (sharp-interface framework).

• المشكلة الجوهرية: تعتمد نظرية رينولدز للنقل (RTT) الكلاسيكية على افتراض أن حقل السرعة $v$ منتظم بما يكفي (عادةً $C^1$ أو على الأقل مستمر ليبشيتز) لضمان القابلية للحل الوحيد للمعادلة التفاضلية العادية (ODE) الكينماتيكية $\dot{x} = v(t, x)$. هذا الحل الوحيد هو ما يحدد "الأحجام المتحركة معاً" (co-moving volumes) المكونة من نفس عناصر المائع.
• الانقطاع: في التدفقات ثنائية الطور مع تغير الطاق (انتقال الكتلة عبر الواجهة) والانزلاق (قفزات السرعة المماسية)، يكون حقل السرعة منقطعاً عند الواجهة $\Sigma(t)$.
  • المكون العمودي: تحدث القفزات بسبب انتقال الكتلة ($\dot{m} \neq 0$).
  • المكون المماسي: تحدث القفزات بسبب الانزلاق بين الأطوار.
• النتيجة: هذا الانقطاع يجعل المعادلة التفاضلية العادية الكينماتيكية غير جيدة التحديد (ill-posed) بالمعنى الكلاسيكي. وتحديداً، قد لا تكون الحلول فريدة (عدم التفرد) أو قد لا توجد بالمعنى القياسي عندما تصل جسيمة مائع إلى الواجهة. وبناءً على ذلك، فإن التعريف القياسي للأحجام المتحركة معاً واشتقاق نظرية رينولدز للنقل يفشلان.

2. المنهجية

يستخدم المؤلفون إطاراً رياضياً صارماً يجمع بين الاحتواءات التفاضلية (differential inclusions) والنظرية القياسية الهندسية (geometric measure theory) لحل مشكلة عدم التحديد في المعادلات الكينماتيكية.

• الاحتواءات التفاضلية: بدلاً من معاملة حقل السرعة كدالة أحادية القيمة $v(t,x)$، يستبدله المؤلفون بـ تنظيم متعدد القيم (convexification) عند الواجهة.
  • بالنسبة لـ $x \in \Sigma(t)$، تُعرف السرعة بأنها الغلاف المحدب للحدود من جانب واحد: $\hat{v}(t, x) = \text{conv}\{v^+(t, x), v^-(t, x)\}$.
  • يحول هذا المعادلة التفاضلية العادية $\dot{x} = v$ إلى احتواء تفاضلي $\dot{x} \in \hat{v}(t, x)$.
• خرائط التدفق للاحتواءات: باستخدام نظرية الاحتواءات التفاضلية (بالإشارة إلى فيليبوف وأوبين/تشيلينا)، يعرّفون خريطة تدفق ثنائية الطور $\hat{\Phi}^t_{t_0}$. وخلافاً للخريطة الكلاسيكية، فإن هذه الخريطة هي خريطة متعددة القيم حيث تمثل $\hat{\Phi}^t_{t_0}(x_0)$ مجموعة كل المسارات الممكنة (المجموعات المتاحة) التي تبدأ من $x_0$.
• مبرهنة التباعد في الزمان والمكان: لإثبات نظرية النقل، يتجنب المؤلفون طريقة "السحب العكسي" (pull-back) الكلاسيكية (التي تتطلب خريطة تدفق سلسة). بدلاً من ذلك، يستخدمون مبرهنة التباعد في الزمان والمكان. يقومون ببناء مجال زمان-مكان $\Omega_0$ (رسم بياني للحجم المتحرك معاً) ويطبقون مبرهنة التباعد في $\mathbb{R}^4$ (الزمن + الفضاء ثلاثي الأبعاد). هذا النهج قوي بما يكفي للتعامل مع انخفاض انتظام الحدود والطبيعة متعددة القيم للتدفق.
• التحليل الهندسي للواجهات: تُعرّف الورقة الواجهة $\Sigma(t)$ بصرامة كعائلة من الأسطح الفائقة المتحركة من نوع $C^{1,2}$. كما تحلل شرط "التقاطع العرضي" (transversality) حيث يتقاطع حد حجم التحكم الأولي مع الواجهة، مما يضمن أن تقارب التكاملات السطحية يستمر حتى مع تطور الحجم.

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

أ. التعريف الصارم للمجموعات المتحركة معاً

تقدم الورقة أول تعريف صارم لـ "الأحجام المتحركة معاً" للتدفقات ثنائية الطور ذات السرعات المنقطعة. من خلال استخدام الاحتواءات التفاضلية، يعرّف المؤلفون المجموعة المتحركة معاً $G(t)$ ليس كدومين متطور منفرد، بل كمجموع لجميع المسارات الممكنة التي تبدأ من مجموعة أولية $G_0$. هذا يعالج عدم تفرد مسارات الجسيمات عند الواجهة رياضياً.

ب. نظرية رينولدز للنقل المعممة (RTT)

النتيجة المركزية هي النظرية 5.2، وهي نسخة جديدة من RTT للتدفقات ثنائية الطور. بالنسبة لحجم متحرك معاً $G(t)$ وحقل قياسي $\phi$، فإن معدل تغير التكامل يُعطى بـ:
$$\frac{d}{dt} \int_{G(t)} \phi \, dx \bigg|_{t=t_0} = \int_{G_0 \setminus \Sigma_0} (\partial_t \phi + \text{div}(\phi v)) \, dx + \int_{\Sigma_0 \cap G_0} [[\phi (v - v_\Sigma) \cdot n_\Sigma]] \, do$$
حيث أن:

• $v_\Sigma$ هي السرعة الجوهرية للواجهة.
• $[[\cdot]]$ يشير إلى القفزة عبر الواجهة.
• تأخذ الصيغة في الاعتبار كلاً من النقل الحجمي والتدفق عبر الواجهة المتحركة، حتى عندما لا يتم تحديد سرعة الواجهة بشكل فريد بواسطة الأطوار الحجمية وحدها.

ج. الاتساق الديناميكي الحراري لعدم التفرد

يوضح المؤلفون أن عدم التفرد الرياضي لمسارات الجسيمات عند الواجهة (تحديداً في حالات الانزلاق وتغير الطور) هو أمر متسق فيزيائياً.

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

د. شروط حفظ الحجم

تستنتج الورقة الشروط الدقيقة التي يكون عندها التدفق ثنائي الطور ثابت الحجم (isochoric). ويثبتون أن حفظ الحجم يتطلب:

1. أن يكون حقل السرعة عديم التباعد ($\text{div } v = 0$) في كلا الطورين الحجميين.
2. عدم وجود تغير في الطور ($\dot{m} = 0$) أو استمرارية الكثافة عبر الواجهة ($[[\rho]] = 0$).
   وهذا يسلط الضوء على أنه في أنظمة السائل-الغاز النموذجية (حيث تقفز الكثافة بشكل كبير)، فإن الحجم ليس محفوظاً إذا حدث تغير في الطور، حتى لو كان التدفق غير انضغاطي في كل طور على حدة.

4. النتائج

• وجود الحلول: يقبل الاحتواء التفاضلي الكينماتيكي حلولاً قوية محلية لجميع القيم الأولية، حتى عند الواجهة.
• انتظام خريطة التدفق: ثبت أن خريطة التدفق متعددة القيم مستمرة محلياً بانتظام ليبشيتز بالنسبة لمقياس هاوسدورف، مما يضمن تطور الأحجام المتحركة معاً بطريقة منضبطة رغم عدم تفرد مسارات الجسيمات الفردية.
• التحقق من شروط القفزة: تتضمن نظرية النقل المستمدة طبيعياً شروط القفزة القياسية للكتلة والزخم (مثل $[[\rho(v-v_\Sigma)\cdot n_\Sigma]] = 0$) دون الحاجة إلى افتراضات ارتجالية.
• نتيجة غير متوقعة: في وجود الانزلاق وتغير الطور، يمكن لحدود الحجم المتحرك معاً أن تطور حواف أو "أقواس" جديدة تنبع من نقطة واحدة (كما هو موضح في المثال 4.4 والشكل 3)، وهي ظاهرة مستح المستحيل حدوثها في التدفقات أحادية الطور الكلاسيكية.

5. الأهمية

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

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