Équation de langevin

Équation de Langevin

L'équation de Langevin (1908) est une équation stochastique pour le mouvement brownien.

Sommaire

1 Théorie de Langevin du mouvement brownien
2 Équation de Langevin
3 Solution de Langevin (1908)
4 Solution moderne
5 Articles connexes
6 Bibliographie
7 Notes et références

Théorie de Langevin du mouvement brownien

Dans l'approche théorique de Langevin, une grosse particule brownienne de masse m, supposée animée à l'instant t d'une vitesse $\vec{v}(t)$ , est soumise à deux forces bien distinctes :

une force de frottement fluide du type $\vec{f} \, = \, - \, k \, \vec{v}$ , où k est une constante positive. Dans le cas d'une particule sphérique de rayon a, cette constante s'écrit explicitement : $k = 6πη a$ (loi de Stokes).

une force complémentaire, notée $\vec{\eta}(t)$ , qui synthétise la résultante des chocs aléatoires des molécules de fluide environnantes. Langevin écrit à propos de cette force supplémentaire qu' « elle est indifféremment positive et négative, et sa grandeur est telle qu'elle maintient l'agitation de la particule que, sans elle, la résistance visqueuse finirait par arrêter ». On appelle au XXI^e siècle une telle force un bruit blanc gaussien^[1].

Équation de Langevin

On applique le principe fondamental de la dynamique de Newton, ce qui conduit à l'équation stochastique de Langevin :

$m \, \frac{d\vec{v}(t)}{dt} \ = \ - \, k \, \vec{v}(t) \ + \ \vec{\eta}(t)$

Solution de Langevin (1908)

Réécriture de l'équation de Langevin

Prenons le produit scalaire de cette équation avec le vecteur position $\vec{r}(t)$ (en omettant la dépendance en temps pour simplifier les notations) :

$m \ \vec{r} \cdot \frac{d^2\vec{r}}{dt^2} \ = \ - \, k \ \vec{r} \cdot \frac{d\vec{r}}{dt} \ + \ \vec{r} \cdot \vec{\eta}$

Remarquons alors d'une part que :

$\frac{ d || \vec{r}||^2 }{ dt } \ = \ \frac{ d ( \vec{r} \cdot \vec{r} ) }{ dt } \ = \ 2 \ \vec{r} \cdot \frac{ d \vec{r} }{ dt }$

$\Longrightarrow \quad \vec{r} \cdot \frac{d\vec{r}}{dt} \ = \ \frac{1}{2} \ \frac{d ||\vec{r}||^2}{dt}$

et d'autre part que :

$\frac{ d^2 ||\vec{r}||^2 }{ dt^2 } \ = \ \frac{d~}{dt} \ \left[ \frac{ d ||\vec{r}||^2}{dt} \right] \ = \ \frac{d~}{dt} \ \left[ \, 2 \ \vec{r} \cdot \frac{d\vec{r}}{dt} \, \right] \ = \ 2 \, ||\vec{v}||^2 \ + \ 2 \ \vec{r} \cdot \frac{d^2\vec{r}}{dt^2}$

$\Longrightarrow \quad \vec{r} \cdot \frac{d^2\vec{r}}{dt^2} \ = \ \frac{1}{2} \ \frac{d^2 ||\vec{r}||^2}{dt^2} \ - \ ||\vec{v}||^2$

En substituant ces expressions dans le produit scalaire obtenu à partir de l'équation de Langevin, on obtient une nouvelle forme de l'équation différentielle :

$m \ \frac{d^2 ||\vec{r}||^2}{dt^2} \ = \ - \, k \ \frac{d ||\vec{r}||^2}{dt} \ + \ 2 \, m \ ||\vec{v}||^2 \ + \ 2 \, \vec{r} \cdot \vec{\eta}$

Moyenne sur le bruit blanc gaussien

On prend alors la moyenne de l'équation précédente sur toutes les réalisations possibles du bruit blanc gaussien. Il vient :

$m \ \langle \ \frac{d^2 ||\vec{r}||^2}{dt^2} \ \rangle \ = \ - \, k \ \langle \ \frac{d ||\vec{r}||^2}{dt} \ \rangle \ + \ 2 \, m \ \langle \ ||\vec{v}||^2 \ \rangle \ + \ 2 \, \langle \ \vec{r} \cdot \vec{\eta} \ \rangle$

On fait l'hypothèse avec Langevin que la valeur moyenne du terme de bruit est nulle^[2] :

$\langle \ \vec{r} \cdot \vec{\eta} \ \rangle \ = \ 0$

Par ailleurs, le processus de moyenne sur le bruit commute avec la dérivation temporelle :

