Homogénéisation

Modèle:Confusion Dans les mathématiques et la physique, lModèle:'homogénéisationModèle:, est un champ scientifique qui s'est développé à partir des années 1970 et qui a pour objet l'étude de systèmes multi-échelles. Plus précisément, l'homogénéisation s'attache à l'étude d'équations aux dérivées partielles dont un terme oscille fortement. Ces oscillations sont généralement liées à l'étude de milieux présentant des hétérogénéités à l'échelle microscopique (par exemple, des matériaux composites). L'objet de la théorie de l'homogénéisation est de proposer une équation « effective » (ou « homogénéisée ») généralement plus simple, qui décrive le comportement de la solution de l'équation considérée dans la limite où la petite échelle tend vers 0. Un des buts de cette théorie est de simplifier ainsi la simulation numérique de systèmes physiques complexes faisant intervenir plusieurs échelles.

Domaines d'applications

Initialement conceptualisée pour des équations elliptiques, la méthode d'homogénéisation par l'analyse asymptotique s'étend à divers types d'équations stationnaires ou non, à commencer par les équations de transport décrites par une équation de Boltzmann dont la diffusion constitue une approximation qui est retrouvée par cette approche. On trouve ainsi des exemples d'application dans des domaines aussi divers que la diffusion de masse ou de chaleur, la mécanique des fluides ou le transfert radiatif. Elle s'applique également à la mécanique des milieux continusModèle:, ou l'électromagnétisme.

Exemple d'une équation elliptique

On traite ici de la méthode utilisant un développement asymptotique sur l'exemple d'une équation elliptique. L'utilisation de cette technique nécessite que le milieu considéré ait une structure spécifique: périodique (comme ci-dessous), presque-périodique, ou encore aléatoire avec des propriétés de stationnarité et d'ergodicitéModèle:,.

On considère ici une équation elliptique pour la fonction inconnue Modèle:Math dans le domaine $𝒟\subset {ℝ}^d$

$\{\begin{array}{ccc}-\mathrm{\nabla }\cdot (𝖠\left(𝐱\right)\cdot \mathrm{\nabla }u(𝐱))& =& f(𝐱),𝐱\in 𝒟\\ [0.6em]u{|}_{\partial 𝒟}& =& g\end{array}$

où $f$est un terme source et $g$ est la donnée au bord imposée. On suppose que la matrice $𝖠$ est définie positive (éventuellement symétrique).

Le problème est défini sur un milieu comportant une échelle de variation lente Modèle:Mvar et une échelle de variation rapide $y=\frac{𝐱}{\epsilon }$ où Modèle:Mvar mesure l'échelle microscopique

$-\mathrm{\nabla }\cdot (𝖠\left(𝐲\right)\cdot \mathrm{\nabla }{u}^{\epsilon }(𝐲))=f(𝐱)$

Lorsque Modèle:Mvar tend vers 0, cette équation peut être efficacement approximée par une équation — dite équation homogénéisée — faisant intervenir une matrice ${𝖠}^{\star }$ qui s'écrit

$-\mathrm{\nabla }\cdot ({𝖠}^{\star }(𝐱)\cdot \mathrm{\nabla }{u}^{\star }(𝐱))=f(𝐱)$

dans le sens où

$\lim {\mathrm{\backslash limits}}_{\epsilon \to 0}{u}^{\epsilon }(𝐱)=u^{\star }(𝐱)$

Dans le cas où  $𝖠(𝐲)$ est un coefficient périodique, la matrice homogénéisée ${𝖠}^{\star }$ est constante, d'où une simplification substantielle du problème.

Analyse asymptotique

• Le milieu est supposé périodique de cellule $ℰ=[0,1{]}^d$. C'est-à-dire que, pour $({𝐞}_i)$ la base canonique de ${ℝ}^d$, on a

$𝖠(𝐱+{𝐞}_i)=𝖠(𝐱)\text{ pour tout }i\in \{1,\cdots ,d\}$

• Modèle:Math et Modèle:Math sont considérées comme des variables indépendantes. On a donc

$u(𝐱)=u^{\epsilon }(𝐱,𝐲)\Rightarrow \mathrm{\nabla }={\mathrm{\nabla }}_x+\frac{1}{\epsilon }{\mathrm{\nabla }}_y,\mathrm{\nabla }\cdot ={\mathrm{\nabla }}_x\cdot +\frac{1}{\epsilon }{\mathrm{\nabla }}_y\cdot$

La solution est développée sous forme d'une série de Hilbert, où chaque terme est périodique par rapport à la seconde variable $𝐲$

$u^{\epsilon }(𝐱,𝐲)=u_0(𝐱,𝐲)+\epsilon u_1(𝐱,𝐲)+{\epsilon }^2{u}_2(𝐱,𝐲)+\cdots$

On obtient ainsi

$\mathrm{\nabla }{u}^{\epsilon }={\epsilon }^{-1}{\mathrm{\nabla }}_y{u}_0+\sum _{i=0}^{\mathrm{\infty }}{\epsilon }^i({\mathrm{\nabla }}_y{u}_{i+1}+{\mathrm{\nabla }}_x{u}_i)$

Le regroupement des termes du même ordre permet d'obtenir à l'ordre 0 l'équation homogénéisée

$-{\mathrm{\nabla }}_x\cdot \left({𝖠}^{\star }{\mathrm{\nabla }}_x{u}_0\right)=f$

où  ${𝖠}^{\star }$ est une matrice constante obtenu par la résolution d'un problème à l'échelle locale.

Modèle:DémonstrationDans le cas monodimensionnel, on peut même obtenir une expression explicite de la matrice homogénéisée: il s'agit de la moyenne harmonique de la matrice $𝖠$:

${𝖠}^{\star }={\left(⟨{𝖠}^{-1}⟩\right)}^{-1}$Modèle:Démonstration

Voir aussi

• Prise de moyenne volumique
• Théories des milieux effectifs
• Micromécanique