13.07.2015 Views

Assistance au calage de modèles numériques en hydraulique ... - TEL

Assistance au calage de modèles numériques en hydraulique ... - TEL

Assistance au calage de modèles numériques en hydraulique ... - TEL

SHOW MORE
SHOW LESS
  • No tags were found...

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

CHAPITRE 7 CALAGE D’UN MODÈLE DE L’HOGNEAUlocal pour modifier les valeurs <strong>de</strong>s coeffici<strong>en</strong>ts <strong>de</strong> résistance <strong>de</strong> ce tronçon. En revanche,le choix <strong>de</strong> conserver les valeurs initiales va nécessairem<strong>en</strong>t influer sur les résultats obt<strong>en</strong>usà l’amont <strong>de</strong> ce tronçon. Le système d’assistance <strong>au</strong> <strong>calage</strong> ne pr<strong>en</strong>d pour l’instantpas <strong>en</strong> compte ce cas <strong>de</strong> figure, et l’utilisateur est libre <strong>de</strong> fournir <strong>au</strong> système sa visionou non d’un écart sur le tronçon considéré.Nous proposons d’implém<strong>en</strong>ter le raisonnem<strong>en</strong>t suivant pour résoudre cette indétermination: les écarts loc<strong>au</strong>x doiv<strong>en</strong>t être considérés conjointem<strong>en</strong>t sur le tronçonincriminé et sur le tronçon immédiatem<strong>en</strong>t à l’amont sur lequel <strong>de</strong>s mesures ont été effectuées.L’ajustem<strong>en</strong>t <strong>de</strong>s coeffici<strong>en</strong>ts <strong>de</strong> résistance <strong>de</strong>vra alors s’effectuer d’une manièrei<strong>de</strong>ntique sur les <strong>de</strong>ux tronçons.Nous effleurons du doigt à travers cet exemple le nombre <strong>de</strong> situations différ<strong>en</strong>tesqui peuv<strong>en</strong>t se prés<strong>en</strong>ter lors d’un <strong>calage</strong> à l’ai<strong>de</strong> <strong>de</strong> données internes discrètes.7.4.5 Ajustem<strong>en</strong>t dissocié <strong>de</strong>s coeffici<strong>en</strong>ts du lit mineur et du lit majeurLes <strong>de</strong>ux cas d’étu<strong>de</strong>s évoqués dans ces trav<strong>au</strong>x n’ont pas permis d’évoquer un ajustem<strong>en</strong>tséparé <strong>de</strong>s coeffici<strong>en</strong>ts <strong>de</strong> résistance du lit mineur et du lit majeur. Cette situationse r<strong>en</strong>contre lorsque <strong>de</strong>ux jeux <strong>de</strong> nive<strong>au</strong>x d’e<strong>au</strong> sont disponibles sur le bief modélisé,l’un pour lequel l’écoulem<strong>en</strong>t est cont<strong>en</strong>u <strong>en</strong>-<strong>de</strong>çà <strong>de</strong> la transition lit mineur–litmajeur, et l’<strong>au</strong>tre pour lequel l’écoulem<strong>en</strong>t se situe <strong>au</strong>-<strong>de</strong>ssus <strong>de</strong> cette transition. Il f<strong>au</strong>tnoter l’importance <strong>de</strong> la définition <strong>de</strong> ce nive<strong>au</strong> <strong>de</strong> transition, qui ne correspond pasnécessairem<strong>en</strong>t <strong>au</strong> nive<strong>au</strong> <strong>de</strong> débor<strong>de</strong>m<strong>en</strong>t dans une plaine d’inondation. En effet, c<strong>en</strong>ive<strong>au</strong> peut interv<strong>en</strong>ir be<strong>au</strong>coup plus bas, pour distinguer la résistance du lit même<strong>de</strong> celle <strong>de</strong>s berges. Dans le cas où le modèle doit <strong>au</strong>ssi pr<strong>en</strong>dre <strong>en</strong> compte une plained’inondation active – c’est-à-dire dans laquelle les vitesses ne sont pas négligeables – ilserait nécessaire <strong>de</strong> considérer un troisième coeffici<strong>en</strong>t <strong>de</strong> résistance par section, ce qu<strong>en</strong>e permet pas actuellem<strong>en</strong>t le co<strong>de</strong> MAGE.Dans une telle situation, le prototype CARMA-1 réalise tout d’abord l’ajustem<strong>en</strong>t<strong>de</strong>s coeffici<strong>en</strong>ts du lit mineur à l’ai<strong>de</strong> <strong>de</strong>s mesures effectuées à faible débit, puis <strong>en</strong>suitel’ajustem<strong>en</strong>t <strong>de</strong>s coeffici<strong>en</strong>ts du lit majeur sur les nive<strong>au</strong>x mesurés à fort débit,<strong>en</strong> conservant constant les coeffici<strong>en</strong>ts du lit mineur. Cette approche permet ainsi <strong>de</strong>garantir <strong>de</strong> meilleures performances du modèle dans toute la gamme <strong>de</strong>s écoulem<strong>en</strong>ts.7.5 Expérim<strong>en</strong>tation d’une optimisation numériqueAu travers <strong>de</strong> ce cas d’étu<strong>de</strong>, nous avons pu voir que l’ajustem<strong>en</strong>t <strong>de</strong>s coeffici<strong>en</strong>ts<strong>de</strong> résistance sur la base d’une comparaison visuelle pouvait s’avérer long et fastidieux(voir annexe H, section H.2). Nous avons donc décidé <strong>de</strong> comparer notre approchevisuelle à une approche par optimisation numérique. Les sections suivantes prés<strong>en</strong>t<strong>en</strong>tbrièvem<strong>en</strong>t le principe <strong>de</strong> l’algorithme utilisé et les résultats obt<strong>en</strong>us.7.5.1 Algorithme du recuit simuléJean-Baptiste FAURE a implém<strong>en</strong>té un algorithme <strong>de</strong> recuit simulé 2 <strong>au</strong>tonome.Nous avons utilisé cet algorithme pour réaliser une optimisation <strong>de</strong>s coeffici<strong>en</strong>ts <strong>de</strong>résistance du modèle numérique <strong>de</strong> l’Hogne<strong>au</strong>. L’idée principale <strong>de</strong> cet algorithme,2. En anglais : simulated annealing.185

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

Saved successfully!

Ooh no, something went wrong!