You have reached the web page for the book
Computational Methods for Process Simulation
and CHEN 4580 - Numerical Methods for Simulation.

Please click on the following to send me an email.
Fred Ramirez


Instructions for downloading
Click on the filename to view and then choose File-SaveAs to save the file.
Short-cut: Shift-LeftMouseButton goes directly to the download prompt!



Chapter 2: Steady-State Lumped Systems

  • bisec.m - matlab file to determine the root of Equation (2.3.2) using the bisection method.
  • cp1.m - matlab file that defines Equation (2.3.2) in a format suitable for the bisection and tangent method.
  • cp2.m - matlab file that defines Equation (2.3.5) in a format suitable for the wegstein method.
  • cp3.m - matlab file that defines Equation (2.3.11) in a format suitable for the wegstein method and converges using direct substitution.
  • cstr.m - matlab file that defines Equations (2.3.30) and (2.3.31) to be used by newtrap.m
  • fig2-28.mat - input file for the program fun.f implementing the example of Section 2.4.3
  • fun.f - fortran program that implements the equation ordering algorithm and the variable group algorithm
  • lin2_1.m - matlab file to solve the set of linear equations in Example 2.1
  • lin2_2.m - matlab file to solve the set of linear equations in Example 2.2
  • lua.m - matlab file to determine a LU factorization
  • mix1.mat - input file that implements the mixer-exchanger-mixer problem with the three design variables remaining unspecified
  • mixer.mat - input file that inplements the mixer-exchanger-mixer problem with specified design variables
  • newtrap.m - matlab file that solves Equations (2.3.30) and (2.3.31)
  • prob2-7.mat - functionality matrix file for the program fun.f implemeting Problem 2.7
  • regfals.m - matlab file that implements the regula-falsi algorithm
  • tangent.m - matlab file that implements the tangent algorithm
  • wegstein.m - matlab file that implements the wegstein accelerator


    Chapter 3: Unsteady-State Lumped Systems

  • dfun32.m - defines the dimensionless differential equation of Equation (3.2.7)
  • dtank32n.m - matlab file to solve the tank problem of Section (3.2)
  • ex310.m - matlab file to solve the problem of Section (3.10)
  • fun34.m - defines the dimensionless model of the tank problem of Section (3.4)
  • fun35.m - defines the enclosed tank vessel problem of Section (3.5)
  • model310.m - defines the model for the problem of Section (3.10)
  • ntrp113.m - file needed for ode113.m
  • numjac.m - file needed for ode15s.m
  • ode113.m - a multi-step integration routine (needs odeget.m, odeset.m and ntrp113.m)
  • ode15s.m - a stiff integration routine (needs odeget.m, odeset.m and numjac.m)
  • ode23s.m - another stiff integration routine (needs odeget.m, odeset.m and numjac.m)
  • ode45n.m - new ode45 integration routine (needs odeget.m and odeset.m)
  • odeget.m - file needed for the ode45n.m among others
  • odeset.m - file needed for the ode45n.m among others
  • tank34.m - matlab file to solve the tank problem of Section (3.4)
  • tank35.m - matlab file to solve the enclosed tank problem of Section (3.5)


    Chapter 4: Reaction-Kinetic Systems

  • ex41.m - matlab file to solve the problem of Section (4.1)
  • ex42.m - matlab file to solve the problem of Section (4.2)
  • ex43.m - matlab file to solve the problem of Section (4.3) using ode23s
  • ex432.m - matlab file to solve the problem of Section (4.3.2)
  • ex43a.m - matlab file to solve the problem of Section (4.3) using ode15s
  • labmodel.m - defines the differential equations for the laboratory CSTR problem of Section (4.3.1)
  • labreact.m - matlab file to solve the laboratory CSTR problem of Section (4.3.1)
  • model41.m - defines the differential equations for the problem of Section (4.1)
  • model42.m - defines the differential equations for the problem of Section (4.2)
  • model43.m - defines the differential equations for the problem of Section (4.3)
  • model432.m - defines the differential equations for the problem of Section (4.3.2)


    Chapter 5: Vapor-Liquid Equilibrium Systems

  • ex51.m - matlab file to solve the problem of Section (5.1)
  • ex51.m - matlab file to solve the problem of Section (5.6)
  • eq51.m - defines the algebraic equations for the problem of Section (5.1)
  • impT.m - file required for the solution of the problem in Section (5.6)
  • model51.m - defines the differential equations for the problem of Section (5.1)
  • model56.m - defines the differential equations for the problem of Section (5.6)
  • Pb.m - file required for the solution of the problem in Section (5.6)
  • Pt.m - file required for the solution of the problem in Section (5.6)
  • Pz.m - file required for the solution of the problem in Section (5.6)


    Chapter 7: Solution of Split Boundary-Value Problems

  • ex71.m - matlab file to solve the tubular reactor problem in Section 7.1 using the shooting technique.
  • ex72.m - matlab file to solve the tubular reactor problem in Section 7.1 using the generalized shooting technique.
  • ex74.m - matlab file to solve the packed-bed reactor problem in Section 7.4 using the superposition principle.
  • ex76.m - matlab file to solve the nonlinear tubular reactor with dispersion in Section 7.6 using the quasilinearization technique.
  • ex77.m - matlab file to solve the problem in Section 7.7 using the method of adjoints.
  • model71.m - defines the model of Section 7.1. Needed by ex71.m
  • model72.m - defines the model of Section 7.2. Needed by ex72.m
  • model74.m - defines the model of Section 7.4. Needed by ex74.m
  • model76.m - defines the model of Section 7.6. Needed by ex76.m
  • model77.m - defines the model of Section 7.7. Needed by ex77.m
  • model77a.m - defines the adjoint equations of Section 7.7. Needed by ex77.m


    Chapter 8: Solution of Partial Differential Equations

  • ex82.m - matlab file to solve the problem in Section 8.2
  • ex83.m - matlab file to solve the problem in Section 8.3 (needs Model83.m)
  • ex84.m - matlab file to solve the problem in Section 8.5
  • Model83.m - defines the differential equations for the heat exchanger dynamics