تتم خطيّة موازنة الكتلة الديناميكية باستبدال علاقة اتزان البخار والسائل غير الخطية بصيغتها الخطية المحلية، ثم تجميع حدود التركيب. بالنسبة إلى المرحلة (n)، تكون موازنة مكوّن السائل
[ H_n\frac{dx_n}{dt}
Lx_{n-1}-Lx_n+Vy_{n+1}-Vy_n. ]
إذا قُرّبت علاقة الاتزان محليًا بواسطة (y_i=m_i x_i)، مع ميل ثابت (m_i)، فإن معادلة المرحلة تصبح معادلة تفاضلية خطية مترابطة مع المراحل المجاورة. ويشكّل هذا الشكل أساس اشتقاق دوال النقل، ودراسة انتشار الاضطرابات، وتصميم أنظمة التحكم في التركيب لمصانع الامتصاص والتقطير التجريبية.
الخلاصة الأساسية: تحوّل الخطيّة موازنة المرحلة إلى نموذج ديناميكي من الرتبة الأولى من خلال اعتبار الاحتجاز ومعدلات التدفق وميل الاتزان المحلي ثابتة. ويكون النموذج الناتج عمومًا مجموعة مترابطة من المعادلات من الرتبة الأولى، وليس معادلة مستقلة لكل مرحلة، لأن كل مرحلة تتبادل المادة مع المراحل المجاورة.
البدء بموازنة المرحلة الديناميكية
تراكم مكوّن السائل
بالنسبة إلى احتجاز سائل ثابت (H_n)، يكون تراكم المكوّن الممتص أو المقطر في المرحلة (n) كما يلي
[ H_n\frac{dx_n}{dt}. ]
هنا، يمثّل (x_n) الكسر المولي للمكوّن في الطور السائل في المرحلة (n). يجب أن تكون الوحدات متسقة: فإذا كان (H_n) احتجازًا موليًا، فينبغي أن يكون كل من (L) و(V) معدلي تدفق مولي.
مساهمات التدفق
يساهم تيار السائل الداخل إلى المرحلة (n) بالحد (Lx_{n-1})، بينما يساهم تيار السائل الخارج بالحد (Lx_n).
ويساهم تيار البخار الداخل من المرحلة المجاورة بالحد (Vy_{n+1})، في حين يزيل البخار الخارج من المرحلة (n) الحد (Vy_n). لذلك،
[ H_n\frac{dx_n}{dt}
Lx_{n-1}-Lx_n+Vy_{n+1}-Vy_n. ]
تفترض هذه الموازنة ثبات (L)، وثبات (V)، وثبات الاحتجاز، وعدم وجود تفاعل، وإهمال التراكم في طور البخار.
تطبيق تقريب الاتزان المحلي
علاقة الاتزان الخطية
يرتبط تركيبا البخار والسائل بواسطة دالة الاتزان
[ y=f(x). ]
حول نقطة تشغيل محددة، يمكن تقريب هذه العلاقة كما يلي
[ y_i \approx y_{i,0}+m_i(x_i-x_{i,0}), ]
حيث
[ m_i=\left.\frac{dy}{dx}\right|{x{i,0}}. ]
في نموذج مبسط ذي ميل ثابت لخط الاتزان، ومع قياس التراكيب بالنسبة إلى نقطة التشغيل، تختزل هذه العلاقة إلى
[ y_i=m x_i. ]
وقد يمثّل الميل (m) علاقة تطايرية محلية أو علاقة قانون هنري، تبعًا للنظام.
التعويض في موازنة المرحلة
باستخدام (y_n=mx_n) و(y_{n+1}=mx_{n+1})،
[ H_n\frac{dx_n}{dt}
Lx_{n-1}-Lx_n+mVx_{n+1}-mVx_n. ]
وبتجميع الحدود التي تتضمن (x_n) نحصل على
[ \boxed{ H_n\frac{dx_n}{dt}
Lx_{n-1}
(L+mV)x_n + mVx_{n+1} } ]
أو، بعد القسمة على الاحتجاز،
[ \boxed{ \frac{dx_n}{dt}
\frac{L}{H_n}x_{n-1}
\frac{L+mV}{H_n}x_n + \frac{mV}{H_n}x_{n+1} }. ]
هذه هي المعادلة الديناميكية الخطية القياسية لمرحلة اتزان مثالية في ظل افتراضات التدفق الثابت.
ماذا تعني «الخطيّة» في هذا النموذج؟
خطيّة علاقة الاتزان
قد تكون علاقة الاتزان الأصلية غير خطية، ولا سيما في أنظمة التقطير ذات التغير الملحوظ في التطايرية النسبية. تستبدل الخطيّة تلك العلاقة غير الخطية بمماسها بالقرب من نقطة التشغيل المستقرة المحددة.
يكون هذا التقريب موثوقًا به بالنسبة إلى تغيرات التركيب الصغيرة نسبيًا حول تلك النقطة. وقد تتطلب الاضطرابات الكبيرة استخدام علاقة اتزان البخار والسائل غير الخطية بالكامل.
خطيّة موازنة المادة
تكون موازنة المادة نفسها خطية بالفعل في (x) عندما تكون (H) و(L) و(V) ثابتة. وتنشأ اللاخطيات الإضافية عندما تتغير هذه الكميات مع التركيب أو درجة الحرارة أو الضغط أو المنسوب أو التدفق الكلي.
بالنسبة إلى نموذج أكثر عمومية، نكتب
[ H_n(x_n)\frac{dx_n}{dt}
L(x_{n-1}-x_n) + V(y_{n+1}-y_n), ]
ثم نوسّع النموذج حول حالة مستقرة:
[ x_i=x_{i,0}+\delta x_i, \qquad y_i=y_{i,0}+\delta y_i. ]
ويُقرّب اضطراب الاتزان بواسطة
[ \delta y_i=m_i\delta x_i. ]
ويكون لنموذج الاضطراب الناتج البنية نفسها للمعادلة المؤطرة، لكن تُقيّم معاملاته عند نقطة التشغيل.
تصحيح مهم للفهارس
يجب الاحتفاظ بتركيب البخار المجاور
من الأخطاء الشائعة في الفهرسة التعويض بـ (y_{n+1}=mx_{n-1}). وفق ترقيم المراحل المعتاد في الجريان المعاكس، يكون البخار الداخل إلى المرحلة (n) هو (y_{n+1})، ويرتبط بتركيب السائل في المرحلة (n+1):
[ y_{n+1}=mx_{n+1}. ]
لذلك، تكون المعادلة المفهرسة بشكل صحيح
[ H_n\frac{dx_n}{dt}
Lx_{n-1}
(L+mV)x_n + mVx_{n+1}. ]
أما التعبير
[ H\frac{dx_n}{dt}
Lx_{n-1}-(L+mV)x_n+mVx_{n-1} ]
فلا يكون صالحًا إلا إذا عُرّف المتغير المسمى (x_{n-1}) عمدًا على أنه تركيب السائل المكافئ upstream للبخار الداخل. وهذا ليس تعويض الاتزان المعتاد المفهرس حسب المرحلة.
سبب أهمية هذا التمييز
يؤدي الاحتفاظ بالحد (x_{n+1}) إلى إنشاء نموذج مترابط للعمود. إذ يستجيب التركيب في إحدى المراحل لكل من المرحلتين المجاورتين، وهو أمر ضروري لتمثيل انتشار الاضطرابات عبر عمود ذي صحون أو عمود مبطن مقسم إلى خلايا.
قد يؤدي استبدال حد المرحلة المجاورة بالمتغير upstream نفسه إلى إنشاء نموذج أبسط لمرحلة واحدة، لكنه يغيّر الاتصال الفيزيائي وقد يعطي تنبؤات ديناميكية غير صحيحة.
من معادلات المراحل إلى دوال النقل
الصيغة في مجال لابلاس
بافتراض اضطرابات ابتدائية معدومة ومعاملات ثابتة، يؤدي أخذ تحويل لابلاس إلى
[ H_n sX_n(s)
LX_{n-1}(s)
(L+mV)X_n(s) + mVX_{n+1}(s). ]
وبإعادة الترتيب،
[ \boxed{ \left[H_n s+(L+mV)\right]X_n(s)
LX_{n-1}(s)+mVX_{n+1}(s) }. ]
تمثّل هذه المعادلة صفًا واحدًا من نموذج دالة النقل الديناميكية للعمود أو نموذج فضاء الحالة.
تقريب المرحلة الواحدة
إذا عومل تركيب البخار الداخل على أنه دخل معلوم بدلًا من كونه حالة مرتبطة بـ (x_{n+1})، فليكن
[ u_n=Lx_{n-1}+Vy_{n+1}. ]
تصبح الموازنة
[ H_n\frac{dx_n}{dt}
u_n-(L+mV)x_n. ]
وتكون دالة النقل المقابلة
[ \boxed{ \frac{X_n(s)}{U_n(s)}
\frac{1}{H_n s+L+mV} } ]
مع ثابت زمني
[ \boxed{ \tau_n=\frac{H_n}{L+mV} }. ]
إذا كان الدخل هو تركيب السائل upstream فقط، وأُهمل إسهام البخار أو أُدرج بشكل منفصل، فإن
[ \frac{X_n(s)}{X_{n-1}(s)}
\frac{L}{H_n s+L+mV}. ]
التطبيق على الامتصاص والتقطير
امتصاص الغاز والسائل
في الامتصاص، ينتقل المذاب عمومًا من الطور الغازي إلى الطور السائل عندما يتجاوز تركيب الغاز الفعلي قيمة الاتزان الخاصة به:
[ y>y^*(x). ]
ويصف ميل الاتزان (m) كيفية تغير تركيب الغاز عند الاتزان مع تركيب السائل. ويزيد الحد (mV) الأكبر من قوة الاقتران في طور البخار ويغير استجابة المرحلة.
في أعمدة الامتصاص المبطنة، يُكتب المفهوم نفسه غالبًا في صورة نموذج تفاضلي موزع بدلًا من موازنات الصحون المنفصلة. وينتج عن التقسيم الحجمي المحدود معادلات شبيهة بمعادلات المراحل، مع الاحتجاز والاقتران بين الخلايا المجاورة.
التقطير
في التقطير، تُمثّل علاقة الاتزان عادةً بواسطة
[ y_i=y^*(x_i), ]
حيث قد تكون (y^*(x)) غير خطية. وحول تركيب تشغيل معين،
[ \delta y_i \approx \left.\frac{dy^*}{dx}\right|{x{i,0}}\delta x_i. ]
وقد يختلف الميل المحلي من مرحلة إلى أخرى، ولذلك يكون لعمود التقطير عادةً معاملات (m_n) و(L_n) و(V_n) و(H_n) بدلًا من قيمة واحدة شاملة لكل من (m) و(L) و(V).
الشروط الحدية
لا تكفي معادلات المراحل الداخلية بمفردها. إذ يوفر المكثف والغلاية ومرحلة التغذية وتيار الراجع وتيار المقطر وتيار القاع ونهايتا العمود الشروط الحدية وعلاقات المتغيرات المتحكم فيها.
فعلى سبيل المثال، يغير الراجع تدفق السائل الداخل إلى المراحل العلوية، بينما يغير الغليان تدفق البخار الداخل إلى المراحل السفلية. وغالبًا ما تكون هذه التأثيرات الحدية المسارات الرئيسية لاضطرابات التحكم.
توسيع النموذج ليشمل الاحتجاز والتدفق المتغيرين
موازنة الكتلة الكلية
عندما يتغير منسوب السائل أو الاحتجاز، يجب تضمين موازنة كلية. وتكون موازنة احتجاز السائل العامة
[ \frac{dH_n}{dt}
L_{n-1}-L_n+V_{n+1}-V_n, ]
مع تعديل الإشارات وتعريفات التيارات وفقًا لاتفاقية المرحلة المختارة.
بعد ذلك، ينبغي كتابة موازنة المكوّن بدلالة احتجاز المكوّن:
[ \frac{d(H_nx_n)}{dt}
\text{المكوّن الداخل}
\text{المكوّن الخارج}. ]
وبتوسيع المشتقة نحصل على
[ H_n\frac{dx_n}{dt} + x_n\frac{dH_n}{dt}
\text{المكوّن الداخل}
\text{المكوّن الخارج}. ]
ويكون الحد الإضافي (x_n,dH_n/dt) مهمًا عندما تتغير مستويات سائل الصحون أو احتجازات الغلاية والمكثف بشكل ملحوظ.
تأثيرات الطاقة والضغط
في نماذج المصانع التجريبية الدقيقة، تؤثر درجة الحرارة في كل من الاتزان ومعدلات التدفق. ويؤثر الضغط في تركيب البخار وظروف التشبع والتدفق الحجمي.
لذلك، تكون الخطيّة المعتمدة على التركيب فقط مناسبة أساسًا للاضطرابات المعتدلة حول نقطة تشغيل ذات درجة حرارة وضغط ثابتين. أما التغيرات التشغيلية الأكبر فتتطلب معادلات مترابطة لموازنة الكتلة والطاقة والضغط.
فهم المفاضلات
فوائد النموذج الخطي
النموذج الخطي مدمج وفعال حسابيًا ومناسب لتحليل دوال النقل. كما يوفر فهمًا مباشرًا للثوابت الزمنية، والاقتران بين المراحل، وتوهين الاضطرابات، وقابلية التحكم.
وهو مفيد خصوصًا لضبط المتحكمات وتجارب المصانع التجريبية التعليمية التي يعمل فيها العمود بالقرب من حالة اسمية مستقرة.
حدود ميل الاتزان الثابت
قد تكون قيمة واحدة لـ (m) غير كافية عندما تتغير التطايرية النسبية بدرجة كبيرة مع التركيب أو درجة الحرارة. وقد يفشل التقريب أيضًا بالقرب من الأزيوتروبات أو في المخاليط شديدة اللاتثالية أو ضمن نطاقات التشغيل الواسعة.
بالنسبة إلى هذه الأنظمة، استخدم ميولًا محلية خاصة بكل مرحلة أو احتفظ بدالة الاتزان غير الخطية.
حدود التدفق والاحتجاز الثابتين
يفترض اعتبار (L) و(V) و(H) ثابتة استبعاد ديناميكيات المنسوب والتأثيرات الهيدروليكية. وقد يؤدي ذلك إلى تنبؤات متفائلة بشأن سرعة الاستجابة، كما قد يخفي تفاعلات مهمة بين التحكم في التدفق والضغط والتركيب.
ينبغي للمصانع التجريبية استخدام احتجازات مقاسة أو نماذج مستقلة للاحتجاز عندما يكون فيضان الصحون أو تسربها أو تكوّن الرغوة أو التغير الكبير في مستوى السائل ممكنًا.
الخلط بين خطوط التشغيل والموازنات الديناميكية
يصف خط التشغيل علاقة موازنة مستقرة للمادة بين تراكيب الطورين. ولا يصف، بمفرده، كيفية تطور التراكيب مع الزمن.
يتطلب النموذج الديناميكي حدود الاحتجاز مثل (H,dx/dt). ويؤدي حذف هذه الحدود إلى تحويل النموذج إلى علاقة مستقرة بدلًا من علاقة عابرة.
اختيار النموذج المناسب لهدفك
يعتمد النموذج المناسب على ما إذا كان الهدف هو تصميم التحكم أو التحليل المفاهيمي أو التنبؤ عالي الدقة بأداء المصنع التجريبي.
- إذا كان تركيزك الأساسي هو تصميم المتحكم: خطِّ (y^*(x)) حول نقطة التشغيل المقصودة، واحتفظ بالاقتران بين المراحل المجاورة، واجمع المعادلات في صيغة فضاء الحالة.
- إذا كان تركيزك الأساسي هو استجابة بسيطة لمرحلة واحدة: عامل تركيب البخار الداخل على أنه دخل معلوم، واستخدم دالة النقل من الرتبة الأولى (1/[Hs+(L+mV)]).
- إذا كان تركيزك الأساسي هو انتشار الاضطرابات عبر العمود بأكمله: احتفظ بـ (x_{n-1}) و(x_{n+1}) في كل معادلة مرحلة، وطبّق الشروط الحدية الصحيحة للمكثف والغلاية والتغذية والراجع والمنتج.
- إذا كان تركيزك الأساسي هو العابر الكبير في المصنع التجريبي: استخدم الاتزان غير الخطي، والاحتجاز المتغير، ومعدلات التدفق المتغيرة، وموازنات الطاقة والضغط المترابطة.
يحافظ نموذج العمود الخطي السليم على الاتصال الفيزيائي بين المراحل، مع تبسيط سلوك الاتزان المحلي حول نقطة تشغيل محددة بوضوح.
جدول الملخص:
| العنصر الرئيسي | رؤية النموذج الخطي |
|---|---|
| الموازنة الديناميكية | H_n dx_n/dt = Lx_{n-1} - Lx_n + Vy_{n+1} - Vy_n |
| تقريب الاتزان | استبدال y=f(x) بـ delta y=m delta x بالقرب من نقطة تشغيل مستقرة |
| معادلة المرحلة المترابطة | H_n dx_n/dt = Lx_{n-1} - (L+mV)x_n + mVx_{n+1} |
| الافتراضات الرئيسية | ثبات الاحتجاز وتدفقات السائل والبخار، وإهمال التراكم في طور البخار |
| صيغة دالة النقل | [H_n s + (L+mV)]X_n = LX_{n-1} + mVX_{n+1} |
| الثابت الزمني للمرحلة الواحدة | tau_n = H_n/(L+mV) |
| نقطة الفهرسة الحرجة | يرتبط حد البخار الداخل بـ x_{n+1}، مما يحافظ على الاقتران بين المراحل المجاورة |
| أفضل التطبيقات | ضبط المتحكمات، وتحليل الاضطرابات، ودوال النقل، ونمذجة فضاء الحالة |
| متى تكون النماذج غير الخطية ضرورية | الاضطرابات الكبيرة، أو الاحتجاز أو التدفق المتغير، أو اللاتثالية الشديدة، أو تغيرات الضغط، أو تأثيرات الطاقة |
هل تحتاج إلى مصنع تجريبي عملي للتدريس أو البحث أو تطوير العمليات؟ تواصل مع LABPARK لمناقشة مصانع الوحدات التشغيلية التجريبية التعليمية والمهنية لامتصاص الغاز والسائل، والتقطير، والعمليات الحيوية والتقنيات الحيوية، وتطبيقات المعالجة البيئية ومعالجة المياه. يمكن لـ LABPARK مساعدة الجامعات ومعاهد البحث والمؤسسات في اختيار معدات تدعم التدريب العملي، وعروض العمليات، وضبط المتحكمات، ودراسات موازنة الكتلة الديناميكية، والعمل التجريبي القابل للتوسّع. شارك عمليتك وسعتها وظروف التشغيل وأهدافك التعليمية أو البحثية، وسيوصي فريقنا بتكوين مناسب للمصنع التجريبي.
المنتجات ذات الصلة
- وحدة تدريب عمليات التقطير متعددة الوسائط - مصنع تجريبي
- وحدة تجريبية تعليمية للتقطير المستمر بالصواني الغربالية لمختبر عمليات الوحدات
- وحدة الامتصاص بالملء التعليمية لعمليات الوحدة النموذجية
- وحدة تجريبية للعمليات الصناعية بنظام التقطير ثنائي النمط
- مصنع تجريبي متعدد الوسائط للامتصاص والامتصاص العكسي لتدريب العمليات الوحدوية
يسأل الناس أيضًا
- كيف يمكن التحقق من صحة معادلات توازن الطاقة للمكثفات وأجهزة إعادة الغليان تجريبياً؟ دليل
- كيف تقوم محطات التقطير التجريبية بالتحقق من صحة بيانات التوازن بين البخار والسائل؟ وتحقق توسيع نطاق فصل المذيبات بشكل موثوق.
- كيف يؤثر اختيار النموذج الديناميكي الحراري على محطات التقطير التجريبية ذات الأزورتروب الحساس للضغط؟ التصميم باستخدام توازن البخار والسائل (VLE) المعتمد.
- لماذا يعد مخطط أداء العمود حاسماً؟ لتحسين سلامة وكفاءة وحدة التقطير النموذجية.
- كيف يمكن للتنبؤ بالعدد الكربوني في الوقت الفعلي أن يحسن مصانع التقطير التجريبية؟ تحسين التحكم.