Annexes

A.2.3 : Transformation de Laplace inverse

Méthode analytique
∞
La transformée de Laplace θ(p) de la fonction T(t) est donnée par : L[T(t )] = θ(p ) = ∫ exp(− p t ) T(t ) dt
0

Il n’existe pas de formule analytique générale permettant de calculer T(t) connaissant θ(p). On connait
cependant l’expression exacte de T(t) pour certaines fonctions particulières θ(p), on en trouvera des exemples
page suivante (cf. Spiegel pour des tables plus complètes). L’utilisation de ces tables associée aux propriétés
particulières de la transformation de Laplace inverse rappelées en annexe A.2.2 peut permettre de résoudre un
certain nombre de cas. On essaiera toujours de décomposer une fonction complexe en somme, produit, série…
de fonctions simples plus facilement inversibles.

Méthodes numériques
Pour les cas de figure pour lesquels on ne peut pas trouver une solution analytique, on peut employer l’une des
deux méthodes numériques suivantes :
Méthode de Stehfest
La transformée inverse de la fonction θ(p) peut se calculer par :
N
ln(2)  j ln(2) 
T(t) =
t ∑ V θ 
j=1
j i
t


N = 20 (double précision) :
V1 = -5,511463844797178.10-6 V2 = 1,523864638447972.10-1 V3 = -1,174654761904762.102
4 5
V4 = 1,734244933862434.10 V5 = -9,228069289021164.10 V6 = 2,37740877871031810.7
V7 = -3,494211661953704.108 V8 = 3,241369852231879.109 V9 = -2,027694830723779.1010
V10 = 8,946482982379724.1010 V11 = -2,870209211471027.1011 V12 = 6,829920102815115.1011
V13 = -1,219082330054374.1012 V14 = 1,637573800842013.1012 V15 = -1,647177486836117.1012
V16= 1,221924554444226.1012 V17 = -6,488065588175326.1011 V18 = 2,333166532137059.1011
V19 = -5,091380070546738.1010 V20 = 5,091380070546738.109

N = 10 (simple précision):

V1 = 1/12 V2 = -385/12 V3 = 1279 V4 = -46871/3 V5 = 505465/6 V6 = -473915/2

V7 = 1127735/3 V8 = -1020215/3 V9 = 328125/2 V10 = -65625/2

Méthode de Fourier
exp(c t )  θ(c ) ∞ 
T(t ) =  + ∑ (Re[θ(c + j ω k )] cos(ω k t ) − Im[θ(c + j ω k )]sin (ω k t ))
t max  2 k =1 
kπ
Avec ωk =
t max

La somme infinie est dans la pratique calculée pour un nombre de fini N de termes , on prendra en général
N > 100. Cette méthode nécessite de choisir deux paramètres : c et tmax . On doit s’assurer a posteriori que
exp(-2 c tmax) T(2 tmax) ≈ 0.
Choix d’une méthode et vérification des résultats
La méthode de Stehfest est plus simple à mettre en oeuvre car elle ne nécessite pas de choisir certains
paramètres. La méthode de Fourier peut conduire à un meilleur résultat dans le cas d’inversion de certaines
fonctions comme les fonctions périodiques par exemple (cf. Maillet).
L’étude du comportement de la fonction θ(p) aux temps longs (t→∞ soit p→0) et aux temps courts (t→0 soit
p→∞) peut conduire à des formules approchées de θ(p) dont on peut alors trouver la transformée de Laplace
inverse analytiquement. La comparaison de ces solutions analytiques avec les résultats de l’inversion numérique
donne une indication sur la justesse de l’inversion numérique.

Yves Jannot 113

 Transferts et échangeurs de chaleur 114

A.2.3 : Transformation de Laplace inverse

p
q=
a

θ(p ) = L{T(t )} T(t) θ(p ) = L{T(t )} T(t)

1 ln(p )
1 − ln(t ) − γ ; γ = 0,57721
p p
1 1
1 δ(t) Dirac
p πt
1 1 2
e −βt t
p+β p p π
ω ω
sin (ωt ) sh (ωt )
2
p +ω 2
p − ω2
2

p p
cos(ωt ) ch (ωt )
p +ω2 2
p 2 − ω2

( ) ( ) t n −1
b 1
1 − exp b 2 t erfc b t n = 1, 2, 3...
(
p b+ p ) p n
(n − 1)!

θ(p ) = L{T(t )} T(t)

x  x2 
e −q x exp − 

2 π α t3  4α t 
1
e −q x α  2  x2 
  exp − 
q  4α t 
πt  
e −q x  x 
erfc 

p  2 αt 
1
e −q x  αt   x2   
 - x erfc x 
2
2  exp -  
pq  π  4αt 
   2 αt 
1
e −q x  x2     t   x2 
 erfc x
2
t + -x  exp - 
p 2  2α     πα   4αt 
   2 αt   
1
e −q x
q+h
α 
 
2  x2 
exp - (
 
 − h α exp hx + α t h 2 erfc x + h αt 
 4αα   )
 πt     2 αt 
e −q x
q (q + h )
(  x
α exp hx + α t h 2 erfc ) 
+ h αt 
 2 αt 
e −q x
p (q + h )
1
h
 x  1
erfc
 x
 − exp hx + α t h 2 erfc
  ( 
+ h αt  )
 2 αt  h  2 αt 
1
e −q x
p q (q + h )
2α
 
h  πt 
2  x 2  1+ h x
exp -
 4αα
−
2
 x
erfc
 1
+ ( 2  x

 h 2 exp hx + α t h erfc ) 
+ h αt 
  h  2 αt   2 αt 
e −q x  α3 t 
−2h 
 π 
 2
 2α t 
( 


 exp - x  + α 1 + h x + 2 h 2 α t exp h x + α t h 2 erfc x + h αt 
 ) ( )
(q + h )2      2 αt 

114