Schéma de Lax-Friedrichs

Considérons l'équation d'advection linéaire en dimension 1 d'espace et de temps, dont une solution u(x,t) doit vérifier : ${\displaystyle (1)\qquad u_{t}+au_{x}=0\,}$ sur le domaine ${\displaystyle b\leq x\leq c,\;0\leq t\leq d}$ avec la condition initiale ${\displaystyle u(x,0)=u_{0}(x)\,}$ et les conditions de bords ${\displaystyle u(b,t)=u_{b}(t)\,}$ ${\displaystyle u(c,t)=u_{c}(t).\,}$

La solution exacte de l'équation d'advection au temps ${\displaystyle t\geqslant 0}$ est ${\displaystyle (2)\qquad u(x,t)=u_{0}(x-at)}$

La méthode des différences finies consiste à chercher une solution discrète définie sur les points de coordonnées ${\displaystyle (x_{i},t^{n})=(b+i\,\Delta x,n\,\Delta t)\ {\text{ pour }}\ i=0,\ldots ,N,\,n=0,\ldots ,M,}$ avec ${\displaystyle N={\frac {c-b}{\Delta x}},\,M={\frac {d}{\Delta t}}}$

On recherche alors la solution discrète ${\displaystyle u_{i}^{n}\approx u(x_{i},t^{n})}$

Le schéma de Lax–Friedrichs appliqué à l'équation d'advection donne : ${\displaystyle (3)\qquad {\frac {u_{i}^{n+1}-{\frac {1}{2}}(u_{i+1}^{n}+u_{i-1}^{n})}{\Delta t}}+a{\frac {u_{i+1}^{n}-u_{i-1}^{n}}{2\,\Delta x}}=0.}$

On obtient alors de manière explicite l'inconnue u_iⁿ⁺¹ : ${\displaystyle (4)\qquad u_{i}^{n+1}={\frac {1}{2}}(u_{i+1}^{n}+u_{i-1}^{n})-a{\frac {\Delta t}{2\,\Delta x}}(u_{i+1}^{n}-u_{i-1}^{n})\,}$ ce qui constitue ainsi le schéma de Lax-Friedrichs, complété par la condition initiale et les conditions de bord : ${\displaystyle u_{i}^{0}=u_{0}(x_{i})}$ ${\displaystyle u_{0}^{n}=u_{b}(t^{n})}$ ${\displaystyle u_{N}^{n}=u_{c}(t^{n}).}$

Une propriété remarquable de ce schéma est qu'il découple les points « pairs » (i.e. u_iⁿ avec i pair) et « impairs » ce qui n'est pas « naturel » car les valeurs u_iⁿ et u_i+1ⁿ devraient a priori suivre des évolutions voisines (pour une solution u régulière)^[1].

Un système hyperbolique de lois de conservation à une dimension d'espace est défini par ${\displaystyle u_{t}+(f(u))_{x}=0}$ où ${\displaystyle f}$ est appelé fonction de flux.

Dans le cas particulier où ${\displaystyle f(u)=au}$, on retrouve un problème linéaire scalaire. Dans le cas général, ${\displaystyle u}$ est un vecteur ayant ${\displaystyle m}$ composantes. La généralisation du schéma de Lax-Friedrichs aux systèmes non linéaires prend la forme^[2] ${\displaystyle u_{i}^{n+1}={\frac {1}{2}}(u_{i+1}^{n}+u_{i-1}^{n})-{\frac {\Delta t}{2\,\Delta x}}(f(u_{i+1}^{n})-f(u_{i-1}^{n})).}$

Ce schéma conservatif est d'ordre 1 en espace et en temps, et donc assez diffusif. On peut en revanche l'utiliser pour construire des schémas d'ordre supérieur pour résoudre des systèmes d'équations aux dérivées partielles hyperbolique, de la même façon que la méthode d'Euler sert à construire les méthodes de type Runge-Kutta plus précises pour la résolution des équations différentielles ordinaires.

Ce schéma peut être écrit sous sa forme conservative : ${\displaystyle u_{i}^{n+1}=u_{i}^{n}-{\frac {\Delta t}{\Delta x}}\left({\hat {f}}_{i+1/2}^{n}-{\hat {f}}_{i-1/2}^{n}\right),}$ où ${\displaystyle {\hat {f}}_{i-1/2}^{n}={\frac {1}{2}}\left(f_{i-1}+f_{i}\right)-{\frac {\Delta x}{2\Delta t}}\left(u_{i}^{n}-u_{i-1}^{n}\right).}$

En l'absence des termes ${\displaystyle u_{i}^{n}}$ et ${\displaystyle u_{i-1}^{n}}$ dans le flux discret, ${\displaystyle {\hat {f}}_{i-1/2}^{n}}$, on retrouve le schéma explicite centré instable ^[1].