Equation Diffusion Thermique 2012
1. Équation de diffusion Soit une fonction u(x, t) représentant la température dans un problème de diffusion thermique, ou la concentration pour un problème de diffusion de particules. L'équation de diffusion est: où D est le coefficient de diffusion et s(x, t) représente une source, par exemple une source thermique provenant d'un phénomène de dissipation. On cherche une solution numérique de cette équation pour une fonction s(x, t) donnée, sur l'intervalle [0, 1], à partir de l'instant t=0. La condition initiale est u(x, 0). Sur les bords ( x=0 et x=1) la condition limite est soit de type Dirichlet: soit de type Neumann (dérivée imposée): 2. Méthode des différences finies 2. a. Méthode. Définitions Soit N le nombre de points dans l'intervalle [0, 1]. On définit le pas de x par On définit aussi le pas du temps. La discrétisation de u(x, t) est définie par: où j est un indice variant de 0 à N-1 et n un indice positif ou nul représentant le temps. Figure pleine page La discrétisation du terme de source est On pose 2. b. Schéma explicite Pour discrétiser l'équation de diffusion, on peut écrire la différence finie en utilisant les instants n et n+1 pour la dérivée temporelle, et la différence finie à l'instant n pour la dérivée spatiale: Avec ce schéma, on peut calculer les U j n+1 à l'instant n+1 connaissant tous les U j n à l'instant n, de manière explicite.
Equation Diffusion Thermique Rule
Supposons λ = 0. Il existe alors de même des constantes réelles B, C telles que X ( x) = Bx + C. Une fois encore, les conditions aux limites entraînent X nulle, et donc T nulle. Equation diffusion thermique rule. Il reste donc le cas λ > 0. Il existe alors des constantes réelles A, B, C telles que Les conditions aux limites imposent maintenant C = 0 et qu'il existe un entier positif n tel que On obtient ainsi une forme de la solution. Toutefois, l'équation étudiée est linéaire, donc toute combinaison linéaire de solutions est elle-même solution. Ainsi, la forme générale de la solution est donnée par La valeur de la condition initiale donne: On reconnait un développement en série de Fourier, ce qui donne la valeur des coefficients: Généralisation [ modifier | modifier le code] Une autre manière de retrouver ce résultat passe par l'application de théorème de Sturm-Liouville et la décomposition de la solution sur la base des solutions propres de la partie spatiale de l'opérateur différentiel sur un espace vérifiant les conditions aux bords.
On considère le cas simplifié de l'équation en une dimension, qui peut modéliser le comportement de la chaleur dans une tige. L'équation s'écrit alors: avec T = T ( x, t) pour x dans un intervalle [0, L], où L est la longueur de la tige, et t ≥ 0. On se donne une condition initiale: et des conditions aux limites, ici de type Dirichlet homogènes:. L'objectif est de trouver une solution non triviale de l'équation, ce qui exclut la solution nulle. Cours 9: Equation de convection-diffusion de la chaleur: Convection-diffusion thermique. On utilise alors la méthode de séparation des variables en supposant que la solution s'écrit comme le produit de deux fonctions indépendantes: Comme T est solution de l'équation aux dérivées partielles, on a: Deux fonctions égales et ne dépendant pas de la même variable sont nécessairement constantes, égales à une valeur notée ici −λ, soit: On vérifie que les conditions aux limites interdisent le cas λ ≤ 0 pour avoir des solutions non nulles: Supposons λ < 0. Il existe alors des constantes réelles B et C telles que. Or les conditions aux limites imposent X (0) = 0 = X ( L), soit B = 0 = C, et donc T est nulle.