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

[ \frac{dT}{dt}

\frac{F}{V}(T_i-T) + \frac{Q}{\rho c_p V}, ]

حيث إن (T) هي درجة حرارة الخزان جيدة الخلط، و(T_i) هي درجة حرارة الدخول، و(F) هو معدل التدفق الحجمي، و(V) هو حجم السائل، و(Q) هو معدل انتقال الحرارة الصافي. ويؤدي إجراء الخطية حول حالة مستقرة إلى نموذج دخل وخرج من الرتبة الأولى مناسب لتحليل دوال التحويل وتصميم وحدات التحكم بالتغذية الراجعة.

تتمثل الصياغة الأساسية في موازنة طاقة ديناميكية معبَّر عنها بمتغيرات الانحراف. ومع اعتبار (T_i) و(Q) مدخلين و(T) خرجًا، يتصرف CSTH ذو الحجم الثابت تقريبًا كنظام من الرتبة الأولى، مع مساهمات من درجة حرارة الدخول ومعدل التسخين.

تعريف مسألة التحكم

اختيار الخرج المتحكم فيه

المتغير المتحكم فيه المعتاد هو درجة حرارة سائل الخزان، (T(t))، لأنها تحدد الحالة الحرارية لتيار العملية الخارج من السخان.

يُفترض أن الخزان ممزوج تمامًا، ولذلك فإن درجة حرارة الخروج تساوي درجة حرارة الخزان:

[ T_{\text{out}}=T. ]

تحديد مدخلات العملية

المدخل المتلاعب فيه الرئيسي هو معدل التسخين، (Q(t))، المزوَّد عبر البخار أو السخان الكهربائي أو الغلاف الحراري.

تُعامل درجة حرارة الدخول (T_i(t)) عادةً باعتبارها مدخل اضطراب، رغم أنه يمكن أيضًا التحكم فيها في بعض تكوينات الوحدات التجريبية. وقد يعمل معدل التدفق (F(t)) أيضًا كاضطراب أو كمتغير متلاعب فيه، لأنه يغير كلًا من زمن المكوث ومعدل انتقال الطاقة.

تحديد افتراضات النمذجة

يفترض النموذج الأولي المفيد عمومًا ما يلي:

• خلط السائل خلطًا تامًا
• ثبات حجم السائل
• ثبات الكثافة (\rho)
• ثبات السعة الحرارية (c_p)
• ثبات تدفقي الدخول والخروج، (F_i=F)
• إهمال الفواقد الحرارية ما لم يتم نمذجتها صراحةً
• إهمال تأثيرات الطاقة الحركية وطاقة الوضع وشغل العمود وشغل التمدد

تنتج هذه الافتراضات نموذجًا مدمجًا مع الاحتفاظ بالديناميكيات الحرارية السائدة.

صياغة موازِنات الحالة الديناميكية

كتابة موازنة الكتلة الكلية

تكون موازنة حجم السائل العامة

[ \frac{dV}{dt}=F_i-F. ]

وبالنسبة إلى CSTH ذي الحجم الثابت،

[ F_i=F, \qquad \frac{dV}{dt}=0. ]

إذا كانت ديناميكيات المستوى مهمة، فيجب الاحتفاظ بحجم الخزان كمتغير حالة بدلًا من اعتباره ثابتًا.

كتابة موازنة الطاقة

بالنسبة إلى سخان خزان مزود بالتحريك وغير تفاعلي، تكون موازنة الطاقة غير المستقرة

[ \rho c_p V\frac{dT}{dt}

\rho c_p F(T_i-T)+Q. ]

وبالقسمة على (\rho c_p V) نحصل على معادلة حالة درجة الحرارة:

[ \frac{dT}{dt}

\frac{F}{V}(T_i-T) + \frac{Q}{\rho c_p V}. ]

وبشكل مكافئ،

[ \frac{dT}{dt} + \frac{F}{V}T

\frac{F}{V}T_i + \frac{1}{\rho c_p V}Q. ]

يمثل الحد الأول في الطرف الأيمن الطاقة الحملية الداخلة مع التغذية، بينما يمثل (Q) صافي الحرارة المضافة إلى السائل. ويمكن تمثيل معدل تبريد باعتبار (Q<0).

إدراج تأثيرات التفاعل عند الحاجة

إذا كان CSTH يعمل أيضًا كمفاعل CSTR، فقد تتضمن موازنة الطاقة حرارة التفاعل:

[ \frac{dT}{dt}

\frac{F_i}{V}(T_i-T) + \frac{(-\Delta H_r)r_A}{\rho c_p} + \frac{Q}{\rho c_p V}, ]

حيث إن (r_A) هو معدل التفاعل.

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

[ r_A=k_0\exp\left(-\frac{E}{RT}\right)c_A, ]

ويجب أن يتضمن النموذج اللاخطي الكامل أيضًا موازنة المكوّن الخاصة بـ (c_A).

تحويل النموذج إلى متغيرات الانحراف

تعريف نقطة التشغيل في الحالة المستقرة

