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

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

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

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

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

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.

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

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

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

Solution Thermodynamics

Name ID # For relatively dilute aqueous solutions the molality and molarity are approximately equal.

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

CHEMISTRY Midterm #2 answer key October 25, 2005

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

Mass Transfer Processes

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

Chemical Engineering Department University of Washington

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

CHEM 112 Exam 3 Practice Test Solutions

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

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

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

36.1 Why is it important to be able to find roots to systems of equations? Up to this point, we have discussed how to find the solution to

Diffusion Mass Transfer

χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body

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)

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

12. The Hamilton-Jacobi Equation Michael Fowler

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

Exercises of Fundamentals of Chemical Processes

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS

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

Thermodynamics General

Gasometric Determination of NaHCO 3 in a Mixture

Fundamental loop-current method using virtual voltage sources technique for special cases

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Analytical Chemistry Calibration Curve Handout

Polynomial Regression Models

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

DEMO #8 - GAUSSIAN ELIMINATION USING MATHEMATICA. 1. Matrices in Mathematica

Lecture 5.8 Flux Vector Splitting

Chapter 9: Statistical Inference and the Relationship between Two Variables

Supplementary Notes for Chapter 9 Mixture Thermodynamics

Numerical Heat and Mass Transfer

Physics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1

Outline. Unit Eight Calculations with Entropy. The Second Law. Second Law Notes. Uses of Entropy. Entropy is a Property.

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

A quote of the week (or camel of the week): There is no expedience to which a man will not go to avoid the labor of thinking. Thomas A.

Homework Notes Week 7

One Dimensional Axial Deformations

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

Linear, affine, and convex sets and hulls In the sequel, unless otherwise specified, X will denote a real vector space.

6.01: Introduction to EECS I Lecture 7 March 15, 2011

THE SUMMATION NOTATION Ʃ

Chemical Reaction Engineering

Assignment 4. Adsorption Isotherms

PHYS 705: Classical Mechanics. Calculus of Variations II

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

CHAPTER 7 ENERGY BALANCES SYSTEM SYSTEM. * What is energy? * Forms of Energy. - Kinetic energy (KE) - Potential energy (PE) PE = mgz

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

ISQS 6348 Final Open notes, no books. Points out of 100 in parentheses. Y 1 ε 2

Section 3.6 Complex Zeros

j) = 1 (note sigma notation) ii. Continuous random variable (e.g. Normal distribution) 1. density function: f ( x) 0 and f ( x) dx = 1

Affine transformations and convexity

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

Mean Field / Variational Approximations

Unit 5: Quadratic Equations & Functions

THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions

CIS526: Machine Learning Lecture 3 (Sept 16, 2003) Linear Regression. Preparation help: Xiaoying Huang. x 1 θ 1 output... θ M x M

Lecture. Polymer Thermodynamics 0331 L Chemical Potential

), it produces a response (output function g (x)

LNG CARGO TRANSFER CALCULATION METHODS AND ROUNDING-OFFS

4.2 Chemical Driving Force

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Homework Chapter 21 Solutions!!

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.

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

Section 8.3 Polar Form of Complex Numbers

Salmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2

Physics 5153 Classical Mechanics. Principle of Virtual Work-1

MAE140 - Linear Circuits - Fall 13 Midterm, October 31

The Feynman path integral

G4023 Mid-Term Exam #1 Solutions

Difference Equations

CHEM 112 Exam 3 Practice Test Solutions

Experiment 1 Mass, volume and density

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

Chapter 13: Multiple Regression

MAE140 - Linear Circuits - Winter 16 Midterm, February 5

Physics 240: Worksheet 30 Name:

Answers Problem Set 2 Chem 314A Williamsen Spring 2000

...Thermodynamics. If Clausius Clapeyron fails. l T (v 2 v 1 ) = 0/0 Second order phase transition ( S, v = 0)

The Fundamental Theorem of Algebra. Objective To use the Fundamental Theorem of Algebra to solve polynomial equations with complex solutions

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

Lecture 3: Probability Distributions

Chemical Equilibrium. Chapter 6 Spontaneity of Reactive Mixtures (gases) Taking into account there are many types of work that a sysem can perform

6.3.4 Modified Euler s method of integration

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

AS-Level Maths: Statistics 1 for Edexcel

Chapter 8 Balances on Nonreactive Processes 8.1 Elements of Energy Balances Calculations 8.1a Reference States A Review

Transcription:

nhe roblem-solvng Strategy hapter 4 Transformaton rocess onceptual Model formulaton procedure Mathematcal Model The mathematcal model s an abstracton that represents the engneerng phenomena occurrng n the conceptual model. It contans four thngs a set of equatons, a lst of assumptons, a degrees-of-freedom analyss, and a lst of varable defntons. If standard nomenclature has been defned for all problems, then the varable descrptons are optonal. Mathematcal models are developed from frst prncples; that s, any aoms, laws, or abstractons assumed and regarded as representng the hghest possble degree of generalzaton. The frst prncples n ths ntroductory course are Newton's second law of moton, materal balances, phase equlbrum, and energy balance as appled to chemcal process unts. The systematc formulaton procedure s descrbed below. ormulaton rocedure 1. Decde what frst prncple(s) apply for the problem by analyzng the nds quanttes: When asked to fnd flow rates or amounts ( m,, V or m, n, V ) and/or compostons ( w,, y, z, or cˆ ), you must wrte materal balances. When no chemcal reactons occur or when the chemcal reactons are known, use ths general form to wrte mass or mole balances: d( amount) sys flow n flow out ± reacton = dt When chemcal reactons do occur, but the reactons are not known, use ths form to wrte atom balances, one for each chemcal element: d( atom) sys flow n flow out = dt When asked to fnd heat (Q) and/or work (W), you must wrte the energy balance. That s, apply the general form: d( energy) sys flow n flow out ± Q ± W s = dt. Determne the type of materal balances to wrte mass?, mole?, or atom? Use mass balances when most nformaton s gven n terms of mass. Use mole balances when most nformaton s gven n terms of moles. Use atom balances when no chemcal reactons are gven for a reactor; otherwse, use mole balances when the chemcal reactons are known. Go to Step 3 to wrte mass or mole balances; otherwse, go to Step 4 to wrte atom balances. 3. Wrte the total balance and all component balances, one for each chemcal substance, when no chemcal reactons occur or when the chemcal reactons are known. Wrte total flow rates (or amounts) as symbols, even when ther values are known. Wrte component flow rates (or amounts) as symbols, even when ther values are known. v07.07.0 007, Mchael E. Hanyak, Jr., All Rghts Reserved age 4-5

nhe roblem-solvng Strategy hapter 4 Wrte total and component balances n lnear form (each term s a number*varable 1 ). Wrte known compostons as numbers n the balances (e.g., 0.35 or 0.45 n f ). Wrte unknown component quanttes as a varable (e.g., N or n f, N) Lne up common terms n ths set of materal balance equatons. Go to Step 5. onsult hapter 5 n ths nhe manual to learn how to wrte the mass or mole balances. 4. Wrte the total mass balance and the atom balances when chemcal reactons occur, but the chemcal reactons are not known. Wrte total flow rates (or amounts) as symbols, even when ther values are known. Wrte component flow rates (or amounts) as symbols, even when ther values are known. Wrte total mass and atom balances n lnear form (each term s a number*varable 1 ). Wrte known compostons as numbers n the balances (e.g., 0.35 or 0.45 n f ). Wrte unknown component quanttes as a varable (e.g., N or n f, N) By the conservaton of mass, total mass s conserved when chemcal reactons occur; however, component masses, component moles, and total moles are not necessarly conserved. onsult hapter 5 n ths nhe manual to learn how to wrte the total mass and atom balances. 5. Update the Lst of Assumptons based on what terms have been dropped from the materal balances. or eample, add no reactons when t s applcable. Note that the absence of a flow n or flow out term s already mpled n the assumpton for the process type; that s, batch, sem-batch, and sem-contnuous process. Also, the assumpton for a process operaton of steady state already mples that the term to the rght of the equal sgn n any balance s zero. 6. Wrte mture equatons, one for each mture that has a least one unknown composton. or a flowng mture (stream), wrte total flow equals the sum of the component flows. or a statonary mture, wrte total amount equals the sum of the component amounts. Wrte known compostons as numbers n mture equatons (e.g., 0.35 or 0.45 n f ). If all mtures n the conceptual model have known compostons, then no mture equatons are wrtten. See hapter 5 n the nhe manual to wrte the mture and composton equatons for a contnuous, batch, sem-batch, or sem-contnuous process. 7. Select one wrtten lnear equaton as the check equaton, because t s not lnearly ndependent. A lnear equaton s dependent f: 1) t s a multple of one of the other equatons, or ) t can be form by algebracally combnng one or more the other equatons together. or mass or mole balances, proceed as follows: All wrtten component balances and mture equatons can be algebracally combned to form the total balance equaton. Thus, one of the lnear equatons s not ndependent. The last wrtten equaton (usually a mture equaton) s to be chosen as the dependent or check equaton. or atom balances, proceed as follows: The total mass balance s chosen as the dependent or check equaton. The chosen dependent equaton s not part of the math model, but t s labeled wth the symbol or the word check whch s crcled. Ths equaton wll be used later as a comprehensve check for any human errors that mght occur when manually solvng the materal balances n the Numercal Soluton. v07.07.0 007, Mchael E. Hanyak, Jr., All Rghts Reserved age 4-6

