Lattice Boltzmann modeling for conduction heat transfer based on the Qadyan method

ملخص تقني: نمذجة لابتز بولتزمان لانتقال الحرارة بالتوصيل بناءً على طريقة قديان

بيان المشكلة
تتناول الورقة البحثية محاكاة انتقال الحرارة بالتوصيل باستخدام طريقة لابتز بولتزمان (LBM). وبينما تعد طريقة LBM فعالة في ديناميكا الموائع، فإن تطبيقها على توصيل الحرارة يواجه قيوداً تتعلاه تجانس الشبكة والاستقرار العددي. فالمخططات القياسية لـ LBM، وخاصة تلك القائمة على تقريب "باتناجر-غروس-كروك" (BGK)، غالباً ما تكون مقيدة بشبكات منتظمة ذات خطوات مكانية متساوية. علاوة على ذلك، يمكن لطرق التجزئة المكانية القياسية، مثل مخططات الفروق المحدودة من الدرجة الثانية، أن تسبب تذبذبات عددية في الحل. تسعى الدراسة إلى تقييم استراتيجيات بديلة للتجزئة المكانية لتعزيز الدقة والاستقرار دون إدخال تعقيدات إضافية على إطار عمل LBM القياسي أو اللجوء إلى شبكات غير منتظمة.

المنهجية
تنفذ البحث وتقارن بين ثلاثة نهج متميزة لحل معادلة لابتز بولتزمان (LBE) لانتقال الحرارة بالتوصيل:

1. طريقة LBM القياسية: تستخدم نماذج شبكية تقليدية (D2Q9، وD2Q5، وD2Q8) مع عمليات التصادم والتدفق القياسية. وتشير الدراسة إلى أنه بالنسبة لمشكلات الانتشار البحتة، فإن النماذج ذات مجموعات السرعة الأقل (مثل D2Q8 وD2Q5) تكون كافية، حيث يفتقر انتشار الحرارة إلى التفضيل الاتجاهي.
2. طريقة لابتز بولتزمان بالفروق المحدودة (FDLBM): يقوم هذا النهج بتجزئة الحد الحملي لمعادلة LBE باستخدام مخططات الفروق المحدودة. وتحديداً، يستخدم فرقاً أمامياً من الدرجة الأولى للمشتقات الزمنية ومخطط "أب ويند" (upwind) من الدرجة الأولى للمشتقات المكانية لتجنب التذبذبات، رغم أن هذا يحد من الدقة المكانية إلى الدرجة الأولى بالقرب من الحدود.
3. طريقة قديان: نهج مبتكر اقترحه "غدياني وآخرون" يجمع بين المخططات شبه المنفصلة وطريقة التربيع التفاضلي (DQM). في هذه الطريقة، يتم تجزئة المشتقات المكانية في معادلة LBE باستخدام معاملات وزن التربيع التفاضلي (DQM). ولتكرار تأثير الـ "أب ويند" المتأصل في مشكلات الانتشار-الحمل، تقوم الطريقة بتعديل معاملات الوزن لتعيين وزن أكبر للعقد التي تقع في اتجاه المنبع مقارنة بالعقد التي تقع في اتجاه المصب. يسمح هذا بحل معادلة LBE كنظام من المعادلات التفاضلية العادية دون الارتباط بنوع محدد من نماذج السرعة المنفصلة.

تتحقق الدراسة من صحة هذه الطرق باستخدام انتقال الحرارة العابر أحادي البعد في لوح يتضمن كلاً من شروط "ديريكليه" (درجة حرارة ثابتة) وشروط "نيومان" (معزول/تدفق صفري). وتُستخدم الحلول التحليلية المستمدة من معادلة الحرارة كمعايير مرجعية. كما يتم حساب معايير الخطأ ($L_1$ و$L_2$ و$L_\infty$) لتقييم الدقة، ويتم حل مشكلة انتشار لوح ثنائي الأبعاد لإثبات قابلية تطبيق الطريقة في الحالات متعددة الأبعاد.

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

• الدقة مقابل الاستقرار: تشير النتائج إلى أنه بينما توفر طريقة LBM التقليدية دقة عالية، فإن طريقة قديان المقترحة توفر بديلاً قوياً يحسن الدقة مقارنة بنهج FDLBM، لا سيما في شروط نيومان. بالنسبة لشروط ديريكليه، أظهر نموذج D2Q5 LBM القياسي أدنى معايير خطأ، لكن طريقة قديان تفوقت بشكل كبير على صيغ FDLBM.
• استقلال الشبكة: نجحت طريقة قديان في توسيع نطاق دقة واستقرار الحلول، مما زاد من الدقة دون إدخال تعقيدات إضافية على طريقة لابتز بولتزمان القياسية أو اللجوء إلى شبكات غير منتظمة.
• التكلفة الحسابية: تتمثل النتيجة الحاسمة في المقايضة بين الدقة والكفاءة؛ إذ تنطوي طريقة قديان على تكلفة حسابية أعلى ووقت محاكاة أطول مقارنة بكل من LBM القياسية وFDLBM، وذلك بسبب زيادة عدد العمليات الحسابية المطلوبة لمعاملات وزن DQM.
• حساسية المعلمات: تثبت الدراسة أن القيم المرجعية لسرعة الشبكة في LBM وFDLBM تعتمد على نوع شروط الحدود (درجة الحرارة مقابل التدفق). وفي المقابل، تستخدم طريقة قديان قيمة ثابتة لسرعة الشبكة ($c=0.015$) بغض النظر عن نوع شرط الحدود، مما يسهل عملية اختيار المعلمات.
• التحقق ثنائي الأبعاد: في محاكاة لوح مربع ثنائي الأبعاد، أنتجت طريقة قدي "كنتورات" متساوية الحرارة وتوزيعات درجة حرارة المركز التي أظهرت توافقاً جيداً مع نتائج طريقة الفروق المحدودة، مما يثبت قدرتها على التعامل مع مشكلات انتقال الحرارة متعددة الأبعاد.

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