كيف يبسط تقريب بوسينيسك (Boussinesq Approximation) نمذجة ديناميكا الموائع وانتقال الحرارة؟ دليل للمحطات التجريبية (Pilot-Plant)

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

يزيل هذا التقريب معظم تعقيدات الكثافة المتغيرة مع الحفاظ على فروق الكثافة التي تحرك الطفو. وهو يسمح للمهندسين والطلاب بالتركيز على أنماط التدفق، والطبقات الحدودية، وانتقال الحرارة، والمجموعات اللابعدية مثل أرقام غراشوف (Grashof) وأرخميدس (Archimedes).

ما الذي يغيره هذا التقريب

حصر تغير الكثافة في الطفو

تُكتب علاقة الكثافة عادةً كالتالي:

[ \rho = \rho_0[1-\beta(T-T_0)] ]

حيث (\rho_0) هي الكثافة المرجعية، و(\beta) هو معامل التمدد الحراري، و(T_0) هي درجة الحرارة المرجعية.

في المعادلات الحاكمة، يتم استبدال (\rho) بـ (\rho_0) في حدود الاستمرارية، واللزوجة، والطاقة. ويتم الاحتفاظ بفرق الكثافة في حد الجاذبية، تقريباً كالتالي:

[ (\rho-\rho_0)\mathbf{g} ]

وهذا يحافظ على الآلية الفيزيائية المسؤولة عن الحمل الحراري الطبيعي.

تصبح معادلة الاستمرارية غير قابلة للانضغاط

يختزل التقريب عموماً معادلة الاستمرارية إلى:

[ \nabla \cdot \mathbf{v}=0 ]

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

معاملة خصائص المائع على أنها ثابتة

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

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

لماذا يساعد هذا في عمليات الوحدات

تصبح معدات الحمل الحراري الطبيعي قابلة للمعالجة

تحتوي المحطات التجريبية غالباً على جدران مسخنة، وأسطح مبردة، وخزانات، وأوعية، وقنوات رأسية حيث تولد فروق درجات الحرارة دورة تدفق. تشمل الأمثلة تجارب التبريد، وأقسام الأنابيب المسخنة، والمبلورات، وخزانات الخلط، والمفاعلات الحيوية، وبعض عمليات التجفيف.

يسمح التقريب بنمذجة هذه التدفقات باستخدام معادلات نافييه-ستوكس غير القابلة للانضغاط ومعادلات الطاقة مع الاحتفاظ بقوة الدفع الحرارية.

يمكن تحليل الحمل الحراري المختلط بكفاءة

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

هذا يجعل من العملي مقارنة مساهمات الحمل القسري والطفو من خلال معاملات لابعدية مثل رقم أرخميدس. وبالمثل، يقيس رقم غراشوف الأهمية النسبية للطفو والمقاومة اللزوجية، بينما يربط رقم براندتل (Prandtl) بين انتشار الزخم وانتشار الحرارة.

تظل الطبقات الحدودية مرئية

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

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

تصبح قياسات المحطات التجريبية أسهل في التفسير

نظراً لأن النموذج يركز على مجموعة صغيرة من المتغيرات الحاكمة، يمكن تنظيم النتائج التجريبية باستخدام علاقات لابعدية. يمكن بعد ذلك ربط قياسات درجة الحرارة، ومعدل التدفق، وتدفق حرارة الجدار، والهندسة بكميات مثل أرقام غراشوف، وبراندتل، ونوسلت (Nusselt)، وستانتون (Stanton)، وأرخميدس.

هذا ذو قيمة في المحطات التجريبية لأن الهدف غالباً هو تحديد السلوك الحاكم وعلاقات القياس (scale-up) بدلاً من إعادة إنتاج كل تأثير جزيئي أو انضغاطي.

توسيع النموذج ليشمل انتقال الحرارة والكتلة

يمكن للتركيز أيضاً أن يحرك الطفو

