CHEN 320 - Numerical Analysis for Chemical Engineers (Section 502) Fall 2016 COMPUTER PROJECT Steady-state and transient analysis of a diffusion-reaction process Assigned: Thursday, November 3. Due: Thursday, December 1. Process description and modeling The objective of this project is to study the steady-state and transient characteristics of the following reaction scheme: A + 2B D A P B U U B which takes place in the membrane reactor shown below: The reactants A and B enter the reaction Film Reaction Zone C AF C BF 0 L Figure 1: Schematic of membrane reactor. 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 and DU , of the species A, B and the following equations can be U, respectively are constant and equal (i.e. DA = DB = DU = D), 1 derived, describing the change in CA , CB , CU with position and time: 2 CA t CA k1 CA C 2 k2 CA = D B r 2 CB t CB 2k1 CA CB2 k3 CB + k4 CU = D r 2 CU t CU + k3 CB k4 CU = D r 2 2 (1) 2 where: CA : concentration of species A CB : concentration of species B CU : concentration of species U 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 diffusion coefficient D: r: distance t: time Defining the following dimensionless variables and parameters: CA CB CU , yB = , yU = , CAF CAF CAF CBF D 2 = , = k1 CAF t, D = 2 CAF k1 CAF L2 k2 k3 k4 r , = , = x= , = 2 2 2 L k1 CAF k1 CAF k1 CAF yA = (2) the system of Eq.1 can be written as: yA 2 yA = D 2 yA yB2 yA x yB = D 2 yB 2yA yB2 yB + yU x2 yU = D 2 yU + yB yU x2 (3) yB yA (0, t) = (1yA ), (0, t) = (yB ) and yU (0, t) = 0, The boundary conditions at x = 0 are x x yB yU yA (1, t) = yA2 , (1, t) = yB , and (1, t) = yU2 . while at x = 1, x x x Questions 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 and yU vs. x, for each of the following three cases: 2 (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 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 , and yU in time (25 pts.). Use as initial conditions yA (x, 0) = yB (x, 0) = yU (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. If you have any questions regarding the requirements of the project, contact Professor Joseph Kwon or TAs. Good Luck!!! 3