Example

Similar documents
GRAVITATION. Einstein Classes, Unit No. 102, 103, Vardhman Ring Road Plaza, Vikas Puri Extn., New Delhi -18 PG 1

Physics 107 HOMEWORK ASSIGNMENT #15

PH126 Exam I Solutions

Electric Potential and Gauss s Law, Configuration Energy Challenge Problem Solutions

MAE 210B. Homework Solution #6 Winter Quarter, U 2 =r U=r 2 << 1; ) r << U : (1) The boundary conditions written in polar coordinates,

Equilibria of a cylindrical plasma

15. SIMPLE MHD EQUILIBRIA

General Relativity Homework 5

Conservation of Linear Momentum using RTT

CHAPTER 25 ELECTRIC POTENTIAL


Hidden Two-Step Phase Transition and Competing

Solutions to Problems : Chapter 19 Problems appeared on the end of chapter 19 of the Textbook

Numerical solution of diffusion mass transfer model in adsorption systems. Prof. Nina Paula Gonçalves Salau, D.Sc.

Solute Transport In Biological Systems Design of An Artificial Kidney Utilizing Urease in Polymeric Beads. Dialysate. Flow out C(z) U Q,C b UB

( )( )( ) ( ) + ( ) ( ) ( )

Youn-Woo Lee School of Chemical and Biological Engineering Seoul National University , 599 Gwanangro, Gwanak-gu, Seoul, Korea

2. Radiation Field Basics I. Specific Intensity

CHAPTER 2 DERIVATION OF STATE EQUATIONS AND PARAMETER DETERMINATION OF AN IPM MACHINE. 2.1 Derivation of Machine Equations

CSTR - PFR - PBR

A Crash Course in (2 2) Matrices

2.25 Advanced Fluid Mechanics

Lecture 2 - Thermodynamics Overview

One-Dimensional, Steady-State. State Conduction with Thermal Energy Generation

Current, Resistance and

CBE Transport Phenomena I Final Exam. December 19, 2013

CHAPTER 10 ELECTRIC POTENTIAL AND CAPACITANCE

EKT 356 MICROWAVE COMMUNICATIONS CHAPTER 2: PLANAR TRANSMISSION LINES

A 1. EN2210: Continuum Mechanics. Homework 7: Fluid Mechanics Solutions

N igerian Journal of M athematics and Applications V olume 24, (2015),

Quantum Mechanics I - Session 5

If there are multiple rxns, use concentrations not conversions. These might occur in combination or by themselves.

Rigid Body Dynamics 2. CSE169: Computer Animation Instructor: Steve Rotenberg UCSD, Winter 2018

EKT 345 MICROWAVE ENGINEERING CHAPTER 2: PLANAR TRANSMISSION LINES

Physics 122, Fall December 2012

SPH4UI 28/02/2011. Total energy = K + U is constant! Electric Potential Mr. Burns. GMm

of Technology: MIT OpenCourseWare). (accessed MM DD, YYYY). License: Creative Commons Attribution- Noncommercial-Share Alike.

2 E. on each of these two surfaces. r r r r. Q E E ε. 2 2 Qencl encl right left 0

Chemical Reaction Engineering

Physics 107 TUTORIAL ASSIGNMENT #8

That is, the acceleration of the electron is larger than the acceleration of the proton by the same factor the electron is lighter than the proton.

Stellar Structure and Evolution

LECTURER: DR. MAZLAN ABDUL WAHID HEAT TRANSFER

Chemical Reaction Engineering

Physics 2212 GH Quiz #2 Solutions Spring 2016

Black Body Radiation and Radiometric Parameters:

ASTR415: Problem Set #6

Flux. Area Vector. Flux of Electric Field. Gauss s Law

Substances that are liquids or solids under ordinary conditions may also exist as gases. These are often referred to as vapors.

Entropy and Free Energy: Predicting the direction of spontaneous change The approach to Chemical equilibrium

= 4 3 π( m) 3 (5480 kg m 3 ) = kg.

University Physics (PHY 2326)

Applied Aerodynamics

Chapter 22 The Electric Field II: Continuous Charge Distributions

3.6 Applied Optimization

working pages for Paul Richards class notes; do not copy or circulate without permission from PGR 2004/11/3 10:50

