ProblemSet1 - ChE 101 2012 Problem Set 1 1. Computational...

Info iconThis preview shows pages 1–2. Sign up to view the full content.

View Full Document Right Arrow Icon

Info iconThis preview has intentionally blurred sections. Sign up to view the full version.

View Full DocumentRight Arrow Icon
This is the end of the preview. Sign up to access the rest of the 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 non-isothermal 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 first-order, 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

Page1 / 4

ProblemSet1 - ChE 101 2012 Problem Set 1 1. Computational...

This preview shows document pages 1 - 2. Sign up to view the full document.

View Full Document Right Arrow Icon
Ask a homework question - tutors are online