Outil mathématique

pour l’ingénieur

2021
Table des matières
Partie 1 : problème 1 ……………………………………………………………………………………………p3
1. Justification de l’existence
2. Solution numérique
3. La convergence de la solution approchée vers la solution exacte 
4. Cas a=0 et b=0
5. avantages et inconvénients du schéma 3

Partie 2 : problème 2……………………………………………………………………………………………p6

1. Présentation du problème
2. Schéma explicite

Conclusion ……………………………………………..……………………………………………………………p8

2
Problème 1

{
'
( β ( x ) u )'−γ ( x ) u ( x ) =f ( x ) 0< x<1 ,
u(0)=u 0 et α u(1)+b u '(1)=c ,

1. Justification de l’existence

On a la formule d’approximation à trois points d’une dérivée d’ordre 1 est la suivante :

g ( x i +h )−g ( x i−h )
g’( x i ¿ =
2h

Pour h/2 on a : g’( x i ¿ = ( h2 )−g( x − h2 )

g x i+ i

On pose : g( x i ¿ = β (x i) u '

Donc : g’( x i ¿ = ( h2 ) u '−β ( x − h2 )u '

β xi + i

h
Pour g( x i+ ), on approche u’ par un schéma avant à deux points :
2

( ) ( )( ) ( )
h h h h
h u x i + + −u xi − +
g xi + = h 2 2 2 2
2 β x i+
2 h

( h2 )=β ( x + 12 )
u ( xi +1 )−u ( x i )
Alors : g xi + i
h

h
Pour g ( x i− ), on approche u’ par un schéma arrière à deux points :
2

( h
) (
h h h h
u x i− + −u x i− − ( ) ( )
g x i− =
2 β x i−
h
2
2 2
h ) 2 2

Alors : (
g x i−
h
2 ) (
= β x i−
1 u ( x i )−u ( x i−1 )
2 ) h

Donc : g’( x i ¿ = ( 12 )(u ( x

β xi + i+1 (h
)
) −u ( xi ) )−β x i− 2 (u ( xi ) −u ( x i−1) )
2
h

On note : u ( x i )=ui , on trouve :

3
 β 1 (ui+1−u i)− β 1 (ui−u i−1)
' '
( β ( x i ) ui ) =¿g’( x i ¿ = i+
2
i−
2

h2

On remplace dans le problème à traiter, on trouve le schéma centré demandé :

β 1 (ui+1−u i)− β 1 (ui−u i−1)

i+ i−
2 2 − 𝛾( x i )U i = 𝑓( x i ).
h2

2. Solution numérique

On a U ext (x )=e− x ( x −1)2

Définition de l'intervalle [0,1], le nombre de points N=40 et le pas h := (1+0)/(N-1):

h=1/39

Les nœuds de calcul 

dom := 0, 1/39, 2/39, 1/13, 4/39, 5/39, 2/13, 7/39, 8/39, 3/13, 10/39, 11/39, 4/13, 1/3, 14/39, 5/13,
16/39, 17/39, 6/13, 19/39, 20/39, 7/13, 22/39, 23/39, 8/13, 25/39, 2/3, 9/13, 28/39, 29/39, 10/13,
31/39, 32/39, 11/13, 34/39, 35/39, 12/13, 37/39, 38/39, 1

Uapp := [1., 0.925353390687743, 0.855093655627851, 0.789014931244723, 0.726919767606324,

0.668618814121469, 0.613930515836559, 0.562680820431904, 0.514702895560824, 0.469836856194528,
0.427929501632615, 0.388834061865951, 0.352409952970511, 0.318522541230975, 0.287042915698669,
0.257847668895330, 0.230818685384126, 0.205842937940368, 0.182812291056582, 0.161623311532520,
0.142177085901004, 0.124379044527894, 0.108138791850987, 0.933699430150876e-1,
0.799899663543362e-1, 0.679200316029123e-1, 0.570848636359441e-1, 0.474126015468880e-1,
0.388346628743038e-1, 0.312856128009013e-1, 0.247030381500100e-1, 0.190274260124456e-1,
0.142020468432284e-1, 0.101728418712723e-1, 0.688831467250827e-2, 0.429942676034852e-2,
0.235949705374573e-2, 0.102410508703159e-2, 0.250997831383061e-3, 0]

3. La convergence de la solution approchée vers la solution exacte 

Convergence de la solution approchée (graphiquement)

Pour N=40, on trace la courbe de Uapp, Uext et de la solution exacte et approchée

4
Tableau de convergence

N Erreur
40 0.0002769750
140 0.000041118875
180 0.000028186691
200 0.0000240775

4. Cas a=0 et b=0

Si a=0 la condition au limite b.U’(1)=c sera de type Neumann
Si b=0 la condition au limite a.U(1)=c sera de type Dirichlet et le problème sera plus facile à résoudre
5. avantages et inconvénients du schéma 3

𝛽(𝑥)𝑢 ′′(𝑥) + 𝛽 ′ (𝑥)𝑢 ′ (𝑥) − 𝛾(𝑥)𝑢(𝑥) = 𝑓(𝑥)

En approchant les dérivées, on trouve le schéma suivant :

βi (U i−1 −2U i +U i +1) ( βi +1−β i−1) .(U i+1−U i−1 )

+ −¿ 𝛾( x i )U i = 𝑓( x i )
h2 4 h2

Avantage :

Schéma facile à établir par rapport au premier schéma

Inconvénients :

Moins de précision de calcul

5
Problème 2
1. Présentation du problème

L’équation de convection–diffusion–réaction exprime la conservation de la concentration u d’une

certaine substance dont les variations en espace et en temps résultent de processus de diffusion, convection
et réaction, par exemple dans un milieu poreux. Elle prend la forme intégrale suivante :
t
∂u
( ⃗x , t )−∆u ( ⃗x ,t ) + b⃗ . ∇ . ( u ( ⃗x ,t ) ) +cu ( ⃗x ,t )−∫ ∆ u ( ⃗x , s ) ds=f ( ⃗x , t )( ⃗x ,t ) ϵ Ωx ¿ 0 ,T ¿ ¿
∂t 0

u ( ⃗x , t )=0 ( ⃗x ,t ) ϵ ∂ Ωx ¿ 0 ,T ¿ ¿

u ( ⃗x , 0 ) =u0 ( ⃗x ) ⃗x ϵ Ω

*On suppose que u0 ( ⃗x )=0

L’objectif principal de cet exercice est de résoudre numériquement le problème intégro-différentiel et

dynamique en deux dimensions d’espace moyennant la méthode des différences finies.

2. Schéma explicite

On approche la dérivée d’ordre 1 par un schéma centré

∂u un+1 n
i , j −ui , j
( ⃗x , t )=
∂t ∆t

On approche l’opérateur gradient par un schéma décentré arrière:

n n
∂u u −u
( ⃗x , t ) = i , j i−1 , j
∂x 2h
n n
∂u u −u
( ⃗x , t )= i , j i , j−1
∂y 2h

On approche l’opérateur laplacien par le schéma centré :

2 n n n
∂ u ( ⃗x ,t ) ui−1 , j −2u i, j +ui +1 , j
=
∂ x2 h2
2 n n n
∂ u ( ⃗x ,t ) ui , j−1 −2u i, j +ui , j +1
2
= 2
∂y h

On approche l’intégrale par la méthode des trapèzes simple

6
 t t 2 t 2
∂ u ( ⃗x , s ) ∂ u ( ⃗x , s )
∫ ∆ u ( ⃗x , s ) ds=∫ ∂x 2
ds +∫
∂ y2
ds
0 0 0

t 2 0 0 0 n n n
∂ u ( ⃗x , s ) u −2 u +u u −2 u + u
∫ ∂ x 2 ds= 2t ( i−1 , j h2i , j i+1 , j + i−1 , j h2i , j i+1 , j )
0

( )
t n n n
∂2 u ( ⃗x , s ) t ui−1 , j−2 ui , j +ui+1 , j
∫ ∂ x2 ds=
2 h
2
0

Et
t 2 0 0 0 n n n
∂ u ( ⃗x , s ) u −2 u +u u −2 u + u
∫ ∂ y 2 ds= 2t ( i , j−1 h2i , j i , j+1 + i , j−1 h2i , j i , j+1 )
0

t 2 n n n
∂ u ( ⃗x , s ) u −2 u +u
∫ ∂ y 2 ds= 2t ( i , j−1 h2i , j i , j+1 )
0

Ainsi l’équation s’écrit comme suit

un+1
(
n
i , j −u i , j un −2 uni , j+ uni+1 , j u ni , j−1−2uni , j +uni , j+1 uni , j−uni−1 , j uni , j−uin, j−1 n
n n n n
t u i−1, j −2ui , j +ui +1 , j u i , j −1
− i−1 , j − +b x +b y + c U i, j − +
∆t h2 h2 2h 2h 2 h2

On cherche la forme matricielle :

(
uni ,+1j =∆ t × f ni, j +U ni , j 1−
4∆t
h2
∆t ∆t 2t∆t
h
∆t
h )
∆t t ∆t
( ∆t t∆t
h 2h
∆t
)
−b x −b y − 2 + c +U ni−1 , j 2 +b x + 2 +U ni+1 , j 2 + 2 +U ni , j−1 2 +b y
2h 2h 2h 2h h
∆
2 ( ) (
Ainsi le schéma explicite du problème s’écrit sous forme : 𝑈 𝑛+1 = 𝐴𝑈 𝑛 + 𝐹 n

7
Conclusion

la méthode des différences finies est une technique courante de recherche de solutions approchées
d'équations aux dérivées partielles qui consiste à résoudre un système de relations (schéma numérique) liant
les valeurs des fonctions inconnues en certains points suffisamment proches les uns des autres.

Grace à ce projet nous avons pu résoudre numériquement l’équation problème 1suivant les étapes de
la méthode des différences finies puis visualiser la solution obtenue.