Physics 2A Chapter 10 - Moment of Inertia Fall 2018

Chapter 1: Mole Balances

Easy. r p 2 f : r p 2i. r p 1i. r p 1 f. m blood g kg. P8.2 (a) The momentum is p = mv, so v = p/m and the kinetic energy is

4. Some Applications of first order linear differential

INTERACTION OF HYDRODYNAMIC ENVIRONMENT ON PERFORMANCE OF HOMOGENEOUS BIOREACTORS WITH ENZYME KINETIC MODELS

556: MATHEMATICAL STATISTICS I

Synthesis of Epichlorohydrin from Glycerol. Hydrochlorination of Glycerol

Physics 4A Chapter 8: Dynamics II Motion in a Plane

Chapter 2: Basic Physics and Math Supplements

Section 5: Magnetostatics

Review: Electrostatics and Magnetostatics

Chapter 13 Gravitation

ABSTRACT SIMULATION OF DYNAMIC PRESSURE- Professor Timothy A. Barbari Department of Chemical Engineering

In the previous section we considered problems where the

Chapter 28: Magnetic Field and Magnetic Force. Chapter 28: Magnetic Field and Magnetic Force. Chapter 28: Magnetic fields. Chapter 28: Magnetic fields

Module Summary Sheets. Mechanics 3 (Version B: Reference to new book) Topic 4: Volumes of revolution and centres of mass by integration

ELECTROSTATICS::BHSEC MCQ 1. A. B. C. D.

This full text version, available on TeesRep, is the post-print (final version prior to publication) of:

Much that has already been said about changes of variable relates to transformations between different coordinate systems.

Astrophysical Fluid Dynamics Solution Set 5 By Eric Bellm, Jeff Silverman, and Eugene Chiang

Chapter 16 Electrochemical Processes

5.4 Second Law of Thermodynamics Irreversible Flow 5

Practice. Understanding Concepts. Answers J 2. (a) J (b) 2% m/s. Gravitation and Celestial Mechanics 287

PHYSICS NOTES GRAVITATION

Phys102 Second Major-182 Zero Version Monday, March 25, 2019 Page: 1

Physics 2B Chapter 22 Notes - Magnetic Field Spring 2018

Exact Solution for Electro- Thermo- Mechanical Behavior of Composite Cylinder Reinforced by BNNTs under Non- Axisymmetric Thermo- Mechanical Loads

Do not turn over until you are told to do so by the Invigilator.

Electrostatics (Electric Charges and Field) #2 2010

Physics Courseware Physics II Electric Field and Force

Supporting information

Homework 3 MAE 118C Problems 2, 5, 7, 10, 14, 15, 18, 23, 30, 31 from Chapter 5, Lamarsh & Baratta. The flux for a point source is:

The Millikan Experiment: Determining the Elementary Charge

THERMODYNAMICS OF SURFACES AND INTERFACES

, and the curve BC is symmetrical. Find also the horizontal force in x-direction on one side of the body. h C

Geometry of the homogeneous and isotropic spaces

Heat transfer has direction as well as magnitude. The rate of heat conduction

Chapter 2: Conversion and Reactor Sizing

