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)

Similar documents
Ph.D. Qualifying Examination in Kinetics and Reactor Design

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

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

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

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

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

Open Systems: Chemical Potential and Partial Molar Quantities Chemical Potential

Exercises of Fundamentals of Chemical Processes

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

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

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

INFLUENCE OF AXIAL MIXING AND PARTICLE SIZE NON-UNIFORMITY ON PERFORMANCE OF CHROMATOGRAPHIC REACTORS

Mass Transfer Processes

#64. ΔS for Isothermal Mixing of Ideal Gases

Supplementary Notes for Chapter 9 Mixture Thermodynamics

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

modeling of equilibrium and dynamic multi-component adsorption in a two-layered fixed bed for purification of hydrogen from methane reforming products

Lecture 12. Transport in Membranes (2)

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

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

Basic concept of reactive flows. Basic concept of reactive flows Combustion Mixing and reaction in high viscous fluid Application of Chaos

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

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.

Mathematical Modeling of an Industrial Naphtha Reformer with Three Adiabatic Reactors in Series

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

Chemical Reaction Engineering

STAT 511 FINAL EXAM NAME Spring 2001

CHEMICAL ENGINEERING

MATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018

Influence Of Operating Conditions To The Effectiveness Of Extractive Distillation Columns

Gasometric Determination of NaHCO 3 in a Mixture

Lecture 16 Statistical Analysis in Biomaterials Research (Part II)

Example: (13320, 22140) =? Solution #1: The divisors of are 1, 2, 3, 4, 5, 6, 9, 10, 12, 15, 18, 20, 27, 30, 36, 41,

Diffusion Mass Transfer

CHEMISTRY Midterm #2 answer key October 25, 2005

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

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

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

Chemical Engineering Department University of Washington

PART I: MULTIPLE CHOICE (32 questions, each multiple choice question has a 2-point value, 64 points total).

Using Spectrophotometric Methods to Determine an Equilibrium Constant Prelab

Analytical Chemistry Calibration Curve Handout

Digital Signal Processing

ESA modelling and cycle design

On the Multicriteria Integer Network Flow Problem

Solution Thermodynamics

If the solution does not follow a logical thought process, it will be assumed in error.

Simultaneous BOP Selection and Controller Design for the FCC Process

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

The optimal delay of the second test is therefore approximately 210 hours earlier than =2.

Please review the following statement: I certify that I have not given unauthorized aid nor have I received aid in the completion of this exam.

Lecture. Polymer Thermodynamics 0331 L Chemical Potential

Math1110 (Spring 2009) Prelim 3 - Solutions

Module 3: The Whole-Process Perspective for Thermochemical Hydrogen

CHAPTER 13. Exercises. E13.1 The emitter current is given by the Shockley equation:

Introduction to Statistical Methods

2 Finite difference basics

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

Chapter 5 rd Law of Thermodynamics

CHEMICAL REACTIONS AND DIFFUSION

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

The equation of motion of a dynamical system is given by a set of differential equations. That is (1)

Chemical Reaction Engineering

Physics 607 Exam 1. ( ) = 1, Γ( z +1) = zγ( z) x n e x2 dx = 1. e x2

Chapter 11: Simple Linear Regression and Correlation

Negative Binomial Regression

Chapter 18, Part 1. Fundamentals of Atmospheric Modeling

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

November 5, 2002 SE 180: Earthquake Engineering SE 180. Final Project

Chemistry 163B Free Energy and Equilibrium E&R ( ch 6)

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

Amiri s Supply Chain Model. System Engineering b Department of Mathematics and Statistics c Odette School of Business

NUMERICAL DIFFERENTIATION

Note 10. Modeling and Simulation of Dynamic Systems

Number Average Molar Mass. Mass Average Molar Mass. Z-Average Molar Mass

Chemical Reaction Engineering. Dr. Yahia Alhamed

Formulas for the Determinant

APPENDIX A Some Linear Algebra

Thermodynamics General

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

An Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors

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.

ECE 534: Elements of Information Theory. Solutions to Midterm Exam (Spring 2006)

Name: SID: Discussion Session:

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

