ChemE 109 – Numerical and Mathematical Methods in Chemical and Biological Engineering
Fall 2019
PROJECT
Steady-state and transient analysis of a diffusion-reaction process
Assigned: Wednesday, November 20. Due: Monday, December 16.
Process description and modeling
The objective of this project is to study the steady-state and transient characteristics of the following
reaction scheme:
k1
A + 2B −−→ D
k2
A −−→ P
k3
B −−→ U
k4
U −−→ F
which takes place in the reaction zone shown below: The reactants A and B enter the reaction zone through the left membrane. The concentrations of the species A and B outside of the film which is located in the left side of the reaction zone are kept at the fixed values CAF and CBF . In the right side of the reaction zone, all the species deposit with different deposition rates. The boundaries in the reaction zone are at r = 0 and r = L.
Under the assumptions that the molecular diffusion inside the reaction zone is one-dimensional and follows Fick’s law, and that the diffusion coefficients DA, DB, DU and DF, of the species A, B, U and F, respectively are constant and equal (i.e. DA = DB = DU = DF = D ̄), the following equations can be derived, describing the change in CA, CB, CU and CF with position and time:
∂CA ∂t
∂CB ∂t
=
=
D ̄ ∂2CA − k1CACB2 − k2CA ∂r2
D ̄ ∂2CB − 2k1CACB2 − k3CB ∂r2
∂CU
∂t ∂r2 3B 4U
∂CF = D ̄∂2CF +k4CU ∂t ∂r2
where:
CA: concentration of species A
=D ̄∂2CU+kC−kC (1)
1
Figure 1: Schematic of reaction zone.
2
CB: concentration of species B
CU : concentration of species U
CF : concentration of species F
k1: rate constant for the first reaction k2: rate constant for the second reaction k3: rate constant for the third reaction k4: rate constant for the fourth reaction D ̄: diffusion coefficient
r: distance
t: time
Defining the following dimensionless variables and parameters:
yA=CA, yB=CB, yU=CU, yF=CF,
CAF CAF CAF
CBF 2 D ̄ β=C , τ=k1CAFt, D=kC2 L2
AF
k1C2 the system of Eq.1 can be written as:
CAF
(2)
x=r, γ= k2
1 AF
k3 , ζ=
k4 k1C2
, δ=
AF AF AF
L
k1C2
∂yA = D∂2yA −yAyB2 −γyA
∂τ ∂x2
∂yB = D∂2yB −2yAyB2 −δyB
∂τ
∂x2
= D∂2yU +δyB −ζyU
(3)
∂yU ∂τ
∂x2
∂yF = D∂2yF +ζyU
∂τ ∂x2
atx=0: ∂yA(0,t)=−(1−yA), ∂yB(0,t)=−(β−yB),yU(0,t)=0,andyF(0,t)=0,
atx=1: ∂yA(1,t)=−εyA,∂yB(1,t)=−ηyB2,∂yU(1,t)=−θyU,and∂yF(1,t)=0, ∂x ∂x ∂x ∂x
Questions
The boundary conditions are:
∂x ∂x
1. (50 pts.) In the first part of the project, the objective is to analyze the steady-state characteristics of the system. Using the O(∆x2) centered finite difference approximation method compute and plot the steady-state profiles yA vs. x, yB vs. x, yU vs. x, and yF vs. x for each of the following three cases:
a. (8 pts.) δ = 0, ζ = 0, D = 0.1, β = 1.5,γ = 0.05, ε = 0.0, η = 0.0, θ = 0.0.
b. (10 pts.) δ = 0.05, ζ = 0.0, D = 0.1, β = 1.5, γ = 0.02, ε = 0.1, η = 0.05, θ = 0.1.
c. (12 pts.) δ = 0.05, ζ = 0.03, D = 0.1, β = 1.5, γ = 0.02, ε = 0.1, η = 0.05, θ = 0.1.
Make sure you use enough nodal points to accurately calculate the solutions. Include your
rationale in the discussion (5 pts.). Explain the criteria that you employed to judge the
3
accuracy of the computed results (5 pts.). Discuss and explain the nature of the profiles that you obtained for the three cases (10 pts.).
2. (40 pts.) In the second part of the project, the objective is to study the time-dependent behavior of the system. For the case (c) from part 1, compute and plot (separately) the evolution of the spatial profiles of yA, yB, yF and yU in time (25 pts.). Use as initial conditions yA(x,0) = yB(x,0) = yU (x, 0) = yF (x, 0) = 0. Include enough profiles (in the same plot) to show clearly the time evolution until a steady-state is reached. Assume that the steady-state is obtained when you reach 99% of the steady-state values found in the previous question. Explain the method that you used to integrate the system in time and criteria that you employed to judge the accuracy of the computed results (5 pts.). Discuss the relationship between the plots that you obtained and your results from the previous question (10 pts.).
Requirements and Report Format
You are expected to work independently on the project. You may use parts of programs given in class or programs you have written on your own, but you may not share pieces of code among fellow students.
The project should be clearly written (deductions, up to 5 points, will be made for messiness, poor organization and poor writing), according to the following format:
• Title page. Please include the title of the project and your name.
• Table of contents. Indicate page numbers for all report sections and figures.
• Introduction. Briefly describe the problem, the development of the dimensionless mathematical model, and the strategy and methods you used to obtain solutions. Specify the class programs that you used, as well as any major modifications made to the class programs. The introduction will be worth 10 points.
• Results and Discussion. Address the questions posed in the previous section. Point values are indicated for each item. Make sure to include all required figures and tables, with a brief explanation of what they represent. Interpret your results adequately.
• Appendices. Include listings of all programs you used to generate your results. Clearly indicate all sections of the programs which are significant modifications of the original class codes. Please also upload your code to courseweb.
If you have any questions regarding the requirements of the project, contact Professor Panagiotis D. Christofides or TAs.
Good Luck!!!
4
Reviews
There are no reviews yet.