4. Compare the electric force holding the electron in orbit ( r = 0.53

Modeling of trickle bed reactor for hydrotreating of vacuum gas oils: effect of kinetic type on reactor modeling

Chapter Sixteen: Electric Charge and Electric Fields

$ i. !((( dv vol. Physics 8.02 Quiz One Equations Fall q 1 q 2 r 2 C = 2 C! V 2 = Q 2 2C F = 4!" or. r ˆ = points from source q to observer

6.641 Electromagnetic Fields, Forces, and Motion Spring 2005

Universal Gravitation

Transcription:

Chapte.4 iffusion with Chemical eaction Example.4- ------------------------------------------------------------------------------ fluiize coal eacto opeates at 45 K an atm. The pocess will be limite by the iffusion of oxygen countecuent to the cabon ioxie, CO, fome at the paticle suface. ssume that the coal is pue soli cabon with a ensity of 80 g/m 3 an that the paticle is spheical with an initial iamete of.5 0-4 m. i (% O an 79% N ) exists seveal iametes away fom the sphee. The iffusivity of oxygen in the gas mixtue at 45 K is.3 0-4 m /s. If a quasi-steay state pocess is assume, calculate the time necessay to euce the iamete of the cabon paticle to 5.0 0-5 m. (ef. Funamentals of Momentum, Heat, an Mass Tansfe by Welty, Wics, an Wilson, 4 th Eition, 00, pg. 496.) Solution ---------------------------------------------------------------------------------------------- The eaction at the cabon suface is C(s) + O (g) CO (g) We have iffusion of oxygen () towa the suface an iffusion of cabon ioxie (B) away fom the suface. The mola flux of oxygen is given by N, c mix y + y (N, + N B, ) In this equation, is the aial istance fom the cente of the cabon paticle. Since N, N B,, we have N, c mix y The system is not at steay state, the mola flux is not inepenent of since the aea of mass tansfe 4π is not a constant. Using quasi steay state assumption, the mass (mole) tansfe ate, 4π N,, is assume to be inepenent of at any instant of time. W 4π N, 4π c mix y constant y, y,inf -3

t the suface of the coal paticle, the eaction ate is much faste than the iffusion ate to the suface so that the oxygen concentation can be consiee to be zeo: y, 0. Sepaating the vaiables an integating gives y, 4π c mix y 0 W W 4π c mix y, > W 4πc mix y, Since one mole of cabon will isappea fo each mole of oxygen consume at the suface W C W 4πc mix y, Maing a cabon balance gives ρc 4 3 π M C t 3 ρc M C 4π t 4πc mix y, Sepaating the vaiables an integating fom t 0 to t gives t t 0 ρc M c C y mix, f i ρc t M c C y mix, ( i f ) The total gas concentation can be obtaine fom the ieal gas law c P T (0.0806)(45) 0.006 mol/m3 Note: 0.0806 m 3 atm/mol o K The time necessay to euce the iamete of the cabon paticle fom.5 0-4 m to 5.0 0-5 m is then ( ) ( ) 5 5 (80) 7.5 0.5 0 t 4 ()(0.006)(.3 0 )(0.) 0.9 s -4

Example.4- ------------------------------------------------------------------------------ Pulveize coal pellets, which may be appoximate as cabon sphees of aius mm, ae bune in a pue oxygen atmosphee at 450 K an atm. Oxygen is tansfee to the paticle suface by iffusion, whee it is consume in the eaction C(s) + O (g) CO (g). The eaction ate is fist oe an of the fom ɺ " C O whee 0. m/s. This is the eaction ate pe unit suface aea of the cabon pellets. Neglecting change in, etemine the steay-state O mola consumption ate in mol/s. t 450 K, the binay iffusion coefficient fo O an CO is.7 0-4 m /s. (ef. Funamentals of Heat Tansfe by Incopea an ewitt.) Solution ---------------------------------------------------------------------------------------------- We have iffusion of oxygen () towa the suface an iffusion of cabon ioxie (B) away fom the suface. The mola flux of oxygen is given by N, c y + y (N, + N B, ) In this equation, is the aial istance fom the cente of the cabon paticle. Since N, N B,, we have N, c y The system is not at steay state, the mola flux is not inepenent of since the aea of mass tansfe 4π is not a constant. Using quasi steay state assumption, the mass (mole) tansfe ate, 4π N,, is assume to be inepenent of at any instant of time. W 4π N, 4π c y constant y, y,inf The oxygen concentation at the suface of the coal paticle, y,, will be etemine fom the eaction at the suface. The mole faction of oxygen at a location fa fom the pellet is. Sepaating the vaiables an integating gives y, 4π c y y, W W 4π c (y, y, ) > W 4πc ( y, ) -5

The mole of oxygen aive at the cabon suface is equal to the mole of oxygen consume by the chemical eaction W 4π ɺ " 4π C O 4π c y, 4πc ( y, ) 4π c y, ( y, ) y, > y, + " y,.7 0 4 4 3.7 0 + 0. 0.63 The total gas concentation can be obtaine fom the ieal gas law. (Note: 0.0806 m 3 atm/mol K) c P T (0.0806)(450) 0.008405 mol/m3 The steay-state O mola consumption ate is W 4πc ( y, ) 4π(0.008405)(.7 0-4 )( 0.63)(0-3 ) W 6.66 0-9 mol/s Example.4-3 ------------------------------------------------------------------------------ biofilm consists of living cells immobilize in a gelatinous matix. toxic oganic solute (species ) iffuses into the biofilm an is egae to hamless poucts by the cells within the biofilm. We want to teat 0. m 3 pe hou of wastewate containing 0. mole/m 3 of the toxic substance phenol using a system consisting of biofilms on otating is as shown below. Waste wate fee steam biofilm C 0 C (z) Inet soli suface Well-mixe contacto C 0 Teate waste wate biofilm z0 etemine the equie suface aea of the biofilm with mm thicness to euce the phenol concentation in the outlet steam to 0.0 mole/m 3. The ate of isappeaance of phenol (species ) within the biofilm is escibe by the following equation -6

c whee 0.09 s - The iffusivity of phenol in the biofilm at the pocess tempeatue of 5 o C is.0 0-0 m /s. Phenol is equally soluble in both wate an the biofilm. (ef. Funamentals of Momentum, Heat, an Mass Tansfe by Welty, Wics, an Wilson, 4 th Eition, 00, pg. 496.) Solution ---------------------------------------------------------------------------------------------- The ate of phenol pocesse by the biofilms, W, is etemine fom the mateial balance on the pocess unit W 0. m 3 /h(0. 0.0) mol/m 3 8.0 0-3 mol/h W is then equal to the ate of phenol iffuse into the biofilms an can be calculate fom c W S N,z S z 0 In this equation, S is the equie suface aea of the biofilm an N,z is the mola flux of phenol at the suface of the biofilm. The mola flux of (phenol) is given by N,z c y + y (N,z + N B,z ) Since the biofilm is stagnant (o noniffusing), N B,z 0. Solving fo N,z give N,z ( y ) c y The mole faction of phenol in the biofilm, y, is much less than one so that c can be consiee to be constant. Theefoe N,z c y c Biofilm Soli suface N,z Maing a mole balance aoun the contol volume S z gives z -7

S N,z z S N,z z+ z + S z 0 iviing the equation by S z an letting z 0 yiels N, z Substituting N,z c (E-) c into equation (E-) we obtain c c > c c (E-) The solution to the homogeneous equation (E-) has two foms ) c C exp ) c B sinh z z + C exp + B cosh z z The fist exponential fom () is moe convenient if the omain of z is infinite: 0 z while the secon fom using hypebolic functions () is moe convenient if the omain of z is finite: 0 z δ. The constants of integation C, C, B, an B ae to be etemine fom the two bounay conitions. We use the hypebolic functions as the solution to Eq. (E-). c B sinh z + B cosh z (E-3) t z 0, c c s c 0 B t z δ, Theefoe c 0 B cosh δ + B sinh δ sinh B B cosh δ δ sinh c 0 cosh δ δ Equation (E-3) becomes -8

sinh c c 0 cosh δ δ sinh z + c 0 cosh z cosh δcosh z sinh δsinh z c c 0 cosh δ Using the ientity cosh( B) cosh()cosh(b) sinh()sinh(b) we have c c 0 δ cosh δ cosh ( z) c z 0 c 0 sinh ( δ z) cosh δ z 0 c 0 tanh δ The mola flux of phenol at the biofilm suface is given by N,z c c0 δ z 0 tanh δ δ The imensionless paamete δ Fo this poblem we have epesents the atio of eaction ate to iffusion ate. 0.09 0.00 m s 0 m 0 s δ 9.49 This value inicates that the ate of eaction is vey api elative to the ate of iffusion. The flux of phenol into the biofilm is then -9

N,z 0 (0.0)( 0 ) 0.00 (9.49) tanh(9.49) 3.9 0-8 mol/(m s) The equie suface aea of the biofilm is finally S W N, z 8.0 0 3 8 (3.9 0 )(3600) 57.0 m Example.4-4. ---------------------------------------------------------------------------------- Consie a spheical oganism of aius within which espiation occus at a unifom volumetic ate of C. That is, oxygen (species ) consumption is govene by a fist-oe, homogeneous chemical eaction. (a) If a mola concentation of C () C,0 is maintaine at the suface of the oganism, obtain an expession fo the aial istibution of oxygen, C (), within the oganism. (b) Obtain an expession fo the ate of oxygen consumption within the oganism. (c) Consie an oganism of aius 0.0 mm an a iffusion coefficient fo oxygen tansfe of 0-8 m /s. If C,0 5 0-5 mol/m 3 an 0 s -, what is the mola concentation of O at the cente of the oganism? What is the ate of oxygen consumption by the oganism? Solution ------------------------------------------------------------------------------------------ (a) If a mola concentation of C () C,0 is maintaine at the suface of the oganism, obtain an expession fo the aial istibution of oxygen, C (), within the oganism. + Figue E- Illustation of a spheical shell 4π The one-imensional mola flux of is given by the equation " N C (E-) pplying a mole balance on the spheical shell shown in Figue E- yiels fo steay state 4π N " 4π " N + 4π 0 + -30

iviing the equation by the contol volume (4π ) an taing the limit as 0, we obtain ( " N ) + 0 (E-) Fo a fist oe eaction, C an substituting the mola flux fom equation (E-) into the above equation, we have C C 0 C C 0 (E-3) In this equation, an ae constants inepenent of. We want to tansfom this equation into the fom y α y 0 (E-4) Let α, we can tansfom equation (4.6-3) into the fom of equation (E-4) by the following algebaic manipulations C α C 0 C C + α C 0 C C + α C 0 Since becomes ( C ) C C + C + C C +, the above equation ( C ) α C 0 Let y C, the equation has the same fom as equation (E-4) with the solution y B sinh(α) + B cosh(α) o C B sinh(α) + B cosh(α), whee α -3

The two constants of integation B an B can be obtaine fom the bounay conitions t 0, C finite o C 0 t, C C 0 (a nown value) pplying the bounay at 0 yiels 0 B pplying the bounay at yiels C 0 C B sinh(α) B sinh( α) Theefoe the concentation pofile fo species within the oganism is C C 0 sinh( α) sinh( α) (E-5) α t the cente of the oganism, the concentation is given by C ( 0) C 0 sinh( α) (b) Obtain an expession fo the ate of oxygen consumption within the oganism. ate of oxygen consumption within the oganism. 4π ( C ) The oxygen concentation within the oganism is given by equation (E-5) C C 0 sinh( α) sinh( α) (E-5) C C sinh( α) 0 α sinh( α ) + cosh( α) C C sinh( α) 0 α cosh( α) sinh( α) C C 0 [ ] ( α ) coth( α) ) -3

/ Let φ α Thiele moulus fo a fist oe eaction. Ignoing the minus sign, the ate of oxygen consumption within the oganism is then ate of oxygen consumption 4π C 0 (φ cothφ - ) ate of oxygen consumption 4π C 0 (φ cothφ - ) (c) Consie an oganism of aius 0.0 mm an a iffusion coefficient fo oxygen tansfe of 0-8 m /s. If C,0 5 0-5 mol/m 3 an 0 s -, what is the mola concentation of O at the cente of the oganism? What is the ate of oxygen consumption by the oganism? α t the cente of the oganism, the concentation is given by C ( 0) C 0 sinh( α) α / 0 8 0 / 4.47 0 4 m α / 4 ( ) 0 0 8 0 / 4.47 α C ( 0) C 0 sinh( α) 5 0-5 4.47 sinh(4.47) 5. 0-6 mol/m 3 ate of oxygen consumption 4π C 0 (φ cothφ - ) ate 4π(0-4 )(0-8 )( 5 0-5 ) [4.47 coth(4.47) - ].8 0-5 mol/s The following Matlab pogam plots the concentation of oxygen within the oganism as a function of position. % Example.4-4 % alfa4.47e4; % m e-4; % m alfa4.47; C05e-5; % mol/m3 o(:50)/50; *o; -33

CC0*sinh(alfa*)./(o*sinh(alfa)); plot(o,c) gi on xlabel('/');ylabel('c_(mol/m^3)') Figue E.4-4 Oxygen concentation pofile in a spheical oganism. We now consie the iffusion of species into a spheical catalyst paticle whee homogeneous fist oe chemical eaction occus. The concentation pofile fo species within the spheical catalyst paticle is then C C sinh( α) sinh( α) (.4-) In this equation C is the concentation of species at the suface of the catalyst paticle an α is efine by the expession α, whee is the fist oe ate constant an is the iffusivity of in the paticle. t the cente of the bea, the concentation is given by α C ( 0) C sinh( α) (.4-) -34