Simulation of a steady state flash

Similar documents
General Thermodynamics for Process Simulation. Dr. Jungho Cho, Professor Department of Chemical Engineering Dong Yang University

Energy, Entropy, and Availability Balances Phase Equilibria. Nonideal Thermodynamic Property Models. Selecting an Appropriate Model

PETE 310 Lectures # 24 & 25 Chapter 12 Gas Liquid Equilibrium

Equation of State Modeling of Phase Equilibrium in the Low-Density Polyethylene Process

10.34 Numerical Methods Applied to Chemical Engineering Fall Homework #3: Systems of Nonlinear Equations and Optimization

I wish to publish my paper on The International Journal of Thermophysics. A Practical Method to Calculate Partial Properties from Equation of State

A Self-Consistent Gibbs Excess Mixing Rule for Cubic Equations of State: derivation and fugacity coefficients

Introduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:

y i x P vap 10 A T SOLUTION TO HOMEWORK #7 #Problem

Supplementary Notes for Chapter 9 Mixture Thermodynamics

The ChemSep Book. Harry A. Kooijman Consultant. Ross Taylor Clarkson University, Potsdam, New York University of Twente, Enschede, The Netherlands

Implicit Integration Henyey Method

Adiabatic Sorption of Ammonia-Water System and Depicting in p-t-x Diagram

Process Modeling. Improving or understanding chemical process operation is a major objective for developing a dynamic process model

( ) 1/ 2. ( P SO2 )( P O2 ) 1/ 2.

between standard Gibbs free energies of formation for products and reactants, ΔG! R = ν i ΔG f,i, we

LNG CARGO TRANSFER CALCULATION METHODS AND ROUNDING-OFFS

CinChE Problem-Solving Strategy Chapter 4 Development of a Mathematical Model. formulation. procedure

Non-Commercial Use Only

Irreversibility of Processes in Closed System

Lecture. Polymer Thermodynamics 0331 L Chemical Potential

Estimation of the composition of the liquid and vapor streams exiting a flash unit with a supercritical component

3. Be able to derive the chemical equilibrium constants from statistical mechanics.

UNIFAC. Documentation. DDBSP Dortmund Data Bank Software Package

Not at Steady State! Yes! Only if reactions occur! Yes! Ideal Gas, change in temperature or pressure. Yes! Class 15. Is the following possible?

University of Washington Department of Chemistry Chemistry 452/456 Summer Quarter 2014

If two volatile and miscible liquids are combined to form a solution, Raoult s law is not obeyed. Use the experimental data in Table 9.

No! Yes! Only if reactions occur! Yes! Ideal Gas, change in temperature or pressure. Survey Results. Class 15. Is the following possible?

Number of cases Number of factors Number of covariates Number of levels of factor i. Value of the dependent variable for case k

Solution Thermodynamics

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

Appendix II Summary of Important Equations

CHEMICAL REACTIONS AND DIFFUSION

(1) The saturation vapor pressure as a function of temperature, often given by the Antoine equation:

INTRODUCTION TO CHEMICAL PROCESS SIMULATORS

Lecture 8. Chapter 7. - Thermodynamic Web - Departure Functions - Review Equations of state (chapter 4, briefly)

Design Equations. ν ij r i V R. ν ij r i. Q n components. = Q f c jf Qc j + Continuous Stirred Tank Reactor (steady-state and constant phase)

Vapor-Liquid Equilibria for Water+Hydrochloric Acid+Magnesium Chloride and Water+Hydrochloric Acid+Calcium Chloride Systems at Atmospheric Pressure

Experimental and Modeling Studies for a Reactive Batch Distillation Column

EN40: Dynamics and Vibrations. Homework 4: Work, Energy and Linear Momentum Due Friday March 1 st

A Modulated Hydrothermal (MHT) Approach for the Facile. Synthesis of UiO-66-Type MOFs

V T for n & P = constant

9.2 Seismic Loads Using ASCE Standard 7-93

The influence of non-ideal vapor-liquid-equilibrium on vaporization of multicomponent hydrocarbon fuels

Mass Transfer Processes

and Statistical Mechanics Material Properties

DETERMINATION OF CO 2 MINIMUM MISCIBILITY PRESSURE USING SOLUBILITY PARAMETER

Summary Correlations from Feedstocks & Products

is the calculated value of the dependent variable at point i. The best parameters have values that minimize the squares of the errors

Non-Ideality Through Fugacity and Activity

Problem Points Score Total 100

Exercises of Fundamentals of Chemical Processes

Thermodynamics II. Department of Chemical Engineering. Prof. Kim, Jong Hak

Prediction of the flash point of ternary ideal mixtures

Chemical Engineering Department University of Washington

Review of Classical Thermodynamics

Thermodynamics General

Prediction of steady state input multiplicities for the reactive flash separation using reactioninvariant composition variables

ORIGIN 1. PTC_CE_BSD_3.2_us_mp.mcdx. Mathcad Enabled Content 2011 Knovel Corp.

Computation of Phase Equilibrium and Phase Envelopes

