Scribd
Importer
80 vues2 pages

Problème de Diffusion-Convection-Réaction

Cet examen porte sur la résolution numérique d'équations aux dérivées partielles par la méthode des différences finies. Il comporte deux exercices, le premier concerne un problème de diffusion-convection-réaction en une dimension, et le second un problème de conduction de chaleur en deux dimensions.

Transféré par

stevi takang
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 DOCX, PDF, TXT ou lisez en ligne sur Scribd
80 vues2 pages

Problème de Diffusion-Convection-Réaction

Cet examen porte sur la résolution numérique d'équations aux dérivées partielles par la méthode des différences finies. Il comporte deux exercices, le premier concerne un problème de diffusion-convection-réaction en une dimension, et le second un problème de conduction de chaleur en deux dimensions.

Transféré par

stevi takang
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 DOCX, PDF, TXT ou lisez en ligne sur Scribd

EXAMEN DU SECOND SEMESTRE

FILIERE : PREPA NIVEAU : II AVRIL 2023


MATIERE : APPLICATION NUMERIQUE II
ENSEIGNANT: DR KINFACK
DUREE: 2H
JEUTSA/Dr DONFACK HUBERT
2022-2023 Logiciel autorisées : Calculatrice autorisée
OUI  NON 

NB : Cette épreuve est théorique. Les ordinateurs sont interdits


Exercice I (15pts)
A) Description du problème
On considère le problème de Diffusion-Convection-Réaction en dimension 1D suivant:
( P) ¿

On considère un maillage régulier du domaine défini par les nœuds

où un entier strictement positif donné et le pas du maillage.


Dans la suite on suppose que la fonction inconnue est quatre fois continûment dérivable sur

.
B) Méthode de Différences Finies (1+1+2+2+2=8 points)

1) Montrer que avec


u ( x +h ) −u(x−h) h2 ' ''
2) Montrer que u' ( x ) = − u ( θ x ) avec θ x ∈ ¿ x−h , x+ h ¿
2h 6
3) Déterminer le problème discret associé au problème (P) en vous servant des
Différences Finies Centrées.
4) Ecrire le système linéaire obtenu en 3) sous forme la matricielle
Ah U h=F h
5) Montrer que le problème discret PD admet une unique solution.
C) Programmation Matlab (2+2+3=7 points)
6) Créer une fonction Matlab qui retourne la matrice Ah ci-dessus ayant en entrée N.
7) Ecrire une fonction Matlab qui retourne le Vecteur F h ci-dessus ayant en entrée N.
8) Ecrire un programme Matlab qui trace la solution approchée et la solution exacte en
appelant les du fonctions ci-dessus.
1
On prendra : N=20; ue ( x ) =sin ⁡(πx ) , a=−1 ; b= ; c=1; α =0 ; β=0
100
Exercice II (5 points)

Page 1 sur 2
Si la température initiale dans la coupe transversale à x est u(x ) et les extrémités de la barre
sont maintenues à la température zéro, le problème pour la fonction de distribution de la
température est
¿
On considère les maillages suivants :
t n=nτ pour n=0 , 1, ⋯ , M +1.
x i=i hpour i=0 , 1 , ⋯ , N +1.
L L
Où τ = et h= avec N et M deux entiers naturels
M +1 N +1
1) Proposer un problème discret associé au problème continu ci-dessus en vous
servant des différences finies progressive sur la dérivée partielle par rapport à t
et les différences finies centrées sur la dérivée partielle par rapport à x.
2) En déduire la forme matricielle du problème discret en 1) sous la forme
n+1 n n
F = A F +B

Page 2 sur 2

Vous aimerez peut-être aussi