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

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

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

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

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

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

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

Non-Ideality Through Fugacity and Activity

Supplementary Notes for Chapter 9 Mixture Thermodynamics

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

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.

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

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

Exercises of Fundamentals of Chemical Processes

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

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

4.2 Chemical Driving Force

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

Name: SID: Discussion Session:

Chapter 8 Solutions Engineering and Chemical Thermodynamics 2e

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

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

Outlet temperature of a WGS reactor (Stage I) for the conversion of CO, applied for the abatement of CO to a fixed value.

Be true to your work, your word, and your friend.

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

1 Matrix representations of canonical matrices

Solution of Linear System of Equations and Matrix Inversion Gauss Seidel Iteration Method

NAME and Section No.

Complex Numbers. x = B B 2 4AC 2A. or x = x = 2 ± 4 4 (1) (5) 2 (1)

ME 300 Exam 2 November 18, :30 p.m. to 7:30 p.m.

Linear Feature Engineering 11

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

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

Solution Thermodynamics

University of Washington Department of Chemistry Chemistry 453 Winter Quarter 2015

Homework Notes Week 7

McCabe-Thiele Diagrams for Binary Distillation

Chemical Engineering Department University of Washington

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

Computation of Phase Equilibrium and Phase Envelopes

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

LNG CARGO TRANSFER CALCULATION METHODS AND ROUNDING-OFFS

MAE140 - Linear Circuits - Winter 16 Midterm, February 5

Lecture. Polymer Thermodynamics 0331 L Chemical Potential

Appendix II Summary of Important Equations

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

Non-Commercial Use Only

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

Physics 4B. A positive value is obtained, so the current is counterclockwise around the circuit.

MAE140 - Linear Circuits - Fall 10 Midterm, October 28

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.

THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions

A New Thermodynamic Function for Phase-Splitting at Constant Temperature, Moles, and Volume

ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM

MAE140 - Linear Circuits - Fall 13 Midterm, October 31

Homework Assignment 3 Due in class, Thursday October 15

Mass Transfer Processes

(2mn, m 2 n 2, m 2 + n 2 )

( ) Phase equilibrium Some basic principles for phase calculations

Section 8.3 Polar Form of Complex Numbers

Electrochemical Equilibrium Electromotive Force

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

Difference Equations

SUMMARY OF STOICHIOMETRIC RELATIONS AND MEASURE OF REACTIONS' PROGRESS AND COMPOSITION FOR MULTIPLE REACTIONS

UNIT 3 EXPRESSIONS AND EQUATIONS Lesson 4: Fundamental Theorem of Algebra. Instruction. Guided Practice Example 1

Assignment 4. Adsorption Isotherms

Adsorption: A gas or gases from a mixture of gases or a liquid (or liquids) from a mixture of liquids is bound physically to the surface of a solid.

8.6 The Complex Number System

Complex Numbers, Signals, and Circuits

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

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

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

TEST 5 (phy 240) 2. Show that the volume coefficient of thermal expansion for an ideal gas at constant pressure is temperature dependent and given by

