This sort of problem is well documented, but, I can not see how to get this solution;

The problem is:

$$\frac{\partial c}{\partial t}=D\frac{\partial^2 c}{\partial x^2}, 0 < x < \infty, t > 0$$

with boundary conditions $c(0,t)= C_0, c(x,0) = 0$ and possibly $\lim_{x \rightarrow \infty}c(x,t) = 0$ as the capillary that the model is modelling is infinitely long.

whose solution is

$$c(x,t) = 2C_0\bigg( 1-\frac{1}{\sqrt{2\pi}}\int^z_{-\infty}\exp(-\frac{s^2}{2}) ds \bigg), z = \frac{x}{\sqrt{2Dt}}$$

Now I have let $c(x,t) = X(x)T(t)$ and then obtained

$$\frac{T'}{DT} = - s^2 = \frac{X''}{X}$$

we have

$\begin{cases}T(t)=c_3(s)e^{-Dts^2}\\X(x)=\begin{cases}c_1(s)\sin xs+c_2(s)\cos xs&\text{when}~s\neq0\\c_1x+c_2&\text{when}~s=0\end{cases}\end{cases}$


$$ u(x,t)=\int_0^\infty C_1(s)e^{-Dts^2}\sin xs~ds+\int_0^\infty C_2(s)e^{-Dts^2}\cos xs~ds$$

but my initial conditions give, first

$$\int_0^\infty C_2(s)e^{-Dts^2}ds = C_0$$

then changing variables $s = \sqrt{\frac{r}{D}}$

$$\int_0^\infty \frac{1}{2}\sqrt{\frac{D}{r}}C_2(r)e^{-tr}ds = C_0$$

$$ L[ \frac{1}{2}\sqrt{\frac{D}{r}}C_2(r)] = C_0$$

the inverse Laplace transformation of the right and side is a multiple of the delta function.... I feel like I have taken the wrong way.

  • $\begingroup$ The expressions $X(x)$ and $T(t)$ you found are both false. Even if you correct them, the particular solution with respect to the boundary conditions will require not only one function $c(x,t)=X(x)T(t)$ but the sum of an infinity of such functions with an infinity of values of $k$. This approach would be theoretically possible, but very arduous. $\endgroup$ – JJacquelin Sep 23 '16 at 9:21
  • $\begingroup$ Note that the method of separation of variables $X(x)T(t)$ is convenient in cases of some "educational" problems where the boundary conditions are especially chosen to make successful the calculus with this method. $\endgroup$ – JJacquelin Sep 23 '16 at 9:32
  • $\begingroup$ Right. Is there a method you can suggest. Or maybe an initial step to take? $\endgroup$ – user197848 Sep 23 '16 at 16:27
  • $\begingroup$ Since the solution is known I suppose that the calculus behind is published or at least the references about the method used. $\endgroup$ – JJacquelin Sep 23 '16 at 17:16
  • $\begingroup$ The key idea is to remplace a discrete infinite sum of terms $\sum_sC_sX(x,s)T(t,s)$ by an integral $\int C(s)X(x,s)T(t,s)ds$. For example see math.stackexchange.com/questions/1486114/… $\endgroup$ – JJacquelin Sep 23 '16 at 17:47

Of course use separation of variables:

Let $c(x,t)=X(x)T(t)$ ,

Then $X(x)T'(t)=DX''(x)T(t)$



$\begin{cases}T(t)=c_3(s)e^{-Dts^2}\\X(x)=\begin{cases}c_1(s)\sin xs+c_2(s)\cos xs&\text{when}~s\neq0\\c_1x+c_2&\text{when}~s=0\end{cases}\end{cases}$

$\therefore c(x,t)=\int_0^\infty C_1(s)e^{-Dts^2}\sin xs~ds+\int_0^\infty C_2(s)e^{-Dts^2}\cos xs~ds$

$c(0,t)=C_0$ :

$\int_0^\infty C_2(s)e^{-kts^2}~ds=C_0$


$\therefore c(x,t)=\int_0^\infty C_1(s)e^{-Dts^2}\sin xs~ds+\int_0^\infty C_0\delta(s)e^{-Dts^2}\cos xs~ds=\int_0^\infty C_1(s)e^{-Dts^2}\sin xs~ds+C_0$

$c(x,0)=0$ :

$\int_0^\infty C_1(s)\sin xs~ds+C_0=0$

$\mathcal{F}_{s,s\to x}\{C_1(s)\}=-C_0$

$C_1(s)=\mathcal{F}^{-1}_{s,x\to s}\{-C_0\}=-\dfrac{2C_0}{\pi s}$

$\therefore c(x,t)=C_0-\dfrac{2C_0}{\pi}\int_0^\infty\dfrac{e^{-Dts^2}\sin xs}{s}~ds=C_0~\text{erfc}\left(\dfrac{x}{\sqrt{4Dt}}\right)$

It luckily satisflies $\lim\limits_{x\to\infty}c(x,t)=0$ .

$\therefore c(x,t)=C_0~\text{erfc}\left(\dfrac{x}{\sqrt{4Dt}}\right)$

| cite | improve this answer | |
  • $\begingroup$ Very nice : The Fourier sin transform to find $C_1(s)$. $\endgroup$ – JJacquelin Sep 25 '16 at 5:23
  • $\begingroup$ Could you expand a little on the inverse Fourier equality. I understand that you have used the Fourier transform. I am trying to follow what is on wikipedia on the transform. $\endgroup$ – user197848 Sep 26 '16 at 15:50
  • $\begingroup$ I get $C_1(s) = -\frac{2}{\pi} \int_{0}^{\infty}C_0 \sin(xs)dx = \frac{2C_0}{\pi s} \bigg([\frac{\cos(ys)}{s}]_{y \rightarrow \infty}-1 \bigg)$ how do i deal with the first term? $\endgroup$ – user197848 Sep 26 '16 at 17:46
  • $\begingroup$ @HMPARTICLE When $\int_0^\infty f(x)\sin sx~dx$ diverges, we shouldn't simply recognize $\mathcal{F}^{-1}_{s,x\to s}\{f(x)\}=\dfrac{2}{\pi}\int_0^\infty f(x)\sin sx~dx$ and should consider at other suitable approach. $\endgroup$ – doraemonpaul Sep 28 '16 at 12:19
  • $\begingroup$ Could you hint at another approach? $\endgroup$ – user197848 Sep 28 '16 at 21:00

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy