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

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

ما الذي يفترضه تقريب الحالة شبه المستقرة

تتكيف الملفات الحرارية بشكل أسرع من حركة الواجهة

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

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

الحرارة الكامنة تهيمن على توازن الحدود المتحركة

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

لذلك يعتمد توازن الواجهة على

[ \text{الحرارة المزالة}=\text{الحرارة الكامنة المنطلقة من المادة الصلبة الجديدة}. ]

تبقى الواجهة بالقرب من درجة حرارة تغير الطور

عادةً ما يتم تعيين درجة حرارة تغير الطور التوازني (T_m) للواجهة بين الصلب والسائل، أو درجة حرارة تبلور فعالة عندما تكون تأثيرات التركيز أو الضغط مهمة.

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

بناء نموذج المقاومة الحرارية

تحديد هندسة الأنبوب والمادة الصلبة

ليكن

(R_i) نصف القطر الداخلي للأنبوب،
(R_t) نصف القطر الخارجي للأنبوب،
(R_s(t)) نصف قطر الواجهة المتحركة بين الصلب والسائل،
(k_w) الموصلية الحرارية لجدار الأنبوب،
(k_s) الموصلية الحرارية للطبقة الصلبة،
(h_i) معامل انتقال الحرارة لجانب المبرد الداخلي،
(T_c) درجة حرارة المبرد الكلية،
(\rho_s) كثافة المادة الصلبة،
(L) الحرارة الكامنة للتصلب أو التبلور.

بالنسبة لترسب ينمو للخارج من سطح الأنبوب، تشغل المادة الصلبة المنطقة الواقعة بين (R_t) و (R_s)، وسمكها هو

[ \delta(t)=R_s(t)-R_t. ]

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

بالنسبة لأنبوب أسطواني، فإن مقاومة غشاء المبرد، وجدار الأنبوب، والطبقة الصلبة لكل وحدة طول هي

