14.09.2014 Views

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 ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

2. Problème de Cauchy stationnaire<br />

Trouverutel que<br />

⎧<br />

⎪⎨<br />

⎪⎩<br />

−∇. ( k(x)∇u ) =f dans Ω<br />

k(x)∇u·⃗n =φ sur Γ m<br />

u =T sur Γ m .<br />

(2.1)<br />

Ce problème est sévèrement mal posé. En effet, seule l’unicité d’une éventuelle<br />

solution est assurée [ISA 97] et l’instabilité de cette hypothétique solution remet<br />

en cause son existence.<br />

Le problème de complétion de données associé au problème de Cauchy (2.1)<br />

s’écrit<br />

Trouver (η d ,τ d ) sur Γ u tels qu’il existeuvérifiant<br />

⎧<br />

⎪⎨<br />

⎪⎩<br />

−∇. ( k(x)∇u ) =f dans Ω<br />

u =T,k(x)∇u·⃗n =φ sur Γ m<br />

u =τ d ,k(x)∇u·⃗n =η d sur Γ u .<br />

(2.2)<br />

Suivant <strong>la</strong> procédure décrite dans [AND 06], on introduit deux champs<br />

distinctsu 1 etu 2 qui diffèrent par les conditions limites qui leur sont imposées.<br />

Sur Γ m ,u 1 vérifie <strong>la</strong> condition de Dirichlet etu 2 <strong>la</strong> condition de Neumann et<br />

inversement sur Γ u . On considère alors, <strong>pour</strong> tout couple (η,τ) correspondant<br />

aux inconnues du problème de complétion de données, les deux problèmes bien<br />

posés suivants dont les champsu 1 etu 2 sont solutions :<br />

⎧<br />

⎪⎨<br />

⎪⎩<br />

−∇. ( k(x)∇u 1<br />

)<br />

=f dans Ω<br />

u 1 =T sur Γ m<br />

k(x)∇u 1·⃗n =η sur Γ u ,<br />

On considère <strong>la</strong> fonctionnelle<br />

∫<br />

E(η,τ) =<br />

Ω<br />

(2.3)<br />

⎧<br />

⎪⎨<br />

⎪⎩<br />

−∇. ( k(x)∇u 2<br />

)<br />

=f dans Ω<br />

u 2 =τ sur Γ u<br />

k(x)∇u 2·⃗n =φ sur Γ m .<br />

(2.4)<br />

k(x) ( ∇u 1 (η)−∇u 2 (τ) ) 2<br />

dx. (2.5)<br />

La <strong>méthode</strong> <strong>énergétique</strong> consiste alors à résoudre le problème de minimisation<br />

suivant :<br />

(η ∗ ,τ ∗ ) = argmin<br />

(η,τ)<br />

E(η,τ). (2.6)<br />

On a ainsi introduit explicitement les problèmes bien posés (2.3) et (2.4), <strong>la</strong><br />

fonctionnelle (2.5) et le problème de minimisation (2.6) <strong>pour</strong> <strong>la</strong> <strong>méthode</strong> <strong>énergétique</strong>.<br />

Dans <strong>la</strong> section suivante, nous introduirons le cadre fonctionnel et variationnel,<br />

notamment les espaces de traces <strong>pour</strong> les données et les inconnues.<br />

Nous montrerons que <strong>la</strong> <strong>méthode</strong> <strong>énergétique</strong> permet bien d’identifier <strong>la</strong> solution<br />

du problème de complétion de données et nous étudierons le problème de<br />

minimisation.<br />

16<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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!