كيف تتعامل النماذج الديناميكية الحرارية مع اتزان الطور فوق الحرج؟ قانون هنري، وطرق معادلات الحالة، والتحقق

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

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

لماذا تتطلب المكونات فوق الحرجة معالجة خاصة

انهيار مفهوم التشبّع

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

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

يعتمد اتزان الطور على الجهد الكيميائي

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

يسعى حساب البخار–السائل لمكوّن (i) عمومًا إلى تحقيق شرط من الشكل:

[ f_i^{(V)} = f_i^{(L)} ]

حيث إن (f_i) هي فيوغاسية المكوّن في طور البخار أو السائل.

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

كيف يختار برنامج المحاكاة البديل الديناميكي الحراري

الخيار الأول: ثوابت قانون هنري

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

الصيغة المبسطة لقانون هنري هي:

[ f_i = H_i x_i ]

حيث إن (H_i) هو ثابت قانون هنري، و(x_i) هو الكسر المولي للمكوّن في الطور السائل.

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

لماذا يكون قانون هنري أكثر ملاءمة عادةً

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

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

الخيار الثاني: العلاقات التجريبية لضغط البخار

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

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

تحديد حد أقصى لمعادلة أنطوان

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

[ T_r = \frac{T}{T_c} > 2 ]

فقد يقيّم البرنامج معادلة ضغط البخار عند (T_r = 2)، بدلاً من السماح للعلاقة بالاستقراء إلى ما لا نهاية.

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

كيف تمثل نماذج الاتزان الأطوار المتبقية

معادلات الحالة التكعيبية

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

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

نماذج معامل النشاط

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

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

النماذج المدمجة

تستخدم العديد من عمليات المحاكاة العملية نهجًا مدمجًا:

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

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

ما الذي يحدث في الأنظمة متعددة الأطوار والمتفاعلة

يجب أن يتحقق اتزان الطور والاتزان الكيميائي معًا

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

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

[ \sum_i \nu_i \mu_i = 0 ]

حيث إن (\nu_i) هو المعامل الستوكيومتري، و(\mu_i) هو الجهد الكيميائي للنوع الكيميائي (i).

يمكن للأنواع فوق الحرجة التأثير في تنبؤات التفاعل

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

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

التقارب العددي ليس تحققًا فيزيائيًا

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

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

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

قانون هنري محلي وليس شاملًا

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

وقد يصبح غير ملائم عند التركيزات العالية للمذاب، أو بالقرب من انتقالات الطور، أو عندما يغيّر الترابط القوي، أو التفاعل الكيميائي، أو تمدد المذيب الملحوظ سلوك المحلول.

يمكن أن يكون استقراء ضغط البخار مضللًا

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

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

تضيف تقديرات UNIFAC قدرًا من عدم اليقين

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

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

تظل بيانات المصنع التجريبي مهمة

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

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

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

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

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

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

جدول الملخص:

الموضوع	الخلاصة الأساسية
السلوك فوق الحرج	يكون ضغط بخار التشبّع التقليدي غير معرّف فوق النقطة الحرجة.
البديل الأساسي	يُستخدم قانون هنري عادةً لتقدير فيوغاسية المذاب المخفف أو ذوبانيته.
البديل الآخر	قد توفر علاقات ضغط البخار تقديرات عددية، لكنها تتطلب الحذر.
تحديد حد أنطوان	تحد بعض البرامج درجة الحرارة المختزلة إلى 2 لمنع الاستقراء غير المستقر.
نماذج الخصائص	تلائم SRK وبنغ–روبنسون الأنظمة ذات الضغط المرتفع؛ بينما تعالج NRTL وUNIQUAC وUNIFAC عدم مثالية السائل.
الأنظمة التفاعلية	يجب تحقيق اتزان الطور والاتزان الكيميائي في الوقت نفسه.
التحقق	لا يُعد التقارب دليلًا على الدقة؛ بل يجب مقارنة التنبؤات ببيانات الذوبانية والطور وبيانات المصنع التجريبي.

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

القائمة