01.05.2013 Views

Modélisation, analyse mathématique et simulations numériques de ...

Modélisation, analyse mathématique et simulations numériques de ...

Modélisation, analyse mathématique et simulations numériques de ...

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.

tel-00656013, version 1 - 3 Jan 2012<br />

154 Asymptotic analysis of blood flow in stented arteries<br />

5.3.4 First or<strong>de</strong>r estimates<br />

The gain obtained when introducing the microscopic corrector is of or<strong>de</strong>r √ ε. In<strong>de</strong>ed,<br />

the following error estimates hold.<br />

Proposition 2. There exist two positive constants c3 and c4, <strong>de</strong>pending only on the mo<strong>de</strong><br />

Ĉk and not on the frequency k, such that:<br />

ûǫ,k − Ûǫ,k H 1 (Ωε) ≤ c3ε , ûǫ,k − Ûǫ,k L 2 (Ω0) ≤ c4ε 3/2 . (5.3.18)<br />

Proof. Denote Rε := ûǫ,k − Ûǫ,k the error to estimate. It is solution of the problem:<br />

⎧<br />

⎪⎨<br />

⎪⎩<br />

Lk Rε = ĈkχΩε\Ω0 −ikMkx2χ Ωε\Ω0 −ikMkε β( x<br />

ε )−β +βχ <br />

Ωε\Ω0 −εMkβδΓ0 , in Ωε<br />

x1 1<br />

Rε = −εMk β( ε , ε )−β<br />

on Γ1,<br />

Rε = 0 on Γε,<br />

Rε is x1-periodic on Γin ∪Γout.<br />

(5.3.19)<br />

The existence and uniqueness of Rε are standard. We focus again on the a priori estimates:<br />

test the system above by Rε and estimate the lhs from below as in (5.3.8), then estimate<br />

from above the rhs. The last step inclu<strong>de</strong>s new terms wrt the zeroth or<strong>de</strong>r approximation,<br />

listed below :<br />

⎧<br />

⎪⎨<br />

⎪⎩<br />

A1 = Ĉk<br />

<br />

Ωε\Ω0<br />

A2 = −ikMk<br />

<br />

<br />

A3 = −ikMkε<br />

Rε dx,<br />

Ωε\Ω0<br />

Ωε<br />

<br />

A4 = −ikMkεβ<br />

A5 = −εβMk<br />

Then, estimating these terms, one g<strong>et</strong>s<br />

|<br />

5<br />

j=1<br />

<br />

Γ0<br />

x2Rε dx,<br />

(β( x<br />

ε )−β)Rε dx,<br />

Ωε\Ω0<br />

Rε dx1.<br />

Rε dx,<br />

Aj| ≤ ε 3<br />

2 c∇Rε L 2 (Ωε)<br />

(5.3.20)

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

Saved successfully!

Ooh no, something went wrong!