next up previous contents index
suivant: III.D.2 Algorithmes conditionnels de programmation. monter: III.D.1 Formules analytiques des intégrales d'ondes planes. précédent: III.D.1.2 Formules d'intégration sur un triangle: calcul   Table des matières   Index


III.D.1.3 Formules de l'intégrale sur un tétraèdre: calcul de I3.

Nous calculons I3 l'intégrale de (III.D.1) sur le tétraèdre [x1,x2,x3,x4] de volume V et de barycentre $\vec{G}$:

(III.D.6) \begin{displaymath}
I_3 = \displaystyle \int _{[x_1,x_2,x_3,x_4]}{e^{i{{\mathbf k}}{{\mathbf X}}}} \ d \, {{\mathbf X}}\ .
\end{displaymath}

On pose, pour $n-1=0 \ldots 3$ modulo 4

\begin{displaymath}
Z_n= e^{ i {{\mathbf k}}\vec{x}_n}
\end{displaymath}

et

\begin{displaymath}
v_n= \frac{\displaystyle {{\mathbf k}}\vec{x}_n}{\displaystyle 2} \ .
\end{displaymath}

On effectue le changement de variable affine

\begin{displaymath}
\vec{x}=\vec{x}_1+ \theta (\vec{x}_2-\vec{x}_1 )+ (1-\theta)...
...}_3-\vec{x}_1 )+ (1-\phi) \zeta (\vec{x}_4-\vec{x}_1 ) \right)
\end{displaymath}

dans (III.D.6), ce qui donne

\begin{displaymath}
I_3 = 6V e^{i{{\mathbf k}}\vec{x}_1} \int_{[0,1]} (1-\theta)...
...\phi) \zeta (\vec{x}_4-\vec{x}_1 ) \right)} \,d \zeta \,d \phi
\end{displaymath}


\begin{displaymath}
I_3 = 6V e^{i{{\mathbf k}}\vec{x}_1} \int_{[0,1]} (1-\theta)...
...e 2i(v_{3}-v_{1})}}{\displaystyle -2i(v_{1}-v_{4})} \,d \theta
\end{displaymath}


\begin{displaymath}
I_3 = 6V Z_1 e^{2i (v_2-v_1 ) } \int_{[0,1]} (1-\theta) \fra...
...e 2i(v_{3}-v_{1})}}{\displaystyle -2i(v_{1}-v_{4})} \,d \theta
\end{displaymath}


\begin{displaymath}
I_3 = 6V Z_1 e^{2i (v_2-v_1 ) } \frac{\displaystyle \frac{\d...
...laystyle 2i(v_{3}-v_{1})}}{\displaystyle -2i(v_{1}-v_{4})} \ .
\end{displaymath}

Après simplification, on a

\begin{displaymath}
I_3 = 6V Z_2 \frac{\displaystyle \frac{\displaystyle e^{i (v...
...playstyle 2i(v_{3}-v_{1})}}{\displaystyle 2i(v_{4}-v_{1})} \ ,
\end{displaymath}

ou

\begin{displaymath}
I_3 = 6V \frac{\displaystyle \frac{\displaystyle e^{i (v_{4}...
...playstyle 2i(v_{3}-v_{1})}}{\displaystyle 2i(v_{4}-v_{1})} \ .
\end{displaymath}

En posant $\alpha_n = {{\mathbf k}}\frac{\displaystyle (\vec{x}_{n-1}-\vec{x}_n )}{\displaystyle 2}$, $\beta_n = {{\mathbf k}}\frac{\displaystyle (\vec{x}_{n+1}-\vec{x}_n )}{\displaystyle 2}$ et $\gamma_n = {{\mathbf k}}\frac{\displaystyle (\vec{x}_{n+2}-\vec{x}_n )}{\displaystyle 2}$, pour tout n-1 de 0 à 3 modulo 4, on a

\begin{displaymath}
I_3 = 6V Z_2 {\frac{\displaystyle I(\gamma_2 , \beta_2)-I(\beta_2 , \alpha_2) }{\displaystyle 2i(\gamma_2-\alpha_2)} }
\end{displaymath}

avec
(III.D.7) \begin{displaymath}
I(\lambda , \mu) = {\frac{\displaystyle e^{i\lambda} {\frac{...
...u}}{\displaystyle \mu}} }{\displaystyle 2i(\lambda-\mu)} } \ .
\end{displaymath}

Finalement, par symétrie, l'intégrale I3 apparaît comme la fonction de $\alpha_n$, $\beta_n$ et $\gamma_n$ suivante

(III.D.8) \begin{displaymath}
I_3(\alpha_n,\beta_n,\gamma_n) = 6V e^{i{{\mathbf k}}\vec{G}...
...eta_n , \alpha_n) }{\displaystyle 2i(\gamma_n-\alpha_n)} } \ ,
\end{displaymath}

fonction que l'on sait continue en $(\alpha_n,\beta_n,\gamma_n) \in {\mathbb{C}}^3$.


next up previous contents index
suivant: III.D.2 Algorithmes conditionnels de programmation. monter: III.D.1 Formules analytiques des intégrales d'ondes planes. précédent: III.D.1.2 Formules d'intégration sur un triangle: calcul   Table des matières   Index
Cessenat Olivier 2007-04-21