a. (All your answers should be in the letter!

The Dirac Equation for a One-electron atom. In this section we will derive the Dirac equation for a one-electron atom.

EPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski

Chapter 12. Ordinary Differential Equation Boundary Value (BV) Problems

INTRODUCTION TO CHEMICAL PROCESS SIMULATORS

Math1110 (Spring 2009) Prelim 3 - Solutions

Principles of Food and Bioprocess Engineering (FS 231) Solutions to Example Problems on Heat Transfer

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis

PHYS 705: Classical Mechanics. Calculus of Variations II

CHEM 112 Exam 3 Practice Test Solutions

Quantum Mechanics for Scientists and Engineers. David Miller

CHEM 112 Exam 3 Practice Test Solutions

One-sided finite-difference approximations suitable for use with Richardson extrapolation

G4023 Mid-Term Exam #1 Solutions

Analytical Chemistry Calibration Curve Handout

( ) 2 ( ) ( ) Problem Set 4 Suggested Solutions. Problem 1

APPENDIX 2 FITTING A STRAIGHT LINE TO OBSERVATIONS

Managing Capacity Through Reward Programs. on-line companion page. Byung-Do Kim Seoul National University College of Business Administration

Using Spectrophotometric Methods to Determine an Equilibrium Constant Prelab

763622S ADVANCED QUANTUM MECHANICS Solution Set 1 Spring c n a n. c n 2 = 1.

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS

Chapter 5 rd Law of Thermodynamics

2.3 Nilpotent endomorphisms

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

PROBLEM SET 7 GENERAL EQUILIBRIUM

n α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0

Chapter 2 - The Simple Linear Regression Model S =0. e i is a random error. S β2 β. This is a minimization problem. Solution is a calculus exercise.

Special Relativity and Riemannian Geometry. Department of Mathematical Sciences

Transcription:

SOLUTION TO HOMEWORK #7 #roblem 1 10.1-1 a. In order to solve ths problem, we need to know what happens at the bubble pont; at ths pont, the frst bubble s formed, so we can assume that all of the number of moles for the components n the gas phase are close to zero (yet there s an equlbrum). Therefore, we use the compostons gven n the problem statement and we assume they are the compostons of the lqud. (x=z, z beng the composton gven). x ET 0.05 x 0.10 x NB 0.40 x M 0.45 Now, we know, from Raoult s law that: y x vap and from the problem statement we know that the pressure s 5 bar. Therefore, we proceed to calculate the bubble pont temperature. The way the bubble pont temperature s done s by dong the followng. We also know that: y 1 and therefore, substtutng y x vap 1 Therefore: 1 x vap ET ET x vap x NB vap vap NB x M M Now, we do the followng: 1 x ET 10 A T B x 10 A T ET B x NB 10 A T B x M 10 A T NB B M

From ths equaton, we can fnd the bubble pont temperature. Ths equaton can be solved for the temperature, whch comes out to be 294K. You can use bsecton, guessng (better be clever guessng), or Goal Seek n excel. Now, we need to calculate the compostons. Ths s done by usng Raoult s law yeldng: y ET 0.4167 y 0.1730 y NB 0.1601 y M 0.2502 One way of checkng f the answer s correct s to use then somethng was done wrong. y 1. If ths equalty s not correct, b. The dew pont calculatons are smlar, but n ths case the equatons change slghtly. x y vap also, snce we are at the dew pont (the pont n whch the frst drop of lqud s formed), we can say that almost no lqud s present. x y vap 1 whch, when the correct terms are substtuted n, we obtan: 1 y ET 10 A T B ET y 10 A T B y NB 10 A T B NB y M 10 A T B M where we use the composton gven (z) as the composton of the vapor. Solvng for the dew pont temperature yelds that the temperature s 314K. Now, by usng Raoult s law, we get that the compostons of the lqud phase are: x x x x ET NB M 0.0039 0.0037 0.5215 0.4409

c. In order to solve the flash problem, we need the help of other equatons. Frst, we begn by calculatng the K-factors for each component. Ths s done by usng the followng equaton: vap vap y K x Where we set the actvty and fugacty coeffcents to one because (Raoult law assumpton). Then, by calculatng each K-factor, we obtan: K ET 10.185 K 2.238 K NB 0.546 K M 0.743 Now, n a flash, the feed stream s separated nto the lqud and vapor streams. Ths s shown below: Vapor Feed V Lqud Therefore, we need to know what are the compostons of the vapor and lqud streams, so these would be part of our equatons used for solvng ths problem. We have: Vapor Stream y 1 Lqud Stream x 1

Now, by dong a mass balance around the flash, for each component, we can obtan four other equatons that could help us n solvng ths problem. These are (z s represent the component s compostons of the feed): z ET F x ET L y ET V Smlarly for propane z F x L y V for n-butane z NB F x NB L y NB V and for 2-methyl propane z M F x M L y M V Now, we can use a bass of 1 mol for the feed (F=1), whch wll gve us: 1 L V Now, we can replace the V n all equatons: ET ET ET ET 1 1 1 z F z x L y L z F z x L y 1 L z F z x L y L NB NB NB NB z F z x L y L M M M M Now, we can solve ths complex problem by dong the followng: frst, guess L (the lqud flow). Once ths value s guessed, we can proceed to solve the compostons by usng the equatons of the K-factors and the mass balances. As stated before, from the K-factors we obtan: K y x so, for each component, we get: y ET 10.185x ET y 2.238x y NB 0.546x NB y M 0.743x M

Now, these relatonshps can be substtuted nto the mass balances as stated above as follows: z F x L 10.185 x 1 L ET ET ET z F x L 2.238 x 1 L z F x L 0.546 x 1 L NB NB NB z F x L 0.743 x 1 L M M M whch wll allow us to calculate the compostons of each stream. x z F /[ L 10.185 1 L ] ET ET x z F /[ L 2.238 1 L ] x z F /[ L 0.546 1 L ] NB NB x z F /[ L 0.743 1 L ] M M We are now left wth only two equatons, namely y 1 and x 1 Thus we equate them y x Whch can be rewrtten as: K x x 0 or K x K x K x K x ( x x x x ) 0 ET ET NB NB M M ET NB M or ( KET 1) xet ( K 1) x ( KNB 1) xnb ( KM 1) xm 0 Thus substtutng the x s:

(10.185 1) z F /[ L 10.185 1 L ] (2.238 1) z F /[ L 2.238 1 L ] ET (0.546 1) z F /[ L 0.546 1 L ] (0.7431) z F /[ L 0.743 1 L ] 0 NB By solvng ths equaton for L (usng F=1), we obtan that: M Therefore, we have solved the flash problem.

#roblem 2

#roblem 3

#roblem 4

#roblem 5 roblem 6: a) Obtan an expresson relatng the mnmum amount of work needed to separate a mxture nto ts pure components (at constant T and ) as a functon of the fugactes of the components n the mxtures and the fugactes of the pure components. b) Show that the expresson can be wrtten only n terms of temperature and the molar fractons, when deal mxture s assumed and the Lews and Randall rule s used. a) We start by wrtng (Chapter 4) G W rev T, rev Thus, f Gm s the gbbs free energy of the mxture of two components, G1 s the Gbbs free energy of the stream of pure component 1 and G2 s the Gbbs free energy of the stream of pure component 2, we wrte W G G G rev m 1 2 But G x G x G m m,1 1 m,2 2 However, o f m, Gm, G (, T ) RT ln o f o o and 1 1 2 2 G G (, T), G G (, T) Substtutng, we get: W G G G n G n G n G (, T ) n G (, T ) o o rev m 1 2 m,1 m,1 m,2 m,2 1 1 2 2

After realzng that nm,1 n1 and nm,2 n2, we get f m,1 f m,2 Wrev RT n1ln n2ln o o f 1 f 2 art b) We can wrte f, x f (Lews and Randall) and therefore: o m Wrev RT n m,1 ln x1 nm,2 ln x2 or Wrev / n RT x m,1 ln x1 xm,2 ln x2 Recognze ths? It s somehow smlar to the gbbs free energy of mxng of deal mxtures, rght?

#roblem 7

#roblem 8

roblem #9

#roblem 10