Analyse numérique d'une méthode énergétique pour la résolution ...
Analyse numérique d'une méthode énergétique pour la résolution ...
Analyse numérique d'une méthode énergétique pour la résolution ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
2. Problème de Cauchy stationnaire<br />
Proposition 2.2.1 Si les données (φ,T) sont compatibles et si <strong>la</strong> solution<br />
(η ∗ ,τ ∗ ) du problème de minimisation (2.6) existe, alors (η ∗ ,τ ∗ ) est solution du<br />
problème de complétion de données (2.2) à une constante près <strong>pour</strong> l’inconnue de<br />
Dirichlet. En outre, lorsque le minimum est atteint, on a<br />
∇u 1 (η ∗ ) =∇u 2 (τ ∗ ) et u 1 (η ∗ ) =u 2 (τ ∗ ) +c dans ¯Ω oùcest une constante.<br />
(2.8)<br />
Démonstration : Si les données (φ,T) sont compatibles, il existe une unique<br />
solutionu d au problème de Cauchy. L’unique solution du problème de complétion<br />
de données est donc donnée par<br />
(η d ,τ d ) = (u d|Γu ,k(x)∇u d·⃗n |Γu ) sur Γ u (2.9)<br />
On a donc<br />
d’où<br />
Ce qui implique<br />
u 1 (η d ) =u 2 (τ d ) =u d dans ¯Ω, (2.10)<br />
0≤E(η ∗ ,τ ∗ )≤E(η d ,τ d ) = 0. (2.11)<br />
∇u 1 (η ∗ ) =∇u 2 (τ ∗ ) et doncu 1 (η ∗ ) =u 2 (τ ∗ ) +c dans ¯Ω, oùcest une constante.<br />
Par conséquent, <strong>pour</strong> toutx∈Γ,<br />
u 1 (x;η ∗ ) |Γ =u 2 (x;τ ∗ ) |Γ +c =<br />
=½{x∈Γ u}(x) (τ ∗ (x) +c) +½{x∈Γ m}(x)T(x) =u d|Γ (x)∈H 1/2 (Γ) (2.12)<br />
etk(x)∇u 1 (x;η ∗ ) |Γ·⃗n =k(x)∇u 2 (x;τ ∗ ) |Γ·⃗n =<br />
=½{x∈Γ u}(x)η ∗ (x) +½{x∈Γ m}(x)φ(x) =k(x)∇u d·n |Γ ∈H −1/2 (Γ). (2.13)<br />
D’où, (η ∗ ,τ ∗ +c) = (η d ,τ d ).<br />
<br />
La solution du problème de minimisation est donc bien solution du problème<br />
de complétion de données à une constante près <strong>pour</strong> l’inconnue de Dirichlet.<br />
De façon à se ramener à des problèmes homogènes, il est nécessaire d’introduire<br />
<strong>la</strong> notion de relèvement des conditions limites de Dirichlet. Cependant, les<br />
opérateurs de relèvement étant définis dans les espacesH s (Γ), il est nécessaire de<br />
prolongerT etτ sur l’ensemble de <strong>la</strong> frontière. Il existe deux opérateurs linéaires<br />
et continusP m etP u définis par<br />
P m ∈ L ( H 1/2 (Γ m ),H 1/2 (Γ) ) , et P u ∈ L ( H 1/2 (Γ u ),H 1/2 (Γ) ) ,<br />
P m g =g p.p. sur Γ m ,∀g∈H 1/2 (Γ m ), P u g =g p.p. sur Γ u ,∀g∈H 1/2 (Γ u ).<br />
On pose alors<br />
˜T =P m T et ˜τ =P u τ ∈H 1/2 (Γ). (2.14)<br />
18<br />
Cette thèse est accessible à l'adresse : http://theses.insa-lyon.fr/publication/2011ISAL0075/these.pdf<br />
© [R. Rischette], [2011], INSA de Lyon, tous droits réservés