Register to reply

Inverse laplace transform for unique diffusion type problem

by Groundwater
Tags: diffusion, inverse, laplace, transform, type, unique
Share this thread:
Groundwater
#1
Feb20-12, 03:35 PM
P: 2
Hi all

I have been working on some unique solutions to advection-diffusion type problems.

One inverse Laplace transform that I seem to continue to encounter is the following:


Inverse Laplace[F(s)] where F(s)=[(1/(((s-α)^2)+β)*exp(-x*sqrt(s/D))]

In their classic 1959 text, Carslaw and Jaeger gave an inverse solution to F(s) for the case when β=0. It is an erfc based solution (as the exp(-x*sqrt(s/D)) would indicate).

If anyone has seen anything like this, please let me know. I have checked the following:

1. The Laplace Transform by Widder
2. Laplace Transforms and Applications by Watson
3. CRC Handbook of Tables for Mathematics, 4th edition
4. Table of Laplace Transforms by Roberts and Kaufman
5. Conduction of Heat in Solids by Carslaw and Jaeger
6. Analytical Solution to the one-dimensional advective-dispersive solute transport equation by van Genuchten and Alves
7....and a pile of elementary PDE Books

It may have to go back to first principles on this one..
Phys.Org News Partner Science news on Phys.org
Experts defend operational earthquake forecasting, counter critiques
EU urged to convert TV frequencies to mobile broadband
Sierra Nevada freshwater runoff could drop 26 percent by 2100
kai_sikorski
#2
Feb20-12, 08:17 PM
PF Gold
kai_sikorski's Avatar
P: 162
The most difficult part if you try to do the inverse transform with integration in the complex plane will be this integral along the branch cut

[tex] \int_0^{\infty } \frac{e^{-t \rho } \left(e^{-i x \sqrt{\frac{\rho }{d}}}-e^{i x \sqrt{\frac{\rho }{d}}}\right)}{\beta +(-\alpha -\rho )^2} \, d\rho[/tex]

Maybe you can find that in a table somewhere. I couldn't even find/figure out the case where [itex]\beta = 0[/itex]. If you can't find the one with non-zero beta, but perhaps can find

[tex] \int_0^{\infty } \frac{e^{-t \rho } \left(e^{-i x \sqrt{\frac{\rho }{d}}}-e^{i x \sqrt{\frac{\rho }{d}}}\right)}{(-\alpha -\rho )^4} \, d\rho[/tex]

Then if beta is small maybe you could expand as a series in beta?
Groundwater
#3
Feb21-12, 07:40 AM
P: 2
Thanks...I think I am a bit out of my league here. I once took an introductory course in ODE that introduced me to Laplace transforms. I have slowly figured out how to apply the technique in the solution of PDE. However, at my level, I am doing well if I can solve a solution using tabulated inverse Laplace transforms.

I did try to solve the inverse Laplace, but I am sure that I didn't even get as far as you. I may have to adjust my boundary condition.

JJacquelin
#4
Feb21-12, 09:52 AM
P: 759
Inverse laplace transform for unique diffusion type problem

Hi !

May be, you could try this :
Split the function into two terms, thanks to
1/(((s-α)^2)+Bē) = c/(s-a -i B) - c/(s-a +i B)
c = 1/(2 i B)
The inverse Laplace of (1/(s-A))*exp(-C*sqrt(s)) can be found in tables.
But, there is a major difficulty, since A is complex : A=-a-i B or A=-a+i B
The formal result will involves Erfc functions in the complex range. The theoretical proof of validity should be rather ardous. Instead of doing it directly, I suggest to compute the Laplace transform of the result and check if it is consistent with the initial function F(s)

I suppose that beta>0 , so beta=Bē. If beta<0, then beta=-Bē and there is no difficulty since A=-a-B or A=-a+B which are reals.


Register to reply

Related Discussions
IVP Laplace Transform Problem - Tricky Inverse Laplace Transform Calculus & Beyond Homework 5
Problem with finding an Inverse Laplace Transform Differential Equations 2
Inverse Laplace transform problem Differential Equations 0
Finding an inverse Laplace Transform for a function - solving IVPs with Laplace Calculus & Beyond Homework 2
Laplace inverse transform Calculus 7