Solution à télcharger. Dans ce fichier, on reprend l'étude des éléments finis P1 pour deux équations et on détaille le calcul des matrices de masse et de rigidité, ainsi que la matrice obtenue à partir d'un terme de convection.
Réponse à la question f) : programme en Scilab et programme en Python.
Réponse à la question h) : programme en Scilab et programme en Python.
Chapitre 13 - Équations aux dérivées partielles
Solutions d'exercices
par Magali Ribot, maître de conférences en Mathématiques Appliquées
-
Exercice 13-2 : Eléments finis pour différentes équations
-
Exercice 13-3 : Eléments finis et conditions aux bords pour l'équation de Poisson
Solution à télécharger. Dans cet exercice, on reprend l'étude faite pour l'équation de Poisson avec conditions aux bords de Dirichlet pour l'étendre à différentes conditions aux bords.
-
Exercice 13-4 : Conditions de stabilité de différents schémas
Solution à télécharger. Vous trouverez dans ce fichier les solutions de l'exercice 13-4 sans le détail des calculs. Pour chaque schéma, la formule du coefficient d'amortissement est donnée, ainsi que la condition de stabilité du schéma (s'il y en a). Les calculs intermédiaires sont laissés au lecteur.
-
Exercice 13-7 : Résolution de l'équation des ondes
Solution à télécharger.
L'objet de cet exercice est la programmation de l'équation des ondes à l'aide des différences finies. Le cas des conditions aux bords de Dirichlet (question a) est traité à l'aide d’un programme en Scilab ou d’un programme en Python. Une animation vous permet de visualiser le résultat.
Pour des conditions aux bords périodiques (question b), le problème peut être traité à l'aide d’un programme en Scilab ou d’un programme en Python.
Projets
Propagation d’un paquet d’ondes quantiqueComment se comporte un paquet d’ondes lorsqu’il se heurte à une barrière de potentiel ? Retrouvez les résultats de l’animation ci-dessous en réalisant la simulation décrite.
Propagation de l’influx nerveux dans un neurone
Comment se propage l'influx nerveux dans un neurone ? La membrane de l'axone, sensible aux potentiels chimiques, présente une différence de potentiel qui permet la propagation d'un signal électrique. Dans ce projet, vous programmerez le modèle de FitzHugh-Nagumo, une version simplifiée du modèle de Hodgkin-Huxley, modèle historique de propagation de l'influx nerveux. Vous mettrez en évidence un « effet de seuil », c'est-à-dire qu'il faut une valeur suffisamment grande de l'impulsion électrique à un bout de l'axone pour voir apparaître la propagation d'un signal électrique le long de l'axone.
Formation de motifs en biologie
Comment se forment les taches d’un léopard ou les rayures d’un zèbre ? Lors de la formation de motifs sur le pelage d'un animal, plusieurs substances chimiques entrent en jeu. Dans ce projet, vous reproduirez ces motifs par la résolution d'un système d'équations de type réaction-diffusion modélisant le comportement d'un système de substances chimiques de type activateur-inhibiteur.
Compléments
Traduction en Python d'exemples donnés en Scilab dans le livre :- Listing 13.1 : Équation de Poisson 2D avec conditions de Dirichlet ; résolution par différences finies
Programme à télécharger - Listing 13.2 : Équation de Poisson 1D avec conditions de Dirichlet ; résolution par éléments finis
Programme à télécharger - Listing 13.3 : Équation de la chaleur avec conditions de Dirichlet ; résolution par différences finies
Programme à télécharger - Listing 13.4 : Équation de la chaleur avec conditions de Dirichlet ; résolution par éléments finis
Programme à télécharger - Listing 13.5 : Équation de la chaleur avec conditions de Dirichlet ; résolution par méthodes de splitting
Programme à télécharger - Listing 13.6 : Équation de transport ; résolution par différences finies
Programme à télécharger - Listing 13.7 : Équation des ondes ; résolution par différences finies
Programme à télécharger
Liens
Freefem++Ce logiciel gratuit permet une initiation très commode à la pratique de la méthode des éléments finis.
http://www.freefem.org