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 />
Annexe A<br />
Compléments sur l’étu<strong>de</strong> du modèle<br />
multicouche<br />
Dans c<strong>et</strong>te annexe, nous regroupons quelques compléments sur l’étu<strong>de</strong> <strong>de</strong> notre nouveau<br />
modèle multicouche <strong>de</strong> type Saint-Venant. Nous établissons à la Section A.1 une<br />
inégalité d’énergie naturelle associée à notre modèle (en 1D), <strong>et</strong> expliquons les difficultés à<br />
généraliser l’estimation <strong>de</strong> BD-entropie établie par D. Bresch <strong>et</strong> B. Desjardins [43] pour un<br />
modèle classique <strong>de</strong> Saint-Venant. Enfin, nous proposons quelques <strong>simulations</strong> <strong>numériques</strong><br />
supplémentaires en 1D, notamment <strong>de</strong> solutions discontinues, à la Section A.2.<br />
A.1 Sur l’énergie du système multicouche en 1D<br />
Dans c<strong>et</strong>te section, afin <strong>de</strong> simplifier l’écriture, nous revenons à notre modèle multicouche<br />
1D <strong>et</strong> considérons une bathymétrie triviale ainsi que la définition <strong>de</strong>s vitesses u i+1/2<br />
comme moyenne arithmétique <strong>de</strong> ui+1 <strong>et</strong> ui (voir la remarque ??). L’équation « générique »<br />
pour la couche i (1 i N) s’écrit donc<br />
<br />
∂t(hiui)+∂x hiu 2 i +ghi∂xH = µ∂x(hi∂xui)<br />
avec les conventions :<br />
⎧<br />
⎪⎨<br />
⎪⎩<br />
u 1/2 = 0, u i+1/2 = ui+1 +ui<br />
2<br />
+u i−1/2w i−1/2 −u i+1/2w i+1/2<br />
−2µ ui −ui−1<br />
+2µ<br />
hi +hi−1<br />
ui+1 −ui<br />
,<br />
hi +hi+1<br />
pour 1 ileqslantN −1, u N+1/2 = 0,<br />
w1/2 = 0, wi+1/2 = −h <br />
i∂xuk pour 1 ileqslantN −1, wN+1/2 = 0.<br />
k=1<br />
83<br />
(A.1.1)<br />
(A.1.2)