كيف تُطبَّق معادلات الاتزان الطوري الديناميكي الحراري وكثافة البخار في وحدات الفلاش التجريبية؟

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

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

كيف يُنشئ وعاء الفلاش طورين

يؤدي خفض الضغط إلى تبخر جزئي

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

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

يحدد الوعاء حالة تشغيل واحدة

في حالة فلاش اتزان جيدة الخلط، يُفترض أن يكون للطورين نفس درجة الحرارة والضغط:

[ T^L=T^V=T ]

[ P^L=P^V=P ]

تصف هذه المساويات الاتزان الحراري والميكانيكي. لكنها لا تحدد الاتزان الكيميائي بحد ذاتها.

يحدد الاتزان الكيميائي التركيب

لكل مكوّن (i)، يتطلب الاتزان الكيميائي تساوي الكمون الكيميائي، ويُعبَّر عنه عادةً بتساوي الهروبية:

[ f_i^L=f_i^V ]

في النظام المثالي أو المبسط، يُمثَّل ذلك غالبًا باستخدام قيمة (K):

[ K_i=\frac{y_i}{x_i} ]

حيث إن (x_i) هو الكسر المولي في الطور السائل، و(y_i) هو الكسر المولي في الطور البخاري.

أما في النموذج الأكثر دقة، فتُستخدم معاملات الهروب ومعاملات النشاط:

[ f_i^V=y_i\phi_i^V P ]

[ f_i^L=x_i\gamma_i f_i^{\text{reference}} ]

ويعتمد الشكل المناسب على الضغط ودرجة الحرارة ولا مثالية المزيج.

كيفية تطبيق معادلات اتزان البخار والسائل

يحدد حساب الفلاش انقسام الطورين

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

الكسر البخاري
الكسر السائل
تركيب البخار
تركيب السائل
درجة حرارة الاتزان عندما تكون مجهولة

تكون موازين المكونات كما يلي:

[ Fz_i=Vy_i+Lx_i ]

حيث إن (F) هو معدل تدفق التغذية، و(V) هو معدل تدفق البخار، و(L) هو معدل تدفق السائل، و(z_i) هو تركيب التغذية، و(y_i) هو تركيب البخار، و(x_i) هو تركيب السائل.

ويجب أن تحقق تراكيب الطورين أيضًا:

[ \sum_i x_i=1 ]

[ \sum_i y_i=1 ]

وبالاقتران مع العلاقة (y_i=K_i x_i)، تحدد هذه العلاقات مسألة فلاش الاتزان.

تحدد معادلة راشفورد-رايس الكسر البخاري

عند تحديد درجة الحرارة والضغط، يمكن حساب الكسر البخاري (\beta=V/F) من خلال:

[ \sum_i \frac{z_i(K_i-1)} {1+\beta(K_i-1)}=0 ]

ثم تُحسب تراكيب السائل والبخار من خلال:

[ x_i=\frac{z_i}{1+\beta(K_i-1)} ]

[ y_i=K_i x_i ]

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

توضح قيم (K) المكونات التي تتبخر

يدخل المكوّن الذي له (K_i>1) إلى الطور البخاري بشكل تفضيلي. أما المكوّن الذي له (K_i<1) فيبقى بشكل تفضيلي في الطور السائل.

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

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

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

أما بالنسبة إلى المزيجات السائلة غير المثالية، فقد تكون هناك حاجة إلى نماذج مثل NRTL أو UNIQUAC أو UNIFAC. وعند الضغط المرتفع أو في الأنظمة التي تحتوي على غازات مذابة، تكون معادلة حالة مكعبة مثل Peng-Robinson أو Soave-Redlich-Kwong أكثر ملاءمة غالبًا.

كيفية تطبيق معادلات كثافة البخار

العلاقة الأساسية لكثافة البخار

بالنسبة إلى مزيج بخاري من الغاز المثالي، تُحسب الكثافة كما يلي:

[ \rho_V=\frac{P\overline{M}_V}{RT} ]

ويكون الوزن الجزيئي لمزيج البخار:

[ \overline{M}_V=\sum_i y_iM_i ]

وبالنسبة إلى مزيج ثنائي، تصبح المعادلة:

[ \rho_V=\frac{P}{RT}\left(y_A M_A+y_B M_B\right) ]

هنا، (P) هو الضغط المطلق، و(T) هي درجة الحرارة المطلقة، و(R) هو ثابت الغاز، و(M_i) هو الوزن الجزيئي للمكوّن (i).

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

يتطلب سلوك الغاز الحقيقي معامل انضغاط

عند الضغوط المرتفعة، قد تكون معادلة الغاز المثالي غير دقيقة. وتكون المعادلة الأكثر عمومية كما يلي:

