This preview shows pages 1–2. Sign up to view the full content.
This preview has intentionally blurred sections. Sign up to view the full version.
View Full Document
Unformatted text preview: ChE 101 2012 Problem Set 1 1. Computational problem of the week: Every week, we will assign one problem designed to increase your exposure to numerical methods, improve your algorithmic thinking, and get you familiar with Matlab (youll thank us once you get to 103C, 126, and 163). This week we address perhaps the most important topic in the study of chemical kinetics: the solution of nonlinear ordinary differential equations. Stripped of its physical richness, all of kinetics can be reduced to solving the most general initial value problem of the form y ( t ) = f ( t, y ) , y ( t ) = y (1) Here, y is a vector of coupled variables { y 1 ,y 2 ,...,y n } , which can include species concen trations, reactor temperature (for nonisothermal problems), etc. Unfortunately, Equation 1 cannot be solved analytically except for very special cases in which the equation is either separable or can be made separable with an integrating factor, so we would like to develop numerical methods to obtain approximate solutions with arbitrary accuracy. (a) Eulers method: The simplest of all numerical integrators, and probably a method youve seen in high school calculus. If we wish to integrate Equation 1 from t up to t f , then we divide up the interval [ t ,t f ] into n segments with length h t/n , where t = t f t . Along each segment, we approximate the function y ( t ) by a straight line, with slope given by f ( t, y ) evaluated at the preceding mesh point. As such, Eulers method can be viewed as a linear interpolation scheme for ODEs, with the simple updating rule y i + 1 = y i + h f ( t i , y i ) i. Implement Eulers method as a Matlab function that takes in a function handle for the right hand side of Equation 1, the initial and final time points t and t f , the number of steps n , and the vector y of initial conditions. The function should return a matrix containing the values of y ( t ) at every mesh point in the simulation. ii. To see how this method performs on a real problem, consider the firstorder, re versible isomerization A k f GGGGGGB FGGGGGG k b B A. Write out the system of coupled differential equations governing the kinetics of this reaction in a batch reactor....
View
Full
Document
 Winter '11
 ARNOLD

Click to edit the document details