Las ecuaciones dimensionales del problema de difusión reacción son:
1)
sujeta a las condiciones de frontera:
en
2)
en
3)
en
4)
5)
sujeta a las condiciones de frontera:
en
6)
en
7)
en
8)
Después de introducir las siguientes variables adimensionales
se obtiene para la concentración adimensional:
9)
sujeta a las condciones de frontera:
en
10)
en
11)
12)
para la temperatura adimensional:
=
13)
sujeta a las condiciones de frontera
en
14)
en
15)
16)
Las ecuaciones 9) a 16) se resuelve con la instrucción
NDSolve para diferentes parámetros y se presentan los resultados en forma gráfica. Otra manera de resolver el problema es mediante colocación ortogonal.