MATLAB: Modeling Chemical Reactions…….

Question

Hi everyone.I am not so familiar with MATLAB or programming.I am writing a code for Solar Model.As a part of it i need to model some chemical reactions.

Ex:

            p     +    p  -------->  D   +  (e^+)  +  neutrino
            D     +    p  -------->  3He +   photons

Will it be possible? If yes, how?

Answer 1

Modeling chemical reactions is not very difficult.

first, rewrite the chemical reaction as a differential equation.

-

example for reaction A->P (with rate constant 'k1')

-

Rate equation: Ra=-k1*Ca

ODE: dCa/dx = -Ra (x either being time or axial position in a PFR)

-

Now you can solve this with an ordinary ODE-solver.