SIC^2 : Logiciel de Simulation Intégrée des Canaux et de leur Contrôle

Les conditions aux limites

David Dorchies — 2017-09-24T10:20:43Z

Les prises sont obligatoirement pourvues d'une condition aux limites. Les conditions aux limites disponibles sont :

Débit imposé (Pompe) : Q(t)
Cote fixe : Z(t)
Débit fonction de la cote (tabulée) : Q(Z) (ou plutôt Z(Q))
Fonction Q puissance alpha : $Q(t)=Q_{ref} [(Z(t)-Z_{0})/(Z_{ref}-Z_{0})]^{\alpha}$

Les nœuds amont et aval sont considérés également comme des endroits où il faut définir des "prises", qui seront utilisées comme conditions aux limites dans les algorithmes de calcul hydraulique. La terminologie peut sembler abusive mais il faut bien choisir un terme générique. Cela est défendable, dans le sens où ce sont les points qui permettent d'injecter ou retirer du débit dans le réseau hydraulique, aussi bien aux nœuds amont et aval, qu'aux nœuds intermédiaires.

Les conditions aux nœuds amont et aval peuvent être choisis dans cette liste. Les détails donnés pour les noeuds aval sont donc applicables aussi ici pour un noeud amont ou intermédiaire.

Toutes les combinaisons sont possibles en régime transitoire.

Quelques combinaisons sont impossibles en régime permanent : Q(t) à l'amont et à l'aval. En effet dans ce cas soit il n'y a pas de solution, soit il y a une infinité de solutions. La configuration la plus classique est Q(t) à l'amont et Q(Z) à l'aval. Les autres options qui n'indiquent pas le débit à l'amont de manière explicite impliqueront des itérations du programme de calcul Fluvia afin de trouver le débit qui satisfait la conditions amont choisie.

Rappel : les débits entrants ou sortants, dans SIC, suivent une convention de signe qui est que les débits entrants sont positifs et les débits sortants sont négatifs. Cela est valable généralement aux conditions aux limites externes (amont et aval), internes (prises) et les infiltrations ou apports répartis le long des sections. Mais, pour éviter les erreurs classiques qu'on rencontrait souvent avec des utilisateurs, on fait des exceptions à cette convention quand il n'y a pas d'ambiguïté. Par exemple une condition limite aval de type Q(Z) ne peut être que sortante, et on attend donc des débits positifs aux interfaces de saisie des données alors que selon la convention il faudrait les choisir négatifs. La calculette les choisira positifs, mais ce seront bien des débits sortants.

Initialisation du calcul aux diffluences

David Dorchies — 2017-09-24T09:57:24Z

Sur un réseau pourvu de diffluences (cas courant sur les canaux d'irrigation contrairement aux rivières), le calcul de la ligne d'eau en régime permanent effectue des itérations pour distribuer le débit provenant de l'Amont dans chaque Bief Aval diffluent.

Par défaut, la répartition est initialisée à partir de la surface des biefs situés à l'aval de la diffluence. Dans le cas où cette répartition par défaut provoquerait des débits initiaux négatifs dans certaines Branches du réseau, l'utilisateur peut modifier ces paramètres en double-cliquant sur le Nœud posant problème et modifier la répartition du débit entre les Biefs à l'Aval du Nœud.

Définition du mode de calcul au noeud

David Dorchies — 2017-09-24T09:46:46Z

La résolution du calcul de la ligne d'eau aux nœuds peut se faire selon deux modes :

égalité des cotes (mode par défaut)
égalité des charges

Dans le cas de l'égalité des charges, il est possible de définir une perte de charge à l'entrée du nœud (pour les biefs Amont arrivant à ce nœud) et à la sortie du nœud (pour les biefs Aval quittant ce nœud). Il est possible d'avoir pour le même nœud des biefs arrivant à ou partant de ce nœud avec des modes différents (cote pour certains et charge pour d'autres). Dans ce cas d'une perte de charge une équation de type Borda (proportionnelle à la vitesse au carré) est utilisée.

