Afficher cette notice en
1 : $$$$ DARCYTRA NOTICE CHAT 11/09/12 21:15:43 7124 2 : DATE 11/09/12 3 : 4 : Procedure DARCYTRA Voir aussi : PECHE 5 : ------------------ DARCYSAT 6 : DARCYTRA TAB1 ; 7 : 8 : TAB1.'SOUSTYPE'.'MODELE'.'DOMAINE'. 9 : 'CARACTERISTIQUES'.'EMMAGASINEMENT'.'CONVECTION'. 10 : 'TEMPS'.'TRACE_CHARGE'.'CHARGE'.'FLUX'. 11 : 'BLOCAGE'.'TRACE_IMPOSE'.'FLUX_IMPOSE'.'SOURCE'. 12 : 'TEMPS_CALCULES'.'TEMPS_SAUVES'. 13 : 'THETA'.'THETA_CONVECTION' 14 : 15 : ou TAB1.'SOUSTYPE'.'MODELE'.'DOMAINE'.'ORIENTATION'. 16 : 'CARACTERISTIQUES'.'POROSITE'.'DECROISSANCE'. 17 : 'COEF_RETARD'.'LANGMUIR'.'FREUNDLICH'. 18 : 'LIMITE_SOLUBILITE'.'COEF_DISSOLUTION'. 19 : 'CONVECTION'.'TEMPS'.'TRACE_CONC'.'CONCENTRATION'. 20 : 'FLUX'.'PRECIPITE'.'DISSOLUTION'. 21 : 'BLOCAGE'.'TRACE_IMPOSE'.'FLUX_IMPOSE'. 22 : 'DISSOLUTION_IMPOSEE'.'SOURCE'.'TEMPS_CALCULES'. 23 : 'TEMPS_SAUVES'.'THETA_DIFF'.'THETA_CONVECTION'. 24 : 'THETA_DEC'.'THETA_DISS'.'PENALISATION'. 25 : 'EPSI_LIM'.'ITMAX_LIM'.'EPSI_RET'.'EPSI_COR'. 26 : 'ITMAX_RET' 27 : 28 : 29 : 30 : Objet : 31 : _______ 32 : 33 : Cette procedure a deux fonctions. 34 : 35 : 1) En presence de l'indice 'CHARGE', on resoud les equations 36 : de DARCY en transitoire pour l'ecoulement par une methode d'elements 37 : finis mixtes hybrides (EFMH). 38 : Les inconnues du probleme sont la charge ('H'), 39 : la trace de charge ('TH') et le flux diffusif ('FLUX'). 40 : 41 : 2) En l'absence de l'indice 'CHARGE', resoud l'equation 42 : de transport par diffusion-convection d'un champ scalaire actif 43 : par un fluide dont la vitesse est connue. L'espece peut se trouver 44 : sous trois formes : solute, adsorbat et precipite, dont les lois 45 : d'echange doivent etre specifiees. A chaque loi correspondent un ou 46 : plusieurs algorithmes auxquels des parametres numeriques doivent etre 47 : fournis. On utilise la modelisation Darcy EFMH. 48 : Les inconnues du probleme sont la concentration ('H'), 49 : la trace de concentration ('TH') et le flux diffusif ('FLUX'). 50 : 51 : 52 : Commentaire : 53 : _____________ 54 : 55 : En entree, TAB1 sert a definir les options et les parametres du 56 : calcul. Les indices de la table TAB1 sont des mots (a coder tel quel) 57 : dont voici la description : 58 : 59 : ___________________________________________________________________ 60 : | | 61 : | Indice Contenu | 62 : | | 63 : ------------------------------------------------------------------- 64 : | | 65 : |------------------------------------------------ | 66 : |Donnees physiques, geometriques et materielles : | 67 : |------------------------------------------------ | 68 : | | 69 : | ------ Indices communs a l'ecoulement et au transport ------ | 70 : | ------------------------------------------------------------ | 71 : | | 72 : |'SOUSTYPE' 'DARCY' (type MOT) | 73 : | | 74 : |'MODELE' Objet modele (MMODEL cree par MODE) | 75 : | | 76 : |'DOMAINE' References geometriques (TABLE creee par DOMA) | 77 : | | 78 : | ------ 1ere possibilite : Resolution de l'ecoulement ------ | 79 : | ----------------------------------------------------------- | 80 : | | 81 : |'CARACTERISTIQUES' Donnees physiques et materielles : | 82 : | conductivite hydraulique (CHAMELEM cree par MATE) | 83 : | | 84 : |'EMMAGASINEMENT' Valeur du coefficient d'emmagasinement | 85 : | (Type CHPO Centre, Comp 'CK', ou FLOTTANT) | 86 : | - Defaut 1. | 87 : | | 88 : | ------ 2eme possibilite : Resolution du transport ------ | 89 : | ------------------------------------------------------------ | 90 : | | 91 : |'CARACTERISTIQUES' Donnees physiques et materielles : | 92 : | diffusivite effective (CHAMELEM cree par MATE) | 93 : | | 94 : |'POROSITE' Valeur de la porosite (Type CHPO Centre, Comp | 95 : | 'CK', ou FLOTTANT) - Defaut 1. | 96 : | | 97 : |'DECROISSANCE' Valeur du terme de decroissance (Type FLOTTANT) | 98 : | Tel que dC/dt = - Lambda * C - Defaut 0. | 99 : | | 100 : |'COEF_RETARD' Coefficient de retard lineaire dans le cas simple, | 101 : | ou Pente a l'origine de la fonction F(C) dans le | 102 : | cas d'isotherme non lineaire de Langmuir | 103 : | ou Coefficient K de l'isotherme de Freundlich | 104 : | (Type CHPO Centre 'SCAL', ou FLOTTANT) | 105 : | | 106 : |'LANGMUIR' Quantite maximale adsorbee sur le solide | 107 : | rapportee a l'unite de volume du fluide et exprimee| 108 : | dans la meme unite que le solute. | 109 : | (Type CHPO Centre 'SCAL', ou FLOTTANT). | 110 : | F = (R-1) C / [1 + ((R-1) C / Fsat)] | 111 : | Si cet indice et le suivant sont absents, | 112 : | l'equilibre d'adsorption est lineaire. Cet indice a| 113 : | priorite sur l'indice FREUNDLICH. | 114 : | | 115 : |'FREUNDLICH' Exposant de la loi de Freundlich F = K (C ^ 1/n) | 116 : | (Type FLOTTANT). | 117 : | Dans ce cas (et si l'indice LANGMUIR n'existe pas),| 118 : | l'indice 'RETARD' contient le coefficient | 119 : | K ramene a une unite de volume de fluide. | 120 : | - Non disponible pour l'instant - | 121 : | | 122 : |'LIMITE_SOLUBILITE' Limite de solubilite (Type MCHAML, Comp 'H') | 123 : | | 124 : |'COEF_DISSOLUTION' Coef. de dissolution (Type CHPO Centre, Comp | 125 : | 'SCAL'). Tel que dC/dt = Coef * (Csat - C) - Par | 126 : | defaut, la dissolution est instantanee | 127 : | | 128 : |'CONVECTION' Debit integre de la vitesse convective a travers | 129 : | chaque face des elements (Type CHPO Face, comp. | 130 : | 'FLUX') | 131 : | | 132 : |---------------------- | 133 : |Conditions initiales : | 134 : |---------------------- | 135 : | | 136 : | ------ Indices communs a l'ecoulement et au transport ------ | 137 : | ------------------------------------------------------------ | 138 : | | 139 : |'TEMPS' TABLE contenant a l'indice 0 la valeur du temps | 140 : | initial (FLOTTANT) | 141 : | | 142 : | | 143 : | ------ 1ere possibilite : Resolution de l'ecoulement ------ | 144 : | ----------------------------------------------------------- | 145 : | | 146 : |'CHARGE' TABLE contenant a l'indice 0 la charge hydraulique | 147 : | (quantite d'element par unite de volume d'eau) | 148 : | (Type CHPO Centre, Comp 'H') | 149 : | | 150 : |'TRACE_CHARGE' TABLE contenant a l'indice 0 la trace de | 151 : | charge initiale (CHPO, 'TH') | 152 : | | 153 : |'FLUX' TABLE contenant a l'indice 0 le flux hydraulique | 154 : | initial integre sur chaque face (Type CHPO Face, | 155 : | comp. 'FLUX') | 156 : | | 157 : | ------ 2eme possibilite : Resolution du transport ------ | 158 : | ------------------------------------------------------------ | 159 : | | 160 : |'CONCENTRATION' TABLE contenant a l'indice 0 la concentration | 161 : | (quantite d'element par unite de volume d'eau) | 162 : | (Type CHPO Centre, Comp 'H') | 163 : | | 164 : |'TRACE_CONC' TABLE contenant a l'indice 0 la trace de | 165 : | concentration initiale (Type CHPO Face, Comp 'TH') | 166 : | | 167 : |'FLUX' TABLE contenant a l'indice 0 le flux diffusif | 168 : | initial integre sur chaque face (Type CHPO Face, | 169 : | comp. 'FLUX') | 170 : | | 171 : |'PRECIPITE' TABLE contenant a l'indice 0 la quantite initiale | 172 : | de precipite par unite de volume de milieu | 173 : | (Type CHPO Centre, Comp 'H') | 174 : | | 175 : |'DISSOLUTION' TABLE contenant a l'indice 0 la quantite initiale | 176 : | pour estimer la dissolution au premier pas de temps| 177 : | (Type CHPO, Comp 'H'), voir plus loin. | 178 : | | 179 : |-------------------------------------- | 180 : |Conditions aux limites / chargements : | 181 : |-------------------------------------- | 182 : | | 183 : | ------ Indices communs a l'ecoulement et au transport ------ | 184 : | ------------------------------------------------------------ | 185 : | | 186 : |'BLOCAGE' Contient les matrices de blocage (RIGIDITE) | 187 : | | 188 : |'TRACE_IMPOSE' Valeurs des traces imposees (charge ou concentra- | 189 : | -tion) (CHARGEMENT 'TH' - Obligatoire si BLOCAGE) | 190 : | | 191 : |'FLUX_IMPOSE' Valeurs des flux imposes integres par face | 192 : | (Type CHARGEMENT de CHPO Face, comp. 'FLUX'- | 193 : | defaut 0.) | 194 : | | 195 : |'SOURCE' Valeurs du terme source par maille et par unite de | 196 : | temps (ex : puits, filiation) | 197 : | Les valeurs a l'indice i sont les valeurs entre | 198 : | les temps i-1 et i. | 199 : | (CHARGEMENT de CHPO Centre, comp 'SOUR'- defaut 0.)| 200 : | | 201 : | ------ 2eme possibilite : Resolution du transport ------ | 202 : | ------------------------------------------------------------ | 203 : | | 204 : |'DISSOLUTION_IMPOSEE' Valeurs des dissolutions imposees par unite| 205 : | de temps et par maille. (Type CHARGEMENT de CHPO, | 206 : | Comp 'H'). Les valeurs a l'indice i sont les | 207 : | valeurs moyennes de dissolution par unite de temps | 208 : | entre les temps i-1 et i. | 209 : | Priorite de la dissolution imposee sur les | 210 : | cinetiques. | 211 : | | 212 : |-------------------- | 213 : |Donnees numeriques : | 214 : |-------------------- | 215 : | | 216 : | ------ Indices communs a l'ecoulement et au transport ------ | 217 : | ------------------------------------------------------------ | 218 : | | 219 : |'TEMPS_CALCULES' Valeur des temps calcules (LISTREEL) | 220 : | Contient obligatoirement le temps final. | 221 : | | 222 : |'TEMPS_SAUVES' Valeur des temps sauvegardes (LISTREEL - defaut : | 223 : | on sauve tous les pas de temps) | 224 : | | 225 : | ------ 1ere possibilite : Resolution de l'ecoulement ------ | 226 : | ----------------------------------------------------------- | 227 : | | 228 : |'THETA' Coefficient de relaxation compris entre 0. et 1. | 229 : | (theta-methode diffusion) (FLOTTANT - defaut 1.) | 230 : | Possibilite de non-convergence lorsque theta<1/2 | 231 : | Valeurs de theta generalement utilisees : | 232 : | Schema de Euler explicite : 0. | 233 : | Schema de Crank-Nicholson : 1/2 | 234 : | Schema de Galerkin : 2/3 | 235 : | Schema de Euler implicite : 1. | 236 : | | 237 : | ------ 2eme possibilite : Resolution du transport ------ | 238 : | ------------------------------------------------------------ | 239 : | | 240 : |'THETA_DIFF' Coefficient de relaxation compris entre 0. et 1. | 241 : | (theta-methode diffusion) ('FLOTTANT' - defaut 1.) | 242 : | | 243 : |'THETA_CONVECTION' Idem pour la convection | 244 : | ('FLOTTANT', Defaut = THETA_DIFF) | 245 : | | 246 : |'THETA_DEC' Idem mais pour la decroissance | 247 : | ('FLOTTANT' - defaut 1/2) | 248 : | | 249 : |'THETA_DISS' Idem mais pour la dissolution | 250 : | ('FLOTTANT' - defaut 1.) | 251 : | | 252 : |'PENALISATION' Coefficient de penalisation pour la prise en | 253 : | compte de la limite de solubilite. La presence de | 254 : | cet indice ou du suivant indique quel schema a ete | 255 : | choisi. | 256 : | (Type 'FLOTTANT') - Valeur conseillee 1.D7 | 257 : | | 258 : |'EPSI_LIM' Precision relative d'arret pour le shema limite de | 259 : | solubilite predicteur-correcteur iteratif | 260 : | (Type FLOTTANT) - Valeur conseillee 5.D-3 | 261 : | | 262 : |'ITMAX_LIM' Nombre maxi d'iterations correspondant aux modules | 263 : | de dissolution avant d'abandonner | 264 : | (Type 'ENTIER') - Defaut 50 | 265 : | | 266 : |'EPSI_RET' Precision relative d'arret pour la resolution | 267 : | iterative (Picard) de l'adsorption non lineaire | 268 : | (Type FLOTTANT) - Defaut 1.D-4 | 269 : | | 270 : |'EPSI_COR' Petit saut de concentration pour calculer le coef. | 271 : | de retard par la methode de la corde lorsque le | 272 : | retard est non-lineaire. | 273 : | (Type FLOTTANT) - Defaut 1.D-4 | 274 : | | 275 : |'ITMAX_RET' Nombre maxi d'iterations correspondant au retard | 276 : | non lineaire avant d'abandonner. | 277 : | (Type 'ENTIER') - Defaut 20 | 278 : |_________________________________________________________________| 279 : 280 : 281 : 282 : --------------------------------- 283 : Resultats (stockes dans TABLE) : 284 : --------------------------------- 285 : 286 : ___________________________________________________________________ 287 : | | 288 : | Indice Contenu | 289 : | | 290 : ------------------------------------------------------------------- 291 : | | 292 : | ------ Indices communs a l'ecoulement et au transport ------ | 293 : | ------------------------------------------------------------ | 294 : | | 295 : |'TEMPS' TABLE contenant les temps sauvegardes (FLOTTANT) | 296 : | | 297 : | ------ 1ere possibilite : Resolution de l'ecoulement ------ | 298 : | ------------------------------------------------------------ | 299 : | | 300 : |'CHARGE' TABLE contenant les charges | 301 : | (Type CHPO Centre, Comp 'H') | 302 : | | 303 : |'TRACE_CHARGE' TABLE contenant les traces de charge | 304 : | (Type CHPO Face, Comp 'TH') | 305 : | | 306 : |'FLUX' TABLE contenant les debits hydrauliques integres | 307 : | par face (Type CHPO Face, comp. 'FLUX') | 308 : | | 309 : | ------ 2eme possibilite : Resolution du transport ------ | 310 : | ------------------------------------------------------------ | 311 : | | 312 : |'CONCENTRATION' TABLE contenant les concentrations | 313 : | (Type CHPO Centre, Comp 'H') | 314 : | | 315 : |'TRACE_CONC' TABLE contenant les traces de concentration | 316 : | (Type CHPO Face, Comp 'TH') | 317 : | | 318 : |'FLUX' TABLE contenant les debits diffusifs integres | 319 : | par face (Type CHPO Face, comp. 'FLUX') | 320 : | | 321 : |'PRECIPITE' TABLE contenant la quantite de precipite par maille| 322 : | (Type CHPO Centre, Comp 'H') | 323 : | | 324 : |'DISSOLUTION' TABLE contenant la quantite de precipite dissoute | 325 : | entre deux pas de temps par unite de volume et par | 326 : | unite de temps. La valeur stockee a l'indice i, | 327 : | est valable entre les temps i-1 et i | 328 : | (Type CHPO, Comp 'H'). | 329 : | ATTENTION, les valeurs de cette table resultat | 330 : | n'ont aucun sens lorsque les temps sauvegardes ne | 331 : | sont pas les memes que les temps calcules. Toute | 332 : | tentative d'exploitation donnera alors des | 333 : | resultats incoherents (erreurs de bilan) | 334 : | | 335 : |'RETARD' Si cet indice a ete prealablement defini comme une | 336 : | TABLE, alors il contient les valeurs du coefficient| 337 : | de retard (Type 'CHPO' centre, Comp 'SCAL'). Sinon,| 338 : | les valeurs du coefficient de retard ne sont pas | 339 : | sauvegardees. | 340 : |_________________________________________________________________| 341 : 342 : 343 : ___________________________________________________________________ 344 : | | 345 : | Indice Contenu | 346 : | | 347 : ------------------------------------------------------------------- 348 : | | 349 : |'TEMPS' TABLE contenant les temps sauvegardes (FLOTTANT) | 350 : | | 351 : |'TRACE_CONC' TABLE contenant les traces de concentration | 352 : | (Type 'CHPO' Centre, Comp 'TH') | 353 : | | 354 : |'CONCENTRATION' TABLE contenant les concentrations | 355 : | (Type 'CHPO' Centre, Comp 'H') | 356 : | | 357 : |'FLUX' TABLE contenant les debits diffusifs integres par | 358 : | face (Type 'CHPO' Face, Comp 'FLUX') | 359 : | | 360 : |'PRECIPITE' TABLE contenant la quantite de precipite par maille| 361 : | (Type 'CHPO' Centre, Comp 'H') | 362 : | | 363 : |'DISSOLUTION' TABLE contenant la quantite de precipite dissoute | 364 : | entre deux pas de temps par unite de volume et par | 365 : | unite de temps. La valeur stockee a l'indice i, | 366 : | est valable entre les temps i-1 et i | 367 : | (Type 'CHPO' Centre, Comp 'H') | 368 : | ATTENTION, les valeurs de cette table resultat | 369 : | n'ont aucun sens lorsque les temps sauvegardes ne | 370 : | sont pas les memes que les temps calcules. Toute | 371 : | tentative d'exploitation donnera alors des | 372 : | resultats incoherents (erreurs de bilan) | 373 : | | 374 : |'RETARD' Si cet indice a ete prealablement defini comme une | 375 : | TABLE, alors il contient les valeurs du coefficient| 376 : | de retard (Type 'CHPO' centre, Comp 'SCAL'). Sinon,| 377 : | les valeurs du coefficient de retard ne sont pas | 378 : | sauvegardees. | 379 : |_________________________________________________________________| 380 : 381 : ___________________________________________________________________ 382 : | | 383 : | Les tables resultats sont indicees par des entiers variant de 0 | 384 : | a N . | 385 : | A l'indice 0 on stocke les valeurs initiales, aux indices | 386 : | suivants les champs correspondant au temps de sortie TEMPS.I . | 387 : | Les champs servant en cas de reprise sont ceux correpondant au | 388 : | dernier indice. | 389 : |_________________________________________________________________| 390 : 391 : Exemple : pour lister le CHPOINT de charge calcule pour la 392 : valeur du parametre d'evolution 2.5, il faudra coder : 393 : LIST ( PECHE TAB1 CHARGE 2.5 ) ; 394 : ou 395 : 396 : si on sait que l'indice i de la table TEMPS contient 397 : la valeur du parametre d'evolution 2.5, on peut coder 398 : LIST ( TAB1 . CHARGE . i ) ; 399 : 400 : 401 : Remarques : 402 : ___________ 403 : 404 : Les resultats etant reperes dans TAB1 il n'y a pas d'objets nommes 405 : crees par cette procedure. 406 : 407 : Une fois sortie de DARCYTRA on peut y re-entrer en definissant de 408 : nouvelles valeurs pour les LISTREELS TEMPS_CALCULES et TEMPS_SAUVES 409 : et en invoquant de nouveau DARCYTRA avec les memes operandes que 410 : lors du premier appel. 411 : 412 : Les FLUX imposes et les flux 413 : calcules sont les flux diffusifs. Les flux convectifs s'obtiennent en 414 : post-traitement en faisant le produit du flux de vitesse (stocke a 415 : l'indice CONVECTION) par la trace de charge (stockee a l'indice 416 : TRACE_CHARGE). 417 :
© Cast3M 2003 - Tous droits réservés.
Mentions légales