Méthode des caractéristiques

Un article de Wikipédia, l'encyclopédie libre.

Cet article est une ébauche concernant les mathématiques.

Vous pouvez partager vos connaissances en l’améliorant. (Comment ?).

En mathématiques, la méthode des caractéristiques est une technique permettant de résoudre les équations aux dérivées partielles. Particulièrement adaptée aux problèmes de transport, elle est utilisée dans de nombreux domaines tels que la mécanique des fluides ou le transport de particules.

Dans certains cas particuliers, la méthode des caractéristiques peut permettre la résolution purement analytique de l'EDP. Dans les cas plus complexes (rencontrés par exemple en modélisation des systèmes physiques), la méthode des caractéristiques peut être utilisée comme une méthode de résolution numérique du problème.

Sommaire

1 Méthode analytique
- 1.1 Principe général
- 1.2 Exemple : Résolution analytique de l'équation de transport
2 Méthode numérique
- 2.1 Principe général
- 2.2 Exemple

[modifier] Méthode analytique

[modifier] Principe général

Pour une équation aux dérivées partielles (EDP) du premier ordre, la méthode des caractéristiques cherche des courbes (appelées « lignes caractéristiques », ou plus simplement « caractéristiques ») le long desquelles l'EDP se réduit à une simple équation différentielle ordinaire (EDO). La résolution de l'EDO le long d'une caractéristique permet de retrouver la solution du problème original.

[modifier] Exemple : Résolution analytique de l'équation de transport

On cherche une fonction $u$ :

$u: (x,t)\in\mathbb{R}\times\mathbb{R}^+ \mapsto u(x,t)\in\mathbb{R}$

solution du problème suivant:

$\left\{\begin{array}{l l} \displaystyle\frac{\partial u}{\partial t} + c\frac{\partial u}{\partial x} = 0, \quad & (x,t)\in\mathbb{R}\times\mathbb{R}^+\\ u(x,0) = u_0(x), & x\in\mathbb{R} \end{array}\right.$

dans laquelle $\textstyle c$ est une constante.

On cherche une ligne caractéristique $\textstyle\left(x(s), t(s)\right)$ le long de laquelle cette EDP du premier ordre se réduirait à une EDO. Calculons la dérivée de $u$ le long d'une telle courbe :

$\frac{du}{ds} = \frac{dt}{ds}\,\frac{\partial u}{\partial t} + \frac{dx}{ds}\,\frac{\partial u}{\partial x}$

On remarquera aisément qu'en imposant $\textstyle\frac{dt}{ds} = 1$ et $\textstyle\frac{dx}{ds} = c$ , on obtient : $\textstyle\frac{du}{ds} = \frac{\partial u}{\partial t} + c\frac{\partial u}{\partial x} = 0$ . La solution de l'équation reste donc constante le long de la ligne caractéristique.

Nous nous retrouvons donc avec trois équations différentielles ordinaires à résoudre :

$\frac{dt}{ds} = 1$ : en posant $\textstyle t(0)=0$ , on obtient :

$\forall s\in\mathbb{R}^+, \quad t(s)=s$

$\frac{dx}{ds} = c$ : en notant $\textstyle x(0)=x_0$ , on obtient :

$\forall s\in\mathbb{R}^+, \quad x(s)=x_0+c\,s=x_0+c\,t$

$\frac{du}{ds} = 0$ :

$\forall s\in\mathbb{R}^+, \quad u(s) = u(0) = u(x_0,0) = u_0(x_0)$

Dans ce cas, les lignes caractéristiques sont donc des droites de pente $\textstyle c$ , le long desquelles la solution reste constante. La valeur de la solution en un point $\textstyle (x,t)\in\mathbb{R}\times\mathbb{R}^+$ peut donc être retrouvée en cherchant la valeur de la condition initiale $\textstyle u_0$ à l'origine $\textstyle x_0=x-c\,t$ de la ligne caractéristique: