Desorption-Fixed bed- Initial conditions

Dear all,

I am trying to solve an axial dispersion model. I have the following desorption conditions.

c = cFeed at t=0 for all z >0

q = qFeed at t=0 for all z >0

c = 0 at t>0 for z=0

For the adsorption stage I already have the conditions:

c0 = zeros(n,1);          % c = 0 at t = 0 for all z
c0(1) = cFeed;            % c = cFeed at t > 0 for z = 0
q0 = zeros(n,1);          % q = 0 at t = 0 for all z  

Any help is much appreciated.

Answer 1

[Torsten]

Then set

c0 = cFeed*ones(n,1)
c0(1) = 0.0;
q0 = qFeed*ones(n,1)

as initial and boundary conditions of the desorption, if this was your question.

[Anna]

ahh yes that makes sense. Thank you!