في عمليات التجفيف، وامتصاص الغاز، والتبلور، والأنظمة التفاعلية أو البيولوجية، يمكن لفروق التركيز أن تغير الكثافة تماماً كما تفعل فروق درجات الحرارة. قد تتضمن العلاقة الخطية كلاً من التمدد الحراري والتركيزي:

[ \rho \approx \rho_0[1-\beta_T(T-T_0)-\beta_C(C-C_0)] ]

وبالتالي يمكن لقوة الطفو الناتجة أن تتضمن مساهمات حرارية وتركيزية.

يسمح هذا التوسيع بنمذجة الحمل الحراري-التركيزي مع الحفاظ على نفس الهيكل الأساسي غير القابل للانضغاط.

يشترك انتقال الحرارة والكتلة في السمات الرياضية

قانون فورييه لتوصيل الحرارة وقانون فيك لانتشار الكتلة لهما أشكال انتشار متشابهة. يُكتب الانتشار الحراري عادةً كالتالي:

[ \alpha=\frac{\lambda}{\rho c_p} ]

بينما يحكم انتشار الكتلة معامل الانتشار الجزيئي (D).

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

المجموعات اللابعدية تربط بين المشكلتين

رقم شميدت (Schmidt)،

[ Sc=\frac{\nu}{D} ]

يقارن بين انتشار الزخم وانتشار الكتلة. بينما يقوم رقم براندتل،

[ Pr=\frac{\nu}{\alpha} ]

بالمقارنة المقابلة لانتقال الحرارة.

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

ما يمكن للمهندسين نمذجته ببساطة أكبر

دورة التدفق ومجالات السرعة

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

توزيع درجة الحرارة وتدفق الحرارة

تظل معادلة الطاقة مقترنة بمجال الزخم من خلال الطفو. ونتيجة لذلك، يمكن للنموذج التنبؤ بكيفية خلق فروق درجات الحرارة للتدفق وكيف يغير هذا التدفق انتقال الحرارة.

بالنسبة للموائع ذات رقم براندتل المرتفع، يمكن للعلاقات القائمة على الجدار أحياناً تقدير تدفق الحرارة بشكل أساسي من قص الجدار وخصائص المائع، مما يقلل الاعتماد على نمذجة الاضطراب (turbulence) التفصيلية في القلب.

انخفاض الضغط وسلوك القص

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

يجب عدم الخلط بين تقريب بوسينيسك وفرضية بوسينيسك للاضطراب، التي تنمذج إجهادات رينولدز المضطربة باستخدام لزوجة مضطربة فعالة. فهما يعالجان تبسيطات مختلفة: أحدهما يتعلق بتغير الكثافة في التدفق الطافي، بينما يتعلق الآخر بإغلاق زخم الاضطراب.

فهم المقايضات

يتطلب تغيرات كسرية صغيرة في الكثافة

شرط الصلاحية المركزي هو:

[ \beta \Delta T \ll 1 ]

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

فروق درجات الحرارة الكبيرة يمكن أن تبطله

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

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

شذوذ الكثافة يتطلب عناية خاصة

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

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

انتقال الكتلة القوي قد يحتاج إلى صياغة أكثر شمولاً

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

يظل الاضطراب تحدياً نمذجياً منفصلاً

لا يحل تقريب بوسينيسك الاضطراب تلقائياً. قد تظل هناك حاجة إلى RANS، أو محاكاة الدوامات الكبيرة (LES)، أو معالجة أخرى للاضطراب لتدفق المحطة التجريبية المضطرب.

وبالمثل، يمكن أن يكون استخدام نموذج لزوجة مضطربة قياسي غير مناسب عندما يكون الاضطراب متبايناً بشدة (anisotropic)، أو يهيمن عليه الطفو، أو يتأثر بهندسة خلط معقدة.

اتخاذ القرار الصحيح لهدفك

يعتمد النموذج الصحيح على نطاق التشغيل، وخصائص المائع، والدقة المطلوبة.

