Solve a PDE with boundary conditions (chemical adsorption in fixed beds)?

Questions
, ,
1

Dear Wolfram team:

I have been trying for week to solve a system of 2 partial differential equations describing the adsorption of a chemical substance on a fixed bed (for example, a column of activated carbon). The 2 equations are the following, taken from McCabe (1993):

EQ1.png
EQ1.png541×273 55 KB

EQ2.png
EQ2.png540×244 74 KB

Unfortunately I cannot get past the general solution (with arbitrary constants) because when I try to put boundary conditions the Mathematica program fails. Maybe I am using the wrong command or syntax, or maybe there are too much or too few boundary conditions.

I have left attached the program, where I tryed to simplify the problem combining both equations in a third.
Thank you in advance for your help.

Best regards,

Alberto Silva

Fixed bed adsorption.nb (364 KB)

2

Dear Wolfram team:

Any idea about how to solve the boundary value problem?
The result should be an adsorption wave approximately like this:

sorptionwave.png
sorptionwave.png1636×532 28.3 KB

With an exponential profile that moves way from the origin as time passes.

However, all the solutions I can get are static (i.e. independent of time) something that is absurd with waves. How can I get a solution that depends both of space (x coordinate) and time?

Thank you in advance,

Best regards,

Alberto Silva

3

I will suggest posting the equations in the actual not rather than an attachment, and using only simple Latin alphabet ascii characters as variables (no subscripts, no Greek letters, etc.) This will make it much easier for readers who may wish to play around with the system in question.

4

The meaning of terms in the code (I have put a “->” in the legend where I have subtitututed a greek or script letter with a latin one) is this:

  1. epsilon → “e” (void fraction, unitless)
  2. Kappa → “Kd” (mass transfer coefficient, assumed constant, units: 1 / time)
  3. “a” : surface area of adsorbing particles
  4. “c”: aqueous concentration of target substance (units:mass/volume)
  5. “w”: adsorbed concentration of target substance in the solid particles (units:mass/volume)
  6. “u”: bulk liquid velocity (units: lenght/time)

Now I will post the code. I will greatly appreciate your help.

ClearAll["Global`*"]
{x, t, t1, eqT, bcon, c, c0} =.

eq1 := u D[c[t, x], x] + e D[c[t, x], t] + (1 - e) D[w[t, x], x]
eq2 := -a Kd c[t, x] + (1 - e) D[w[t, x], x]
eqT := FullSimplify[eq1 - eq2]

bcon1 := {DirichletCondition[c[t, x] == c0  , {x == 0, t >= t1}], 
  DirichletCondition[c[t, x] == c0 E^(-((x a Kd)/u)) , {t == t1}], 
  DirichletCondition[c[t, x] == c0*cb , {x == L, t == tb}]}

GeneralSolution -> FullSimplify [DSolve[{eqT == 0}, {c[t, x]}, {x, t}]]
RelevantSolution -> 
 FullSimplify [DSolve[{eqT == 0, bcon1}, {c[t, x]}, {x, t}]]
5

MODERATION TEAM NOTE: This comment was moved. New location of this comment is at:

http://community.wolfram.com/groups/-/m/t/1470321

6

Due to incredible work that wen into Tim Laska’s comment, it was moved to a separate stand-alone thread and placed in Staff Picks editorial collection. This is now reflected with a badge on @Tim Laska 's Profile, - thank you!

New location of this comment: http://community.wolfram.com/groups/-/m/t/1470321

7

Could you share the title or DOI of the book from which you took the screenshots? Thanks!

8

Dear Tom Water:

Good Morning.
I now post the full reference of the screenshot:

TITLE: Unit operations of chemical engineering - 5th ed; International
Editions 1993

ISBN: 0-07-112738-0

AUTHORS: Warren L. McCabe, Julian C. Smith, Peter Harriott.

Citation details:

CHAPTER: "Chapter 25: Absorption "

SECTION: “Basic equations for absorption: Rates of mass transfer” (pages 825-826):

N.B.: apologies for the 7-year delay, I should have done so in 2018…