next up previous contents index
suivant: III.E.2 Problème de Helmholtz tridimensionnel. monter: III.E Déterminant de la matrice D du précédent: III.E Déterminant de la matrice D du   Table des matières   Index


III.E.1 Problème de Helmholtz bidimensionnel.

Pour le triangle numéroté par l'indice k, on définit les quantités suivantes:

  1. les trois longueurs notées L1,L2,L3,
  2. les trois normales notées $\vec{\nu}_1$, $\vec{\nu}_2$ et $\vec{\nu}_3$,
  3. le périmètre noté P,
  4. la surface notée S.
On s'intéresse au cas asymptotique $h \rightarrow 0$ et p=3 pour un problème bidimensionnel.

Théorème 19  

On introduit la matrice limite Dr définie par

(III.E.1) \begin{displaymath}%
D^r=\lim_{h \to 0}{\frac{\displaystyle D}{\displaystyle h\omega^2}} \ ,
\end{displaymath}

ainsi que les quantités réduites indexées (Pr,Sr,L1r,L2r,L3r) qui sont obtenues des précédentes en divisant par $h\omega^2$ sauf Sr où l'on divise par $(h\omega^2)^2$. Alors:
  1. La matrice limite Dr existe et est inversible.
  2. Le déterminant de Dr se calcule par le produit de deux quantités, l'une faisant intervenir la géométrie du triangle, l'autre la répartition des fonctions de base. Des ondes planes équiréparties maximisent le déterminant de Dr et il vaut alors
    (III.E.2) \begin{displaymath}
\mbox{\fbox{$\mbox{det} D^r =
27.\frac{\displaystyle (P^rS^r)^2}{\displaystyle L_1^rL_2^rL_3^r}
$}}
\end{displaymath}

  3. La majoration du conditionnement à partir du déterminant est maximale pour des ondes planes équiréparties et on majore le conditionnement de D par
    (III.E.3) \begin{displaymath}%
\frac{\displaystyle \lambda_{max}}{\displaystyle \lambda_{m...
...\le \frac{\displaystyle 12}{\displaystyle \pi} (4\sigma)^4 \ .
\end{displaymath}

    dont le terme majorant $\sigma=h/\rho$ (I.2.4) est minimal dans le cas d'un triangle équilatéral.

La preuve de ce théorème fait l'objet de toute cette section. La preuve est longue, technique et calculatoire. Elle est effectuée en plusieurs étapes et suit l'ordre de présentation des résultats effectuée théorème 19:

a)
simplification du déterminant de la matrice D (où l'on ne garde que les termes d'ordre 1 en h) en un produit dépendant de l'unité d'échelle du triangle h et un terme indépendant de h,
b)
découplage du déterminant (toujours de la matrice D où l'on n'a gardé que les termes du premier ordre en h) entre les termes géométriques et les fonctions de base,
c)
optimisation du déterminant en fonction du choix des directions de propagation des ondes planes,
d)
récapitulatif des points précédents et fin du calcul du déterminant (ceci montre le point 2 du théorème 19),
e)
majoration du conditionnement (point 3 du théorème 19).

Le lemme suivant donne la forme de la matrice réduite Dr.

Lemme 29   Les termes de la matrice Dk (calculés chapitre III.A) sont équivalents au premier ordre à la matrice, encore notée D (par abus) dans toute cette section, dont les termes Dl,m pour $(l,m) \in \{1,2,3 \}^2$ sontIII.E1
(III.E.4) \begin{displaymath}%
D^{l,m}=\omega^2 \sum_{n=1}^{3}{L_n(1-\vec{\nu}_n \, \vec{v}_{m})(1-\vec{\nu}_n.\vec{v}_{l})} \ .
\end{displaymath}

Les termes de la matrice réduite sont donc
(III.E.5) \begin{displaymath}%
D^r_{l,m}=\sum_{n=1}^{3}{L^r_n(1-\vec{\nu}_n \, \vec{v}_{m})(1-\vec{\nu}_n.\vec{v}_{l})} \ .
\end{displaymath}


\begin{proof}
% latex2html id marker 36824Rappelons que les termes de la matri...
... \buildl \rightarrow_{h\to 0} 1
\crcr}}\,
\right.
\end{displaymath}
\end{proof}
On pose:

(III.E.6) \begin{displaymath}%
\vec{w}_n = \left(
\begin{array}{l}
1-\vec{\nu}_n.\vec{v}_...
...}_{2} \\
1-\vec{\nu}_n.\vec{v}_{3} \\
\end{array}\right) \ ,
\end{displaymath}

et D sera la matrice définie par
(III.E.7) \begin{displaymath}%
\null\,\vcenter{\openup1\jot \let\\ =\@
\ialign{\strut\hfi...
...+L_3D_3) \cr
& D_n = [\vec{w}_n\otimes\vec{w}_n] \ .
\crcr}}\,
\end{displaymath}

Remarque 54   En toute rigueur, on a

\begin{displaymath}
D_k = \omega^2(L_1D_1+L_2D_2+L_3D_3) + O(h^2) \ .
\end{displaymath}

On pose en outre, dans toute cette section,
(III.E.8) \begin{displaymath}
\xi=(\vec{w}_1,\vec{w}_{2},\vec{w}_{3})_{mixte} \ .
\end{displaymath}

a) Simplification du déterminant de D.

Lemme 30   La matrice D est inversible si et seulement si les trois vecteurs $\vec{w}_n$ sont indépendants dans ${\mathbb{R}}^3$:

\begin{displaymath}
\mbox{\fbox{$\mbox{det} D =
\omega^6L_1L_2L_3.(\vec{w}_{1},\vec{w}_{2},\vec{w}_{3})_{mixte}^2$}}
\end{displaymath}