إذا كان تركيزك الأساسي هو التحليل التعليمي: استخدم تقريب بوسينيسك لدراسة الحمل الحراري الطبيعي والمختلط باستخدام معادلات التدفق غير القابل للانضغاط والمجموعات اللابعدية مثل (Gr)، و(Ar)، و(Pr)، و(Nu).
إذا كان تركيزك الأساسي هو التنبؤ بانتقال الحرارة في المحطة التجريبية: تحقق من أن (\beta\Delta T\ll1)، ثم استخدم الخصائص المرجعية مع الاحتفاظ بتغير الكثافة في حد الطفو.
إذا كان تركيزك الأساسي هو التجفيف أو امتصاص الغاز: وسّع علاقة الطفو لتشمل تأثيرات التركيز واستخدم (Sc)، و(Pr)، وعلاقات نوع ستانتون لمقارنة انتقال الكتلة والحرارة.
إذا كان تركيزك الأساسي هو التبلور أو المعالجة الحيوية منخفضة الحرارة: تحقق من سلوك الكثافة غير الخطي والشذوذ المرتبط بالتجمد قبل تطبيق نموذج التمدد الحراري الخطي.
إذا كان تركيزك الأساسي هو المعالجة عالية اللزوجة: ضع في اعتبارك أرقام براندتل الكبيرة وفكر في علاقات انتقال الحرارة القائمة على طبقة الجدار أو قص الجدار حيث يتم استيفاء افتراضاتها.
إذا كان تركيزك الأساسي هو CFD عالي الدقة أو تغيرات كبيرة في الخصائص: استخدم خصائص تعتمد على درجة الحرارة والتركيب، وفكر في صياغة قابلة للانضغاط بالكامل عندما لا تعود تغيرات الكثافة صغيرة.

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

جدول الملخص:

الموضوع	الرؤية الرئيسية
الافتراض الأساسي	يعامل التدفق كغير قابل للانضغاط بخصائص ثابتة تقريباً.
الطفو	يحتفظ بتغير الكثافة فقط في حد قوة الجاذبية.
التطبيقات الرئيسية	الحمل الحراري الطبيعي والمختلط في الخزانات، والأوعية، والقنوات، والمفاعلات، والمفاعلات الحيوية، ومعدات انتقال الحرارة.
المعلمات المفيدة	أرقام غراشوف، وأرخميدس، وبراندتل، ونوسلت، وشميدت، وستانتون.
انتقال الحرارة والكتلة	يمكن أن يشمل الطفو المدفوع بالحرارة والتركيز باستخدام علاقات كثافة خطية.
شرط الصلاحية	يتطلب تغيرات كسرية صغيرة في الكثافة، يتم التحقق منها عادةً بـ (\beta\Delta T \ll 1).
القيود الرئيسية	أقل موثوقية لتغيرات الحرارة أو التركيز الكبيرة، وشذوذ الكثافة، والانضغاطية القوية، أو تغير الخصائص الكبير.
الاضطراب	يتطلب معالجة منفصلة مثل RANS أو LES عندما يكون التدفق مضطرباً.
القيمة الهندسية	يدعم التحليل الفعال للدورة، والطبقات الحدودية، ومجالات الحرارة، وتدفق الحرارة، وانخفاض الضغط، وسلوك القياس (scale-up).

طبق نمذجة تدفق الطفو مع المحطة التجريبية المناسبة

هل تحتاج إلى دراسة الحمل الحراري الطبيعي، أو الحمل الحراري المختلط، أو انتقال الحرارة، أو انتقال الكتلة في ظل ظروف تشغيل واقعية؟ اتصل بـ LABPARK لمناقشة المحطات التجريبية لعمليات الوحدات التعليمية والمهنية للهندسة الكيميائية، والمعالجة الحيوية والتكنولوجيا الحيوية، وتطبيقات المعالجة البيئية ومعالجة المياه.

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