24.02.2014 Views

Partial Differential Equations - Modelling and ... - ResearchGate

Partial Differential Equations - Modelling and ... - ResearchGate

Partial Differential Equations - Modelling and ... - ResearchGate

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

176 B. Maury<br />

Proposition 8. The problem (8) admits a unique weak solution u ∈ V =<br />

H0 1 (Ω), which is characterized as the solution to the minimization problem<br />

⎧<br />

Find u ∈ K such that<br />

⎪⎨<br />

∫<br />

J(u) = inf J(v), with J(v) =1<br />

v∈K 2 Ω<br />

⎪⎩<br />

K = {v ∈ H0 1 (Ω) |∇v =0a.e. in O},<br />

∫<br />

|∇u| 2 − fv<br />

Ω<br />

where f has been extended by 0 inside O.<br />

Furthermore, the restriction of u to Ω \ O is in H 2 (Ω \ O).<br />

Proof. Existence <strong>and</strong> uniqueness are direct consequences of the Lax–Milgram<br />

theorem applied in K = {v ∈ V |∇v =0a.e.inO}, which gives, in addition,<br />

the characterization of u as the solution to (9). Now the restriction of u to<br />

Ω \ O satisfies −△u = f, with regular Dirichlet boundary conditions on the<br />

boundary of Ω \ O which decomposes as ∂O∪∂Ω. AsΩ is a convex polygon<br />

<strong>and</strong> ∂O is smooth, st<strong>and</strong>ard theory ensures u| Ω\O<br />

∈ H 2 (Ω \ O).<br />

(9)<br />

We introduce the penalized version of the problem (9)<br />

⎧<br />

⎪⎨<br />

Find u ε ∈ V such that J ε (u ε ) = inf J ε (v),<br />

v∈V<br />

⎪⎩ J ε (v) = 1 ∫<br />

|∇v| 2 + 1 ∫ ∫<br />

|∇v| 2 − fv,<br />

2<br />

2ε<br />

for which linear convergence can be expected:<br />

Ω<br />

O<br />

Ω<br />

(10)<br />

Proposition 9. Let u be the solution to the problem (9), u ε the solution to<br />

the problem (10). It holds ‖u − u ε ‖ H1 (Ω) = O(ε).<br />

Proof. Let us show that<br />

B : v ∈ H 1 0 (Ω) ↦−→ ∇v ∈ L 2 (O) 2<br />

has a closed range. Consider µ ∈ Λ with µ = ∇v. We define w ∈ H 1 0 (O) as<br />

w = v − m(v), where m(v) isthemeanvalueofv over O. BythePoincaré–<br />

Wirtinger inequality, one has<br />

‖w‖ H1 (O) ≤ C ‖µ‖ L 2 (O) 2 .<br />

Now, as O⊂⊂Ω, there exists a continuous extension operator from H 1 (O)<br />

to H0 1 (Ω), so that we can extend w to obtain ˜w ∈ H0 1 (Ω) with a norm controlled<br />

by ‖µ‖ L 2 (O) 2, which proves the closed character of B(V ). The linear<br />

convergence is then given by Corollary 1.

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

Saved successfully!

Ooh no, something went wrong!