nhe roblem-solvng Strategy hapter 4 8. Label the ndependent equatons, wrtten thus far n the mathematcal model, wth numbers begnnng wth one. Also crcle those numbers, n order to dentfy those equatons as beng part of the mathematcal model. or eample,,,,,, etc. 9. Mentally, do a degrees-of-freedom analyss on ths lnear set of ndependent equatons; that s, subtract the number of equatons from the number of varables (DO = # vars - # eqns). If the DO <= zero, then these equatons were not properly wrtten. If the DO >= one, then you have properly wrtten these equatons. The DO must equal the number of flow rates (or amounts) known n the conceptual model, before you can solve for all of the unknown varables n the lnear set. If no flow rates (or amounts) are known n the conceptual model, then you must assume one. That s, you pck a bass for a flow rate (or amount) on whch to do the calculatons later n the Numercal Soluton. Add ths bass to the Lst of Assumptons; for eample, = 100 mol / h or m 100 kg. = Symbols that represent constants lke R for the gas constant and M for molecular weght are not counted as varables n the mathematcal model. 10. If the DO s greater than the number of flow rates (or amounts) known n the conceptual model (mass, molar, and/or volumetrc quanttes), you must add more equatons to the mathematcal model untl ts DO becomes equal to number of flow rates (or amounts) known n the conceptual model. heck the conceptual model to see f any nformaton has not been used yet. Wrte an equaton for each pece of nformaton that has not be used. or eample, molar reacton converson, reacton yeld, ecess reactants, ecess ar, etc. These equatons should contan varables that appear n the materal balances. Thus, they usually do not ncrease the varable count but ncrease the equaton count n the DO analyss. Re-calculate the DO usng all equatons that now appear n the mathematcal model. If t does not equal the number of flow rates (or amounts) known n the conceptual model, then add more equatons, as descrbed below, untl t does. 11. Based on what s known (mass, molar, or volumetrc quanttes) and what varables are used n the materal balances (mass or molar quanttes), you can proceed one of two ways: If the known and balance quanttes are the same (.e., mass or moles), then suffcent equatons est n the lnear set to solve for all unknown flow rates (or amounts). Thus, proceed to Step 1 below. If the known and balance quanttes are not the same, then more equatons must be added to the mathematcal model to relate known quanttes to the unknown quanttes used n the materal balance. Wrte equatons to go from V m n or V m n, and then re-calculate the DO for all equatons that are now n the mathematcal model. onsult Table gotche n hapter 3 of the nhe manual for the nterrelatonshps between mass, molar, and volumetrc quanttes. 1. Eamne the nds quanttes and add more equatons to the mathematcal model, f necessary. Wrte equatons to relate the mass or molar quanttes n the materal balances to any mass, molar, or volumetrc flow rates (or amounts) n the nds lst. That s, wrte equatons to go from m V or n m V. Wrte composton equatons when mass fractons, mole fractons, volume fractons and/or concentratons are to be found. or eample, wrte for mass or mole fractons: v07.07.0 007, Mchael E. Hanyak, Jr., All Rghts Reserved age 4-7

nhe roblem-solvng Strategy hapter 4 m j= m w j or, j= j wrtten for each component j or m j= m w j or n j= n j wrtten for each component j onsult Table gotche n hapter 3 of the nhe manual for the nterrelatonshps between mass, molar, and volumetrc quanttes. 13. Show a degrees-of-freedom (DO) analyss for all equatons wrtten thus far n the mathematcal model. 14. heck f the DO agrees wth the number of known varables that appear n the conceptual model. Remember that symbols for constants lke R and M are not counted as varables. Numerous constants lke molecular weghts can be found n lterature tables. or eample, Table B.1 n Elementary rncples of hemcal rocesses by elder and Rousseau [005]. If t does not, then you may need to wrte more equatons. If the problem statement drectly or ndrectly ndcates that multple phases (gas, lqud, and/or sold) co-est n the process unt or set of process unts, you wll need to apply approprate frst prncples for phase equlbrum, such as functonal equatons lke Tdp = vlet[, V f = 1, Z ] for vapor-lqud equlbrum and dstrbuton coeffcent relatonshps for water solublty and envronmental parttonng coeffcents. onsult hapter 6 n the nhe manual to fnd the approprate relatonshps for phase equlbrum. After addng the approprate equatons, re-check your DO. 15. If applcable because of what s left n the nds lst (.e., a temperature, pressure, heat, and/or work) that has not been accounted for n the mathematcal model so far, wrte the energy balance wth ts mture enthalpy relatonshps. or eample, the energy to heat process Stream so that ts temperature s hgher n Stream s gven by: ˆ ˆ H H + Q Hˆ = hm[ T, X ] Hˆ = hm[ T,, X ] where Ĥ s the specfc enthalpy of the mture n the -th stream ( or ) n unts of energy/mass or energy/mole. Usng the above procedure, you should be able to develop a workable mathematcal model. The net three pages present three eamples for mathematcal models that were developed usng the above procedure. or these three eamples, hydrogen and oygen react to form water, and 50 mol% of the oygen gets converted (.e., the amount reacted over the amount fed equals 0.50). The frst eample s for a contnuous process where the chemcal reacton s known and mole balances are wrtten to count the molecules. In these mole balances, an etent-of-reacton term s wrtten to account for those molecules beng consumed or produced by chemcal reacton. or eample, term R I accounts for two hydrogen molecules beng consumed per reacton event. The second eample s also for a contnuous process, but the chemcal reacton s presumed not to be known and atom balances are wrtten to count the atoms. The thrd eample s a batch process where the chemcal reacton s known and mole balances are wrtten to count the molecules between the ntal and fnal tmes. Note that no materal flows nto or out of the system boundary for a batch process. v07.07.0 007, Mchael E. Hanyak, Jr., All Rghts Reserved age 4-8

nhe roblem-solvng Strategy hapter 4 ontnuous Reactor Eample { when the chemcal reactons are gven } Dagram Assumptons: T = 5 = 3 atm h = 100 mol / s, H,, N.50.10.40 reactor T 0 = 3 atm h = gas H, O, N, 1. contnuous process. steady state 3. M s deal gas Gvens: nds: Rn I: H + O H O V n L / h 50 mol% converson of O Mathematcal Model Total: 1 R I H : 0.50 H R I O : 0.10 1 R I N : 0.40 N H O: n, + R I # vars = 7 M : n = H + + N + # eqns = 5 0.10 onverson:. 50 0.10 deal gas: DO = # vars = 10 V = R T # eqns = 7 DO = 3 Snce three varables [, T, ] and a constant [R] are known, the math model can therefore be solved! The symbol R wthout a dot and subscrpt wll always represent the deal gas constant. The symbol R I wth a dot and/or a subscrpt wll always represent the etent of reacton. Eample net unts for the R g mol of H g rns I term n the H balance are: R I s. What does R I mean? g rns s If no reactons est n a contnuous problem, then the materal balances would not contan R terms. v07.07.0 007, Mchael E. Hanyak, Jr., All Rghts Reserved age 4-9

nhe roblem-solvng Strategy hapter 4 ontnuous Reactor Eample { when the chemcal reactons are not gven } Dagram Assumptons: T = 5 = 3 atm h = 100 mol / s, H,, N.50.10.40 reactor T 0 = 3 atm h = gas H, O, N, 1. contnuous process. steady state 3. M s deal gas Gvens: nds: 50 mol% converson of O V n L / h Mathematcal Model Total: m m H atom: (0.50 ) H n, O atom: (0.10 ) 1 n, N atom: (0.40 ) N # vars = 6 M : n = H + + N + # eqns = 4 0.10 onverson:. 50 0.10 deal gas: DO = # vars = 9 V = R T # eqns = 6 DO = 3 Snce three varables [, T, ] and a constant [R] are known, the math model can therefore be solved! In any atom balance, each term has net unts of number of atoms per tme. or eample, the number of H atoms n the feed stream s gven by the epresson: {note a mol and g-mol are synonymous} g atoms of H g mol of H g mol of (0.50 ) s g mol of H g mol of m s m When the flow rate s kg-mol/h or lb-mol/mn, the net unts are kg-atoms/h or lb-atoms/mn, respectvely. v07.07.0 007, Mchael E. Hanyak, Jr., All Rghts Reserved age 4-10

nhe roblem-solvng Strategy hapter 4 Batch Reactor Eample Dagram ntal tme, t closed contaner Gvens: T = 5 = 1 atm h n, H,, N.50.10.40 Rn I: H + O H O 50 mol% converson of O fnal tme, t f closed contaner nds: Tf 0 f = 1.5 atm hf = gas n f f, H f, =? f, N f, V f n L Assumptons: 1. batch process. unsteady state 3. deal gas behavor 4. n = 100 g-mol Mathematcal Model Total: 1 RI = n f n H : RI = n f,h 0.50 n O : 1 RI = n f, 0.10 n N : 0 = n f, N 0.40 n H O: RI = n f, # vars = 7 M : n f = n f, H + n f, + n f, N + n f, # eqns = 5 0.10 n n f, onverson:. 50 0.10 n V deal gas: f f f f Snce three varables [ n, T, f f DO = v07.07.0 007, Mchael E. Hanyak, Jr., All Rghts Reserved age 4-11 # vars = 10 = n RT # eqns = 7 DO = 3 ] and a constant [R] are known, the math model can therefore be solved! The symbol R wthout a dot and subscrpt wll always represent the deal gas constant. The symbol R I wth a subscrpt wll always represent the etent of reacton. Eample net unts for the RI term n g mol of H the H balance are: R I s ( g rns). What does RI mean, physcally? g rns If no reactons est n a batch problem, then the materal balances would not contan RI terms.