Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Finding the Inverse Laplace Transform (Bromwich Contour)

  1. Aug 13, 2011 #1
    1. The problem statement, all variables and given/known data

    Find the inverse laplace transform: [tex] {\cal L}^{-1}\{e^{-\sqrt{s-a}c}\} [/tex]

    where [itex] a [/itex] and [itex] c [/itex] are constants. [itex] s [/itex] is the complex variable.

    2. Relevant equations

    Bromwich integral: [tex] \frac{1}{2\pi{}i}\lim_{T\rightarrow \infty}\int_{\gamma-Ti}^
    {\gamma+Ti}e^{st-\sqrt{s-a}c}ds [/tex]

    Residue theorem: [tex] \oint f(s) ds=\sum Res(f(s)) [/tex]

    3. The attempt at a solution

    There is a branch point at [itex] s=a [/itex]. Also, there are no poles. Therefore, the sum of the residues is equal to zero. So I have a keyhole contour centered at [itex] a [/itex] with an outer radius of [itex] R [/itex] and inner radius of [itex] \varepsilon [/itex]. I chose [itex] \gamma=2a [/itex]. I attached a general image of the contour that I used. The contour is centered at [itex] a [/itex].

    Therefore, [itex] s=Re^{i\theta}+a [/itex] along the the outer circle,
    [itex] s=xe^{\pi i}+a [/itex] going from left to right along the branch cut,
    [itex] s=\varepsilon e^{i\theta}+a [/itex] along the the inner circle, and
    [itex] s=xe^{-\pi i}+a [/itex] going from right to left along the branch cut.

    From Jordan's lemma the integral along the outer radius goes to zero as [itex] R\rightarrow \infty [/itex]. Similarly, the integral along the inner circle (keyhole) goes to zero as [itex] \varepsilon\rightarrow 0 [/itex]


    [tex] \frac{1}{2\pi{}i}\lim_{T\rightarrow \infty}\int_{\gamma-Ti}^
    {\gamma+Ti}e^{st-\sqrt{s-a}c}ds + \frac{1}{2\pi{}i}\int_{branch cut}=0 [/tex]

    [tex] \frac{1}{2\pi{}i}\lim_{T\rightarrow \infty}\int_{\gamma-Ti}^
    {\gamma+Ti}e^{st-\sqrt{s-a}c}ds =

    \frac{-1}{2\pi{}i}\int_{branch cut}

    =\frac{-1}{2\pi{}i}\lim_{\substack{R\rightarrow \infty\\\varepsilon\rightarrow 0}}


    \int_{-\varepsilon+a}^{-R+a}e^{st-\sqrt{s-a}c}ds) [/tex]

    [tex] = \frac{-1}{2\pi{}i}(-\int_{a}^{\infty}e^{(a-x)t-\sqrt{x}ic}e^{\pi{}i}dx+\int_{a}^{\infty}e^{(a-x)t+\sqrt{x}ic}e^{-\pi{}i}dx) [/tex]

    [tex] = \frac{-1}{2\pi{}i}(\int_{a}^{\infty}e^{(a-x)t-\sqrt{x}ic}dx-\int_{a}^{\infty}e^{(a-x)t+\sqrt{x}ic}dx)

    =\frac{e^{ta}}{\pi{}}\int_{a}^{\infty}e^{-xt}\sin{(\sqrt{x}c)}dx [/tex]

    And this is where I'm stuck. I think this integral doesn't have a closed form. I can expand the integral into an infinite series but I don't know how to integrate the terms. I tried doing the Inverse Laplace transform in mathematica. The answer it gave me is:

    [tex] \frac{ce^{\frac{-c^2}{4t}+at}}{2\sqrt{\pi{}}}t^{\frac{3}{2}} [/tex]

    So what I ended up getting doesn't agree with Mathematica. Also I know from tables that the inverse laplace transform of [tex] e^{-\sqrt{s}} [/tex] is [tex] \frac{a}{2\sqrt{\pi{}t^3}}e^{\frac{-a^2}{4t}} [/tex] I really don't know what I did wrong. I don't know what to do from here. Thanks in advance.

    Attached Files:

  2. jcsd
  3. Aug 13, 2011 #2

    Ray Vickson

    User Avatar
    Science Advisor
    Homework Helper

    If you can find the inverse transform of exp(-c*sqrt(s)), you can use standard theorems to get the inverse transform of exp(-c*sqrt(s-a)). Now c*sqrt(s) = sqrt(r*s), where r = c^2, so if f(t) = inverse of exp(-sqrt(s)), we have exp(-sqrt(r*s)) = int_{t=0..inf} f(t) exp(-t*(r*s)) dt = int_{t=0..inf} f(t/r)*exp(-t*s) dt/r (r*t <---t), so inv of exp(-sqrt(r*s)) is f(t/r)/r = f(t/c^2)/c^2.

  4. Aug 13, 2011 #3


    User Avatar
    Staff Emeritus
    Science Advisor
    Homework Helper
    Education Advisor

    Try using the substitution [itex]u=s-a[/itex] first to get[tex]\int e^{st}e^{-c\sqrt{s-a}}\,ds = e^{at}\int e^{ut}e^{-c\sqrt{u}}\,du[/tex]Then you could complete the square
    [tex]ut-c\sqrt{u} = t\left(\sqrt{u}-\frac{c}{2t}\right)^2-t\left(\frac{c}{2t}\right)^2[/tex]to get
    [tex]e^{at}\int e^{ut}e^{-c\sqrt{u}}\,du = e^{at-\frac{c^2}{4t}}\int e^{t(\sqrt{u}-\frac{c}{2t})^2}\,du[/tex]
  5. Aug 13, 2011 #4

    I follow everything you wrote except "exp(-sqrt(r*s)) = int_{t=0..inf} f(t) exp(-t*(r*s)) dt"

    [tex] f(t)= {\cal L}^{-1}\{e^{-\sqrt{s}}\}[/tex]

    [tex] e^{-\sqrt{rs}}=e^{-c\sqrt{s}}=(\int_{0}^{\infty}f(t)e^{-t{s}}dt)^c [/tex]

    I don't see how to get that to [tex] \int_{0}^{\infty}f(t)e^{-trs}dt [/tex] :confused:

    Vela, that looks interesting. So with that substitution I end up with a branch point at u=0 right? I'll see what i get with the contour and get back to you. Thanks for the help guys. ^^
  6. Aug 13, 2011 #5
    Your analysis is correct as I see it except you're not integrating over the correct bounds. If [itex]z=a+re^{\pi i}[/itex], what must the lower limit on r be in order to reach a?
    Last edited: Aug 13, 2011
  7. Aug 13, 2011 #6

    Ray Vickson

    User Avatar
    Science Advisor
    Homework Helper

    It's just the statement exp(-sqrt(S)) = int_{0..inf} f(t) exp(-S*t) dt, with S replaced by r*s!

  8. Aug 13, 2011 #7
    From similar steps as in the first post and letting [itex] u=s-a [/itex] I get the final integral
    [tex] \frac{1}{2\pi{}i}\lim_{T\rightarrow \infty}\int_{\gamma-Ti}^
    {\gamma+Ti}e^{st-\sqrt{s-a}c}ds =\frac{e^{ta}}{\pi{}}\int_{0}^{\infty}e^{-xt}\sin{(\sqrt{x}c)}dx [/tex]

    I evaluated the integral with mathematica and it agrees with what i got by just taking the inverse directly with mathematica. The only thing it is a conditional expression valid with t>0. But doing the inverse directly it doesn't give that condition... thoughts?

    jackmell, is this what you were referring to? that the limits should be from zero to infinity? I'm not sure how u figured that out without making the substitution. Without making the substitution [itex] z=a+Re^{\pi i} [/itex]...[itex] z [/itex] goes from [itex] -R+a [/itex] to [itex] -\varepsilon+a [/itex]. Then taking the limit as R goes to infinity and epsilon goes to zero aren't the limits from a to infinity?

    ray, so what ur suggesting is to do a substitution...I saw the constant outside the square root and had immediately thought it didn't match with the table...lol guess i made this problem more complicated than it needed to be. -_-;

    So out of curiosity how would I evaluate this final integral? I can do a series expansion....but i don't know how to integrate the individual terms.
  9. Aug 13, 2011 #8
    Assume a is on the real axis for simplicity, then you're integrating over a branch-cut [itex](-\infty,a)[/itex]. Ok then, if I let [itex]z=a+re^{\pi i}[/itex], and [itex]z=a+re^{-\pi i}[/itex] like you and I both did for the analysis, and let r go from [itex](\infty,0)[/itex], does not z traverse over that entire branch-cut except for the indentation around the branch-point?
    Last edited: Aug 13, 2011
  10. Aug 13, 2011 #9


    User Avatar
    Staff Emeritus
    Science Advisor
    Homework Helper
    Education Advisor

    By default, Mathematica uses the one-sided Laplace transform, so it's understood the inverse transform is valid only for t>0.

    When you evaluate the Bromwich integral, you get to choose which way to close the contour. When t>0, the integral over the outer circle vanishes only when you close the contour in the left half-plane. When t<0, you have to close it in the right half-plane. There's nothing going on the right half-plane, so you get f(t)=0 for t<0.
  11. Aug 14, 2011 #10
    I don't think we have to Vela. I'll conjecture we can still close the contour over the left half-plane when t<0 if we take the limit over all the contours simultaneously because in the limit, Cauchy's Theorem still applies. I can't at the moment however, show this.
    Last edited: Aug 14, 2011
  12. Aug 14, 2011 #11
    ohhh, i get it. i didn't change my bounds when I did the substitution! thanks.

    So what happens at t=0? Is the solution not valid? :confused:
  13. Aug 14, 2011 #12
    What happens to the expression for t>0 as t tends to zero?
  14. Aug 15, 2011 #13
    [tex] \frac{ce^{\frac{-c^2}{4t}+at}}{2\sqrt{\pi{}}t^{\frac{3}{2}}} [/tex] goes to zero as t goes to zero...well i'm almost 100% sure its valid for t=0, since other inverse laplace transforms are valid for t=0. I was wondering how that translates to the bromwich contour since mathematica gave the conditional existence of t>0 for the final integral i got from using residue theory.
Share this great discussion with others via Reddit, Google+, Twitter, or Facebook