,..., k N. , k 2. ,..., k i. The derivative with respect to temperature T is calculated by using the chain rule: & ( (5) dj j dt = "J j. k i.

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

IMPROVING THE MATHEMATICAL MODEL OF A HETEROCATALYTIC TUBE REACTOR

Notes on Analytical Dynamics

BAE 820 Physical Principles of Environmental Systems

FUNDAMENTALS OF CHEMISTRY Vol. II - Microkinetics Versus Macrokinetics: The Role of Transport Phenomena in Determining Reaction Rates - Rota R.

Airflow and Contaminant Simulation with CONTAM

Simultaneous Optimization of Berth Allocation, Quay Crane Assignment and Quay Crane Scheduling Problems in Container Terminals

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

EEE 241: Linear Systems

Problem Set 9 Solutions

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

Review of Classical Thermodynamics

Chapter 7 Channel Capacity and Coding

Transcription:

Desgn Equatons Batch Reactor d(v R c j ) dt = ν j r V R n dt dt = UA(T a T) r H R V R ncomponents V R c j C pj j Plug Flow Reactor d(qc j ) dv = ν j r 2 dt dv = R U(T a T) n r H R Q n components j c j C pj Strred Tank Reactor d(v R c j ) dt = Q f c jf Qc j + ν j r V R dt dt = UA(T a T) n ncomponents r H R V R + Q f j c jf (H jf H j ) ncomponents V R j c j C pj Contnuous Strred Tank Reactor (steady-state and constant phase) 0 = Q f c jf Qc j + ν j r V R n n components Tf 0 = UA(T a T) r H R V R + Q f c jf C pj dt j T 2

General Informaton on Mass Transfer wth Reacton For heterogeneous reactons one must determne the rate per unt volume of pellet. R jp = 1 R j dv = S p dc j D j V p V p V p dr r=rp The dmensonless steady-state concentraton wthn a symmetrcal pellet s found wth [ 2 a 2 ] R js c + R = 0 R = R j c js D j R js For a frst-order reacton (A B) wth a sphercal pellet surface concentraton of c As ( ) R c A = c As r snh Φ( 3r R ) snh(3φ) R Ap = 1 Φ [ 1 tanh3φ 1 ] ( k 1 c As ) 3Φ Φ = R 3 k 1 For some heterogeneous cases (A B) t s approprate to use D e η = R j R js η 1 [ 1 Φ tanh3φ 1 ] 3Φ Φ = V [ p n + 1 S p 2 knc ] s n 1 0.5 D e where n s the reacton order n component A and the overall reacton order. 3

Useful Integrals 1 a + bx dx = 1 ln(a + bx) b (a + bx) n dx = (a + bx)n+1 (n + 1)b ;n 1 ( dx (a + bx)(a + b x) = 1 a ab a b ln + b ) x a + bx a + bx bx a + b dx = x b + ab a b b 2 ln ( a + b x ) (a + bx) m (a + b x) n dx = [ 1 (a + bx) m ] (a + bx) m 1 (n 1)b (a + b n 1 mb x) (a + b x) n 1dx x m (a + bx) n dx = xm+1 (a + bx) n m + n + 1 + an m + n + 1 x m (a + bx) n 1 dx where: Smpson s three-eghth s rule for numercal ntegraton: X3 X 0 f(x)dx = 3h 8 [f(x 0) + 3f(X 1 ) + 3f(X 2 ) + f(x 3 )] h = X 3 X 0 3 X 1 = X 0 + h X 2 = X 0 + 2h Integraton of N + 1 ponts, where N s even: XN ] [f 0 + 4f 1 + 2f 2 + 4f 3 + 2f 4 +... + 4f (N 1) + f N where: X 0 f(x)dx = h 3 h = X N X 0 N 4

Problem 1 The elementary, lqud-phase, rreversble reacton A + B C s carred out n a CSTR operated at 300 K. Reactant B also somerzes to form undesred product D. The rates of reacton are gven by: A + B C B D r 1 = k 1 c A c B r 2 = k 2 c B There are two reactant streams that are mxed before enterng the reactor. The frst reactant stream contans A (but no B) at a concentraton of 1 mol/lter and a flow rate of 1 lter/mn. The second reactant stream contans B (but no A) at a concentraton of 1 mol/lter and a flow rate of 1 lter/mn. The two reactant streams mx to form a sngle feed stream that s equal molar n A and B, wth a total volumetrc flow rate of 2 lters/mn. Addtonal Informaton: At 300 K, k 1 = 1 lters/mol-mn, k 2 = 0.1 mn -1 Reactor volume V R = 100 lters a) What s the resdence tme n the CSTR? b) What fracton of ncomng A s converted n the CSTR? c) What fracton of ncomng B s converted n the CSTR? d) What are the concentratons of products C (desred) and D (undesred) n the output stream of the CSTR? Problem 2 Sulfurc acd s produced by the exothermc, reversble reacton of SO 2 and ar over a vanadum pentoxde catalyst at 800-900 o F. The reacton s exothermc and has no sde-reactons. The feed gas to the catalytc reactors s produced by burnng sulfur n ar to produce SO 2 then blendng n more ar. a) One early reactor desgn uses 3 or 4 stages of catalyst wth ntercoolng between the stages to obtan an overall converson of SO 2 to SO 3 of 95 to 98%. Ths s followed by a scrubber at lower T wth addton of water to condense H 2 SO 4 from the gas. The effluent gas contanng 1000-5000 ppm SO 2 was vented. 1) What tradeoffs determne the T of the reactor? 2) Why are multple reactor stages useful to ncrease converson? b) A later system desgn addresses the emssons of SO 2 by usng two scrubbers. Gas from the frst scrubber s reheated to 800 o F then passes over addtonal catalyst stages wth ntercoolng to get 99-99.9% converson to SO 3. The gas s then cooled and passed through another scrubber to remove the SO 3 as sulfurc acd. 1) Why does ths work? c) What would be requred to desgn the process wth a classc recycle of unreacted SO 2 so that 99.99% converson would be acheved? 5

Problem 3 A gas phase thermal decomposton reacton takes place n a plug-flow reactor. The reactant A s decomposed nto products B and C. The reacton has the followng rate constant: log 16010 ] = 2 T 1 10( k[ s ) 15.1 The system operates sothermally at temperature T = 1000Kand pressure P = 1 atm. The feed stream conssts of pure reactant A, and has molar flow rate N Af = 5 kmol/hr and volumetrc flow rate Q f. The feed stream s at the same temperature and pressure as the reactor. The gases follow deal gas behavor. The converson of A at the effluent s x A = 0.80. Please compute: a) Q f and Q, whch are the volumetrc flow rates at the reactor nlet and outlet, respectvely. b) The reactor volume (V R ) requred to attan a converson of A of x A = 0.80. R = 0.082 lter atm mol -1 K -1 Problem 4 The followng frst-order catalytc reacton s conducted n a fxed-bed reactor. The reactant and product are gases. A 2B r = kc A The ntrnsc rate constant s k = 2.35 10 13 exp( 19,124/T) s -1, where T s n Kelvn. The feed conssts of pure A at 2 atm and a total molar flow rate of N Af = 12 mol/s. At ths feed rate there are no external mass transfer lmtatons. The reactor s sothermal at a temperature of 578 K. The catalyst pellets are sphercal wth a radus of 0.3 cm. The effectve dffusvty s D eff = 0.004 cm 2 /s. The catalyst bed has a vod fracton of 0.412. Under these condtons, the fxed bed reactor was desgned to yeld a converson of A correspondng to x A = 0.9. There was a problem wth the purty of the feed for several days, the catalyst was posoned by an mpurty and now the catalyst s ten tmes less actve. For all practcal purposes the ntrnsc rate constant s now k posoned = 2.35 10 12 exp( 19,124/T) s -1. Assumng the dffusvty wll not change sgnfcantly wth temperature, what new reactor temperature wll you need to ensure the exstng reactor contnues to provde x A = 0.9? You should assume you can adjust the feed and reactor to ths new temperature so the system s stll sothermal. R = 82 cm 3 atm K 1 mol 1 6

Problem 5 The lqud-phase parallel reactons A B A C r 1 = k 1 c A r 2 = k 2 c A are to be carred out n a CSTR. The feed enters at 65 C and conssts of pure A at a concentraton of c Af = 5 mol/lter. Varable Value Unts ρq 93,000 g/mn H R1-3,200 cal/mn H R2 2,000 cal/mn V R 100 lter Q f 100 lter/mn T a 50 C C pf 0.22 cal/g-k U 0.65 cal/mn-cm 2 -K A 15,400 cm 2 k 1 3.16 10 14 exp(-12,500/t) mn -1 ; T n K k 2 2.52 10 9 exp(-8,500/t) mn -1 ; T n K Under these condtons the reactor has three possble operatng condtons (two stable and one unstable) as llustrated n Fgure 1 on Page 8. What process parameters would you change to ensure operaton at the upper branch of stable solutons? You do not need to have the same operatng temperature as llustrated n Fgure 1. You do need to realze hgh converson of reactant A. Explan how changng the parameter(s) you select wll work to make the system have a stable soluton. 7

8