لتكن قيم الحالة المستقرة هي (T_s) و(T_{i,s) و(Q_s). في الحالة المستقرة،

[ \frac{dT}{dt}=0, ]

ولذلك تحقق نقطة التشغيل

[ 0= \frac{F}{V}(T_{i,s}-T_s) + \frac{Q_s}{\rho c_p V}. ]

تمثل نقطة التشغيل هذه المرجع الذي تُحلَّل حوله التغيرات الديناميكية.

تعريف الانحرافات

عرّف

[ T'=T-T_s, ]

[ T_i'=T_i-T_{i,s}, ]

[ Q'=Q-Q_s. ]

بالنسبة إلى (F) و(V) و(\rho) و(c_p) الثابتة، فإن طرح معادلة الحالة المستقرة من المعادلة الديناميكية يعطي

[ \frac{dT'}{dt} + \frac{F}{V}T'

\frac{F}{V}T_i' + \frac{1}{\rho c_p V}Q'. ]

وباستخدام تعريف زمن المكوث

[ \tau=\frac{V}{F}, ]

يصبح النموذج

[ \boxed{ \frac{dT'}{dt} + \frac{1}{\tau}T'

\frac{1}{\tau}T_i' + \frac{1}{\rho c_p V}Q' } ]

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

ربط النموذج بالمعلمات العامة

يمكن أيضًا كتابة النموذج بالصورة العامة

[ \frac{dT'}{dt}+aT'

\frac{1}{\tau}T_i' + K_QQ', ]

حيث

[ a=\frac{F}{V}=\frac{1}{\tau}, ]

و

[ K_Q=\frac{1}{\rho c_p V}. ]

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

اشتقاق دوال التحويل

تطبيق تحويل لابلاس

بتطبيق تحويل لابلاس، وبافتراض انحراف ابتدائي صفري لدرجة الحرارة، نحصل على

[ sT'(s)+\frac{1}{\tau}T'(s)

\frac{1}{\tau}T_i'(s) + K_QQ'(s). ]

ومن ثم،

[ T'(s)

\frac{1/\tau}{s+1/\tau}T_i'(s) + \frac{K_Q}{s+1/\tau}Q'(s). ]

الحصول على دالة تحويل درجة حرارة الدخول

تكون دالة التحويل من اضطراب درجة حرارة الدخول إلى درجة حرارة الخزان

[ \boxed{ \frac{T'(s)}{T_i'(s)}

\frac{1}{\tau s+1} } ]

وبالتالي تؤثر درجة حرارة الدخول في درجة حرارة الخزان من خلال استجابة ديناميكية من الرتبة الأولى ذات ثابت زمني (\tau=V/F).

الحصول على دالة تحويل معدل التسخين

تكون دالة التحويل من معدل التسخين إلى درجة حرارة الخزان

[ \boxed{ \frac{T'(s)}{Q'(s)}

\frac{K_Q}{s+1/\tau}

\frac{K_T}{\tau s+1} } ]

حيث

[ K_T=\frac{\tau}{\rho c_p V} =\frac{1}{\rho c_p F}. ]

وبالتالي، بالنسبة إلى مدخل معدل حرارة مباشر،

[ \boxed{ \frac{T'(s)}{Q'(s)}

\frac{1}{\rho c_p F}\frac{1}{\tau s+1} } ]

تكون الاستجابة الكلية لدرجة الحرارة مجموع الاستجابتين للاضطراب والمدخل المتلاعب فيه:

[ T'(s)

\frac{1}{\tau s+1}T_i'(s) + \frac{K_T}{\tau s+1}Q'(s). ]

توسيع النموذج ليشمل ديناميكيات المستوى

إضافة مستوى السائل كمتغير حالة

إذا لم يكن مستوى الخزان ثابتًا، فعرّف (h(t)) بأنه ارتفاع السائل و(A) بأنها مساحة المقطع العرضي للخزان. تصبح موازنة الكتلة

[ A\frac{dh}{dt}=F_i-F. ]

ويمكن كتابة موازنة الطاقة كما يلي:

[ Ah\frac{dT}{dt}

F_i(T_i-T) + \frac{Q}{\rho c_p}. ]

يمثل العامل (Ah) حجم السائل. وبالتالي، فإن التغيرات في المستوى تعدّل السعة الحرارية، ومن ثم ديناميكيات درجة الحرارة.

تحديد ضرورة اقتران المستوى

يكون نموذج درجة الحرارة وحدها مناسبًا عندما يحافظ التحكم في المستوى على (V) شبه ثابت، وتكون تغيرات المستوى صغيرة.

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

فهم المعنى التحكمي للنموذج

تفسير الثابت الزمني

الثابت الزمني

[ \tau=\frac{V}{F} ]

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

بعد تغير خطوة، تصل عملية من الرتبة الأولى إلى نحو 63% من تغير درجة الحرارة النهائي بعد ثابت زمني واحد، وتقترب من قيمتها النهائية بصورة تقاربية.

تفسير كسب العملية

بالنسبة إلى تغيرات معدل التسخين، يكون الكسب في الحالة المستقرة

[ K_T=\frac{1}{\rho c_p F}. ]

ينخفض ارتفاع درجة الحرارة الناتج عن مدخل حراري معين مع زيادة معدل التدفق أو الكثافة أو السعة الحرارية.

أما بالنسبة إلى تغيرات درجة حرارة الدخول، فإن كسب الحالة المستقرة يساوي الواحد في نموذج الخواص الثابتة المثالي:

[ \lim_{s\to 0}\frac{T'(s)}{T_i'(s)}=1. ]

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

استخدام النموذج لتصميم التحكم بالتغذية الراجعة

يحدد النموذج العناصر الأساسية اللازمة لحلقة التحكم في درجة الحرارة:

• المتغير المتحكم فيه: درجة حرارة الخزان (T)
• المتغير المتلاعب فيه: معدل التسخين (Q)
• الاضطراب الرئيسي: درجة حرارة الدخول (T_i)
• الثابت الزمني للعملية: (\tau=V/F)
• كسب العملية: (K_T=1/(\rho c_p F))
• القياس: مرسل درجة الحرارة
• عنصر التحكم النهائي: صمام البخار أو السخان أو مشغل معدل تسخين الغلاف

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

فهم المفاضلات

لا تتعامل مع المفاعل اللاخطي باعتباره من الرتبة الأولى دون تحقق

يكون نموذج الرتبة الأولى صالحًا لنظام حراري غير تفاعلي أو ضعيف الاقتران بالقرب من نقطة تشغيل ثابتة.

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

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

إذا فُقدت الحرارة عبر جدران الخزان أو العزل أو الحيز البخاري، فقد يلزم إضافة حد آخر:

[ Q_{\text{loss}}=UA(T-T_{\text{amb}}). ]

وتتضمن موازنة الطاقة عندئذٍ صافي المدخل الحراري

[ Q_{\text{net}}=Q-UA(T-T_{\text{amb}}). ]

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

ميّز بين معدل الحرارة وأمر السخان

ينبغي أن يمثل المتغير (Q) معدل انتقال حرارة فعليًا إذا استُخدم الكسب (1/(\rho c_p V)) مباشرةً.

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

مراعاة ديناميكيات القياس والمشغل

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

إذا كانت هذه التأثيرات مهمة، فقد يأخذ نموذج العملية الصورة

[ G(s)= \frac{K_T e^{-\theta s}} {(\tau s+1)(\tau_m s+1)}, ]

حيث إن (\theta) هو التأخير و(\tau_m) هو الثابت الزمني للقياس أو المشغل.

اختيار النموذج المناسب لهدفك

استخدم تعقيد النموذج الذي يتوافق مع سؤال التحكم وظروف التشغيل.

• إذا كان تركيزك الأساسي هو التحكم الأساسي في درجة الحرارة: استخدم نموذج الرتبة الأولى ذي الحجم الثابت، مع اعتبار (T_i) اضطرابًا و(Q) مدخلًا متلاعبًا فيه.
• إذا كان تركيزك الأساسي هو التحكم في المستوى ودرجة الحرارة: احتفظ بكل من موازنة الكتلة وموازنة الطاقة، مع اعتبار (h) أو (V) و(T) متغيري حالة.
• إذا كان تركيزك الأساسي هو سلوك المفاعل: أضف موازنة المكوّن وحركية التفاعل المعتمدة على درجة الحرارة، ثم أجرِ الخطية فقط حول حالة مستقرة محددة بوضوح.
• إذا كان تركيزك الأساسي هو ضبط وحدة التحكم في الوحدة التجريبية: أدرج تأثيرات المشغل والمستشعر والفقد الحراري والتأخير المقاسة عندما تكون كبيرة بما يكفي للتأثير في الاستجابة الملحوظة.
• إذا كان تركيزك الأساسي هو المحاكاة أو اكتشاف الأعطال: تحقّق من النموذج مقابل بيانات اختبار الخطوة قبل استخدامه لتدريب المشغلين أو التحسين أو التشخيص.

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

جدول الملخص:

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

المنتجات ذات الصلة

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

يسأل الناس أيضًا

• تعدد حالات CSTR: أي المتغيرات الفيزيائية تتوافق مع التوسع الشامل؟
• كيف يؤثر تكوين المروحة على توزيع أوقات الإقامة (RTD) وأداء المفاعل؟ قم بتحسين المزج والتحويل في وحدات التجارب النموذجية.
• كيف تؤثر عدة مفاعلات خلط مستمر في سلسلة على الاستقرار؟ أتقن ديناميكيات المصنع النموذجي المعقد
• كيف يمثل نموذج MRF دوران المحرك المائي (الإمبلر) في ديناميكا الموائع الحسابية (CFD): شرح شروط الحدود الرئيسية
• كيف يمكن لمصنع تجريبي CSTR إثبات التقسيم المرحلي للمفاعل؟ تصور وفورات الحجم.