Ce document présente un compte rendu d'un TP sur la résolution d'équations aux dérivées partielles par la méthode des différences finies. Il décrit les méthodes pour les équations elliptiques et paraboliques avec différentes conditions aux limites et montre des exemples de résolution numérique.
0 évaluation0% ont trouvé ce document utile (0 vote)
34 vues7 pages
Ce document présente un compte rendu d'un TP sur la résolution d'équations aux dérivées partielles par la méthode des différences finies. Il décrit les méthodes pour les équations elliptiques et paraboliques avec différentes conditions aux limites et montre des exemples de résolution numérique.
Ce document présente un compte rendu d'un TP sur la résolution d'équations aux dérivées partielles par la méthode des différences finies. Il décrit les méthodes pour les équations elliptiques et paraboliques avec différentes conditions aux limites et montre des exemples de résolution numérique.
Ce document présente un compte rendu d'un TP sur la résolution d'équations aux dérivées partielles par la méthode des différences finies. Il décrit les méthodes pour les équations elliptiques et paraboliques avec différentes conditions aux limites et montre des exemples de résolution numérique.
Téléchargez comme DOCX, PDF, TXT ou lisez en ligne sur Scribd
Télécharger au format docx, pdf ou txt
Vous êtes sur la page 1/ 7
NOM PRENOM MATRICULE
LAMARI Omar 202031062952
HOCEINE Chaima 202031088115
COMPTE RENDU DU TP 3 : RESOLUTION DES EQUATIONS AUX
DERIVEES PARTIELLES
I. Le but du TP
L’objectif de ce TP est de résoudre les équations aux dérivées partielles du 2éme ordre linéaire en appliquant les méthodes numériques étudiées au cours. Les méthodes de résolution à programmées dans ce TP sont la méthode des différences finies.
II. Différences finies pour les équations elliptiques
II.1. Conditions aux limites Dirichlet
On applique sur une plaque les températures suivantes : coté supérieur100°C, côté droit 50°C, côté gauche 75°C, coté inférieur 0°C. La formule utilisée dans la méthode des différences finis pour les EDPs elliptiques :
clc clear all x=10; y=10; n=5; T=zeros(n,n); T(end,:)=150; T(:,1)=75; tol=1e-5; iter=0; nr=inf; while nr>tol Tpr=T for j=2:n-1 T(1,j)=0.25*(T(1,j+1)+T(1,j-1)+2*T(2,j)); T(1,n)=0.25*(2*T(1,n-1)+2*T(2,n)); for i=2:n-1 T(i,n)=0.25*(T(i+1,n)+T(i-1,n)+2*T(i,n-1)); T(i,j)=0.25*(Tpr(i+1,j)+Tpr(i-1,j)+Tpr(i,j+1)+Tpr(i,j-1)); end end iter=iter+1 nr=norm(Tpr-T) end xvec=linspace(0,x,n); yvec=linspace(0,y,n); [yy,xx]=meshgrid(yvec,xvec); mesh(xx,yy,T); xlabel('x'); ylabel('y'); zlabel('temperature');
Si le nombre de points de maillage est élevée, donc on pourra obtenir une solution plus précise, « cela influe sur la représentation de la géométrie obtenue. Ainsi qu’on pourra observer une meilleure convergence de la solution. Dans le cas de bord isolé inférieur : on remarque que les températures dans les différents points changent par rapport à la première méthode (bord inférieur =0)
III. Différences finies pour les équations paraboliques
Utiliser la méthode des différences finies pour résoudre la distribution de température à t=3s d'une tige longue et mince d'une longueur de 10 cm et les valeurs suivantes : Ax=1cm et Δt=0,1s. A t=0, la température de la tige est nulle et les conditions aux limites sont fixées à tout moment à T(0)=100°C et T(10)=75°C. Notez que la tige est en aluminium. Par conséquent, k=0,835 cm2/s La formule utilisée dans la méthode des différences finis pour les EDPs paraboliques :
i =T i + λ ( T i+1 −2T i +T i−1 ) Tels que : λ=k ∆ t /(∆ x)
T l+1 l l l l 2 Remarque
La première figure représente la distribution de température dans la tige de t=0 à t=3s avec un pas de 0,1s. La deuxième figure représente la distribution de température dans la tige à t =0s et à t=3. On remarque que la distribution de température s’équilibre au fil du temps.
IV. Conclusion
La précision de la température est directement proportionnelle à l’augmentation du nombre de
points de maillage (n) dans la méthode des différences finies.