\begin{proof}
On d\'efinit les trois vecteurs $\vec{e}_n=\vec{w}_{n+1} \wedge \v...
...L_3.(\vec{w}_{1},\vec{w}_{2},\vec{w}_{3})_{mixte}^2
\end{displaymath}\end{proof}

b) Découplage du déterminant entre les termes géométriques et les fonctions de base.

Lemme 31   $L_1\vec{\nu}_1+L_2\vec{\nu}_2+L_3\vec{\nu}_3 = 0$ .


\begin{proof}
On introduit $R$\ la rotation d'angle $\pi/2$\ dans le plan. Alors...
...ec{x}_{3}-\vec{x}_{1}+\vec{x}_{1}-\vec{x}_{2})=\vec{0}
\end{equation}\end{proof}

Lemme 32   Pour $\xi = (\vec{w}_{1},\vec{w}_{2},\vec{w}_{3})_{mixte}$ (cf (III.E.9)) et P=(L1+L2+L3), on a:
(III.E.9) \begin{displaymath}
\xi = \frac{\displaystyle P}{\displaystyle L_3} (\vec{\nu}_1...
...1 \\
v_{3}^1 &
v_{3}^2 &
1 \\
\end{array}\right\vert \ .
\end{displaymath}


\begin{proof}
% latex2html id marker 37180Par d\'efinition,
\begin{displaymat...
...uation.h2adetddecou.013}) on tire (\ref{equation.h2adetddecou.007}).
\end{proof}

c) Optimisation du déterminant.

Le calcul du produit mixte de trois vecteurs donne le volume du simplexe déterminé par ces trois vecteurs: nous avons pensé que le déterminant est de module maximal pour des vecteurs vn équirépartis dans le plan. C'est ce que nous allons montrer dans le lemme suivant.

Lemme 33   On montre que $\vert\xi\vert$ est maximal pour des vecteurs vn équirépartis dans le plan.


\begin{proof}
On pose, pour $(\alpha,\beta) \in ]-\pi,\pi]^2$\ et $\alpha \ne \b...
... = \frac{\displaystyle 3\sqrt{3}}{\displaystyle 2} \ .
\end{equation}\end{proof}

d) Récapitulatif et fin du calcul du déterminant.

Lemme 34   $L_1L_2(\vec{\nu}_{1},\vec{\nu}_{2})_{mixte}=$ 2 fois l'aire S du triangle.


\begin{proof}
% latex2html id marker 37493Le raisonnement est identique \\lq a ce...
..._{2}-\vec{x}_{3},\vec{x}_{3}-\vec{x}_{1})_{mixte}\vert
\end{equation}\end{proof}


\begin{proof}
% latex2html id marker 37518Nous effectuons la preuve des deux p...
...displaystyle 27P^2}{\displaystyle L_1L_2L_3}S^2
$}}
\end{displaymath}\end{proof}

e) Majoration du conditionnement. Nous allons montrer le dernier point du théorème 19.

Lemme 35   Des ondes planes équiréparties assurent la majoration suivante du conditionnement des matrices D et Dr:
(III.E.10) \begin{displaymath}%
\mbox{\fbox{$\frac{\displaystyle \lambda_{max}}{\displaysty...
...leq \frac{\displaystyle 12}{\displaystyle \pi} (4\sigma)^4
$}}
\end{displaymath}

$\sigma=h/\rho$ est défini en (I.2.4). Cette majoration est obtenue pour un élément régulier, c'est-à-dire vérifiant les hypothèses de régularité H1, H2 et H3 de la section I.2.1.2. Dans le cas d'un triangle équilatéral, nous maximisons le terme majorant $\sigma$.


\begin{proof}
% latex2html id marker 37586On note $\lambda_{min}$\ et $\lambda...
... triangle \'equilat\'eral,
ce que l'on comprend bien intuitivement.
\end{proof}

f) Remarques.

Remarque 55   Ce résultat est généralisable au problème de Helmholtz scalaire dans un espace tridimensionnel avec un tétraèdre et 4 fonctions de base (cf Annexe III.E.2). En revanche, on ne peut pas le généraliser au cas des n-polygones avec n fonctions de base $\vec{v}_n$. On peut en donner un contre-exemple pour n=4. Prenons $(\vec{x}_1,\vec{x}_2,\vec{x}_3,\vec{x}_4)$ un carré, $(\lambda_1,\lambda_2,\lambda_3,\lambda_4)=(1,-1,1,-1)$. Les vecteurs wn sont alors liés.

Remarque 56   Il est trivial de remarquer que la matrice limite Dr ( $h \rightarrow 0$) avec quatre fonctions de base ou plus ($p \geq 4$) n'est pas inversible. En effet, la nouvelle matrice Dr obtenue s'écrira toujours sous la forme (en gardant toujours le même sens à la notation en indice r):

\begin{displaymath}
\null\,\vcenter{\openup1\jot \let\\ =\@
\ialign{\strut\hfil...
..._3^rD_3) \cr
& D_n = [\vec{w}_n\otimes\vec{w}_n] \ ,
\crcr}}\,
\end{displaymath}

avec

\begin{displaymath}
\vec{w}_n = \left(
\begin{array}{l}
1-\vec{\nu}_n.\vec{v}_{...
...ec{v}_{3} \\
1-\vec{\nu}_n.\vec{v}_{4} \\
\end{array}\right)
\end{displaymath}

Le rang de la matrice Dr est au plus de trois, pour une matrice de taille p x p et $p \geq 4$.


next up previous contents index
suivant: III.E.2 Problème de Helmholtz tridimensionnel. monter: III.E Déterminant de la matrice D du précédent: III.E Déterminant de la matrice D du   Table des matières   Index
Cessenat Olivier 2007-04-21