Accueil | Cast3M

   1 : $$$$ TRAJ     NOTICE  CHAT      11/09/12    21:18:23     7124           
   2 :                                              DATE     11/09/12
   3 :   
   4 :    Operateur TRAJ                            Voir aussi :
   5 :    --------------  
   6 : 
   7 : 
   8 :     Cas d'une formulation elements finis :
   9 :     MODL4 MCH4 = TRAJ MOT1  |CHPO1|   |TAB1 |  ('PORO' MCH1) TAB2 ;
  10 :                             |TAB4 |   |MODL1|
  11 : 
  12 :     Cas d'une formulation mixte hybride(modele DARCY)
  13 :     MODL4 MCH4 = TRAJ MOT1 MODL1 |CHPO2|  ('PORO' MCH1)  TAB2  ;
  14 :                                  |TAB5 |  ('DISP' MCH2)
  15 :                                           ('DIFF' MCH3)
  16 : 
  17 :       Objet
  18 :      -------
  19 :      L'operateur TRAJ permet de calculer les trajectoires de particules
  20 :      lachees dans un domaine maille pour lequel on connait :
  21 :      soit un champs de vitesses ou de flux, constant au cours du temps,
  22 :      soit des champs de vitesses ou de flux, donnes pour differentes
  23 :      valeurs du temps.
  24 :      
  25 :      Commentaire
  26 :      -----------
  27 :  
  28 :      TAB1   est un objet de type TABLE et de sous type DOMAINE.
  29 :             C'est le resultat de l'operateur DOMA applique au maillage
  30 :             sur lequel on fait le calcul.
  31 :             
  32 :      MOT1   indique le type de calcul que l'on veut faire :
  33 :             'CONVECTION_EXPLICITE' la position des particules est
  34 :                  calculee de proche en proche en fonction de la
  35 :                  vitesse locale (c'est l'option par defaut).
  36 :             'CONVECTION_ANALYTIQUE' Calcul des lignes de courant par 
  37 :                  integration analytique (uniquement en formulation EFMH)
  38 :             'CONVECTION_DIFFUSION' Calcul des trajectoires par
  39 :                  iterations successives en prenant en compte les
  40 :                  phenomenes de convection-dipersion-diffusion. (Cette
  41 :                  option n'est pour l'instant developpee que pour le
  42 :                  modele DARCY)
  43 :                  
  44 :      MODL1  Objet modele (type MMODEL) decrivant la formulation
  45 :             utilisee (cf. MODE). Les formulations actuellement
  46 :             prevues sont DARCY et NAVIER_STOKES.
  47 : 
  48 :      CHPO1  champ de vitesse defini aux noeuds du maillage TAB1.MAILLAGE
  49 :             Les composantes ont pour noms VX VY (VZ).
  50 : 
  51 :      TAB4   est le nom d'un objet de type TABLE et de sous type
  52 :             TRANSITOIRE. Cette table contient a l'indice 'TEMPS'
  53 :             une table de flottants, et a l'indice 'VITESSE',une
  54 :             table de CHPOINT. Ces deux tables sont indicees par des
  55 :             entiers 0 1 2 ...N :
  56 :               TAB4.VITESSE.N est un CHPOINT ayant les memes
  57 :                   caracteristiques que CHPO1.
  58 :               TAB4.TEMPS.N est le temps correspondant.
  59 : 
  60 :      CHP02  Objet de type CHPOINT contenant le debit a travers
  61 :             chaque face. Le support geometrique de ce champ est
  62 :             le maillage contitue les centres des FACES. Le nom de
  63 :             la composante du CHPOINT est FLUX (cf. HDEB).
  64 : 
  65 :      TAB5   est un objet de type TABLE et de sous-type
  66 :             DARCY_TRANSITOIRE. Elle contient a l'indice 'TEMPS'
  67 :             une table de flottants, et a l'indice 'FLUX', une table de
  68 :             CHPOINT. Ces deux tables sont indicees par des
  69 :             entiers 0 1 2 ...N (cf DARCYTRA) :
  70 :               TAB5.FLUX.N est un CHPOINT ayant les memes
  71 :                 caracteristiques que CHPO2.
  72 :               TAB5.TEMPS.N est le temps correspondant.
  73 :                
  74 :      MCH1   objet MCHAML contenant la porosite au centre de gravite
  75 :             de chaque element. Au moment du calcul la vitesse sera
  76 :             divisee par la porosite. Si cette valeur  est absente la
  77 :             porosite est supposee egale a 1.
  78 : 
  79 :      MCH2   objet MCHAML a deux composantes contenant respectivement
  80 :             la dispersivite longitudinale et la dispersivite
  81 :             transversale au centre de gravite de chaque element. Les
  82 :             dispersivites sont imposees constantes par element.
  83 :             Par defaut, la dispersivite est nulle. 
  84 :             (utilise uniquement dans le calcul 'CONVECTION_DIFFUSION')
  85 : 
  86 :      MCH3   objet MCHAML a une composante contenant la diffusion
  87 :             isotrope effective au centre de gravite de chaque
  88 :             element. Par defaut, la diffusion est nulle.
  89 :             (utilise uniquement dans le calcul 'CONVECTION_DIFFUSION')
  90 : 
  91 :      TAB2   table a plusieurs indices contenant la description
  92 :             du lacher de particules :
  93 : 
  94 :      TAB2.'TEMPS_LIMITE' contient un reel : le temps maximal de calcul
  95 : 
  96 :      TAB2.'CFL' contient un reel : le nombre de Courant a respecter.
  97 :             Le pas de temps de calcul en depend. En moyenne, il y aura
  98 :             1/CFL sauts de particule par maille.
  99 :             Ce nombre doit etre compris entre 1.E-8 et 1.(Defaut 0.05)
 100 :             (utilise uniquement dans le calcul 'CONVECTION_EXPLICITE').
 101 : 
 102 :      TAB2.'DELTAT_SAUVE' contient un reel : le pas de temps avec lequel
 103 :             on conserve les resultats pour un post-traitement. Si
 104 :             cette valeur est nulle, tous les temps de calcul seront
 105 :             sauvegardes. 
 106 : 
 107 :      TAB2.'TEMPS_SAUVES' contient un 'LISTREEL' : liste des temps
 108 :             auxquels les resultats sont sauvegardes pour un
 109 :             post-traitement.
 110 :             Cet indice 'ET' le precedent, 'DELTAT_SAUVE', ne peuvent
 111 :             pas co-exister. Si aucun des deux n'existe, tous les temps
 112 :             de calcul seront sauvegardes. 
 113 :             Disponible uniquement dans le calcul 'CONVECTION_DIFFUSION'.
 114 : 
 115 : 
 116 :      TAB2.'IMPERMEABLE' : objet MAILLAGE de points faces de type 'POI1'
 117 :             determinant les faces infranchissables par les particules.
 118 :             Disponible uniquement dans un calcul 'CONVECTION_DIFFUSION'.
 119 : 
 120 :      TAB2.'TEMPS_LACHER' : LISTREEL contenant les temps initiaux des
 121 :             differents lachers. Pour chacun des reels de cette liste,
 122 :             il doit y avoir un maillage de points de lacher (voir
 123 :             plus bas).
 124 : 
 125 :      TAB2.1 : maillage (elements 'POI1') representant les positions des
 126 :             particules qui seront lachees au premier temps de la liste
 127 :             TAB2.'TEMPS_LACHER' 
 128 :        "           "
 129 :        "           "
 130 :        "           "
 131 :      TAB2.i : maillage de POI1 representant les positions des
 132 :             particules qui seront lachees au ieme temps de la liste
 133 :             TAB2.'TEMPS_LACHER' 
 134 : 
 135 :      MODL4  est un modele, concatenation d'un modele elementaire par
 136 :             trajectoire. La formulation est 'TRAJECTOIRE'.
 137 :             Le maillage associe est un maillage de SEG2 appuyes sur
 138 :             les points des trajectoires (voir 'EXTR' 'ZONE' et 'COOR'
 139 :             pour en obtenir les coordonnees).
 140 : 
 141 :      MCH4   est un MCHAML base sur le modele MOD4, a valeurs aux
 142 :             noeuds des SEG2 des trajectoires. La composante
 143 :             'TMPS' contient le temps en chacun des points de chaque
 144 :             trajectoire. 
 145 : 
 146 :             
 147 :      Exemple:   (voir valitraj.dgibi et trajec.dgibi)
 148 :      --------
 149 : 
 150 :      si CARRE est un maillage et VN un champ de vitesse aux noeuds,
 151 :      et MODCAR un modele de type NAVIER_STOKES, on pourra coder la
 152 :      sequence suivante. 
 153 : 
 154 : *    DEFINITION DE LA TABLE DE LACHER
 155 :      LACHER = TABLE ;
 156 :      LACHER.TEMPS_LACHER = PROG 0.  2. ;
 157 :      LACHER.TEMPS_LIMITE =  180.;
 158 :      LACHER.CFL =  0.05 ;
 159 :      LACHER.DELTAT_SAUVE =  0. ;
 160 : *    points du premier lacher a t=0. :      
 161 :      LACHER.1 = (0.0 0.0) 'ET' (0.0 1.) 'ET' (0.0 2.0) ;
 162 : *    points du second lacher a t= 2. :      
 163 :      LACHER.2 = 'MANU' 'POI1' (0.0 1.)  ;
 164 :     
 165 : *    CALCUL DES TRAJECTOIRES
 166 :      MODL4 MCH4 = 'TRAJ' MODCAR VN LACHER ;
 167 : *    TRACE DES TRAJECTOIRES
 168 :      CROB = 'EXTR' MODL4 'MAIL' ;
 169 :      'TRACER' CROB  ;
 170 : *    TRACE DES ISOVALEURS DU TEMPS
 171 :      'TRACER' MCH4 MODL4 ;
 172 :       
 173 : *    IMPRESSION DU TEMPS DE SORTIE DE LA 1ERE PARTICULE
 174 :      TABZON = 'EXTR' MODL4 'ZONE' ;
 175 :      NBB    = 'NBEL' TABZON . 2 ;
 176 :      TT1    = 'EXTR' MCH4 'TMPS' 1 NBB 2 ;
 177 :      'LISTE' TT1 ;
 178 : 
 179 :         
 180 :      Remarques :
 181 :      -----------
 182 :      
 183 :     1. Les elements de maillages actuellement utilisables sont,
 184 :        en 2D : TRI3  QUA4
 185 :        en 3D : PRI6  CUB8  TET4
 186 :     2. Le nombre de Courant est utilise pour calculer le pas de
 187 :        temps dans chaque element du maillage par la formule
 188 :        dt = (L * C) / Vm
 189 :        ou L est 2 fois la plus petite distance entre le point centre
 190 :               et le milieu de faces de l'element.
 191 :           Vm est la vitesse moyenne sur l'element.
 192 :           dt est le pas de temps du calcul.
 193 :           C  est le nombre de Courant.
 194 :        C doit etre suffisament petit pour permettre de faire plusieurs
 195 :        pas par element. 
 196 :        Nous avons fait nos tests de validation avec C=0.05 .
 197 :      3. Lorsque TAB2.DELTAT_SAUVE est donne different de 0.,
 198 :         on conserve les positions des particules avec le pas de temps
 199 :         donne. Toutefois si la particule sort du domaine entre deux pas
 200 :         de temps on conserve la position et le temps du point de sortie.
 201 :      4. Dans le cas d'un calcul transitoire, la vitesse ou le flux sont
 202 :         interpoles lineairement en fonction du temps.
 203 : 
 204 : 
 205 : 
 206 :  
 207 : 
 208 :  
 209 : 
 210 :  
 211 :  
 212 :  
 213 :  
 214 :