هل تغيّر التجزئة المتوازية النتيجة؟

شغّلتم المسألة نفسها على 16 نواة ثم على 64، والنتيجتان لا تتطابقان. وهذا شائع، وشائع أن يُشخَّص خطأً.

السلوك المتوقع

تجزئة المجال لا تغيّر المعادلات؛ بل تقرر فقط أي معالج يحسب أي جزء. والحل المتقارب ينبغي أن يكون مستقلاً عن عدد الأجزاء.

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

إذا كان الفرق كبيراً

1. الحل لم يتقارب. وهذا أشيع الأسباب بفارق كبير. ففي حلّ نصف متقارب يكبر الفرق الناتج عن ترتيب الجمع، ويتباعد التشغيلان. انظروا ملاحظة التقارب.

2. الفيزياء عابرة. ففي جريان فوضوي أو متذبذب لا تتطابق الحقول اللحظية أصلاً؛ وما ينبغي مقارنته ليس قيمة لحظية بل متوسطاً زمنياً.

3. حدّ تجزئة يقطع منطقة حرجة. بعض النماذج (دوال الجدار خصوصاً، وبعض إغلاقات الأطوار المتعددة) تحتاج معلومات الخلايا المجاورة. فإن قطع حدّ تجزئة طبقةً حديةً رقيقة أو سطحاً بينياً للطور، نشأ تنافر موضعي.

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

ماذا تفعلون

شدّدوا التقارب أولاً؛ فالفرق يزول عادةً. فإن لم يزل، فثبّتوا عدد الأجزاء واذكروا ذلك في التقرير — ففي دراسة تحسين يكون الاتساق بين الحسابات أهم من الدقة المطلقة. انظروا ملاحظة النماذج البديلة.

المراجع

1. Ferziger, J. H., Perić, M. & Street, R. L. (2020). Computational Methods for Fluid Dynamics, 4th ed. Springer. DOI: 10.1007/978-3-319-99693-6
2. Celik, I. B., Ghia, U., Roache, P. J., Freitas, C. J., Coleman, H. & Raad, P. E. (2008). Procedure for estimation and reporting of uncertainty due to discretization in CFD applications. Journal of Fluids Engineering, 130(7), 078001. DOI: 10.1115/1.2960953