$\langle \ \frac{d ||\vec{r}||^2}{dt} \ \rangle \ = \ \frac{d~}{dt} \ \langle \ ||\vec{r}||^2 \ \rangle \quad \mathrm{et} \quad \langle \ \frac{d^2 ||\vec{r}||^2}{dt^2} \ \rangle \ = \ \frac{d^2~}{dt^2} \ \langle \ ||\vec{r}||^2\ \rangle$

ce qui conduit à l'équation différentielle pour les moyennes :

$m \ \frac{d^2~}{dt^2} \ \langle \ ||\vec{r}||^2\ \rangle \ = \ - \, k \ \frac{d~}{dt} \ \langle \ ||\vec{r}||^2 \ \rangle \ + \ 2 \, m \ \langle \ ||\vec{v}||^2 \ \rangle$

On pose alors :

$u(t) \ = \ \frac{1}{2} \ \frac{d~}{dt} \ \langle \ ||\vec{r}(t)||^2 \ \rangle$

de telle sorte que l'équation différentielle se réécrive sous la forme simple :

$m \ \frac{du(t)}{dt} \ = \ - \, k \ u(t) \ + \ m \ \langle \ || \vec{v} ||^2 \ \rangle$

Équipartition de l'énergie

On obtient une estimation du dernier terme de vitesse quadratique moyenne en utilisant le théorème d'équipartition de l'énergie de la physique statistique classique^[3]. Pour le mouvement d'une particule dans un espace à d dimensions, on obtient :

$\frac{1}{2} \ m \ \langle \ ||\vec{v}||^2 \ \rangle \ = \ \frac{d}{2} \ k_B \ T$

où $k B$ est la constante de Boltzmann, et $T$ la température absolue en kelvins. L'énergie thermique moyenne $k B T$ par particule peut se réécrire :

$k_B \ T \ = \ \frac{R \, T}{\mathcal{N}_A}$

où $R$ est la constante des gaz parfaits, et $\mathcal{N}_A$ le nombre d'Avogadro. L'équation différentielle se met donc finalement sous la forme :

$m \ \frac{du(t)}{dt} \ + \ k \ u(t) \ = \ \frac{d \, RT}{\mathcal{N}_A}$

Cette équation différentielle linéaire du premier ordre à coefficients constants avec second membre admet la solution exacte :

$u(t) \ = \ \frac{ d \, RT}{k \, \mathcal{N}_A} \ + \ \lambda \ \mathrm{e}^{ - \, t / \tau}$

où $λ$ est une constante, et $τ$ le temps caractéristique de relaxation, qui vaut :

$\tau \ = \ \frac{m}{k} \ = \ \frac{m}{6 \pi \eta a} \ \simeq$ 10^-8 secondes

dans les conditions d'observations expérimentales usuelles du mouvement brownien.

Coefficient de diffusion d'Einstein

Dans les conditions expérimentales usuelles, on est toujours dans le régime où : $t \gg \tau$ , et on observe alors :

$u(t) \ \sim \ \frac{d \, RT}{k \, \mathcal{N}_A} \ = \ \frac{d \, RT}{6 \pi \eta a \, \mathcal{N}_A}$

Compte-tenu de la définition de $u (t)$ , on a :

$u(t) \ = \ \frac{1}{2} \ \frac{d~}{dt} \ \langle \ || \vec{r}(t) ||^2 \ \rangle \ \sim \ \frac{d \, RT}{6 \pi \eta a \, \mathcal{N}_A}$

ce qui donne par intégration par rapport au temps t la loi de la diffusion classique :

$\langle \ ||\vec{r}(t)||^2 \ \rangle \ \sim \ \frac{2d \, RT}{6 \pi \eta a \, \mathcal{N}_A} \ t \ = \ 2d \ D \ t$

où le coefficient de diffusion $D$ s'écrit explicitement :

$D \ = \ \frac{RT}{6 \pi \eta a \, \mathcal{N}_A}$

On retrouve bien le résultat d'Einstein (1905).

Solution moderne

Cf. l'article de Bertrand Duplantier, paragraphe 1.5.3 , p. 177.

Bibliographie

Paul Langevin ; Sur la théorie du mouvement brownien, Comptes-Rendus de l'Académie des Sciences 146 (1908), 530-532. Possibilité de consulter et de télécharger le texte complet depuis le site Gallica de la Bibliothèque nationale de France : format pdf .

Bertrand Duplantier ; Le mouvement brownien, (2005), dans : Séminaire Poincaré Einstein, 1905-2005, (Paris, 8 avril 2005). Texte complet disponible ici.

Notes et références

↑ En termes modernes, un bruit blanc gaussien est un processus stochastique de moyenne nulle :

$\langle \, \vec{\eta}(t) \, \rangle \ = \ \vec{0}$

et totalement décorrélé dans le temps ; sa fonction de corrélation à deux points vaut en effet :

