Scribd
Importer
11 vues3 pages

Simulation de l'Équation des Ondes 2D

Le projet vise à résoudre numériquement l'équation des ondes en 2D sur un domaine défini, en utilisant la méthode des volumes finis et des conditions aux limites de Neumann ou périodiques. Il comprend une partie théorique sur la discrétisation de l'équation et la gestion des conditions aux limites, ainsi qu'une mise en œuvre pratique pour écrire un programme qui simule l'évolution de l'onde. Les résultats seront visualisés à l'aide d'isocontours.

Transféré par

ameneallal
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
11 vues3 pages

Simulation de l'Équation des Ondes 2D

Le projet vise à résoudre numériquement l'équation des ondes en 2D sur un domaine défini, en utilisant la méthode des volumes finis et des conditions aux limites de Neumann ou périodiques. Il comprend une partie théorique sur la discrétisation de l'équation et la gestion des conditions aux limites, ainsi qu'une mise en œuvre pratique pour écrire un programme qui simule l'évolution de l'onde. Les résultats seront visualisés à l'aide d'isocontours.

Transféré par

ameneallal
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

Université de Poitiers

Master Gestion de l'Energie


Année 2023

Simulation Numérique Multiphysique


Projet n°3
Equation des ondes

On veut résoudre l'équation des ondes :

 2
 c 2  (1) sur le domaine  :  0,1  0,1
t 2

c est la célérité de l’onde.

On se propose de résoudre numériquement l'équation (1) avec la méthode des volumes finis.
On se place en 2D dans la configuration Cell centred.

En termes de conditions aux limites on utilisera soit des conditions de Neumann, soit des
conditions périodiques.


0
y

 
0 0
x x


j
0 
i 
0
y

Conditions aux limites

On adoptera les notations suivantes :


Le nombre de volumes de control dans la direction x est NICV, et dans la direction y: NJCV.
On numérote toujours à partir de 0. Et nous avons les relations suivantes :
NI =NICV+1
NIM= NI-1

NJ =NJCV+1
NJM= NJ-1
Pour la discrétisation spatiale de la dérivée première, nous prendrons : Un schéma centré
d'ordre 2.

n
 2   pn1  2 pn   pn 1
On utilisera le schéma suivant pour la dérivée temporelle : 2  
t  P t 2

I Partie théorique :
A) Discrétisez l'équation (1) pour obtenir l'équation générique
APP  AEE  AW W  ANN  ASS  SMB

Donnez les coefficients AE, AW, AN, AS, AP.


Donnez le second membre SMB.

B) Gérez la prise en compte des conditions aux limites de Dirichlet et de Neumann.

II Mise en œuvre et écriture du programme

A) Nous utiliserons les tableaux suivants :

Métrique :

xc(i), yc(j) : les coordonnées du centre du volume de contrôle (i,j).


dx , dy : les dimension des volumes de contrôle (i,j)
Se, Sw, Sn, Ss les surfaces du volume de contrôle (i,j)
Vp: le volume du volume de contrôle (i,j).
DXw(i), DXe(i), DYs(j), DYn(j) : les distances.

Coefficients :

GAM(i,j) : coefficient Gama.


de(i,j), dw(I,j), dn(i,j), ds(i,j) : coefficients diffusifs.
AE(i,j), AW(i,j), AN(i,j), AS(i,j), AP(i,j), SMB(i,j) : coefficients de l’équations générique.
Paramètres du programme

LX, LY la taille du domaine (ce sont des nombres réels.


NICV, NJCV le nb de volumes de contrôle dans chaque direction.
NTM : le nombre maximal de pas de temps
dt : la valeur du pas de temps
c : la valeur de la célérité de l’onde.

Conditions initiales

Nous testerons différentes conditions initiales qui seront fournies.

Organisation du travail

1) Ecrire le programme qui permet de résoudre l’équation des ondes


avec des conditions aux limites de Neumann ou des conditions aux limites périodiques.

Test du programme

2) On visualisera l’évolution de  en traçant les isocontours.

Vous aimerez peut-être aussi