Lab 4: Two-level Random Intercept Model

Assignment 4. Adsorption Isotherms

Thermodynamics. Section 4

Amplification and Relaxation of Electron Spin Polarization in Semiconductor Devices

Chapter 13: Multiple Regression

Physics 4C. Chapter 19: Conceptual Questions: 6, 8, 10 Problems: 3, 13, 24, 31, 35, 48, 53, 63, 65, 78, 87

Interval Regression with Sample Selection

Determination of Structure and Formation Conditions of Gas Hydrate by Using TPD Method and Flash Calculations

Chapter 18, Part 1. Fundamentals of Atmospheric Modeling

Games and Market Imperfections

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Title: Radiative transitions and spectral broadening

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Solution Thermodynamics

NAME and Section No. it is found that 0.6 mol of O

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

Influence Of Operating Conditions To The Effectiveness Of Extractive Distillation Columns

Parts Manual. EPIC II Critical Care Bed REF 2031

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

Department of Physical Chemistry, Faculty of Chemistry, University of Bucharest, Bd. Regina Elisabeta 4-12, , Bucharest, Romania

Optimization of the thermodynamic model of a solar driven Aqua - ammonia absorption refrigeration system

Three-Phase Distillation in Packed Towers: Short-Cut Modelling and Parameter Tuning

Second Order Analysis

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

Physical Chemistry I for Biochemists. Chem340. Lecture 16 (2/18/11)

Particle Deposition in AERMOD: Overview

THE IGNITION PARAMETER - A quantification of the probability of ignition

The International Association for the Properties of Water and Steam

Numerical Investigation of Electroosmotic Flow. in Convergent/Divergent Micronozzle

Module 3: The Whole-Process Perspective for Thermochemical Hydrogen

McCabe-Thiele Diagrams for Binary Distillation

Problem Set #6 solution, Chem 340, Fall 2013 Due Friday, Oct 11, 2013 Please show all work for credit

Chemical Reaction Engineering

Structure and Property Prediction of Sub- and Super-Critical Water

Lecture 12. Transport in Membranes (2)

Quantitative Genetic Models Least Squares Genetic Model. Hardy-Weinberg (1908) Principle. change of allele & genotype frequency over generations

Modelli Clamfim Equazione del Calore Lezione ottobre 2014

Linear Approximation with Regularization and Moving Least Squares

PROBABILITY PRIMER. Exercise Solutions

Transcription:

Smulaton of a steady state flash Descrpton: Statonary flash smulaton of an Ethanol(1) - Water(2) - mxture Wth followng assumptons: Apart from heater and mass flows, no energy s transferred across the system boundary. Lqud and vapor phase are each deally mxed. Lqud and vapor phase are n equlbrum. The vapor phase behaves as an deal gas and the lqud as an ncompressble flud. VLE descrbed by extended Raoult s law. Vapor pressure calculaton by Antone. Evaporaton enthalpy calculaton by PPDS12. Real flud behavor of the lqud phase s expressed through Wlson s ge model. Excess varables neglected. Wlson and Antone parameter are taken from: Gmehlng, Kolbe; PPDS12 parameter are taken from: www.ddbst.de; Equaton System: 86383: statflash.moseqs IndexSpecfcaton: e[0]86383.nc = 2 Varable Specfcaton: 86386: IV statflash.mosvar Parameter Specfcaton: none. Results Specfcaton: 86388: RE statflash CBzzMath Org.mosvar Herarchcal vew of equatons: Equaton System: 86383: statflash.moseqs Descrpton: Statonary flash wth followng addtonal assumpton to bascflash: 1. Apart from heater and mass flows, no energy s transferred across the system boundary. 2. flash geometry s cylndrc. 1