Remarque : Dans le cas où l'écoulement à un nœud est en régime torrentiel sur un ou plusieurs des biefs arrivant ou partant de ce nœud, le mode de calcul sera basculé automatiquement en mode "égalité des charges", même si c'était le mode "égalité des cotes" qui avait été choisi ici. C'est en effet plus réaliste physiquement. Pour que les calculs en régime torrentiel soient autorisés il faut activer l'option correspondante dans les paramètres de calcul permanent.

Format du fichier .par

David Dorchies — 2016-10-14T16:56:32Z

Un fichier des variables paramétrées est par créé par scénario (et contenant les informations également pour ses variantes éventuelles) et est par défaut mis à jour automatiquement lors de l'enregistrement du scénario dans EdiSic. Le fichier se nomme [Nom du projet XML]_[N° de scénario].par.

Dans le nom du fichier ci-dessus, nous avons fait le choix de conserver le nom du projet initial [Nom du projet XML] dans le nom, même si le nom du projet est ensuite modifié, que ce soit par la commande "enregistrer sous" d'EdiSic, ou à la main sous Windows ou sous Dos. Si vous avez une routine générant ce fichier pour des exécutions automatiques, vous n'aurez donc pas à la modifier même si le nom du fichier projet est ensuite modifié.

C'est un fichier texte contenant des lignes et des colonnes. Les séparateurs acceptés pour les colonnes sont les espaces, la tabulation, le point-virgule et la virgule. Le format est le suivant :

Les lignes ne commençant ni par « L » ni par « X » sont ignorées. Par convention, les lignes de commentaires commencent par le double symbole « // ».
Les lignes commençant par X1, X2, X3, ... Xn servent à définir la valeur de la variable paramétrée
Les lignes commençant par L1, L2, L3, ... Ln servent à définir la localisation d'une variable paramétrée dans le réseau

EdiSic permet de définir la localisation des variables paramétrées à partir des interfaces de saisie du mode hydraulique (cf. gestion des variables paramétrées). Mais il est aussi possible de définir de nouvelles variables paramétrées directement dans le fichier .par en définissant leur localisation dans le réseau à l'aide des lignes commençant par L.

Une fois que des variables ont été définies comme variables paramétrées dans l'interface EdiSic, des noms (X1, X2, etc) sont affectés automatiquement en prenant le premier nom disponible non utilisé, par ordre croissant. Il n'est pas possible ensuite de les modifier dans la version actuelle d' EdiSic (sauf éventuellement en retouchant directement le fichier xml à la main).

Syntaxe des lignes commençant par X

Pour une valeur fixe, la syntaxe sera :
X[n] [Valeur]
Exemple :
X1 0.25

Par une valeur évoluant en fonction du temps (Loi fonction du temps), la syntaxe sera

X[n] LOI [E ou R]

X[n] [Temps en secondes] [Valeur]

X[n] [Temps en secondes] [Valeur]

...

Le paramètre [E ou R] permet de choisir entre le mode échelon [E] ou le mode rampe [R].

Exemple pour un tour d'eau évoluant à 0, 1 et 2 heures :

// Loi fonction du temps en mode "E" Echelon

X1 LOI E

X1 0 -0.1

X1 3600 -0.2

X1 7200 -0.1

Syntaxe des lignes commençant par L

