Scribd
Importer
26 vues2 pages

Stabilité des schémas numériques

Ce document présente cinq exercices sur la modélisation numérique des équations aux dérivées partielles. Les exercices portent sur la stabilité des schémas numériques en fonction du pas de temps et du type de schéma, explicite ou implicite.

Transféré par

Akkad Ahmed
Copyright
© All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd
26 vues2 pages

Stabilité des schémas numériques

Ce document présente cinq exercices sur la modélisation numérique des équations aux dérivées partielles. Les exercices portent sur la stabilité des schémas numériques en fonction du pas de temps et du type de schéma, explicite ou implicite.

Transféré par

Akkad Ahmed
Copyright
© All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd

Modélisation numérique des produits et des structures

Travaux Dirigés 3: Temporel et stabilité

Exercice 1: Stabilité en temps ordre 1


On se propose de résoudre le problème scalaire instationnaire du premier ordre suivant :
du
 u  0, avec u(t=0)=1
dt
1. Calculer la solution exacte du problème.

2. Définir la notion de calcul explicite, implicite et Crank-Nicholson.

3. A l'aide d'un schéma explicite en temps et après avoir determiné un critère de

stabilité en temps pour t, calculer 5 solutions successives. On choisira un t = 0.1 s.

4. Même question que la précédente pour un schéma implicite en temps.

5. Développer un programme Matlab pour résoudre de l'équation instationnaire au

premier ordre

Exercice 2: Stabilité en temps ordre 2


On considère l'équation scalaire du second ordre en temps :

d2u du
 4 u  1, avec u(t=0)= (t=0)=0
dt 2 dt
1. Calculer la solution exacte du problème.

2. A l'aide d'un schéma explicite en temps et après avoir determiné un critère de stabilité

en temps pour t (d'après le concept de positivité), calculer 5 solutions successives. On

choisira un = 1/4s.

3. Même question que la précédente pour un schéma implicite en temps

4. Développer un programme Matlab pour résoudre de l'équation instationnaire au

deuxième ordre

1/2
Exercice 3: Stabilité et décentrement
Soit l'équation de convection suivante :
C(x, t) C(x, t)
 0
t x
Le schéma numérique pour obtenir la solution à l'instant n+1 est le suivant:
t n
C nj 1  C nj 
2x

C j1  C nj1 
1. Ce schéma est-il explicite ou implicite ?
2. Montrer que le schéma est toujours instable en utilisant le concept de positivité ?
3. Même question avec un décentrement droite puis gauche du second terme.

Exercice 4: Stabilité-viscosité numérique


Soit l'équation d'advection suivante :
C(x, t) C(x, t)
 0
t x
Le schéma numérique pour obtenir la solution à l'instant n+1 est le suivant:
t n  2 t 2 n
C n 1
j C 
n
j
2x

C j1  C j1 
n

2x 2 
C j1  2Cnj  Cnj1 
1. Ce schéma est-il explicite ou implicite ?
2. Etudier la positivité du schéma.
3. Etudier sa stabilité par la méthode de décomposition de Neumann.
4. Comparez les deux critères et concluez.

Exercice 5: Stabilité ordre 2-système


Soit le système d'équations différentielles à deux variables suivant:
 MÜ(t)  K U(t)  F,  t>0
avec
1 1  1 1 1  0
M  1 , K  1 , F    et U(t  0)   
U(t  0) 0

 2   3 2   
1. Déterminer le pas de temps critique donné pour un schéma aux différences finies centré
Tmin
par : t crit  où est la plus petite période de vibration de la structure.

2. Calculer les solutions correspondants aux instants ,2 et 3 après discrétisation
du terme de dérivée temporelle par un schéma aux différences finies centré. On prendra :
t  0.5t crit

2/2

Vous aimerez peut-être aussi