I need a Matlab code or a Polymath code to solve this differential equation and get a graph of CA (of some species) vs radial direction.
Catalyst: Nafion 100g
T = 400K
P = 1340KPa
Details:
For the first-order reaction A -> B with rate law r = kC, the differential equation that describes diffusion and reaction in the catalyst pellet is:
d^2CA/dr^2 + 2dCA/dr - kCA = 0
Where:
r is the radius of the catalyst pellet (cm)
C is the concentration in the catalyst (mmol/cm^3)
k is the first-order rate constant (1/s)
De is the effective diffusivity (cm^2/s)
DAB is the pellet porosity
a is the constriction factor
ε is tortuosity
The boundary conditions for the differential equation in spherical coordinates are:
C = finite at r = 0
CA = CAs at r = R
Where R is the pellet radius (cm) and CAs is the concentration (mmol/cm^3) at r = R.