Connected Equatons: Eq: 86378: StatEnergyBalance.mosequ (usng Nota: 8639: not - Desc.: Energy balance (statonary) 0 = F F,n h F,n h B,n F B,n h D,n F D,n + Q (1) Eq: 86379: StatComponentBalance.mosequ (usng Nota: 8639: not Desc.: Component balance (statonary) 0 = x F F F,n x B F B,n y D F D,n (2) Connected EQ-Systems: 86375: bascflash.moseqs 8638: tankgeometry.moseqs Connecton Level (1) EQ-Systems connected to 86383: statflash.moseqs: Equaton System: 86375: bascflash.moseqs Descrpton: basc model of a flash wth followng assumptons: 1. lqud and vapor phase are each deally mxed. Lqud and vapor phase are n equlbrum. 2. The vapor phase behaves as an deal gas and the lqud as an ncompressble flud. 3. VLE descrpton by extended Raoult s law. Vapor pressure calculaton by Antone 5. Actvty calculatons n the lqud phase by Wlson s g E -Model Connected Equatons: 2

Eq: 86350: sumrelatonvapor.mosequ (usng Nota: 8639: not - Desc.: Molar summaton relaton (vapor) NC 1 = =1 y D (3) Eq: 86353: molarenthalpyvapor.mosequ (usng Nota: 8639: not - Desc.: Molar enthalpy (vapor) NC h D,n = =1 y D h D,n + h E,D,n () Eq: 86356: molarenthalpyfeed.mosequ (usng Nota: 8639: not - Desc.: Molar enthalpy (feed) NC h F,n = =1 x F h F,n + h E,F,n (5) Eq: 8636: molarevapenthalpylqud.mosequ (usng Nota: 8639: not Desc.: Molar evaporaton enthalpy (lqud, pure component), source: www.ddbonlne.ddbst.de h B,LV,n = R T c (A P P DS12, (1 T T c ) 1 3 + BP P DS12, (1 T T c + C P P DS12, (1 T T c ) + D P P DS12, (1 T T c ) 2 3 ) 2 + E P P DS12, (1 T T c ) 6 ) Eq: 86368: molarenthalpycomponentvapor.mosequ (usng Nota: 8639: not Desc.: Molar enthalpy (vapor, pure component), assumptons: deal gas (6) h D,n = h o,n + C cp, 3 + A cp, (T T o ) + B cp, 2 ((T ) 3 (T o ) 3 ) + D cp, ((T ) 2 (T o ) 2 ) ((T ) (T o ) ) (7) 3

Eq: 86367: molarenthalpycomponentfeed.mosequ (usng Nota: 8639: not Desc.: Molar enthalpy (feed, pure component), assumptons: ncompressble flud deal gas behavor of correspondng vapor phase (pressure ndependency of h n V-phase) h F,n = h o,n + A cp, (T F T o ) + B cp, ((T F ) 2 (T o ) 2 ) 2 + C cp, ((T F ) 3 (T o ) 3 ) + D cp, 3 ((T F ) (T o ) ) + v L,n (p F p F,LV ) h F,LV,n (8) Eq: 8637: actvtycoeffwlsonparameterlqud.mosequ (usng Nota: 8639: not Desc.: Parameter calculaton (lqud) based on Wlson s g E -model (lqud), source: Gmehlng, Kolbe S.239 α B = NC =1 vl,n v L,n v L,n exp( λ T ) (9) Eq: 86372: volumeflowlqud.mosequ (usng Nota: 8639: not - Desc.: Volume flow (lqud) F B,n = F B,v ρ B,n (10) Eq: 86363: pressuredrop.mosequ (usng Nota: 8639: not - Desc.: Pressure drop (Feed > flash) p = p F p (11) Eq: 86361: volumelqud.mosequ (usng Nota: 8639: not - Desc.: Volume n flash (lqud) NC HU L,v = ( v L,n x B + v L,n,E ) HU L,n 1000 =1 (12)

Eq: 86351: sumrelatonlqud.mosequ (usng Nota: 8639: not - Desc.: Molar summaton relaton (lqud) NC 1 = =1 x B (13) Eq: 86352: molarenthalpylqud.mosequ (usng Nota: 8639: not - Desc.: Molar enthalpy (lqud) NC h B,n = =1 x B h B,n + h B,n,E (1) Eq: 86359: molarholdupcomponent.mosequ (usng Nota: 8639: not Desc.: Molar component holdup wthn flash HU n = x B HU L,n + y D HU V,n (15) Eq: 86357: levellqud.mosequ (usng Nota: 8639: not Desc.: Lqud level calculaton L L = HU L,v π (d) 2 (16) Eq: 86369: vaporpressureantonelqud.mosequ (usng Nota: 8639: not Desc.: Vapor pressure by Antone equaton (lqud, pure component), source: Kolbe, Mehlng p.59 p LV,B = (10) 3 (10) A Antone, B Antone, C Antone, +(T 273.15) (17) Eq: 86371: molardenstylqud.mosequ (usng Nota: 8639: not - Desc.: Molar densty wth no excess volume (lqud) NC 1 = ρ B,n =1 x B v L,n (18) 5

Eq: 86370: vaporpressureantonefeed.mosequ (usng Nota: 8639: not Desc.: Vapor pressure by Antone equaton (feed, pure component), source: Kolbe, Mehlng p.59 p LV,F = (10) 3 (10) A Antone, B Antone, C Antone, +(T F 273.15) (19) Eq: 86362: volume.mosequ (usng Nota: 8639: not Desc.: Volume of flash HU v = HU L,v + HU V,v (20) Eq: 86365: molarevapenthalpyfeed.mosequ (usng Nota: 8639: not Desc.: Molar evaporaton enthalpy (feed, pure component), source: www.ddbonlne.ddbst.de h F,LV,n = R T c (A P P DS12, (1 T F ) 1 3 + BP P DS12, (1 T F T c T c ) 2 3 + CP P DS12, (1 T F ) + D P P DS12, (1 T F T c T c ) 2 + E P P DS12, (1 T F ) 6 ) T c (21) U Eq: 86358: nternalenergy.mosequ (usng Nota: 8639: not - Desc.: Internal energy wthn flash (22) = HU L,n (h B,n p ( NC =1 xb v L,n + v L,n,E )) + HU V,n (h D,n R T z) (10) 6 Eq: 86360: volumevapor.mosequ (usng Nota: 8639: not - Desc.: Volume n flash (vapor) HU V,v = HU V,n R T z p 1000 (23) 6

Eq: 8635: sumrelatonfeed.mosequ (usng Nota: 8639: not - Desc.: Molar summaton relaton (feed) NC 1 = =1 x F (2) Eq: 86373: actvtycoeffwlsonlqud.mosequ (usng Nota: 8639: not Desc.: Actvty coeffcent calculatons by Wlson s g E -model (lqud) γ B = 1 x B + α B (1 x B ) α B exp((1 x B ) ( x B + α B (1 x B ) NC =1 αb α B ( NC =1 αb α B) xb + (1 x B ))) (25) Eq: 86355: VLEextendedRaoult.mosequ (usng Nota: 8639: not - Desc.: Extended Raoult s law for VLE n flash y D = x B pb,lv p γ B (26) Eq: 86366: molarenthalpycomponentlqud.mosequ (usng Nota: 8639: not Desc.: Molar enthalpy (lqud, pure component), assumptons: ncompressble flud deal gas behavor of correspondng vapor phase (pressure ndependency of h n V-phase) h B,n = h o,n + A cp, (T T o ) + B cp, ((T ) 2 (T o ) 2 ) 2 + C cp, ((T ) 3 (T o ) 3 ) + D cp, 3 ((T ) (T o ) ) + v L,n (p p B,LV ) h B,LV,n (27) 7

Equaton System: 8638: tankgeometry.moseqs Descrpton: Tank geometry Connected Equatons: Eq: 86381: volumetank.mosequ (usng Nota: 8639: not - Desc.: Flash s volume HU v = A L (28) Eq: 86382: ratodameterlengthtank.mosequ (usng Nota: 8639: not Desc.: dameter to heght rato of tank r D,L = d L (29) Eq: 86380: crosssectonarea.mosequ (usng Nota: 8639: not - Desc.: Flash s cross secton area A = π (d)2 (30) Equaton nstances: Eq: 86378: StatEnergyBalance.mosequ (usng Nota: 8639: not. Descrpton: Energy balance (statonary). 0 = e0.f F,n e0.h F,n e0.h B,n e0.f B,n e0.h D,n e0.f D,n + e0.q (31) Eq: 86379: StatComponentBalance.mosequ (usng Nota: 8639: not. Descrpton: Component balance (statonary). 0 = e0.x F =1 e0.f F,n e0.x B =1 e0.f B,n e0.y D =1 e0.f D,n (32) 0 = e0.x F =2 e0.f F,n e0.x B =2 e0.f B,n e0.y D =2 e0.f D,n (33) Eq: 86350: sumrelatonvapor.mosequ (usng Nota: 8639: not. Descrpton: Molar summaton relaton (vapor). 1 = (e0.y D =1 + e0.y D =2) (3) 8

Eq: 86351: sumrelatonlqud.mosequ (usng Nota: 8639: not. Descrpton: Molar summaton relaton (lqud). 1 = (e0.x B =1 + e0.x B =2) (35) Eq: 8635: sumrelatonfeed.mosequ (usng Nota: 8639: not. Descrpton: Molar summaton relaton (feed). 1 = (e0.x F =1 + e0.x F =2) (36) Eq: 86355: VLEextendedRaoult.mosequ (usng Nota: 8639: not. Descrpton: Extended Raoult s law for VLE n flash. e0.y D =1 = e0.x B =1 e0.pb,lv =1 e0.γ B =1 e0.p (37) e0.y D =2 = e0.x B =2 e0.pb,lv =2 e0.γ B =2 e0.p (38) Eq: 86352: molarenthalpylqud.mosequ (usng Nota: 8639: not. Descrpton: Molar enthalpy (lqud). e0.h B,n = (e0.x B =1 e0.h B,n =1 + e0.xb =2 e0.h B,n =2 ) + e0.hb,e,n (39) Eq: 86353: molarenthalpyvapor.mosequ (usng Nota: 8639: not. Descrpton: Molar enthalpy (vapor). e0.h D,n = (e0.y D =1 e0.h D,n =1 + e0.yd =2 e0.h D,n =2 ) + e0.hd,e,n (0) Eq: 86356: molarenthalpyfeed.mosequ (usng Nota: 8639: not. Descrpton: Molar enthalpy (feed). e0.h F,n = (e0.x F =1 e0.h F,n =1 + e0.xf =2 e0.h F,n =2 ) + e0.he,f,n (1) Eq: 86357: levellqud.mosequ (usng Nota: 8639: not. Descrpton: Lqud level calculaton. e0.l L = e0.hu L,v e0.π (e0.d) (2) (2) Eq: 86358: nternalenergy.mosequ (usng Nota: 8639: not. Descrpton: Internal energy wthn flash. e0.u = e0.hu L,n (e0.h B,n e0.p ((e0.x B =1 e0.vl,n =1 + e0.xb =2 e0.vl,n =2 ) + e0.ve,l,n )) + e0.hu V,n (e0.h D,n (10) (6) Eq: 86359: molarholdupcomponent.mosequ (usng Nota: 8639: not. Descrpton: Molar component holdup wthn flash. (3) e0.hu n =1 = e0.x B =1 e0.hu L,n + e0.y D =1 e0.hu V,n () 9

e0.hu n =2 = e0.x B =2 e0.hu L,n + e0.y D =2 e0.hu V,n (5) Eq: 86360: volumevapor.mosequ (usng Nota: 8639: not. Descrpton: Volume n flash (vapor). e0.hu V,v = e0.hu V,n e0.r e0.t e0.z e0.p 1000 (6) Eq: 86361: volumelqud.mosequ (usng Nota: 8639: not. Descrpton: Volume n flash (lqud). e0.hu L,v = ((e0.v L,n =1 e0.xb =1 + e0.v L,n =2 e0.xb =2) + e0.v E,L,n ) e0.hu L,n 1000 (7) Eq: 86362: volume.mosequ (usng Nota: 8639: not. Descrpton: Volume of flash. e0.hu v = e0.hu L,v + e0.hu V,v (8) Eq: 86363: pressuredrop.mosequ (usng Nota: 8639: not. Descrpton: Pressure drop (Feed > flash). e0. p = e0.p F e0.p (9) Eq: 8636: molarevapenthalpylqud.mosequ (usng Nota: 8639: not. Descrpton: Molar evaporaton enthalpy (lqud, pure component), source: www.ddbonlne.ddbst.de. e0.h B,LV,n =1 = e0.r e0.t=1 c (e0.a P P DS12,=1 (1 e0.t e0.t=1 c ) ( 3) 1 + e0.b P P DS12,=1 (1 e0.t e0.t=1 c ) ( 2 3) + e0.cp P DS12,=1 (1 e0.t e0.t=1 c ) + e0.d P P DS12,=1 (1 e0.t e0.t=1 c ) (2) + e0.e P P DS12,=1 (1 e0.t e0.t=1 c ) (6) ) (50) e0.h B,LV,n =2 = e0.r e0.t=2 c (e0.a P P DS12,=2 (1 e0.t e0.t=2 c ) ( 3) 1 + e0.b P P DS12,=2 (1 e0.t e0.t=2 c ) ( 2 3) + e0.cp P DS12,=2 (1 e0.t e0.t=2 c ) + e0.d P P DS12,=2 (1 e0.t e0.t=2 c ) (2) + e0.e P P DS12,=2 (1 e0.t e0.t=2 c ) (6) ) (51) 10

Eq: 86365: molarevapenthalpyfeed.mosequ (usng Nota: 8639: not. Descrpton: Molar evaporaton enthalpy (feed, pure component), source: www.ddbonlne.ddbst.de. e0.h F,LV,n =1 = e0.r e0.t c =1 (e0.a P P DS12,=1 (1 e0.t F + e0.b P P DS12,=1 (1 e0.t F e0.t c =1 e0.t c =1 ) ( 1 3) ) ( 2 3) + e0.cp P DS12,=1 (1 e0.t F e0.t=1 c ) + e0.d P P DS12,=1 (1 e0.t F e0.t=1 c ) (2) + e0.e P P DS12,=1 (1 e0.t F ) (6) ) e0.t c =1 (52) e0.h F,LV,n =2 = e0.r e0.t c =2 (e0.a P P DS12,=2 (1 e0.t F + e0.b P P DS12,=2 (1 e0.t F e0.t c =2 e0.t c =2 ) ( 1 3) ) ( 2 3) + e0.cp P DS12,=2 (1 e0.t F e0.t=2 c ) + e0.d P P DS12,=2 (1 e0.t F e0.t=2 c ) (2) + e0.e P P DS12,=2 (1 e0.t F ) (6) ) e0.t c =2 (53) Eq: 86366: molarenthalpycomponentlqud.mosequ (usng Nota: 8639: not -. Descrpton: Molar enthalpy (lqud, pure component), assumptons:. ncompressble flud deal gas behavor of correspondng vapor phase (pressure ndependency of h n V-phase) e0.h B,n =1 = e0.hn,o =1 + e0.a cp,=1 (e0.t e0.t o =1) + e0.b cp,=1 2 ((e0.t ) (2) (e0.t o =1) (2) ) + e0.c cp,=1 3 + e0.d cp,=1 ((e0.t ) () (e0.t o =1) () ) + e0.v L,n =1 (e0.p e0.pb,lv =1 ) e0.h B,LV,n =1 ((e0.t ) (3) (e0.t o =1) (3) ) (5) 11

e0.h B,n =2 = e0.hn,o =2 + e0.a cp,=2 (e0.t e0.t o =2) + e0.b cp,=2 2 ((e0.t ) (2) (e0.t o =2) (2) ) + e0.c cp,=2 3 + e0.d cp,=2 ((e0.t ) () (e0.t o =2) () ) + e0.v L,n =2 (e0.p e0.pb,lv =2 ) e0.h B,LV,n =2 ((e0.t ) (3) (e0.t o =2) (3) ) (55) Eq: 86368: molarenthalpycomponentvapor.mosequ (usng Nota: 8639: not -. Descrpton: Molar enthalpy (vapor, pure component), assumptons:. deal gas e0.h D,n =1 = e0.hn,o =1 + e0.a cp,=1 (e0.t e0.t o =1) + e0.b cp,=1 2 ((e0.t ) (2) (e0.t o =1) (2) ) + e0.c cp,=1 3 + e0.d cp,=1 ((e0.t ) () (e0.t o =1) () ) e0.h D,n =2 = e0.hn,o =2 + e0.a cp,=2 (e0.t e0.t o =2) + e0.b cp,=2 2 ((e0.t ) (2) (e0.t o =2) (2) ) + e0.c cp,=2 3 + e0.d cp,=2 ((e0.t ) () (e0.t o =2) () ) ((e0.t ) (3) (e0.t o =1) (3) ) (56) ((e0.t ) (3) (e0.t o =2) (3) ) (57) Eq: 86367: molarenthalpycomponentfeed.mosequ (usng Nota: 8639: not -. Descrpton: Molar enthalpy (feed, pure component), assumptons:. ncompressble flud deal gas behavor of correspondng vapor phase (pressure ndependency of h n V-phase) e0.h F,n =1 = e0.hn,o =1 + e0.a cp,=1 (e0.t F e0.t o =1) + e0.b cp,=1 2 ((e0.t F ) (2) (e0.t o =1) (2) ) + e0.c cp,=1 3 + e0.d cp,=1 ((e0.t F ) () (e0.t o =1) () ) + e0.v L,n =1 (e0.pf e0.p F,LV =1 ) e0.hf,lv,n =1 ((e0.t F ) (3) (e0.t o =1) (3) ) (58) 12

e0.h F,n =2 = e0.hn,o =2 + e0.a cp,=2 (e0.t F e0.t o =2) + e0.b cp,=2 2 ((e0.t F ) (2) (e0.t o =2) (2) ) + e0.c cp,=2 3 + e0.d cp,=2 ((e0.t F ) () (e0.t o =2) () ) + e0.v L,n =2 (e0.pf e0.p F,LV =2 ) e0.hf,lv,n =2 ((e0.t F ) (3) (e0.t o =2) (3) ) (59) Eq: 86369: vaporpressureantonelqud.mosequ (usng Nota: 8639: not -. Descrpton: Vapor pressure by Antone equaton (lqud, pure component), source: Kolbe, Mehlng p.59. ( ) e0.p B,LV =1 = (10) (3) (10) e0.p B,LV =2 = (10) (3) (10) ( e0.b e0.a Antone,=1 Antone,=1 e0.c Antone,=1 +(e0.t 273.15) e0.b e0.a Antone,=2 Antone,=2 e0.c Antone,=2 +(e0.t 273.15) Eq: 86370: vaporpressureantonefeed.mosequ (usng Nota: 8639: not. Descrpton: Vapor pressure by Antone equaton (feed, pure component), source: Kolbe, Mehlng p.59. ( ) e0.p F,LV =1 = (10)(3) (10) e0.p F,LV =2 = (10)(3) (10) ( e0.b e0.a Antone,=1 Antone,=1 e0.c Antone,=1 +(e0.t F 273.15) e0.b e0.a Antone,=2 Antone,=2 e0.c Antone,=2 +(e0.t F 273.15) ) ) (60) (61) (62) (63) Eq: 86371: molardenstylqud.mosequ (usng Nota: 8639: not. Descrpton: Molar densty wth no excess volume (lqud). 1 = e0.ρ B,n (e0.x B =1 e0.v L,n =1 + e0.xb =2 e0.v L,n =2 ) (6) Eq: 86372: volumeflowlqud.mosequ (usng Nota: 8639: not. Descrpton: Volume flow (lqud). e0.f B,n = e0.f B,v e0.ρ B,n (65) Eq: 86373: actvtycoeffwlsonlqud.mosequ (usng Nota: 8639: not. Descrpton: Actvty coeffcent calculatons by Wlson s g E -model (lqud). e0.γ B 1 =1 = e0.x B =1 + e0.αb =1 (1 e0.xb =1 ( ) exp (1 e0.x B e0.α=1 B =1) ( e0.x B =1 + e0.αb =1 (1 e0.xb =1 ) (e0.α=1 B + e0.αb =2 ) ) e0.αb =1 ((e0.α=1 B + e0.αb =2 ) e0.αb =1 ) e0.xb =1 + (1 e0.xb =1 )) (66) 13

e0.γ B 1 =2 = e0.x B =2 + e0.αb =2 (1 e0.xb =2 ( ) exp (1 e0.x B e0.α=2 B =2) ( e0.x B =2 + e0.αb =2 (1 e0.xb =2 ) (e0.α=1 B + e0.αb =2 ) ) e0.αb =2 ((e0.α=1 B + e0.αb =2 ) e0.αb =2 ) e0.xb =2 + (1 e0.xb =2 )) (67) Eq: 8637: actvtycoeffwlsonparameterlqud.mosequ (usng Nota: 8639: not. Descrpton: Parameter calculaton (lqud) based on Wlson s g E -model (lqud), source: Gmehlng, Kolbe S.239. e0.α B =1 = (e0.vl,n =1 + e0.vl,n =2 ) ( ) e0.vl,n =1 e0.λ=1 e0.v L,n exp e0.t =1 e0.α B =2 = (e0.vl,n =1 + e0.vl,n =2 ) ( ) e0.vl,n =2 e0.λ=2 e0.v L,n exp e0.t =2 (68) (69) Eq: 86380: crosssectonarea.mosequ (usng Nota: 8639: not. Descrpton: Flash s cross secton area. e0.a = e0.π (e0.d) (2) (70) Eq: 86381: volumetank.mosequ (usng Nota: 8639: not. Descrpton: Flash s volume. e0.hu v = e0.a e0.l (71) Eq: 86382: ratodameterlengthtank.mosequ (usng Nota: 8639: not. Descrpton: dameter to heght rato of tank. e0.r D,L = e0.d e0.l (72) Varable Specs 86386: IV statflash.mosvar Desgn varables e0.a Antone,=1 = 7.2371 e0.a Antone,=2 = 7.19621 e0.a P P DS12,=1 = 9.1919 e0.a P P DS12,=2 = 5.6297 e0.a cp,=1 = 9.008 e0.a cp,=2 = 32.22 e0.b Antone,=1 = 1592.86 e0.b Antone,=2 = 1730.63 1

e0.b P P DS12,=1 = 2.8118 e0.b P P DS12,=2 = 13.962 e0.b cp,=1 = 0.2139 e0.b cp,=2 = 0.0019225 e0.c Antone,=1 = 226.18 e0.c Antone,=2 = 233.26 e0.c P P DS12,=1 = 8.6931 e0.c P P DS12,=2 = 11.673 e0.c cp,=1 = 8.386E 5 e0.c cp,=2 = 1.058E 5 e0.d P P DS12,=1 = 11.776 e0.d P P DS12,=2 = 2.178 e0.d cp,=1 = 1.3723E 9 e0.d cp,=2 = 3.59E 9 e0.e P P DS12,=1 = 31.75 e0.e P P DS12,=2 = 0.31666 e0.f F,n = 1.75 e0.l = 0.5 e0.l L = 0.25 e0.r = 8.31 e0.t = 353.15 e0.t F = 353.15 e0.t=1 c = 516.2 e0.t=1 o = 298.15 e0.t=2 c = 67.3 e0.t=2 o = 298.15 e0. p = 25000.0 e0.λ =1 = 95.68 e0.λ =2 = 506.7 e0.π = 3.1159265359 e0.d = 0.16 e0.h B,E,n = 0.0 e0.h D,E,n = 0.0 e0.h E,F,n = 0.0 e0.h n,o =1 = 23800.0 e0.h n,o =2 = 21820.0 e0.p = 75000.0 e0.v E,L,n = 0.0 e0.v L,n =1 = 5.869E 5 e0.v L,n =2 = 1.807E 5 e0.x F =1 = 0.15 e0.z = 1.0 Iteraton varables e0.a = 1.0 15

e0.f B,n = 1.0 e0.f B,v = 1.0 e0.f D,n = 1.0 e0.hu L,n = 1.0 e0.hu L,v = 1.0 e0.hu V,n = 1.0 e0.hu V,v = 1.0 e0.hu v = 1.0 e0.hu=1 n = 1.0 e0.hu=2 n = 1.0 e0.q = 1.0 e0.u = 1.0 e0.α=1 B = 1.0 e0.α=2 B = 1.0 e0.γ=1 B = 1.0 e0.γ=2 B = 1.0 e0.ρ B,n = 1.0 e0.h B,n = 1.0 e0.h D,n = 1.0 e0.h F,n = 1.0 e0.h B,LV,n =1 = 1.0 e0.h B,n =1 = 1.0 e0.h D,n =1 = 1.0 e0.h F,LV,n =1 = 1.0 e0.h F,n =1 = 1.0 e0.h B,LV,n =2 = 1.0 e0.h B,n =2 = 1.0 e0.h D,n =2 = 1.0 e0.h F,LV,n =2 = 1.0 e0.h F,n =2 = 1.0 e0.p F = 1.0 e0.p B,LV =1 = 1.0 e0.p F,LV =1 = 1.0 e0.p B,LV =2 = 1.0 e0.p F,LV =2 = 1.0 e0.r D,L = 1.0 e0.x B =1 = 1.0 e0.x B =2 = 1.0 e0.x F =2 = 1.0 e0.y=1 D = 1.0 e0.y=2 D = 1.0 16

Results 86388: RE statflash CBzzMath Org.mosvar Desgn varables e0.a Antone,=1 = 7.2371 e0.a Antone,=2 = 7.19621 e0.a P P DS12,=1 = 9.1919 e0.a P P DS12,=2 = 5.6297 e0.a cp,=1 = 9.008 e0.a cp,=2 = 32.22 e0.b Antone,=1 = 1592.86 e0.b Antone,=2 = 1730.63 e0.b P P DS12,=1 = 2.8118 e0.b P P DS12,=2 = 13.962 e0.b cp,=1 = 0.2139 e0.b cp,=2 = 0.0019225 e0.c Antone,=1 = 226.18 e0.c Antone,=2 = 233.26 e0.c P P DS12,=1 = 8.6931 e0.c P P DS12,=2 = 11.673 e0.c cp,=1 = 8.386E 5 e0.c cp,=2 = 1.058E 5 e0.d P P DS12,=1 = 11.776 e0.d P P DS12,=2 = 2.178 e0.d cp,=1 = 1.3723E 9 e0.d cp,=2 = 3.59E 9 e0.e P P DS12,=1 = 31.75 e0.e P P DS12,=2 = 0.31666 e0.f F,n = 1.75 e0.l = 0.5 e0.l L = 0.25 e0.r = 8.31 e0.t = 353.15 e0.t F = 353.15 e0.t=1 c = 516.2 e0.t=1 o = 298.15 e0.t=2 c = 67.3 e0.t=2 o = 298.15 e0. p = 25000.0 e0.λ =1 = 95.68 e0.λ =2 = 506.7 e0.π = 3.1159265359 e0.d = 0.16 e0.h B,E,n = 0.0 e0.h D,E,n = 0.0 e0.h E,F,n = 0.0 17

e0.h n,o =1 = 23800.0 e0.h n,o =2 = 21820.0 e0.p = 75000.0 e0.v E,L,n = 0.0 e0.v L,n =1 = 5.869E 5 e0.v L,n =2 = 1.807E 5 e0.x F =1 = 0.15 e0.z = 1.0 Iteraton varables e0.a = 0.020106192982976 e0.f B,n = 1.100135560068 e0.f B,v = 3.0600035E 5 e0.f D,n = 0.33998565399321 e0.hu L,n = 0.232789065957856 e0.hu L,v = 0.00502658257 e0.hu V,n = 1.283990729133E e0.hu V,v = 0.00502658257 e0.hu v = 0.0100530969188 e0.hu=1 n = 0.02021282023757 e0.hu=2 n = 0.2126761816263232 e0.q = 13780.5013337676 e0.u = 0.065335268665 e0.α=1 B = 0.238168850216206 e0.α=2 B = 0.773536035272 e0.γ=1 B = 3.2920391080781 e0.γ=2 B = 1.020862183793210 e0.ρ B,n = 6311.912910192 e0.h B,n = 280525.9122197867 e0.h D,n = 23627.8265326519 e0.h F,n = 279798.270272677 e0.h B,LV,n =1 = 39069.79076511685 e0.h B,n =1 = 270032.671230127 e0.h D,n =1 = 230961.03110951 e0.h F,LV,n =1 = 39069.79076511685 e0.h F,n =1 = 270030.999873013 e0.h B,LV,n =2 = 156.18630081196 e0.h B,n =2 = 281522.3579753125 e0.h D,n =2 = 239958.6728128122 e0.h F,LV,n =2 = 156.18630081196 e0.h F,n =2 = 281521.906225312 e0.p F = 100000.0 e0.p B,LV =1 = 10833.6198873657 18

e0.p F,LV e0.p B,LV e0.p F,LV =1 = 10833.6198873657 =2 = 7266.7526108215 =2 = 7266.7526108215 e0.r D,L = 0.32 e0.x B =1 = 0.0867236920119577 e0.x B =2 = 0.913276307988023 e0.x F =2 = 0.85 e0.y=1 D = 0.122322212912 e0.y=2 D = 0.5875756777870588 Notaton 8639: not flash.mosnot Base lne symbols A Area [m]; Parameter B Parameter C Parameter D Parameter E Parameter F Flow[mol/s] HU Holdup I Integral (by tme) K Equlbrum constant [1]; Control constant L Heght [m] M Molar mass [g/mol] Q Heat flow [W] R Ideal gas constant [J/mol/K] T Temperature [K]; Tme constant [s] U Internal Energy [J] F Flow dfference L Level dfference [m] Q Heat flow dfference [J/s] T Temperature dfference [K] p Pressure dfference [Pa] α Wlson Parameter [1] γ Actvty coeffcent [1] λ Interacton parameter (Wlson Modell) π Natural constant π [1] ρ Densty [mol/m] d Dameter [m] h Specfc enthalpy[j/mol]; p Pressure [Pa] r rato [1] t Tme [s] 19

v Specfc molar volume [mol/m] x Mol fracton (lqud) [1] y Mol fracton (vapor) [1] z Compressblty factor [1] Superscrpts B D E F L LV V c max n o v Bottom Dstllate Excess Feed Lqud; Heght Lqud/Vapor Vapor Crtcal Maxmum Molar Reference volumetrc Subscrpts Antone P P DS12 cp Antone equaton PPDS12 equaton (http://www.ddbst.de) molar heat capacty Indces 1..N C Komponentenndex 20

2024 © sciencedocbox.com Privacy Policy | Terms of Service | Feedback