استقلال الشبكة: الدراسة التي تمنح حساب الموائع مصداقيته

المشكلة

تقسيم نطاق متصل إلى خلايا يُدخل خطأً يعتمد على حجم الخلية. فإذا نعّمت الشبكة تغيّرت الإجابة. وإن ظلّت تتغيّر فأنت تُبلّغ عن خاصية لشبكتك لا للجريان.

ودراسة التقارب الشبكي تبرهن أن مزيداً من التنعيم لن يغيّر النتيجة تغييراً ذا معنى. وبدونها تكون نتيجة حساب الموائع غير قابلة للتفنيد.

الإجراء

1. ابنِ ثلاث شبكات على الأقل منعّمة منهجياً. ينبغي أن يكون التنعيم منتظماً عبر النطاق، بنسبة تنعيم r = h_خشنة / h_ناعمة لا تقلّ عن 1.3. والتنعيم فقط حيث تبدو النتائج مثيرة يُبطل التحليل، لأنه يغيّر توزيع الخطأ بدل تقليله انتظاماً.

2. تابع كميّة تهمّك فعلاً. معامل الرفع، أو معامل السحب، أو هبوط الضغط، أو أقصى إجهاد قصّ عند الجدار، أو التدفق الكتلي — قيمة قياسية مرتبطة مباشرة بالسؤال الهندسي. ورسوم الحقول ليست مقياس تقارب؛ فقد يحمل رسمان كنتوريان متطابقان بصرياً قوى متكاملة مختلفة اختلافاً ذا معنى.

3. احسب رتبة الدقة المرصودة:

p = ln((f₃ − f₂)/(f₂ − f₁)) / ln(r)

حيث f₁ هو الحل على أنعم شبكة. وقارن p برتبة مخطّط التقسيم الرسمية. فمخطّط من الرتبة الثانية يُظهر p ≈ 2 يسلك كما هو متوقع. أما p البعيد عن الرتبة الرسمية فيشير إلى أن الحل لم يبلغ بعدُ النطاق المقارب، أو أن مصدر خطأ آخر يهيمن.

4. قدّر خطأ التقسيم باستقراء ريتشاردسون، وأبلغ عن مؤشر التقارب الشبكي (GCI) — نطاق عدم يقين متحفّظ حول نتيجة أنعم شبكة. والمؤشر هو ما يحوّل «بدت الشبكة جيدة» إلى رقم يستطيع المراجع وزنه.

استبانة الطبقة الحدّية

في الجريانات المحدودة بجدار، تكون الاستبانة قرب الجدار عادةً القيد الحاكم. والبارامتر الحاكم هو y⁺، المسافة اللابعدية من الجدار إلى مركز أول خلية.

• محلولة عند الجدار: y⁺ ≈ 1، مع 10–15 خلية عبر الطبقة الحدّية. لازمة للتنبؤ الدقيق بالانفصال والانتقال الحراري عند الجدار واحتكاك السطح.
• دوال الجدار: y⁺ بين 30 و300 تقريباً، بنمذجة الطبقة القريبة من الجدار بدل حلّها. أرخص كثيراً، ومقبولة حيث لا يكون الانفصال هو السؤال المطروح.

والمصيدة هي المنطقة العازلة، نحو 5 < y⁺ < 30، حيث لا تصحّ أي من المعالجتين. والوقوع فيها يُنتج إجابات واثقة وخاطئة بلا عَرَض ظاهر.

وافحص y⁺ بعد الحل لا قبله. فهو يعتمد على قصّ الجدار المحسوب، وأي تقدير سابق للحل تخمينٌ ابتدائي يجب التحقق منه.

أخطاء لن تصلحها شبكة أنعم أبداً

استقلال الشبكة يحدّ خطأ التقسيم وحده. وهو صامت عمّا سواه:

• خطأ النموذج — نموذج اضطراب غير ملائم، أو تطبيق RANS المستقر على جريان غير مستقر فعلاً
• خطأ الشروط الحدّية — نطاق أصغر مما ينبغي، أو ملف دخول لا يطابق الواقع، أو مستوى تماثل مفروض على جريان غير متماثل
• خطأ الهندسة — نموذج CAD مُبسَّط أزال عنصراً ذا دلالة أيروديناميكية

فالحل المتقارب شبكياً تماماً للمسألة الخاطئة يبقى خاطئاً. التقارب الشبكي ضروري لا كافٍ.

ما الذي يُبلَّغ عنه

اذكر أعداد الخلايا ونسبة التنعيم، والكميّة المتابَعة عند كل مستوى، والرتبة المرصودة، ومؤشر GCI على أنعم شبكة، ومدى y⁺ المتحقق، ونموذج الاضطراب، وأبعاد النطاق.

هذه الفقرة وحدها هي الفارق بين نتيجة يستطيع المراجع الاعتماد عليها وأخرى عليه أن يأخذها على الثقة.