Rank: 3796 based on 26 downloads (last 30 days) and 1 file submitted
photo

Wahyu Srigutomo

E-mail

Personal Profile:
Professional Interests:

 

Watch this Author's files

 

Files Posted by Wahyu
Updated   File Tags Downloads
(last 30 days)
Comments Rating
15 Feb 2006 Gaver-Stehfest algorithm for inverse Laplace transform Gaver-Stehfest algorithm for inverse Laplace transform Author: Wahyu Srigutomo optimization, inverse laplace trans..., gaverstehfest, algorithm 26 10
  • 3.375
3.4 | 8 ratings
Comments and Ratings on Wahyu's Files View all
Updated File Comment by Comments Rating
12 Apr 2011 Gaver-Stehfest algorithm for inverse Laplace transform Gaver-Stehfest algorithm for inverse Laplace transform Author: Wahyu Srigutomo Saskia

I am trying to solve the heat equation with both conduction and advective terms:

alpha*(d^2)T/(d^2)z-v*dT/dz-dT/dt (1)

where v=u*Cw/Cv, alpha is a constant. According to the literature the temperature frequency cariation solution of equation 1 is:

That=T0*exp(gamma*z)*exp(i*omega*t) (2)

where T0 and omega are constants. To get the actual temperature (T) value I am trying to use the Gaver-Stehfest numerical inversion. However, the values I get are not realistic, while my input parameters are fine. What am I doing wrong?

My input function is:

function f=fun_s(p)

%Define constants u=3.08*10^-8; %ms^-1 Cw=4.18*10^6; %Jm^-3K^-1 Cv=0.7*10^6; %Jm^-3K^-1 k=0.3; %Wm^-1K^-1 L=1; %m alpha=0.4*10^-6; %m^2 s^-1 omega=(2*pi())/(365.25*86400); theta0=11.5112; z=1;

v=u*Cw/Cv; gamma=(v-sqrt(v^2+4*i*omega*alpha))/(2*alpha);

f=theta0*(exp(gamma*z))*(exp(i*omega*p));

which I then run with the following to get a daily temperature value:

t=1:2*pi/365:2*pi; l=length(t);

for i=1:l T(i)=gavsteh('fun_s',t(i),18); end

Any help on this will be greatly appreciated!

28 May 2010 Gaver-Stehfest algorithm for inverse Laplace transform Gaver-Stehfest algorithm for inverse Laplace transform Author: Wahyu Srigutomo Peng

this one totally helps a lot

09 Jan 2010 Gaver-Stehfest algorithm for inverse Laplace transform Gaver-Stehfest algorithm for inverse Laplace transform Author: Wahyu Srigutomo Keramat, Alireza

Thank you. It was interesting because of its simplicity.

26 Aug 2009 Gaver-Stehfest algorithm for inverse Laplace transform Gaver-Stehfest algorithm for inverse Laplace transform Author: Wahyu Srigutomo Mark

Worked great for me. Nice job.

26 Feb 2009 Gaver-Stehfest algorithm for inverse Laplace transform Gaver-Stehfest algorithm for inverse Laplace transform Author: Wahyu Srigutomo Maxime, Montaru

You need to use multi-precision toolbox to enhance results. Very good results for one staircase simulation with a big L. Beware when you simulate multiple staircase.

Contact us