Afficher cette notice en
$$$$ DARCYTRA NOTICE CHAT 11/09/12 21:15:43 7124 DATE 11/09/12 Procedure DARCYTRA Voir aussi : PECHE ------------------ DARCYSAT DARCYTRA TAB1 ; TAB1.'SOUSTYPE'.'MODELE'.'DOMAINE'. 'CARACTERISTIQUES'.'EMMAGASINEMENT'.'CONVECTION'. 'TEMPS'.'TRACE_CHARGE'.'CHARGE'.'FLUX'. 'BLOCAGE'.'TRACE_IMPOSE'.'FLUX_IMPOSE'.'SOURCE'. 'TEMPS_CALCULES'.'TEMPS_SAUVES'. 'THETA'.'THETA_CONVECTION' ou TAB1.'SOUSTYPE'.'MODELE'.'DOMAINE'.'ORIENTATION'. 'CARACTERISTIQUES'.'POROSITE'.'DECROISSANCE'. 'COEF_RETARD'.'LANGMUIR'.'FREUNDLICH'. 'LIMITE_SOLUBILITE'.'COEF_DISSOLUTION'. 'CONVECTION'.'TEMPS'.'TRACE_CONC'.'CONCENTRATION'. 'FLUX'.'PRECIPITE'.'DISSOLUTION'. 'BLOCAGE'.'TRACE_IMPOSE'.'FLUX_IMPOSE'. 'DISSOLUTION_IMPOSEE'.'SOURCE'.'TEMPS_CALCULES'. 'TEMPS_SAUVES'.'THETA_DIFF'.'THETA_CONVECTION'. 'THETA_DEC'.'THETA_DISS'.'PENALISATION'. 'EPSI_LIM'.'ITMAX_LIM'.'EPSI_RET'.'EPSI_COR'. 'ITMAX_RET' Objet : _______ Cette procedure a deux fonctions. 1) En presence de l'indice 'CHARGE', on resoud les equations de DARCY en transitoire pour l'ecoulement par une methode d'elements finis mixtes hybrides (EFMH). Les inconnues du probleme sont la charge ('H'), la trace de charge ('TH') et le flux diffusif ('FLUX'). 2) En l'absence de l'indice 'CHARGE', resoud l'equation de transport par diffusion-convection d'un champ scalaire actif par un fluide dont la vitesse est connue. L'espece peut se trouver sous trois formes : solute, adsorbat et precipite, dont les lois d'echange doivent etre specifiees. A chaque loi correspondent un ou plusieurs algorithmes auxquels des parametres numeriques doivent etre fournis. On utilise la modelisation Darcy EFMH. Les inconnues du probleme sont la concentration ('H'), la trace de concentration ('TH') et le flux diffusif ('FLUX'). Commentaire : _____________ En entree, TAB1 sert a definir les options et les parametres du calcul. Les indices de la table TAB1 sont des mots (a coder tel quel) dont voici la description : ___________________________________________________________________ | | | Indice Contenu | | | ------------------------------------------------------------------- | | |------------------------------------------------ | |Donnees physiques, geometriques et materielles : | |------------------------------------------------ | | | | ------ Indices communs a l'ecoulement et au transport ------ | | ------------------------------------------------------------ | | | |'SOUSTYPE' 'DARCY' (type MOT) | | | |'MODELE' Objet modele (MMODEL cree par MODE) | | | |'DOMAINE' References geometriques (TABLE creee par DOMA) | | | | ------ 1ere possibilite : Resolution de l'ecoulement ------ | | ----------------------------------------------------------- | | | |'CARACTERISTIQUES' Donnees physiques et materielles : | | conductivite hydraulique (CHAMELEM cree par MATE) | | | |'EMMAGASINEMENT' Valeur du coefficient d'emmagasinement | | (Type CHPO Centre, Comp 'CK', ou FLOTTANT) | | - Defaut 1. | | | | ------ 2eme possibilite : Resolution du transport ------ | | ------------------------------------------------------------ | | | |'CARACTERISTIQUES' Donnees physiques et materielles : | | diffusivite effective (CHAMELEM cree par MATE) | | | |'POROSITE' Valeur de la porosite (Type CHPO Centre, Comp | | 'CK', ou FLOTTANT) - Defaut 1. | | | |'DECROISSANCE' Valeur du terme de decroissance (Type FLOTTANT) | | Tel que dC/dt = - Lambda * C - Defaut 0. | | | |'COEF_RETARD' Coefficient de retard lineaire dans le cas simple, | | ou Pente a l'origine de la fonction F(C) dans le | | cas d'isotherme non lineaire de Langmuir | | ou Coefficient K de l'isotherme de Freundlich | | (Type CHPO Centre 'SCAL', ou FLOTTANT) | | | |'LANGMUIR' Quantite maximale adsorbee sur le solide | | rapportee a l'unite de volume du fluide et exprimee| | dans la meme unite que le solute. | | (Type CHPO Centre 'SCAL', ou FLOTTANT). | | F = (R-1) C / [1 + ((R-1) C / Fsat)] | | Si cet indice et le suivant sont absents, | | l'equilibre d'adsorption est lineaire. Cet indice a| | priorite sur l'indice FREUNDLICH. | | | |'FREUNDLICH' Exposant de la loi de Freundlich F = K (C ^ 1/n) | | (Type FLOTTANT). | | Dans ce cas (et si l'indice LANGMUIR n'existe pas),| | l'indice 'RETARD' contient le coefficient | | K ramene a une unite de volume de fluide. | | - Non disponible pour l'instant - | | | |'LIMITE_SOLUBILITE' Limite de solubilite (Type MCHAML, Comp 'H') | | | |'COEF_DISSOLUTION' Coef. de dissolution (Type CHPO Centre, Comp | | 'SCAL'). Tel que dC/dt = Coef * (Csat - C) - Par | | defaut, la dissolution est instantanee | | | |'CONVECTION' Debit integre de la vitesse convective a travers | | chaque face des elements (Type CHPO Face, comp. | | 'FLUX') | | | |---------------------- | |Conditions initiales : | |---------------------- | | | | ------ Indices communs a l'ecoulement et au transport ------ | | ------------------------------------------------------------ | | | |'TEMPS' TABLE contenant a l'indice 0 la valeur du temps | | initial (FLOTTANT) | | | | | | ------ 1ere possibilite : Resolution de l'ecoulement ------ | | ----------------------------------------------------------- | | | |'CHARGE' TABLE contenant a l'indice 0 la charge hydraulique | | (quantite d'element par unite de volume d'eau) | | (Type CHPO Centre, Comp 'H') | | | |'TRACE_CHARGE' TABLE contenant a l'indice 0 la trace de | | charge initiale (CHPO, 'TH') | | | |'FLUX' TABLE contenant a l'indice 0 le flux hydraulique | | initial integre sur chaque face (Type CHPO Face, | | comp. 'FLUX') | | | | ------ 2eme possibilite : Resolution du transport ------ | | ------------------------------------------------------------ | | | |'CONCENTRATION' TABLE contenant a l'indice 0 la concentration | | (quantite d'element par unite de volume d'eau) | | (Type CHPO Centre, Comp 'H') | | | |'TRACE_CONC' TABLE contenant a l'indice 0 la trace de | | concentration initiale (Type CHPO Face, Comp 'TH') | | | |'FLUX' TABLE contenant a l'indice 0 le flux diffusif | | initial integre sur chaque face (Type CHPO Face, | | comp. 'FLUX') | | | |'PRECIPITE' TABLE contenant a l'indice 0 la quantite initiale | | de precipite par unite de volume de milieu | | (Type CHPO Centre, Comp 'H') | | | |'DISSOLUTION' TABLE contenant a l'indice 0 la quantite initiale | | pour estimer la dissolution au premier pas de temps| | (Type CHPO, Comp 'H'), voir plus loin. | | | |-------------------------------------- | |Conditions aux limites / chargements : | |-------------------------------------- | | | | ------ Indices communs a l'ecoulement et au transport ------ | | ------------------------------------------------------------ | | | |'BLOCAGE' Contient les matrices de blocage (RIGIDITE) | | | |'TRACE_IMPOSE' Valeurs des traces imposees (charge ou concentra- | | -tion) (CHARGEMENT 'TH' - Obligatoire si BLOCAGE) | | | |'FLUX_IMPOSE' Valeurs des flux imposes integres par face | | (Type CHARGEMENT de CHPO Face, comp. 'FLUX'- | | defaut 0.) | | | |'SOURCE' Valeurs du terme source par maille et par unite de | | temps (ex : puits, filiation) | | Les valeurs a l'indice i sont les valeurs entre | | les temps i-1 et i. | | (CHARGEMENT de CHPO Centre, comp 'SOUR'- defaut 0.)| | | | ------ 2eme possibilite : Resolution du transport ------ | | ------------------------------------------------------------ | | | |'DISSOLUTION_IMPOSEE' Valeurs des dissolutions imposees par unite| | de temps et par maille. (Type CHARGEMENT de CHPO, | | Comp 'H'). Les valeurs a l'indice i sont les | | valeurs moyennes de dissolution par unite de temps | | entre les temps i-1 et i. | | Priorite de la dissolution imposee sur les | | cinetiques. | | | |-------------------- | |Donnees numeriques : | |-------------------- | | | | ------ Indices communs a l'ecoulement et au transport ------ | | ------------------------------------------------------------ | | | |'TEMPS_CALCULES' Valeur des temps calcules (LISTREEL) | | Contient obligatoirement le temps final. | | | |'TEMPS_SAUVES' Valeur des temps sauvegardes (LISTREEL - defaut : | | on sauve tous les pas de temps) | | | | ------ 1ere possibilite : Resolution de l'ecoulement ------ | | ----------------------------------------------------------- | | | |'THETA' Coefficient de relaxation compris entre 0. et 1. | | (theta-methode diffusion) (FLOTTANT - defaut 1.) | | Possibilite de non-convergence lorsque theta<1/2 | | Valeurs de theta generalement utilisees : | | Schema de Euler explicite : 0. | | Schema de Crank-Nicholson : 1/2 | | Schema de Galerkin : 2/3 | | Schema de Euler implicite : 1. | | | | ------ 2eme possibilite : Resolution du transport ------ | | ------------------------------------------------------------ | | | |'THETA_DIFF' Coefficient de relaxation compris entre 0. et 1. | | (theta-methode diffusion) ('FLOTTANT' - defaut 1.) | | | |'THETA_CONVECTION' Idem pour la convection | | ('FLOTTANT', Defaut = THETA_DIFF) | | | |'THETA_DEC' Idem mais pour la decroissance | | ('FLOTTANT' - defaut 1/2) | | | |'THETA_DISS' Idem mais pour la dissolution | | ('FLOTTANT' - defaut 1.) | | | |'PENALISATION' Coefficient de penalisation pour la prise en | | compte de la limite de solubilite. La presence de | | cet indice ou du suivant indique quel schema a ete | | choisi. | | (Type 'FLOTTANT') - Valeur conseillee 1.D7 | | | |'EPSI_LIM' Precision relative d'arret pour le shema limite de | | solubilite predicteur-correcteur iteratif | | (Type FLOTTANT) - Valeur conseillee 5.D-3 | | | |'ITMAX_LIM' Nombre maxi d'iterations correspondant aux modules | | de dissolution avant d'abandonner | | (Type 'ENTIER') - Defaut 50 | | | |'EPSI_RET' Precision relative d'arret pour la resolution | | iterative (Picard) de l'adsorption non lineaire | | (Type FLOTTANT) - Defaut 1.D-4 | | | |'EPSI_COR' Petit saut de concentration pour calculer le coef. | | de retard par la methode de la corde lorsque le | | retard est non-lineaire. | | (Type FLOTTANT) - Defaut 1.D-4 | | | |'ITMAX_RET' Nombre maxi d'iterations correspondant au retard | | non lineaire avant d'abandonner. | | (Type 'ENTIER') - Defaut 20 | |_________________________________________________________________| --------------------------------- Resultats (stockes dans TABLE) : --------------------------------- ___________________________________________________________________ | | | Indice Contenu | | | ------------------------------------------------------------------- | | | ------ Indices communs a l'ecoulement et au transport ------ | | ------------------------------------------------------------ | | | |'TEMPS' TABLE contenant les temps sauvegardes (FLOTTANT) | | | | ------ 1ere possibilite : Resolution de l'ecoulement ------ | | ------------------------------------------------------------ | | | |'CHARGE' TABLE contenant les charges | | (Type CHPO Centre, Comp 'H') | | | |'TRACE_CHARGE' TABLE contenant les traces de charge | | (Type CHPO Face, Comp 'TH') | | | |'FLUX' TABLE contenant les debits hydrauliques integres | | par face (Type CHPO Face, comp. 'FLUX') | | | | ------ 2eme possibilite : Resolution du transport ------ | | ------------------------------------------------------------ | | | |'CONCENTRATION' TABLE contenant les concentrations | | (Type CHPO Centre, Comp 'H') | | | |'TRACE_CONC' TABLE contenant les traces de concentration | | (Type CHPO Face, Comp 'TH') | | | |'FLUX' TABLE contenant les debits diffusifs integres | | par face (Type CHPO Face, comp. 'FLUX') | | | |'PRECIPITE' TABLE contenant la quantite de precipite par maille| | (Type CHPO Centre, Comp 'H') | | | |'DISSOLUTION' TABLE contenant la quantite de precipite dissoute | | entre deux pas de temps par unite de volume et par | | unite de temps. La valeur stockee a l'indice i, | | est valable entre les temps i-1 et i | | (Type CHPO, Comp 'H'). | | ATTENTION, les valeurs de cette table resultat | | n'ont aucun sens lorsque les temps sauvegardes ne | | sont pas les memes que les temps calcules. Toute | | tentative d'exploitation donnera alors des | | resultats incoherents (erreurs de bilan) | | | |'RETARD' Si cet indice a ete prealablement defini comme une | | TABLE, alors il contient les valeurs du coefficient| | de retard (Type 'CHPO' centre, Comp 'SCAL'). Sinon,| | les valeurs du coefficient de retard ne sont pas | | sauvegardees. | |_________________________________________________________________| ___________________________________________________________________ | | | Indice Contenu | | | ------------------------------------------------------------------- | | |'TEMPS' TABLE contenant les temps sauvegardes (FLOTTANT) | | | |'TRACE_CONC' TABLE contenant les traces de concentration | | (Type 'CHPO' Centre, Comp 'TH') | | | |'CONCENTRATION' TABLE contenant les concentrations | | (Type 'CHPO' Centre, Comp 'H') | | | |'FLUX' TABLE contenant les debits diffusifs integres par | | face (Type 'CHPO' Face, Comp 'FLUX') | | | |'PRECIPITE' TABLE contenant la quantite de precipite par maille| | (Type 'CHPO' Centre, Comp 'H') | | | |'DISSOLUTION' TABLE contenant la quantite de precipite dissoute | | entre deux pas de temps par unite de volume et par | | unite de temps. La valeur stockee a l'indice i, | | est valable entre les temps i-1 et i | | (Type 'CHPO' Centre, Comp 'H') | | ATTENTION, les valeurs de cette table resultat | | n'ont aucun sens lorsque les temps sauvegardes ne | | sont pas les memes que les temps calcules. Toute | | tentative d'exploitation donnera alors des | | resultats incoherents (erreurs de bilan) | | | |'RETARD' Si cet indice a ete prealablement defini comme une | | TABLE, alors il contient les valeurs du coefficient| | de retard (Type 'CHPO' centre, Comp 'SCAL'). Sinon,| | les valeurs du coefficient de retard ne sont pas | | sauvegardees. | |_________________________________________________________________| ___________________________________________________________________ | | | Les tables resultats sont indicees par des entiers variant de 0 | | a N . | | A l'indice 0 on stocke les valeurs initiales, aux indices | | suivants les champs correspondant au temps de sortie TEMPS.I . | | Les champs servant en cas de reprise sont ceux correpondant au | | dernier indice. | |_________________________________________________________________| Exemple : pour lister le CHPOINT de charge calcule pour la valeur du parametre d'evolution 2.5, il faudra coder : LIST ( PECHE TAB1 CHARGE 2.5 ) ; ou si on sait que l'indice i de la table TEMPS contient la valeur du parametre d'evolution 2.5, on peut coder LIST ( TAB1 . CHARGE . i ) ; Remarques : ___________ Les resultats etant reperes dans TAB1 il n'y a pas d'objets nommes crees par cette procedure. Une fois sortie de DARCYTRA on peut y re-entrer en definissant de nouvelles valeurs pour les LISTREELS TEMPS_CALCULES et TEMPS_SAUVES et en invoquant de nouveau DARCYTRA avec les memes operandes que lors du premier appel. Les FLUX imposes et les flux calcules sont les flux diffusifs. Les flux convectifs s'obtiennent en post-traitement en faisant le produit du flux de vitesse (stocke a l'indice CONVECTION) par la trace de charge (stockee a l'indice TRACE_CHARGE).
© Cast3M 2003 - Tous droits réservés.
Mentions légales