[ \rho_V=\frac{P\overline{M}_V}{ZRT} ]

حيث إن (Z) هو معامل انضغاط البخار المستخرج من معادلة الحالة.

وتأخذ قيمة (Z) المختلفة عن الواحد في الاعتبار التأثيرات بين الجزيئية. ويصبح هذا التصحيح مهمًا عندما يعمل الوعاء التجريبي عند ضغط مرتفع أو بالقرب من نقطة الندى للمزيج أو من المنطقة الحرجة.

تربط الكثافة بين الديناميكا الحرارية وسلوك الوعاء

تتيح كثافة البخار تحويل المخزون المولي أو الكتلي المحسوب للبخار إلى حجم:

[ m_V=\rho_V V_V ]

حيث إن (V_V) هو حجم حيز البخار.

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

تدعم كثافة البخار تفسير قراءات الأجهزة

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

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

كيف يُكمل اتزان الطاقة النموذج

لا يتحدد التبخر بالضغط وحده

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

بالنسبة إلى فلاش أديباتي، يُكتب اتزان الطاقة الكلي عادةً كما يلي:

[ Fh_F=Vh_V+Lh_L ]

حيث إن (h_F) و(h_V) و(h_L) هي الإنثالبيات النوعية أو المولية للتغذية والبخار والسائل.

إذا كان الوعاء يحتوي على سخان أو مبرد أو فقد حراري ملحوظ، فإن الاتزان يتضمن الحمل الحراري المقابل:

[ Fh_F+Q=Vh_V+Lh_L ]

تقع درجة حرارة الفلاش ضمن منطقة تغير الطور

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

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

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

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

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

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

بدء التشغيل وتغيرات الضغط

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

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

حسابات مستوى السائل والاحتجاز

يغير تدفق خروج السائل مخزون السائل:

[ \frac{dN_L}{dt}=F_L-L_{\text{out}} ]

وبالمثل، يتغير مخزون البخار وفقًا لتبخر التغذية وتدفق خروج البخار:

[ \frac{dN_V}{dt}=V_{\text{generated}}-V_{\text{out}} ]

وتربط كثافة البخار بين (N_V) وحجم البخار المشغول وضغط التشغيل. وتؤدي كثافة السائل وهندسة الوعاء دورًا مماثلًا في حسابات مستوى السائل.

تغيرات التركيب

يغير التغير في ضغط الوعاء قيم (K_i). ويغير ذلك انقسام الطورين والتركيبين (x_i) و(y_i).

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

قياسات أداء الفصل

تشمل مؤشرات الأداء الرئيسية ما يلي:

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

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

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

قد تبالغ افتراضات الاتزان في تقدير الفصل

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

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

قد تكون كثافة الغاز المثالي غير كافية

تُعد العلاقة

[ \rho_V=\frac{P\overline{M}_V}{RT} ]

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

تنتقل أخطاء نموذج الخواص عبر الحساب بأكمله

تؤثر قيمة (K) غير الدقيقة في الكسر البخاري وفي تركيبي الطورين. ويغير التركيب الناتج الوزن الجزيئي للبخار، مما يؤثر بدوره في تقدير كثافة البخار.

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

قد يؤدي التبخر الأكبر إلى تقليل استرداد السائل

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

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

تحد الأجهزة من جودة التحقق

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

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

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

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

إذا كان تركيزك الأساسي هو التنبؤ بانقسام الطورين: استخدم موازين المكونات، وقيم (K) أو تساوي الهروبية، ومعادلة راشفورد-رايس عند ضغط ودرجة حرارة الوعاء المقاسين أو المحددين.
إذا كان تركيزك الأساسي هو مخزون البخار وديناميكيات الضغط: احسب الوزن الجزيئي للبخار من تركيب البخار عند الاتزان، واستخدم (\rho_V=P\overline{M}_V/(ZRT)) لتقدير احتجاز البخار.
إذا كان تركيزك الأساسي هو استهلاك الطاقة: اجمع بين حساب فلاش الاتزان وإنثالبيات التغذية والبخار والسائل لتحديد الحمل المطلوب للسخان أو المبرد.
إذا كان تركيزك الأساسي هو أداء التقطير اللاحق: استخدم تراكيب ومعدلات تدفق البخار والسائل المتنبأ بها بوصفها ظروف التغذية لمرحلة الفصل التالية.
إذا كان تركيزك الأساسي هو التحقق من النموذج: قارن درجات الحرارة والضغوط ومعدلات تدفق الطورين والكثافات والتراكيب المتنبأ بها بقياسات الوحدة التجريبية عبر ظروف تشغيل متعددة.

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

جدول الملخص:

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

صمّم وحدة الفلاش التجريبية وتحقق من أدائها

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

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