$\langle \, \eta_i(t_1) \ \eta_j(t_2) \, \rangle \ = \ \Gamma \ \delta_{ij} \ \delta(t_1-t_2)$

Dans cette formule, $Γ$ est une constante positive, $δ i j$ le symbole de Kronecker, et $δ(t)$ la distribution de Dirac, qui est identiquement nulle lorsque $t_1 \ne t_2$ Dans ces deux formules, la moyenne est prise sur toutes les réalisations possibles du bruit blanc gaussien.
↑ Langevin écrit :

« La valeur moyenne du terme $\vec{r} \cdot \vec{\eta}$ est évidemment nulle à cause des irrégularités des actions complémentaires $\vec{\eta}$ . »

Ce n'est en réalité pas si évident que cela ; lire par exemple l'article de Bertrand Duplantier, page 176, note 52. Cet auteur donne un peu plus loin dans le même article la dérivation moderne de la solution de l'équation stochastique de Langevin (paragraphe 1.5.3 , p. 177).
↑ Le théorème d'équipartition de l'énergie de la mécanique statistique classique dit que la valeur moyenne de l'énergie associée à un degré de liberté \emph{quadratique} d'un système mécanique en équilibre thermique avec un thermostat à la température $T$ est égale à $k B T / 2$ . Pour une particule ponctuelle qui n'est soumise à aucune force dans un espace à d dimensions, il y a exactement d degrés de liberté \emph{quadratiques}, qui correspondent aux d contributions à l'énergie cinétique : $mv_1^2/2, \, mv_2^2/2, \, \dots, mv_d^2/ 2$ , d'où le résultat utilisé ici.

Portail des mathématiques
Portail de la physique

Ce document provient de « %C3%89quation de Langevin ».

Catégories : Processus stochastique | Physique statistique

Wikimedia Foundation. 2010.

Contenu soumis à la licence CC-BY-SA. Source : Article Équation de langevin de Wikipédia en français (auteurs)

Игры ⚽ Нужна курсовая?

Regardez d'autres dictionnaires:

Equation de Langevin — Équation de Langevin L équation de Langevin (1908) est une équation stochastique pour le mouvement brownien. Sommaire 1 Théorie de Langevin du mouvement brownien 2 Équation de Langevin 3 Solution de Langevin (1908) … Wikipédia en Français
Équation de Langevin — L équation de Langevin (1908) est une équation stochastique pour le mouvement brownien. Sommaire 1 Théorie de Langevin du mouvement brownien 2 Équation de Langevin 3 Solution de Langevin (1908) … Wikipédia en Français
Equation differentielle stochastique — Équation différentielle stochastique Une équation différentielle stochastique (EDS) est une généralisation de la notion d équation différentielle prenant en compte un terme de bruit blanc. Les EDS permettent de modéliser des trajectoires… … Wikipédia en Français
Langevin equation — In statistical physics, a Langevin equation (Paul Langevin, 1908) is a stochastic differential equation describing the time evolution of a subset of the degrees of freedom. These degrees of freedom typically are collective (macroscopic) variables … Wikipedia
Équation différentielle stochastique — Une équation différentielle stochastique (EDS) est une généralisation de la notion d équation différentielle prenant en compte un terme de bruit blanc. Les EDS permettent de modéliser des trajectoires aléatoires, tels des cours de bourse ou les… … Wikipédia en Français
Langevin — Cette page d’homonymie répertorie les différents sujets et articles partageant un même nom. Toponymes Bloc Langevin, immeuble d Ottawa abritant le cabinet du Premier ministre canadien Rivière Langevin, fleuve français situé dans le département d… … Wikipédia en Français
Langevin — may refer to:People* Adélard Langevin (1855 1915), Canadian priest * Charles Langevin (1789 1869), Businessman * Chris Langevin (1959 ), Canadian ice hockey player * Dave Langevin (1954 ), American ice hockey player * Hector Louis Langevin (1826… … Wikipedia
Langevin dynamics — is an approach to mechanics using simplified models and using stochastic differential equations to account for omitted degrees of freedom.A molecular system in the real world is unlikely to be present in vacuum. Jostling of solvent or air… … Wikipedia
Paul Langevin — Pour les articles homonymes, voir Langevin. Paul Langevin Naissance 23 janvier 1872 … Wikipédia en Français
Stochastic differential equation — A stochastic differential equation (SDE) is a differential equation in which one or more of the terms is a stochastic process, thus resulting in a solution which is itself a stochastic process. SDE are used to model diverse phenomena such as… … Wikipedia

Dictionnaires et Encyclopédies sur 'Academic'

Équation de langevin

Équation de Langevin

Sommaire

Théorie de Langevin du mouvement brownien

Équation de Langevin