How to prove solution to Fick's second law by substitution

Join the discussion
Ask a follow-up here, or get your own question answered by working scientists, mathematicians and engineers — people, not an autocomplete.
Real named experts · corrections over time · the nuance an AI answer skips
4 replies · 4K views
neromax
Messages
6
Reaction score
0

Homework Statement


I am to prove that a solution to the differential equation Fick's second law is valid by substitution.

Homework Equations


Fick's second law:
[tex]\frac{\partial C}{\partial t} = \frac{\partial}{\partial x} \left( D \frac{C}{\partial x} \right)[/tex]
Solution to Fick's second law:
[tex]C(x,t) = \left( \frac{C_1 + C_2}{2} \right) - \left( \frac{C_1 - C_2}{2} \right) \text{erf} \left( \frac{x}{2 \sqrt{Dt}} \right)[/tex]
The (Gauss) error function (erf) which I found online:
[tex]\text{erf}(z) = \frac{2}{\sqrt{\pi}} \int_o^z e^{-y^2} dy[/tex]

The Attempt at a Solution


I presume I am supposed to differentiate the proposed solution with respect to t once and compare it to the proposed solution differentiated twice with respect to x? But I am not sure how i shall handle the integral.

Can someone help me/point me to literature or give me some pointers on how to proceed?
 
on Phys.org
@HallsofIvy - Thanks for a reply. The problem bugging me though, is that i do not have [itex]x[/itex], but [itex]\frac{x}{2\sqrt{Dt}}[/itex], such that:

[tex]\frac{2}{\sqrt{\pi}} \frac{\partial}{\partial x} \int_0^{\frac{x}{2\sqrt{Dt}}} e^{-y^2}dy[/tex]
 
Is this the correct change of limits and function of the integral?
[tex] <br /> \frac{\partial}{\partial x} C(x,t) = -B \frac{\partial}{\partial x}\int_0^{\frac{x}{2\sqrt{Dt}}} e^{-y^2}dy = -B \frac{\partial}{\partial x} \int_0^x e^{-y^2\left(4Dt \right)^{-1}}dy = -B e^{-x^2\left(4Dt \right)^{-1}}\\<br /> \frac{\partial^2}{\partial^2 x}C(x,t) = -B \frac{\partial}{\partial x} e^{-x^2\left(4Dt \right)^{-1}} = -B e^{-x^2\left(4Dt \right)^{-1}} \frac{\partial}{\partial x}\left( -x^2\left(4Dt \right)^{-1} \right) = B \frac{2}{Dt}\cdot x e^{-x^2\left(4Dt \right)^{-1}}[/tex]
 
Last edited:
Solved it!

I finally understood how it worked :wink:

[tex] <br /> D\frac{\partial}{\partial x} C(x,t) = -BD \frac{\partial}{\partial x}\int_0^{\frac{x}{2\sqrt{Dt}}} e^{-y^2}dy = -BD e^{-\left(\frac{x}{2\sqrt{Dt}}\right)^2} \frac{\partial}{\partial x} \left(\frac{x}{2\sqrt{Dt}}\right) = -BD e^{-\frac{x^2}{4Dt}} \left(\frac{1}{2\sqrt{Dt}}\right) \\= \frac{-BD}{2\sqrt{Dt}} e^{-\frac{x^2}{4Dt}} \\<br /> \frac{\partial}{\partial x} \left( D \frac{C(x,t)}{\partial x} \right) = \frac{\partial}{\partial x} \frac{-BD}{2\sqrt{Dt}} e^{-\frac{x^2}{4Dt}} = \frac{-BD}{2\sqrt{Dt}} e^{-\frac{x^2}{4Dt}} \frac{\partial}{\partial x} \left( - \frac{x^2}{4Dt} \right) \\= \frac{BDx}{4}\left( Dt \right)^{-\frac{3}{2}}e^{-\frac{x^2}{4Dt}}\\<br /> <br /> \frac{\partial}{\partial t} = -B \frac{\partial}{\partial t}\int_0^{\frac{x}{2\sqrt{Dt}}} e^{-y^2}dy = -B e^{-\left(\frac{x}{2\sqrt{Dt}}\right)^2} \frac{\partial}{\partial t} \left(\frac{x}{2\sqrt{Dt}}\right) = -B e^{-\left(\frac{x}{2\sqrt{Dt}}\right)^2} \left( - \frac{Dx}{4(Dt)^{3/2}} \right) \\= \frac{BDx}{4}\left( Dt \right)^{-\frac{3}{2}}e^{-\frac{x^2}{4Dt}} \qquad \square<br /> [/tex]
 
Last edited: