Pdepe set up for Langmuir diffusion process

Hi all. Is it possible to solve langmuir diffusion process with matlab pdepe

dc(x,t)/dt=dif*d2c(x,t)/dx2

dgama(t)/dt=k1*c(0,t)*(gama0-gama(t))-k2*gama(t)

boundary c

dif*dc(x,t)/dx=k1*c(x,t)*(gama0-gama(t))-k2*gama(t) at x=0

c(d,t)=1e-6

ic c(x,0)=1e-6 gama(0)=0 Atemp in attached files and model recalls on article in attached files

Thank You

Answer 1

[Torsten]

pdepe can not be used to solve mixtures of partial differential equations (your equation for c) and ordinary differential equations (your equation for gama).

You will have to discretize your equation for c in space on your own and solve the resulting complete system of ordinary differential equations using ODE15S, e.g.

Look up "Method-of-Lines" for more details.

Best wishes

Torsten.