©
Chemical reaction parameters
Computes reaction rate constants.
2024.Jul.03 18:16:26
File to upload File* (optional), a numerical matrix (as next). •
Experimental
data
Kinetic data (each row: time, 3 concentrations).
Ignored iff a file is uploaded. •
Action   Verify with parameters below; or Optimize. •
Parameters Rate constants (initial guess). •
Scale factors Scale factors ("steps"). •
dt, tol, max iter. Integr. step, tolerance, max. iterations for adjustment. •
Show values ? Shows the graph coordinates.

Proceeds to verify or optimize the kinetics rate constants: verify, using the given values to compare computed vs. experimental concentrations;  or optimize, to get the values for the best fit, through sequential simplex minimization.

The problem is based on a particular embedded reaction mechanism, A → B → C, through a system of 3 ODEs, ordinary differential equations.

Equivalent Matlab and Mathematica scripts are given.

*An uploadable file with data identical to the above is supplied for simplicity here: ABC.dat.

Plots the experimental and computed data vs time.

The data are from Winkel [2011].

Acknowledgement — MSc Engineer Ricardo Nunes (Mathematica).

References: Plate: ChemReacOpt

• Google: "chemical kinetics" • Wikipedia: Chemical kinetics

• Winkel, B..pdf, 2011, International Journal of Mathematical Education in Science and Engineering, 42(1), 37–51.

• Burkardt, J., 2011, ASA047 (Nelder-Mead minimization).

• Wikipedia: Nelder-Mead method

• 1765-12-22: Pfaff, Johann Friedrich (†1825-04-21, 59 yrs.).

 
 
Valid HTML 4.01! IST http://web.ist.utl.pt/~ist11038/compute/com/Fx-ChemReacOpt.php
Created: 2019-12-22 — Last modified: 2024-05-30