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

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)

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

Saved successfully!

Ooh no, something went wrong!