Les lignes commençant par L contiennent un ensemble de colonnes avec des couples [Type d'objet]=[Valeur]. La signification des types d'objet est la suivante :

BF : n° du Bief
ND : n° du Noeud
PR : n° de la Prise
ST : n° de la Structure d'ouvrage
OUV : n° d'Ouvrage
SN : n° de Section
PBF : n° de la Portion de bief (pour les Strickler et les infiltrations)
CAR : Variable à paramétrer sur l'objet (voir la liste des variables ci-dessous)

Le numéro des biefs, des noeuds et des sections est visible dans l'infobulle qui apparaît lorsqu'on passe la souris sur les différents objets du réseau dans la vue arborescente d'EdiSic. Le fichier .par écrit par EdiSic comporte la ligne « L » en commentaire pour chaque variable définie à titre d'exemple de syntaxe.

Il est possible de définir plusieurs localisations pour un seule variable. Un avertissement sans conséquence s'affichera lors de la simulation et la même valeur sera appliquée aux différentes localisation demandée.

Voici différents exemple de localisation pour différents objets du réseau :

Type d'objet	Syntaxe de la ligne
Débit de la prise n°1 du noeud n°1	`L1 ND=1 PR=1 CAR=Q`
Ouverture de la vanne de l'ouvrage du premier niveau de structure de la prise n°1 du noeud n°2	`L2 ND=2 PR=1 ST=1 OUV=1 CAR=W`
Coefficient de débit de la vanne de l'ouvrage n°1 de la structure en travers située dans la section n°4 du bief n°1	`L3 BF=1 SN=4 OUV=1 CAR=CoefQR`
Strickler de la première portion de Strickler du bief n°1	`L4 BF=1 PBF=1 CAR=KMin`

La liste des variables disponibles par le paramètre CAR est la suivante (attention les minuscules et majuscules sont importants) :

Variable	Signification	Localisation possible
Q	Débit ou débit objectif	Prise avec une condition limite en débit ou avec un débit objectif contrôlant un ouvrage (`ND=... PR=... CAR=Q`)
Z	Cote de l'eau	Prise avec une condition limite en cote (`ND=... PR=... CAR=Z1`)
KMin	Strickler du lit mineur	Portion de bief (`BF=... PBF=... CAR=KMin`)
KMoy	Strickler du lit moyen	Portion de bief (`BF=... PBF=... CAR=KMoy`)
Inf	Infiltration	Portion de bief (`BF=... PBF=... CAR=Inf`)

Aux interfaces d'EdiSic les frottements sont indiqués en Strickler en Français et en Manning en Anglais et en Espagnol. Mais de manière interne, ces frottements sont écrits dans les fichiers .xml ainsi que dans les fichiers .par en Strickler. Si vous êtes plutôt anglo-saxon et avez l'habitude de manipuler des Manning il faudra cependant penser à les convertir (Ks = 1/n) en Strickler si vous les écrivez dans un fichier .par.

On pourra aussi piloter des paramètres d'ouvrages en travers (BF=... SN=... OUV=...) et les ouvrages de prises (ND=... PR=... ST=... OUV=...). Les variables disponibles sont les suivantes :

Variable	Signification	Type d'ouvrage possible
CoteRadier	Cote du radier	Tous les ouvrages
Largeur	Largeur	Toutes les vannes et seuils sauf vanne circulaire et tube
Ouverture	Ouverture de vanne	Toutes les vannes sauf GEC-Alsthom
CoefQR	Coefficient de débit rectangulaire	Toutes les équations de vannes et seuil sauf Goussard
SurverseHauteur	Hauteur de la pelle	Toutes les vannes
CoefQSurverse	Coefficient de débit de la surverse	Toutes les vannes avec surverse au dessus de la pelle
TanAl	Fruit des berges	Vannes et seuils trapézoïdaux, vannes GEC-Alsthom
CoefQT	Coefficient de débit triangulaire	Vannes et seuils trapézoïdaux et vannes GEC avec équation CEM02
CoteAxe	Cote de l'axe	Vannes GEC-Alsthom
Rayon	Rayon de la vanne	Vannes GEC-Alsthom
D	Caractéristique "D" de la loi de Goussard	Vannes AMIL
JMax	Perte de charge maximale	Vannes GEC-Alsthom
S1S2	Rapport des surfaces pour le filtrage	Vannes AVIS, AVIO, Mixtes
Decal	Décalage de prise	Vannes AVIS, AVIO, Mixtes
Decrement	Décrément	Vannes AMIL, AVIS, AVIO
CoteAmont	Cote amont max.	Vannes AMIL

Exemple de fichier .par

// Localisation définie par EdiSic

//Nd : OFF 3 > Pr : Prise1 > CL Débit 

// L1 Nd=3 Pr=1 Car=Q

// Loi fonction du temps en mode "E" Echelon

X1 LOI E

X1 0 -0.1

X1 3600 -0.2

X1 7200 -0.1

// Ajout par l'utilisateur de Nd : DAM > Pr : Prise 1 > CL Débit

L2 Nd=1 pr=1 car=Q

// Loi fonction du temps en mode "R" Rampe

X2 LOI R

X2 0 2.5

X2 3600 3

X2 7200 2.5 

// Ajout par l'utilisateur des Nd 3 et 5 Prise 1 avec débit fixe de 0.4 m3/s

L3 Nd=3 Pr=1 Car=Q

L3 Nd=5 Pr=1 Car=Q

X3 -0.4

Les vannes GEC-Alsthom (AMIL, AVIS, AVIO, Mixtes)

David Dorchies — 2015-04-20T15:48:25Z

Les vannes de régulation AMIL, AVIS et AVIO et mixtes peuvent être modélisées directement sans faire appel à un module de régulation. Cela présente l'avantage de modéliser ce type d'ouvrage dès le calcul et régime permanent et d'avoir une modélisation beaucoup plus précise.

Lors de la création de ces vannes (à chaque modification de modèle ou à chaque modification de rayon de la vanne) on indique, à titre d'information, la perte de charge maximale admissible Jmax, le débit correspondant Q pour Jmax, le débit maximum Qmax qui peut passer à travers cette vanne ainsi que la charge correspondante J pour Qmax. Pour les vannes Amil cette perte de charge maximum Jmax n'est pas pertinente (pas d'indication du constructeur) et est donc remplacée par la perte de charge Jlimite correspondant à la transition noyée-dénoyée. Cette indication permet de choisir la vanne adaptée aux conditions hydrauliques locales du canal modélisé.

Les abaques des vannes, que nous avons régénérées à partir des équations que nous avons programmées sont fournies en Annexe. On peut vérifier qu'elles sont identiques à celles d'origine fournies par le constructeur Gec-Alsthom (maintenant Hydrostec).

Type

Le "type" correspond à la gamme des vannes choisies (Amil, Avio Haute Chute, Avio Basse Chute, Avis Haute Chute, Avis Basse Chute, Mixte sans masque d'ouverture max 45°, Mixte sans masque d'ouverture max 55°, Mixte avec masque).

Référence

La "référence" correspond à celle du catalogue constructeur Gec-Alsthom. Si vous voulez modéliser une vanne différente vous pouvez sélectionner l'option "autre". Dans ce cas nous utilisons certaines dimensions nécessaires pour le calcul déterminées à partir du rayon R indiqué en s'inspirant des valeurs proches des vannes classiques :

Amil :
- D=1.6*R
- E=0.95*R (hauteur du tablier pour calcul des débordements éventuels)
- b=0.91*R (largeur du tablier en bas)
Avis Haute Chute :
- E=0.81*R (hauteur du tablier pour calcul des débordements éventuels)
- b=1.19*R (largeur du tablier en bas)
Avis Basse Chute :
- E=0.95*R (hauteur du tablier pour calcul des débordements éventuels)
- b=1.19*R
Avio Haute Chute :
- L=0.5*R (largeur buse)
Avio Basse Chute :
- L=R (largeur buse)

Décrément

Le "décrément" est la différence entre les niveaux minimum et maximum régulés. Pour des raisons de stabilité ce décrément ne peut être nul. Pour une vanne AMIL la valeur par défaut proposée par le constructeur (et par le logiciel SIC) est de d=0.02*D. Pour une vanne AVIS, la valeur par défaut est d=0.028*R (basse chute) et d=0.032*R (haute chute). Pour une vanne AVIO, la valeur par défaut est d=0.028*R (basse chute et haute chute).

Cote de l'axe

La "cote de l'axe" est la cote de l'axe de rotation de la vanne. Elle correspond à la cote de consigne pour Q=0 pour les vannes AVIS, AVIO et les vannes AMIL en réglage alternatif (ou pour Q=Qmax pour une AMIL en réglage classique). La cote de l'axe par défaut proposée par le constructeur est calculée et proposée par le logiciel.

Cote de radier

La "cote de radier" est la cote du fond de l'ouvrage de génie civil à l'endroit où la vanne bute en fermeture complète.

Rayon

Le "rayon" est le rayon du tablier de la vanne. Il est proposé automatiquement en fonction de la référence de la vanne choisie, mais il peut être modifié. C'est un paramètre essentiel pour le calcul du fonctionnement de la vanne.

Équations de débit

L'équation "GOU 93" correspond à l'équation déduite d'un document fourni par Goussard, ancien ingénieur de Gec-Alsthom, et rédigé en 1993.

L'équation "CEM 02" correspond à une équation mise au point par le Cemagref en 2002.

L'ouverture maximum

La variable "ouverture maximum" correspond à l'ouverture réelle bloquée de la vanne si l'équation choisie est "CEM 02 (V)" et correspond à l'ouverture maximum atteinte pour le réglage lié au décrément de la vanne si l'équation choisie est "GOU 93". Une valeur par défaut correspondant à la butée mécanique est proposée par défaut en fonction du type et de la référence de la vanne. Cette valeur peut être modifiée si une butée différente limite ou augmente l'ouverture maximum atteignable.

Pour une vanne Amil en réglage classique la vanne sera à cette ouverture max dès que le niveau de l'eau amont sera >= à la cote de l'axe. Pour une vanne Amil en réglage alternatif la vanne sera à cette ouverture max dès que le niveau de l'eau amont sera >= à la cote de l'axe + décrément. Pour une vanne Amil en réglage intermédiaire la vanne sera à cette ouverture max dès que le niveau de l'eau amont sera >= à la cote de l'axe + Z amont max (cf. ci dessous pour la définition de cette variable).

Pour les vannes Avis et Avio, la vanne sera à cette ouverture max dès que le niveau de l'eau aval sera <= à la cote de l'axe - décrément. Modifier cette valeur revient à modifier la butée contrôlant l'ouverture maximum de la vanne et à modifier l'équilibrage de la vanne.

Ouverture minimum

Par défaut les vannes GEC ont par défaut une fermeture minimum égale à 0.5% du rayon de la vanne. Cela permet de mimer le fait que ces vannes ne sont pas parfaitement étanches et facilite le calcul en régime permanent pour les vannes AVIS et AVIO car une vanne fermée complètement générera une perte de charge très importante qui se répercutera d'aval en amont.

Z amont max

La variable "Z amont max" correspond à la hauteur d'eau max par rapport à la cote de l'axe (correspond au réglage dans le cas de la vanne AMIL) : 0 <= Zmax <= d. Zmax=0 pour le réglage classique et Zmax=d pour le réglage alternatif. Les réglages intermédiaires sont possibles. Pour le réglage classique les niveaux amont vont donc varier entre Z=cote de l'axe - d pour Q=0 et Z=cote de l'axe pour Q=Qmax. Pour le réglage alternatif les niveaux amont vont donc varier entre Z=cote de l'axe pour Q=0 et Z=cote de l'axe + d pour Q=Qmax. Cette variable n'est pas utilisée pour les vannes autres que les vannes Amil.

Filtrage S1/S2

La variable "filtrage S1/S2" permet de filtrer la cote influençant le déplacement du flotteur de la vanne. Nous n'utilisons pas cette variable pour l'instant pour les vannes Amil. On pourrait imaginer l'utiliser pour mimer l'amortissement lié à l'amortisseur à piston installé sur ces vannes. Pour les vannes Avis et Avio, cette valeur permet de mimer le filtrage lié au bac dans lequel le flotteur est situé. Théoriquement, le coefficient S1/S2 à renseigner ici correspond au ratio entre la surface libre de l'eau dans le bac autour du flotteur divisé par la surface de l'orifice qui fait communiquer le bac au bief de canal à l'aval. Mais dans la pratique ce coefficient doit être augmenté pour tenir compte des pertes de charges supplémentaires liées aux tuyauteries éventuelles de cette prise d'eau.

Décalage prise bac

La variable "décalage prise bac" correspond à déporter la prise d'eau du bac vers le bief aval. Dans la pratique, en particulier lorsqu'il y a plusieurs vannes en parallèle, la prise d'eau du ou des bacs est reliée à un tuyau communiquant au bief aval plusieurs mètres à l'aval de la vanne dans une zone tranquillisée. Ce décalage est indiqué en mètres (SIC cherche alors la valeur de la cote de l'eau dans la section de calcul la plus proche).

Limitation pour les vannes mixtes

Il n'est possible de définir qu'une seule vanne mixte par structure (Le modèle ne permet pas de placer plusieurs vannes mixtes en parallèle).

La surverse d'une vanne (débordement au dessus d'une vanne)

David Dorchies — 2015-04-13T12:01:16Z

Pour les différentes vannes disponibles parmi les types d'ouvrages, il y a la possibilité de modéliser la surverse éventuelle se faisant par dessus l'ouvrage en cas de hautes eaux.

Pour se faire, il suffit de :

cocher la case "Surverse" au bas de l'écran,
définir la hauteur de la vanne
choisir l'équation de seuil (les même que pour un seuil rectangulaire)
définir le coefficient de débit pour l'équation de seuil

La surverse se déclenche à partir du moment où la cote de l'eau à l'amont est supérieure à la cote de radier + l'ouverture de la vanne + la hauteur de la vanne.

La largeur prise en compte pour la surverse correspond :

à la largeur de la vanne pour une vanne rectangulaire ;
à la largeur au niveau des berges pour une vanne trapézoïdale.

Si la surverse n'est pas cochée, la hauteur de la vanne est considérée comme infinie et aucune surverse n'est possible.

Les équations de débit aux ouvrages

David Dorchies — 2015-04-13T11:57:58Z

"CEM 88(V)" (Déversoir / Vanne de fond) est l'équation proposée par défaut pour une vanne rectangulaire ou circulaire . Elle prend en compte les différents écoulements : noyé, dénoyé, à surface libre et en charge. De plus, elle assure une continuité du débit entre ces différents écoulements. Elle peut aussi être utilisée dans le cas d'un seuil, c'est à dire sans possibilité de mise en charge. Elle peut être utile dans le cas d'un seuil avec une pelle faible ou nulle.
"CEM 88(D)" (Déversoir / Orifice) est l'équation proposée par défaut pour un seuil rectangulaire. Elle prend aussi en compte les différents écoulements (noyé, dénoyé, à surface libre et en charge) et assure une continuité du débit entre ces différents écoulements. Elle peut être utile dans le cas d'une vanne avec une pelle importante.
"CUN 80" sont les équations décrites par Cunge dans son livre [1], ou avec plus de détails dans un article de Mahmood et Yevjevich [2]. Les équations "CUN 80" prennent aussi en considération les différentes conditions d'écoulement : noyé, dénoyé, à surface libre et en charge (comme les équations Cemagref "CEM 88"), mais il n'y a aucune continuité entre les conditions d'écoulement à surface libre et en charge. Cela peut entraîner des problèmes de calcul quand nous sommes près de cette transition.
"Noyée" correspond à l'équation classique d'une vanne noyée en racine de cote amont - cote aval.
"Dénoyé" correspond à l'équation classique d'un seuil dénoyé en cote amont à la puissance 1.5.
"GOU 93 Amil" et "GOU 93 Avi" sont les équations proposées par défaut pour les vannes GEC-Alsthom AVIS/AVIO et AMIL. Elles sont issues de documents internes de GEC-Alsthom et permettent de reproduire les débits fournis dans les abaques du constructeur.
"CEM 02(V)" est l'équation proposée par défaut pour les vannes trapézoïdales. C'est la somme de l'équation déversoir / vanne de fond "CEM 88(V)" pour la partie rectangulaire et de l'équation de Gourley pour la partie triangulaire.
"CEM 02(D)" est l'équation proposée par défaut pour les seuils trapézoïdaux. C'est la somme de l'équation déversoir / orifice "CEM 88(D)" pour la partie rectangulaire et de l'équation de Gourley pour la partie triangulaire.

[1] Cunge, Holly, Verwey, 1980, "Practical aspects of computational river hydraulics", Pitman, p. 169 pour les seuils et p. 266 pour les vannes

[2] Mahmood K., Yevjevich V., 1975, "Unsteady flow in open channels, Volume 1 and 2", Water resources publications, Fort Collins, USA, 923 p

Les ouvrages disponibles dans SIC

David Dorchies — 2015-04-13T11:57:54Z

Les ouvrages sont définis au niveau de la structure en travers d'une section singulière ou au niveau d'une structure d'ouvrages de 1^er ou 2^ème niveau d'une prise.

Au sein d'une structure, on peut définir un ou plusieurs ouvrages placés en parallèle.

Différents types d'ouvrages peuvent être combinés au sein d'une même structure. Pour chacun de ces types d'ouvrage, SIC propose un choix d'équation de débit défini ci-dessous entre parenthèses :

Seuil rectangulaire (CEM 88 (D), Dénoyé, CUN 80, CEM 88 (V)) ;
Vanne rectangulaire (CEM 88 (V), Noyé, CUN 80, CEM 88 (D)) ;
Vanne circulaire (CEM 88 (V), Noyé, CUN 80, CEM 88 (D)) ;
Vannes auto-régulées GEC-Alsthom, AMIL, AVIO-AVIS et MIXTES (GOU 93, CEM 02 (V)) ;
Seuil trapézoïdal (CEM 02 (V), CEM 02 (D)) ;
Vanne trapézoïdale (CEM 02 (V), CEM 02 (D)) ;
Tube
APM

Exécution d'un calcul en régime transitoire (SIRENE)

David Dorchies — 2012-03-09T14:46:07Z

Pour lancer une simulation en régime transitoire à partir de l'édition d'un scénario, il suffit de cliquer sur le bouton . Cela peut aussi se faire à partir de l'explorateur de projet en sélectionnant la ligne du scénario et en cliquant sur le bouton « Calcul transitoire ». Avant le calcul, le logiciel propose d'enregistrer les modifications effectuées et effectue des tests de cohérence.

Les test de cohérence sont actuellement conçus pour le régime permanent, si les conditions aux limites et le sens de l'écoulement ne correspondent pas à celles acceptées en régime permanent, il ne faut pas tenir compte de ces avertissements.

Durant son exécution, le programme de calcul en régime transitoire (SIRENE) va :

Calculer la ligne d'eau à chaque pas de temps et enregistrer, à la fréquence définie dans les paramètres de temps, les résultats demandés dans les paramètres de calcul > Écriture des Résultats
Afficher à l'écran les avertissements et les erreurs en cours de calcul

Les lois fonction du temps

David Dorchies — 2012-03-09T14:22:04Z

Sur les conditions aux limites (débit ou cote fixé) et sur les caractéristiques d'ouvrage (largeur, cote de seuil, ouverture...), il est possible de définir des valeurs fixes ou des valeurs qui évoluent au cours du temps.

Partout où c'est possible, pour passer d'une valeur fixe à une lois fonction du temps, il suffit de cliquer sur le symbole présent à droite du champ à saisir.

On définit ensuite dans un tableau les évènements en définissant la date et la nouvelle valeur à cette date.

Il existe deux méthodes d'interpolation des données entre deux évènements :

le mode échelon
le mode rampe

Le mode échelon correspond à une manœuvre ponctuelle à un pas de temps donné. La manœuvre est effective un pas de temps après la date de l'évènement.

Le mode rampe correspond à une évolution linéaire entre deux évènements. La valeur calculée à chaque pas de temps est une interpolation entre l'évènement précédent et l'évènement suivant.

Il est possible de revenir à une valeur fixe ne dépendant pas du temps en cliquant sur le bouton « Valeur fixe ».