reaction_twoway_ode

reaction_twoway_ode, a MATLAB code which defines the ordinary differential equations (ODE) which model a two-way chemical reaction between species W1 and W2.

The reference suggests an experiment in which the Forward Euler method is used, with a timestep of 1/50. The problem is to be solved three times, with the second reaction coefficient, k2, set to 10 (the default), 100, and 1000.

Licensing:

The computer code and data files described and made available on this web page are distributed under the GNU LGPL license.

Reference:

1. Willem Hundsdorfer, Jan Verwer,
   Numerical solution of time-dependent advection-diffusion-reaction equations,
   Springer, 2003
   ISBN: 978-3-662-09017-6

Source Code:

• reaction_twoway_conserved.m, returns the value of a quantity that should be conserved.
• reaction_twoway_deriv.m, returns the right hand side of the ODE.
• reaction_twoway_exact.m, returns the exact solution of the ODE.
• reaction_twoway_parameters.m, returns parameters of the ODE.
• timestamp.m, prints the YMDHMS date as a timestamp.

Tests:

• reaction_twoway_ode_test.m, calls all the tests.
• reaction_twoway_ode_test.sh, runs all the tests.
• reaction_twoway_ode_test.txt, the output file.

Plots:

• reaction_twoway_conserved.png, displays the conserved quantity.
• reaction_twoway_solution.png, displays the solution.