Matlab Review.

You will find Matlab solutions and m-files you can download and have that illustrate the use of matrix operations, solving linear ODEs, and solving sets of nonlinear algebraic equations. The problem statement for the matrix operations file and the ODE file is in the pdf file titled Matlab Review. Problem 1 is solved using the m-file titled MatLabLec1.m. Problem 2 is solved using two approaches. In the first you need two files, MatLabLec2.m and m1p2.m; MatLabLec2.m is the script file and m1p2.m is the function file (subroutine) that contains the list of differential equations to be solved. In the second approach, both files are merged into one function file, MatLabLec2_Complete.m. Study both approaches. The third problem is Text Exercise 3.3 and is solved with the file exer3_3.m.

Read about how to solve ODEs using explicit methods.
Runge Kutta Primer

Class problem illustrating variable density in a CSTR.
Using the data from Exercise 4.8 we consider three cases: (1) variable density with the volume specified as constant, (2) variable density with the mass specified as constant and set equal to the mass present initially if charged with unreacted feed, and (3) constant density. The Problem Statement and the Matlab program script file that solves it (cstr_variable_rho.m).

Class Example of 3 CSTRs in Series

Matlab Script File for the Benzene Pyrolysis Text Example

Nonisothermal Batch Reactor Problem Statement & the MatLab m-file that solves it.

Nonisothermal CSTR Problem Statement and the MatLab m-file that solves it.

CSTR Stability Problem Statement and the MatLab m-file that solves it.

Slides illustrating safety considerations in the operation of a batch reactor and how a temperature run away contributed to an accident. This is case study prepared by Professor Ronald J. Willey of Northeastern Univerisity for SACHE. The slides are for personal study use only.

There are two different pdf files that describe parameter estimation and three Matlab m-files. One set is (Parameter Estimation Using Solved Reactor Models and parameter_estimation.m). The second set is (Parameter Estimation Using the Material Balances and the two programs parameter_1_est.m and parameter_2_est.m).

Read about Effective Diffusivity

Programs that illustrate collocation and solve mass transfer with reaction in a porous catalyst pellet

Read a document in the Learning Module "Numerical Approximation to Reaction with Diffusion" titled "Collocation Primer" that describes the collocation approach and then run the Octave program or download the Matlab program (fstord_pellet_Master_MATLAB.m) to explore the solution as the Thiele modulus and number of collocation points are changes.

Read pdf file for the Learning Module "Reaction with Diffusion - Examining External Gradients" that describes the numerical approach to handling external gradients and then run a program (fstord_external_MASTER_MATLAB.m) to explore the impact of an external gradient on the concentration profile and the rate of reaction. Solve Example 7.1 and explore the effect of external mass transfer limitations.

The collocation function (put this in the same directory as the preceding programs)

Solve Example 7.5 using collocation to illustrate how the pellet profile solution is merged with an ODE solver for a fixed-bed

Revisit Example 7.5 to explore how an external gradient is incorporated into a fixed bed calculation (Example7_5_external.m) and to see the effect on bed volume


Mixing Problem Statement and the Matlab file that solves the mixing problem (mixing). You can also download the Excel file that was illustrated in class (MixingExcel).