A family of second order, linear, unconditionally stable methods for the Cahn-Hilliard-Navier-Stokes equations

ملخص تقني: عائلة من الطرق من الدرجة الثانية، الخطية، والمستقرة غير مشروطة لمعادلات "كاهن-هيليارد-نافييه-ستوكس"

بيان المشكلة
تمثل المحاكاة العددية الدقيقة لديناميكيات الموائع متعددة الأطوار المعقدة، وتحديداً تلك التي تُنمذج بواسطة معادلات "كاهن-هيليارد-نافييه-ستوكس" (CHNS)، تحديات كبيرة بسبب الصلابة الرياضية لمؤثر "كاهن-هيليارد" من الدرجة الرابعة واقترانه غير الخطي بمعادلات "نافييه-ستوكس" غير القابلة للانضغاط. تواجه مخططات الخطوات الزمنية الصريحة القياسية قيوداً شديدة على الخطوة الزمنية (على سبيل المثال، $\Delta t \sim O(\Delta x^4)$) للحفاظ على الاستقرار. وبينما توفر طرق التقسيم المحدب أو الطرق الضمنية الكاملة استقراراً طاقياً غير مشروط، إلا أنها غالباً ما تتطلب حل أنظمة غير خطية ضخمة ومقترنة عند كل خطوة زمنية، مما يزيد من التكلفة الحسابية. وهناك حاجة لخوارزميات تجمع بين دقة زمنية من الدرجة الثانية، واستقرار غير مشروط، وقابلية للحل الخطي دون التضحية بخصائص الحفظ الفيزيائية للنظام.

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

1. استقراء الحدود غير الخطية: يتم التعامل مع حدود الحمل والاقتران غير الخطية صراحةً باستخدام كميات مستقرأة (يُرمز لها بـ $H_{n+\theta}$)، مما يتجنب الحاجة إلى التكرارات غير الخطية.
2. صياغة المتغير المساعد: يتم التعامل مع حد الطاقة الحرة غير الخطي عبر متغير مساعد $q$ (مشتق من جهد البئر المزدوج $F(\phi)$)، مما يعيد صياغة المسألة للحفاظ على خطيتها.
3. تنظيم الانحناء الزمني: يتم إدخال آلية تثبيت محكومة بمعلمة $\epsilon \geq 0$. يتضمن ذلك مؤثر استكمال محدد $J^\epsilon_{n+\theta}$ يدمج الانحناء الزمني المنفصل للحل، وهو مفهوم مقتبس من العمل السابق على تنظيم "نافييه-ستوكس".

تتطلب الخوارزمية الناتجة (الخوارزمية 2.1) حل أنظمة خطية فقط ذات معاملات ثابتة عند كل خطوة زمنية. وتكون هذه الطريقة معلمة بـ $\theta \in (1/2, 1]$ و $\epsilon \geq 0$.

المساهمات الرئيسية والنتائج النظرية

• الاستقرار غير المشروط: يثبت البحث تقدير طاقة منفصل صارم يوضح أن الطريقة المقترحة مستقرة زمنياً على المدى الطويل بشكل غير مشروط لـ $\theta \in (1/2, 1]$ و $\epsilon \geq 0$. يستخدم الإثبات مصفوفات موجبة محددة لتعريف المعايير المنفصلة، ويوضح أن الطريقة تحقق قانون تبدد الطاقة المنفصل دون قيود على الخطوة الزمنية.
• القابلية للحل الخطي: على عكس طرق التقسيم المحدب أو الطرق الضمنية الكاملة التي قد تتطلب حلولاً غير خطية، توفر هذه الطريقة أنظمة خطية عند كل خطوة زمنية، مما يقلل التعقيد الحسابي بشكل كبير.
• الدقة من الدرجة الثانية: تؤكد التحليلات النظرية والتجارب العددية أن الطريقة تحافظ على دقة زمنية من الدرجة الثانية.

النتائج العددية
يتحقق المؤلفون من صحة الطريقة من خلال سلسلة من الحسابات المرجعية:

• تحليل التقارب: باستخدام طريقة الحلول المصنعة (MMS)، تظهر الطريقة تقارباً زمنياً من الدرجة الثانية تقريباً لمتغيرات السرعة، والضغط، وحقول الطور عبر تكوينات مختلفة من العناصر المنتهية ($P2-P1-P2-P2$ و $P3-P2-P3-P3$).
• المتانة والتنظيم: في مسألة تجويف "ليد-دريفن" (Lid-driven cavity) ثنائية الطور، تدرس الدراسة التفاعل بين $\theta$ و $\epsilon$. لوحظ أنه بالنسبة لقيم $\theta$ القريبة من حد "كرانك-نيكولسون" (على سبيل المثال، $\theta=0.51$) مع $\epsilon=0$، يمكن أن تصبح الطريقة غير مستقرة بسبب التذبذبات الزمنية الضعيفة التخميد. إن إدخال تنظيم الانحناء الموجب ($\epsilon > 0$) يعمل بفعالية على تخميد هذه التذبذبات، مما يستعيد المحدودية والمتانة حتى مع خطوات زمنية أكبر.
• المعايير الفيزيائية:
  • التفكك الإسفنجي (Spinodal Decomposition): تلتقط الطريقة بدقة تبدد الطاقة وحفظ الكتلة عبر محاكاة طويلة الأمد ($T=100$)، مما يظهر انفصالاً طورياً وتكتلاً فعالاً.
  • استرخاء القطرة: تنجح الخوارزمية في حل حركة الواجهة المدفوعة بالتوتر السطحي، حيث تنتقل من قطرة مربعة إلى حالة توازن دائرية مع الحفاظ على الكتلة.
  • تجويف "ليد-دريفن" وعدم استقرار "رايلي-تيلور": تتعامل الطريقة مع القص القوي، وتشوه الواجهة، وعدم الاستقرار المدفوع بالكثافة (بما في ذلك لفات "كيلفن-هلمهولتز" في الأنظمة منخفضة اللزوجة) دون ظهور تذبذبات زائفة مرئية.

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