Can I use the PDEPE solver to solve the following dimensionless Transient PFR PDE? If not, what method do you reconmend?

How do I expresss the last term of the PDE in the code?

Answer 1

[David Goodmanson]

Hello Rogelio,

For psi(t,eta) I assume that the boundary conditions are

psi(0,eta)                     = 0     for all eta
and
d psi(t,eta) / deta |          = 0     for all t
                    | eta = 0

and a solution is intended for t >= 0. In that case a solution that fits the bc's is

  psi(t,eta) = 1 -exp(-t/lambda)

i.e. no eta dependence at all.

[Roger Vegeta]

psi(0,eta) = 0 ; for all eta.

d(psi)/d(eta) = 0 ; for eta = 0 for all time.

[David Goodmanson]

yes. so the solution from my previous comment works.