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 ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
tel-00656013, version 1 - 3 Jan 2012<br />
134 Système continu<br />
cependant obtenir <strong>de</strong>s bornes L ∞ sur les solutions, ce qui suffira. Pour cela, nous préférons<br />
la formulation en u <strong>et</strong> v <strong>de</strong> ce système, à savoir :<br />
avec les conditions initiales<br />
⎧<br />
⎪⎨<br />
⎪⎩<br />
⎧<br />
⎪⎨<br />
⎪⎩<br />
u ′ = 0,<br />
v ′ = − 1<br />
ε R(u,v),<br />
u(0) = − 1<br />
√ a R0 := U0,<br />
v(0) = 0.<br />
(B.2.9)<br />
(B.2.10)<br />
Ce système différentiel ordinaire possè<strong>de</strong> une unique solution globale que nous allons « expliciter<br />
» grâce au développement <strong>de</strong> Taylor <strong>de</strong> R <strong>et</strong> à la formule <strong>de</strong> Duhamel. Tout d’abord<br />
la fonction u est i<strong>de</strong>ntiquement égale à U0 pour tout temps. Ensuite, la formule <strong>de</strong> Taylor<br />
appliquée à R entre les points (U0,A(U0)) <strong>et</strong> (U0,v(t)) donne, pour tout t > 0 :<br />
où la fonction α s’écrit :<br />
R(U0,v(t)) = α(t) (v(t)−A(U0)) ,<br />
α(t) := ∂vR U0; σv(t)+(1−σ)A(U0) ,<br />
avec le paramètre σ compris entre 0 <strong>et</strong> 1. Nous pouvons maintenant écrire la formule <strong>de</strong><br />
Duhamel pour le système (B.2.9) -(B.2.10); elle donne, pour tout t > 0 : pour tout t > 0 :<br />
⎧<br />
⎪⎨<br />
⎪⎩<br />
u(t) = U0,<br />
v(t) = 1<br />
ε A(U0)<br />
t<br />
e<br />
0<br />
− t<br />
τ<br />
α(s)<br />
ε ds α(τ)dτ .<br />
(B.2.11)<br />
Ainsi, grâce à l’hypothèse (B.1.6) sur le gradient <strong>de</strong> R, la fonction α reste bornée pour<br />
tout t > 0 (quel que soit le paramètre σ) comme suit :<br />
β0(U0) ≤ α(t) ≤ h(U0). (B.2.12)<br />
Choisissons donc la valeur <strong>de</strong> R0, donc celle <strong>de</strong> U0 = 1<br />
√ a R0. Comme les conditions initiales<br />
(B.1.2) <strong>de</strong> notre problème <strong>de</strong> Cauchy vérifient<br />
il vient pour les variables diagonales :<br />
|u0(x)|, |v0(x)| N0,<br />
|w0(x)|, |z0(x)| (1+ √ a)N0.