An Entropy Stable Formulation of Two-equation Turbulence Models with Particular Reference to the k-epsilon Model
تقدم هذه الورقة إطاراً عددياً متيناً لنماذج الاضطراب ثنائية المعادلات، ولا سيما نموذج k-epsilon، من خلال دمج مبادئ إنتاج الإنتروبيا والمتوسط الزماني-المكاني لاشتقاق معادلات متناظرة ومستقرة إنتروبياً تضمن الاتساق والإيجابية دون الاعتماد على اللزوجة الاصطناعية.
المؤلفون الأصليون:Guillermo Hauke, Thomas J. R. Hughes
تخيل أنك تحاول التنبؤ بكيفية تدفق نهر حول صخرة، أو كيف يندفع الهواء فوق سيارة مسرعة. في عالم الفيزياء، تُعرف هذه الدراسة بديناميكا السوائل. لكن السوائل مخادعة؛ فهي لا تتدفق بسلاسة مثل الماء في مجرى هادئ، بل تلتف وتدور وتصطدم بنفسها في أنماط فوضوية تسمى الاضطراب (Turbulence). لفهم هذه الحركات الجامحة، يستخدم العلماء مجموعة من القواعد تسمى معادلات نافييه-ستوكس. فكر في هذه المعادلات كدليل تعليمات نهائي لكيفية حركة السوائل. ومع ذلك، عندما تصبح السوائل مضطربة، يصبح الدليل معقداً للغاية لدرجة أن أسرع الحواسيب الفائقة في العالم لا تستطيع حله بشكل مثالي للتصاميم الواقعية.
وللتغلب على ذلك، يستخدم المهندسون طريقاً مختصراً. فبدلاً من تتبع كل دوامة صغيرة، يأخذون "متوسطاً" للتدفق، مما يؤدي إلى تنعيم الفوضى وتحويلها إلى صورة يمكن التحكم بها. وهذا ما يسمى بـ "متوسط رينولدز" (Reynolds-averaging). ولكن تكمن المشكلة هنا في أنه عندما تقوم بتنعيم الفوضى، فإنك تفقد بعض القواعد الأساسية التي تحافظ على استقرار الكون. إحدى أهم هذه القواعد هي القانون الثاني للديناميكا الحرارية، والذي ينص باختاً على أن الاضطراب (أو "الإنتروبيا") في أي نظام مغلق يميل دائماً إلى الزيادة. الأمر يشبه كيف أن الغرفة الفوضوية لا تنظف نفسها أبداً من تلقه؛ فالطاقة تتبدد دائماً على شكل حرارة. إذا تجاهل محاكاة الحاسوب هذه القاعدة، فقد يصبح غير مستقر، وينتج نتائج جامحة وغير منطقية تؤدي إلى تعطل البرنامج. السؤال الكبير الذي طرحه العلماء هو: هل يمكننا بناء نموذج للاضطراب يحترم "قاعدة الفوضى" هذه لتبقى محاكاتنا الحاسوبية مستقرة وواقعية فيزيائياً؟
هذه الورقة البحثية، التي كتبها غييرمو هوكي وتوماس جيه آر هيوز، تتناول هذه المشكلة تحديداً. يقترح المؤلفان طريقة جديدة لكتابة معادلات نماذج الاضطراب ذات المعادلتين (وتحديداً نموذج k-ϵ الشهير) بحيث تلتزم بطبيعتها بالقانون الثاني للديناميكا الحرارية. لقد فعلوا ذلك من خلال تقديم خدعة رياضية ذكية تسمى "المتوسط الزماني-المكاني" (space-time averaging). تخيل أنك لا تنظر إلى تدفق مضطرب من خلال كاميرا تصوير سريع فحسب، بل أيضاً من خلال عدسة ضبابية قليلاً تعمل على متوسط التفاصيل المكانية الصغيرة. ومن خلال القيام بذلك، يكتشفون صلة خفية بين طاقة الدوامات المضطربة ومعدل تلاشي تلك الطاقة (التبدد).
النتيجة الرئيسية للورقة هي أنه إذا قمت بتغيير المتغيرات في معادلاتك لتمثل هذه الكميات المتوسطة، فإن النظام بأكمله يصبح "متماثلاً". في عالم الرياضيات، التماثل هو قوة خارقة؛ فهو يعني أن المعادلات أكثر استقراراً وأقل عرضة للانهيار أثناء المحاكاة الحاسوبية. يوضح المؤلفان أنه لكي ينجح هذا، يجب أن تتبع الطريقة التي ينتشر بها الاضطراب عبر الحرارة والزخم (معاملات الانتشار) قيوداً محددة. لقد اختبروا هذا النسخة "المستقرة إنتروبياً" من نموذج k-ϵ باستخدام الحاسوب عبر محاكاة تدفق الهواء فوق لوح مسطح. كانت النتائج واعدة: أنتج النموذج نتائج مستقرة وواقعية تطابقت جيداً مع البيانات التجريبية، وفعلت ذلك دون الحاجة إلى "الإصلاحات" الاصطناعية (مثل إضافة لزوجة وهمية) التي غالباً ما تتطلبها النماذج الأقدم. ويقترح المؤلفان أن هذا النهج يوفر طريقة أكثر اتساقاً من الناحية الفيزيائية لتصميم نماذج الاضطراب، مما يضمن أن تحترم المحاكاة الحاسوبية القوانين الأساسية للطبيعة.
ملخص تقني: صياغة مستقرة إنتروبياً لنماذج الاضطراب ثنائية المعادلة مع إشارة خاصة إلى نموذج k-ϵ
بيان المشكلة يعتمد تصميم الخوارزميات العددية القوية لديناميكا السوائل بشكل كبير على الاتساق والاستقرار. وبينما يوفر القانون الثاني للديناميكا الحرارية (إنتاج الإنتروبيا) مبدأ استقرار فيزيائي أساسي لمعادلات نافيير-ستوكس، فإن تطبيقه على معادلات رينولدز-أفيريجد نافيير-ستوكس (RANS) ليس بالأمر الهين. في صياغات RANS القياسية، لا سيما تلك التي تستخدم نماذج اضطراب ثنائية المعادلة مثل k-ϵ، تضيع الصلة المباشرة بين الإنتروبيا المتوسطة والمتغيرات المعتمدة المتوسطة بسبب عدم خطية علاقة غيبس. وبناءً على ذلك، غالباً ما تفتقر الصياغات القياسية إلى دالة إنتروبيا معممة، مما يجعل من الصعب دمج مبادئ الاستقرار الفيزيائي مباشرة في المخطط العددي دون اللجوء إلى اللزوجة الاصطناعية. لقد كانت المحاولات السابقة لاستعادة استقرار الإنتروبيا، مثل تلك التي قام بها Jansen وآخرون، مقتصرة على النماذج أحادية المعادلة أو الأنظمة التي تتضمن ست معادلات لموتر إجهاد رينولدز، مما ترك نموذج k-ϵ واسع الاستخدام دون صياغة مستقرة إنتروبياً بشكل صارم.
المنهجية يقترح المؤلفون إطار عمل لاشتقاق صياغة مستقرة إنتروبياً لنماذج الاضطراب ثنائية المعادلة، وتحديداً نموذج k-ϵ، من خلال إدخال المتوسط الزماني-المكاني.
المتوسط الزماني-المكاني وتحويل المتغيرات: يقدم البحث مؤثر متوسط زماني-مكاني لسد الفجوة بين معدل التبديد (ϵ) والإنتروبيا المعممة. من خلال إجراء متوسط مكاني للطاقة الحركية المضطربة (k) عبر حجم صغير، يوضح المؤلفون أن ϵ يأخذ في الاعتبار التقلبات المكانية لـ k. يؤدي هذا إلى تعريف مقياس سرعة جديد، q1=(νϵ)1/4، بالإضافة إلى مقياس السرعة المضطربة القياسي q0=2k. يتم إعادة تعريف نظام متغيرات الحفظ ليشمل الكثافة، والزخم، والطاقة الكلية، و q0، و q1. كما يتم تعديل الطاقة الكلية لتشمل مساهمات من كل من q0 و q1.
دالة الإنتروبيا المعممة والتماثل: باستخدام هذه المتغيرات الجديدة، يبني المؤلفون دالة إنتروبيا معممة H(U)=−ρs^، حيث s^ هي إنتروبيا معدلة بناءً على المتغيرات المتوسطة زمانياً-مكانياً. ويتم اشتقاق متغيرات الإنتروبيا V المقابلة كتدرج لـ H بالنسبة لمتغيرات الحفظ U. هذا التغيير في المتغيرات يحول الصيغة شبه الخطية لمعادلات k-ϵ إلى نظام زائدي متماثل (تحديداً نظام انتشار-حمل متماثل). تصبح مصفوفات المعاملات للاشتقاق الزمني وحدود الحمل متماثلة، كما ثبت أن مصفوفة الانتشار متماثلة وموجبة شبه محددة تحت قيود معينة.
القيود من أجل إنتاج الإنتروبيا: لضمان استيفاء النظام لمتراجحة كلاوزيوس-دوبليه (إنتاج إنتروبيا غير سالب)، يشتق المؤلفون قيوداً محددة على معاملات نموذج الاضطراب:
يجب أن يتساوى الانتشار المضطرب لمعدل التبديد (ϵ) مع الانتشار المضطرب للطاقة الحركية المضطربة (k). ومن حيث أعداد براندتل، يتطلب هذا أن يكون Prϵ=Prk.
يجب أن تحقق حدود المصدر في معادلة q1 متباينات محددة لضمان بقاء إنتاج الإنتروبيا المعممة غير سالب. يتضمن ذلك تعديل حد التبديد في دالة المصدر ليتغلب على حدود الإنتاج تحت ظروف معينة.
التنفيذ العددي: تم تنفيذ المعادلات المصاغة باستخدام طريقة العناصر المحدودة "غاليركين غير المتصل زمنياً/المربعات الصغرى". يتم حل النظام باستخدام خوارزمية المتنبئ-المصحح المتعدد (predictor-multicorrector) مع محلل خطي من نوع GMRES. تعالج عملية التنفيذ مشكلة الصلابة (stiffness) لمعادلات k-ϵ، لا سيما بالقرب من الجدران، من خلال استخدام شرط حدودي محدد لـ q1 مشتق من توازن التبديد والانتشار المتقاطع.
النتائج الرئيسية تعرض الورقة محاكاة عددية لطبقة حدودية مسطحة قريبة من حالة عدم الانضغاط باستخدام نموذج Lam-Bremhorst k-ϵ المعدل ليتوافق مع هذا الإطار المستقر إنتروبياً.
التحقق من الصحة: تظهر النتائج الخاصة بمعامل احتكاك الجلد وملفات السرعة (في إحداثيات الجدار) توافقاً جيداً مع الارتباطات التجريبية (Karman-Schönherr) والحسابات الراسخة الأخرى (Patel et al., Wilcox).
الاستقرار: يظهر المخطط العددي متانة ملحوظة. ومن الجدير بالذكر أن معادلة q1 (معادلة التبديد المحولة) تظهر استقراراً ممتازاً، على النقيض من الصعوبات العددية المعتادة المرتبطة بمعادلة ϵ القياسية.
الاستقلال عن الشبكة: تشير دراسة التكرار إلى أن الحل غير حساس لحجم عنصر القرب من الجدار (y+)، حتى لقيم تصل إلى y+=5.0. ويعزى ذلك إلى تنفيذ الشرط الحدي المحدد لـ q1، والذي يسمح بشبكات أقرب من الجدار أكثر خشونة مقارنة بنماذج k-ϵ القياسية ذات رقم رينولدز المنخفض التي تتطلب عادةً y+<1.
إنتاج الإنتروبيا: تؤكد عمليات المحاكاة أن إجمالي إنتاج الإنتروبيا (التبديد اللزج بالإضافة إلى مساهمات الاضطراب) يظل غير سالب عبر الطبقة الحدودية، مما يؤكد صحة الاشتقاق النظري.
الأهمية والادعاءات يزعم المؤلفون أن هذا العمل نجح في توسيع مفهوم الاستقرار الإنتروبي ليشمل نموذج الاضطراب k-ϵ ثنائي المعادلة واسع الاستخدام. فمن خلال إدخال المتوسط الزماني-المكاني والمتغيرات المرتبطة به q0 و q1، تمكنوا من استعادة دالة إنتروبيا معممة تسمح بتماثل النظام.
تكمن الأهمية الأساسية في القدرة على تصميم خوارزميات عددية ترث مبدأ الاستقرار الفيزيائي للقانون الثاني للديناميكا الحرارية مباشرة في الصياغة المنفصلة. وهذا يوفر بديلاً صارماً لإضافة اللزوجة الاصطناعية من أجل الاستقرار. ويؤكد المؤلفون أن القيود المشتقة (تحديداً Prϵ=Prk) ضرورية لكل من إنتاج الإنتروبيا وإيجابية الحل، وهو ما يتوافق مع النتائج السابقة لـ Mohammadi و Pironneau. ويخلص المؤلفون إلى أن هذا الإطار يوفر مساراً نحو تصميم نماذج اضطراب أكثر تماسكاً من الناحية الفيزيائية ومخططات عددية أكثر استقراراً، كما هو موضح في التطبيق الناجح لنموذج Lam-Bremhorst.