لا يوجد رقم واحد لدقة CFD. تعرّف على مصادر الخطأ وما تذكره دراسة تحقق منشورة، وكيف تحدد الخطأ الذي يحتمله قرارك.

رسم توضيحي من DEEB-FLOW لمنحنى عدد نسلت يقارن التجربة بالمحاكاة
رسم تصوّري من مكتبة حالات DEEB-FLOW. وهو ليس نتيجة من هذا المقال.

الجواب المختصر

انتهى العرض. بدت الحقول مقنعة والاستنتاجات واضحة، ثم يأتي من آخر القاعة السؤال الذي ينتظره كل مهندس محاكاة: ما مدى دقة هذا، زائد أو ناقص كم؟

هناك جوابان سيئان. الأول «نحو 5%» يُقال بثقة ومن دون أساس. والثاني محاضرة من عشر دقائق عن إغلاق معادلات الاضطراب تنتهي من دون رقم واحد.

الجواب الصادق أن CFD ليست له دقة واحدة. المحاكاة دقيقة لكمية معينة، في جريان معين، وبدرجة يمكن تقديرها. وكفاية هذه الدرجة تتوقف على القرار الذي ستبنيه عليها.

لماذا لا يوجد رقم واحد؟

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

لذلك حين يذكر أحدهم دقة CFD بصورة عامة، اسأل عن أي كمية وأي جريان يتحدث. إن لم يستطع الإجابة فالرقم لا يعني شيئاً.

من أين يأتي الخطأ؟

من ثلاثة مصادر. والتمييز بينها هو معظم المهارة.

  • الخطأ العددي. تُحل المعادلات على شبكة محدودة وتتوقف بعد عدد محدود من التكرارات. يمكن تقدير هذا الجزء من دون أي تجربة، بتنعيم الشبكة والتأكد من أن الحل توقف عن التغير.
  • خطأ المدخلات. معدل التدفق عند المدخل قيس بمقياس له تفاوته. الخواص أُخذت من جدول. الهندسة هي نموذج CAD، أما القطعة الحقيقية ففيها لحامات وترسبات وحشية تبرز مليمترين في مجرى الجريان.
  • خطأ النموذج. الاضطراب والاحتراق والجريان متعدد الأطوار تُنمذج بعلاقات إغلاق عويرت على جريانات بسيطة. هذا هو الخطأ الذي يقصده الناس عادة، وهو الوحيد الذي يحتاج تقديره إلى قياسات.

هناك رقم رابع يُنسى غالباً: عدم اليقين في القياس الذي تقارن به. التجربة تقدير آخر له حدود خطئه.

كيف تبدو مقارنة منشورة؟

القواذف فوق الصوتية اختبار صعب، لأن الجهاز كله يعتمد على الخلط المضطرب. قارنت دراسة منشورة عام 2021 في مجلة Applied Thermal Engineering عدة نماذج اضطراب بالقياسات. كان نموذج k-omega SST الأفضل توافقاً، ومع ذلك اختلفت نسبة السحب التي تنبأ بها عن التجربة بمقدار 14% في المتوسط وحتى 20% كحد أقصى.

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

مقارنة بالأرقام

هكذا كنت سأجيب عن سؤال آخر القاعة. الأرقام مثال مفترض لتوضيح الحساب.

قستَ هبوط ضغط مقداره 12.4 kPa. الحساس وطريقة تركيبه يعطيان عدم يقين مقداره 0.5 kPa. أعطت المحاكاة 13.1 kPa، فخطأ المقارنة 0.7 kPa، أي نحو 5.6% من القيمة المقيسة.

E = 13.1 − 12.4 = 0.7 kPa

هذا لا يعني بعدُ أن النموذج مخطئ بنسبة 5.6%. دراسة تنعيم الشبكة تضع عدم اليقين العددي عند 2% من القيمة المحسوبة، أي 0.26 kPa. ومعدل التدفق عند المدخل معروف بدقة 3%، وهبوط الضغط في هذا النظام يتناسب مع معدل التدفق مرفوعاً إلى أس قريب من 1.8، فيصبح ذلك نحو 5.4% في هبوط الضغط، أي 0.71 kPa.

u_val = √(0.5² + 0.26² + 0.71²) ≈ 0.91 kPa

خطأ المقارنة البالغ 0.7 kPa أصغر من 0.91 kPa. الفرق بين المحاكاة والاختبار يقع ضمن ما يفسره عدم اليقين في القياس والشبكة وتدفق المدخل. خطأ النموذج يقع في مكان ما دون 0.9 kPa تقريباً، ولا يمكن تحديده بدقة أكبر بهذه البيانات. هذا هو منطق معيار ASME V&V 20 مختصراً في سطر واحد.

وهو يبيّن أيضاً الخطوة التالية. أكبر حد هو معدل تدفق المدخل، فقياس أدق للتدفق يحسّن المقارنة أكثر بكثير من شبكة أنعم.

ما الخطأ المقبول إذن؟

هو الخطأ الذي لا يغيّر قرارك.

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

المقارنات هي المجال الذي يثبت فيه CFD قيمته. تصميمان متشابهان يُحسبان باستراتيجية الشبكة نفسها والنماذج نفسها يشتركان في معظم أخطائهما، ولذلك يكون الفرق بينهما أوثق من أي من القيمتين منفردة.

أخطاء أراها مراراً

  • ذكر دقة عامة. «محاكاتنا دقيقة حتى 5%» جملة بلا فاعل.
  • المقارنة ببيانات مجهولة عدم اليقين. إن لم تعرف جودة القياس، فتوافق بنسبة 3% واختلاف بنسبة 10% يخبرانك بالقدر الضئيل نفسه.
  • الضبط حتى تتطابق نقطة واحدة. عدّل ثابت اضطراب حتى تتوافق حرارة المخرج فتكون قد طابقت رقماً واحداً. كل حقل آخر صار أقل موثوقية.
  • إهمال الجزء العددي. من دون تقدير لحساسية الشبكة لا يمكن فصل الخطأ العددي عن خطأ النموذج.

ما الذي كنت سأفعله

قبل الدراسة، اكتب القرار والدقة التي يتطلبها. وأثناءها، قدّر عدم اليقين العددي للكمية المهمة. افعل ذلك على الرقم نفسه، فرسم الحقول لن يظهره.

إذا توفرت بيانات اختبار فقارن مع وضع كل مصادر عدم اليقين على الطاولة. وإن لم تتوفر فاذكر ذلك في التقرير، وبيّن أي النتائج مقارنة وأيها مطلقة، وسمِّ القياس الذي يحسم المسألة.

عندها يصبح لسؤال آخر القاعة جواب من جملتين. في هبوط الضغط، عدم اليقين العددي 2% ونتفق مع اختبارك ضمن تفاوته. أما حرارة الجدران فلا بيانات لدينا عنها، فاعتبرها ترتيباً بين البدائل.

المصادر والمطالعة

خدمة ذات صلة

مراجعة النماذج والدراسات المخصصة

معرفة ما يدعمه النموذج، وأين يبقى عدم اليقين، وما يحتاج إلى تحقق.

استعرض الخدمات
كل المقالات