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.

174 B. Maury<br />

Proof. By the definition of u ε ,<br />

J ε (u ε )= 1 2 a(uε ,u ε ) −〈ϕ, u ε 〉 + 1 2ε b(uε ,u ε ) ≤ J ε (u) = 1 a(u, u) −〈ϕ, u〉,<br />

2<br />

so that<br />

0 ≤ 1 2ε b(u ε,u ε ) ≤ 1 2 a(u, u) − 1 2 a(uε ,u ε )+〈ϕ, u ε − u〉<br />

≤ 1 2 a(u + uε ,u− u ε )+〈ϕ, u ε − u〉,<br />

which yields the estimate by continuity of a(·, ·) <strong>and</strong>ϕ.<br />

Proposition 6. Under the assumptions (1), we denote by u ε h<br />

the solution to<br />

the problem (6). Then<br />

(<br />

|u ε h − u| ≤C min |v h − u| + √ )<br />

|u ε − u| .<br />

v h ∈V h ∩K<br />

Proof. As u ε h minimizes a(v − uε ,v− u ε )+b(v − u ε ,v− u ε )/ε over V h ,<br />

α |u ε h − u ε | 2 ≤ a(u ε h − u ε ,u ε h − u ε )<br />

≤ a(u ε h − u ε ,u ε h − u ε )+ 1 ε b(uε h − u ε ,u ε h − u ε )<br />

(<br />

≤ min a(v h − u ε ,v h − u ε )+ 1<br />

v h ∈V h ε b(v h − u ε ,v h − u ε )<br />

≤<br />

min<br />

v h ∈V h ∩K<br />

(<br />

a(v h − u ε ,v h − u ε )+ 1 ε b(v h − u ε ,v h − u ε )<br />

As v h is in K, the second term is b(u ε ,u ε )/ε, which is bounded by C |u ε − u|<br />

(by Lemma 1). Finally, we get<br />

(<br />

|u ε h − u ε |≤C min |v h − u ε | + √ )<br />

|u ε − u| ,<br />

v h ∈V h ∩K<br />

from which we conclude.<br />

Proposition 7. Under the assumptions (1), it holds<br />

|u ε h − u| ≤ C √ ε<br />

inf<br />

v h ∈V h<br />

|u ε − v h | + |u ε − u| ,<br />

)<br />

)<br />

.<br />

where u ε h<br />

is the solution to (6).<br />

Proof. One has<br />

|u ε h − u| ≤|u ε h − u ε | + |u ε − u| ,<br />

<strong>and</strong> we control the first term by Céa’s lemma applied to the bilinear form<br />

a + b/ε, whose ellipticity constant behaves like 1/ε.

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

Saved successfully!

Ooh no, something went wrong!