Visualiser les détails du résultat
| Identifiant | Projet | Catégorie | Visibilité | Date de soumission | Dernière mise à jour |
|---|---|---|---|---|---|
| 0000453 | SIC | Exe Sirene | public | 2011-07-21 09:46 | 2011-11-22 10:11 |
| Rapporteur | dorch | Assigné à | utilisateur6 | ||
| Priorité | normale | Impact | critique | Reproductibilité | toujours |
| Statut | fermé | Résolution | corrigé | ||
| Version du produit | 5.20h | ||||
| Version ciblée | 5.23a | Résolu dans la version | 5.22e | ||
| Résumé | 0000453: CL Q(Z) : Oscillations du débit et de la cote en calcul linéaire | ||||
| Description | Sur un régime hydraulique permanent, sur une condition limite Q(Z) basée sur la hauteur normale, le débit et la cote à l'aval du modèle se mettent à osciller jusqu'à provoquer des écoulements torrentiels ou une mise à sec dans le bief et arrêter la simulation. Le problème est résolu lorsqu'on passe en non linéaire comme il est conseillé de le faire dans le cas de l'utilisation d'une Q(Z). Dans l'ancienne version de SIRENE (<5.12), ces oscillations n'avaient pas lieu. | ||||
| Balises | Aucune balise n'est associée. | ||||
| Temps projeté (jours) | |||||
|
|
|
|
|
Les problèmes rencontrés sont en partie résolus par 0000496. Cependant, on revient à la situation existante avant 0000303 (SIC 5.21a), avec un calcul de dérivée simple entre deux points et les problèmes d'instabilités possibles autour d'un point avec une forte rupture de pente. L'exemple joint semble cependant démontrer que ce n'est pas une oscillation de la dérivée qui serait en cause mais plutôt un problème de RST sur la CL aval. JPB conseille de vérifier les variables en cause (SNOO, TNOO...) et de les comparer avec le calcul analytique à réaliser sur la Q(Z) aval : R*dQ+S*dZ=T dQ=df1*dz dZ=T/(R*df1+S) |
|
|
JPB : Mauvaise utilisation des équations sur les noeuds sans casier. La résolution de la Q(Z) revient à chercher la valeur nulle d'une équation. Avec les casiers, c'est la valeur nulle de l'intégrale de cette équation que l'on recherche et cette intégrales est approximée par un trapèze. Le problème est que l'on résoud l'intégrale à chaque fois y compris pour les noeuds sans casier et que la somme des deux termes de l'intégrale s'annulent mais pas les termes eux-même. Visiblement les deux termes changent de signe à chaque pas de temps ce qui crée les oscillations constantes. |
|
|
Modification du calcul de TNOO dans SIR38 : IF(tNd(N)%CasZS%nbrPt > 0) THEN !JPB le 21/11/2011 !Casier (intégrale par trapèze) TNOO(N)=TNOA(N)+TNO(N)+TNO1A(N)/DT-TNO1(N)/DT ELSE !Pas de casier (fonction) TNOO(N)=TNO(N) ENDIF |
| Date de modification | Nom d’utilisateur | Champ | Changement |
|---|---|---|---|
| 2011-07-21 09:46 | dorch | Nouveau bogue | |
| 2011-07-21 09:46 | dorch | Statut | nouveau => affecté |
| 2011-07-21 09:46 | dorch | Assigné à | => utilisateur6 |
| 2011-07-21 11:17 | dorch | Description mise à jour | Voir les révisions |
| 2011-07-30 23:48 | dorch | Version ciblée | => 5.22e |
| 2011-10-21 17:14 | dorch | Relation ajoutée | dépendant de 0000496 |
| 2011-10-21 17:14 | dorch | Fichier ajouté: 1Bief.xml | |
| 2011-10-21 17:47 | dorch | Note ajoutée: 0000565 | |
| 2011-11-09 15:13 | dorch | Version ciblée | 5.22e => 5.23a |
| 2011-11-21 11:02 | dorch | Note ajoutée: 0000584 | |
| 2011-11-22 10:11 | dorch | Note ajoutée: 0000585 | |
| 2011-11-22 10:11 | dorch | Statut | affecté => fermé |
| 2011-11-22 10:11 | dorch | Résolution | ouvert => corrigé |
| 2011-11-22 10:11 | dorch | Résolu dans la version | => 5.22e |
| 2012-07-11 08:25 | dorch | Catégorie | Sirene => Exe Sirene |