[ \mathcal{R}_i'=\frac{1}{h_i,2\pi R_i}, ]

[ \mathcal{R}_w'=\frac{\ln(R_t/R_i)}{2\pi k_w}, ]

[ \mathcal{R}_s'=\frac{\ln(R_s/R_t)}{2\pi k_s}. ]

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

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

[ \mathcal{R}_{\text{tot}}'

\frac{1}{h_i,2\pi R_i} + \frac{\ln(R_t/R_i)}{2\pi k_w} + \frac{\ln(R_s/R_t)}{2\pi k_s}. ]

حساب معدل إزالة الحرارة اللحظي

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

[ q'(R_s)

\frac{T_m-T_c}{\mathcal{R}_{\text{tot}}'}. ]

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

اشتقاق معدل نمو المادة الصلبة

كتابة توازن الحرارة الكامنة

عندما تتحرك الواجهة للخارج بمقدار (dR_s)، يكون حجم المادة الصلبة المتكون حديثاً لكل وحدة طول من الأنبوب هو

[ dV'=2\pi R_s,dR_s. ]

الحرارة الكامنة المقابلة المنطلقة هي

[ dQ_{\text{latent}}'

\rho_s L,2\pi R_s,dR_s. ]

القسمة على الوقت ومساواة هذا المعدل بالحرارة المزالة يعطي

[ q'(R_s)

2\pi R_s\rho_s L\frac{dR_s}{dt}. ]

التعويض بتعبير المقاومة ينتج معادلة النمو الحاكمة:

[ \boxed{ \frac{dR_s}{dt}

\frac{T_m-T_c} { 2\pi R_s\rho_s L \left[ \frac{1}{h_i,2\pi R_i} + \frac{\ln(R_t/R_i)}{2\pi k_w} + \frac{\ln(R_s/R_t)}{2\pi k_s} \right] } } ]

تعطي هذه المعادلة معدل النمو الشعاعي اللحظي للطبقة الصلبة.

تفسير انخفاض معدل النمو

مع زيادة (R_s)، تزداد أيضاً مقاومة المادة الصلبة

[ \frac{\ln(R_s/R_t)}{2\pi k_s} ]

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

العامل (2\pi R_s) في مصطلح الحرارة الكامنة مهم أيضاً: الزيادة الشعاعية المعينة تخلق حجماً أكبر من المادة الصلبة لكل وحدة طول من الأنبوب مع نمو محيط الواجهة.

تبسيطات مفيدة للتقديرات الهندسية

استخدام تقريب مستوٍ لطبقة رقيقة

إذا كان سمك المادة الصلبة صغيراً مقارنة بنصف قطر الأنبوب، (\delta \ll R_t)، يمكن تقريب المقاومة الأسطوانية بمقاومة مستوية. تصبح معادلة النمو

[ \frac{d\delta}{dt}

\frac{T_m-T_c} {\rho_s L \left[ \frac{1}{h_i} + \frac{t_w}{k_w} + \frac{\delta}{k_s} \right]}, ]

حيث (t_w) هو سمك جدار الأنبوب.

إذا تم تجميع مقاومة غشاء المبرد والجدار في مقاومة ثابتة (R_0)، فإن

[ \frac{d\delta}{dt}

\frac{T_m-T_c} {\rho_s L\left(R_0+\delta/k_s\right)}. ]

التكامل لتقدير السمك مقابل الوقت

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

[ t

\frac{\rho_s L}{T_m-T_c} \left[ R_0\delta+\frac{\delta^2}{2k_s} \right]. ]

تظهر هذه النتيجة نظامين محدودين.

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

[ \delta \propto \sqrt{t}. ]

استخدام الشكل الأسطواني للترسبات الكبيرة

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

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

تطبيق النموذج على عملية وحدة تعتمد على الأنابيب

تقدير معامل انتقال الحرارة الفعال

يحدد نموذج المقاومة أي جزء من النظام يحد من النمو. يتم تمثيل غشاء المبرد الداخلي بـ (1/h_i)، بينما تساهم المادة الصلبة النامية بـ (\delta/k_s) في الحد المستوي.

يمكن تعريف معامل إجمالي فعال من المقاومة الكلية، مع مراعاة منطقة مرجعية متسقة:

[ U=\frac{1}{R_{\text{total}}}. ]

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

تقييم ظروف تشغيل المبرد

تؤثر درجة حرارة المبرد ومعدل التدفق على القوة الدافعة الحرارية والمعامل الداخلي (h_i). يؤدي خفض (T_c) عموماً إلى زيادة القوة الدافعة (T_m-T_c)، بينما يمكن أن تؤدي زيادة سرعة المبرد إلى تقليل مقاومة الغشاء الداخلي.

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

مراعاة تأثيرات درجة الحرارة والتركيب

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

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

ربط النموذج بسلوك الأنبوب المحوري

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

[ \Delta T(z,t)=T_m(z,t)-T_c(z,t), ]

ويجب حل معادلة النمو على طول الأنبوب وكذلك في الوقت.

ينتج عن هذا سمك محلي (R_s(z,t))، والذي يمكنه تحديد المناطق التي يحتمل فيها حدوث ترسب مفرط أو انسداد أو تبلور غير منتظم.

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

التقريب يغفل التسخين المحسوس العابر

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

عندما تكون التأثيرات المحسوسة مهمة، يكون نموذج الحدود المتحركة العابر الكامل أكثر ملاءمة.

الخصائص الثابتة قد تكون غير واقعية

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

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

ظروف السطح يمكن أن تغير سلوك الواجهة

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

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

النمو يمكن أن يقلل من سعة العملية المفيدة

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

لذلك، فإن تعظيم معدل النمو اللحظي ليس دائماً هو نفسه تحسين عملية الوحدة الشاملة.

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

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

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

يحول توازن المقاومة في الحالة شبه المستقرة حدود التصلب المتحركة إلى حساب عملي لانتقال الحرارة مع توضيح سبب تباطؤ النمو مع زيادة سمك الطبقة الصلبة.

جدول الملخص:

العنصر الأساسي	الرؤية الهندسية
الافتراض الجوهري	تتكيف الملفات الحرارية بشكل أسرع من حركة الواجهة بين الصلب والسائل.
القوة الدافعة	درجة حرارة تغير الطور ناقص درجة حرارة المبرد الكلية، ΔT = Tm − Tc.
مسار انتقال الحرارة	غشاء المبرد الداخلي، وجدار الأنبوب، والطبقة الصلبة النامية.
توازن النمو	الحرارة المزالة تساوي الحرارة الكامنة المنطلقة من المادة الصلبة المتكونة حديثاً.
النمو الأسطواني	استخدم نصف قطر الواجهة الحالي والمقاومة الصلبة اللوغاريتمية للترسبات السميكة.
تقدير الطبقة الرقيقة	يكون النمو في البداية خطياً تقريباً، ثم يقترب من δ ∝ √t مع زيادة مقاومة المادة الصلبة.
عوامل التحكم الرئيسية	درجة حرارة المبرد، وسرعة المبرد، والخصائص الحرارية، والحرارة الكامنة، وسمك الترسب.
قيود النموذج	قد يتطلب التسخين المحسوس العابر، والخصائص المتغيرة، وتأثيرات التركيب، والتحجيم، والترسبات المسامية نمذجة ذات دقة أعلى.

تنبأ وحسن عمليتك القائمة على الأنابيب

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

توفر LABPARK وحدات تجريبية للعمليات التعليمية والمهنية للجامعات ومعاهد البحوث والمؤسسات، بما في ذلك حلول للهندسة الكيميائية، والمعالجة الحيوية والتكنولوجيا الحيوية، ومعالجة المياه والبيئة. يمكننا مساعدتك في تقييم ظروف المبرد، ومقاومة انتقال الحرارة، ونمو الواجهة، وتغير درجة الحرارة المحوري، ومتطلبات التوسع حتى تتمكن من:

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

اتصل بـ LABPARK لتحديد مادتك المستهدفة، وهندسة الأنبوب، وظروف المبرد، وهدف العملية، وتطوير نهج للمصنع التجريبي يناسب احتياجاتك.