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

Posted 10 months ago
1737 Views
|
5 Replies
|
3 Total Likes
|
 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):  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 Attachments: Answer
5 Replies
Sort By:
Posted 10 months ago
 Dear Wolfram team:Any idea about how to solve the boundary value problem? The result should be an adsorption wave approximately like this: 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 Answer
Posted 10 months ago
 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. Answer
Posted 10 months ago
 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: epsilon -> "e" (void fraction, unitless) Kappa -> "Kd" (mass transfer coefficient, assumed constant, units: 1 / time) "a" : surface area of adsorbing particles "c": aqueous concentration of target substance (units:mass/volume) "w": adsorbed concentration of target substance in the solid particles (units:mass/volume) "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}]] ` Answer
Posted 9 months ago
 MODERATION TEAM NOTE: This comment was moved. New location of this comment is at:http://community.wolfram.com/groups/-/m/t/1470321 Answer
Posted 8 months ago
 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 Answer
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Reply Preview
Attachments