Accueil | Cast3M

$$$$ TRAJ     NOTICE  CHAT      11/09/12    21:18:23     7124           
                                             DATE     11/09/12
  
   Operateur TRAJ                            Voir aussi :
   --------------  


    Cas d'une formulation elements finis :
    MODL4 MCH4 = TRAJ MOT1  |CHPO1|   |TAB1 |  ('PORO' MCH1) TAB2 ;
                            |TAB4 |   |MODL1|

    Cas d'une formulation mixte hybride(modele DARCY)
    MODL4 MCH4 = TRAJ MOT1 MODL1 |CHPO2|  ('PORO' MCH1)  TAB2  ;
                                 |TAB5 |  ('DISP' MCH2)
                                          ('DIFF' MCH3)

      Objet
     -------
     L'operateur TRAJ permet de calculer les trajectoires de particules
     lachees dans un domaine maille pour lequel on connait :
     soit un champs de vitesses ou de flux, constant au cours du temps,
     soit des champs de vitesses ou de flux, donnes pour differentes
     valeurs du temps.
     
     Commentaire
     -----------
 
     TAB1   est un objet de type TABLE et de sous type DOMAINE.
            C'est le resultat de l'operateur DOMA applique au maillage
            sur lequel on fait le calcul.
            
     MOT1   indique le type de calcul que l'on veut faire :
            'CONVECTION_EXPLICITE' la position des particules est
                 calculee de proche en proche en fonction de la
                 vitesse locale (c'est l'option par defaut).
            'CONVECTION_ANALYTIQUE' Calcul des lignes de courant par 
                 integration analytique (uniquement en formulation EFMH)
            'CONVECTION_DIFFUSION' Calcul des trajectoires par
                 iterations successives en prenant en compte les
                 phenomenes de convection-dipersion-diffusion. (Cette
                 option n'est pour l'instant developpee que pour le
                 modele DARCY)
                 
     MODL1  Objet modele (type MMODEL) decrivant la formulation
            utilisee (cf. MODE). Les formulations actuellement
            prevues sont DARCY et NAVIER_STOKES.

     CHPO1  champ de vitesse defini aux noeuds du maillage TAB1.MAILLAGE
            Les composantes ont pour noms VX VY (VZ).

     TAB4   est le nom d'un objet de type TABLE et de sous type
            TRANSITOIRE. Cette table contient a l'indice 'TEMPS'
            une table de flottants, et a l'indice 'VITESSE',une
            table de CHPOINT. Ces deux tables sont indicees par des
            entiers 0 1 2 ...N :
              TAB4.VITESSE.N est un CHPOINT ayant les memes
                  caracteristiques que CHPO1.
              TAB4.TEMPS.N est le temps correspondant.

     CHP02  Objet de type CHPOINT contenant le debit a travers
            chaque face. Le support geometrique de ce champ est
            le maillage contitue les centres des FACES. Le nom de
            la composante du CHPOINT est FLUX (cf. HDEB).

     TAB5   est un objet de type TABLE et de sous-type
            DARCY_TRANSITOIRE. Elle contient a l'indice 'TEMPS'
            une table de flottants, et a l'indice 'FLUX', une table de
            CHPOINT. Ces deux tables sont indicees par des
            entiers 0 1 2 ...N (cf DARCYTRA) :
              TAB5.FLUX.N est un CHPOINT ayant les memes
                caracteristiques que CHPO2.
              TAB5.TEMPS.N est le temps correspondant.
               
     MCH1   objet MCHAML contenant la porosite au centre de gravite
            de chaque element. Au moment du calcul la vitesse sera
            divisee par la porosite. Si cette valeur  est absente la
            porosite est supposee egale a 1.

     MCH2   objet MCHAML a deux composantes contenant respectivement
            la dispersivite longitudinale et la dispersivite
            transversale au centre de gravite de chaque element. Les
            dispersivites sont imposees constantes par element.
            Par defaut, la dispersivite est nulle. 
            (utilise uniquement dans le calcul 'CONVECTION_DIFFUSION')

     MCH3   objet MCHAML a une composante contenant la diffusion
            isotrope effective au centre de gravite de chaque
            element. Par defaut, la diffusion est nulle.
            (utilise uniquement dans le calcul 'CONVECTION_DIFFUSION')

     TAB2   table a plusieurs indices contenant la description
            du lacher de particules :

     TAB2.'TEMPS_LIMITE' contient un reel : le temps maximal de calcul

     TAB2.'CFL' contient un reel : le nombre de Courant a respecter.
            Le pas de temps de calcul en depend. En moyenne, il y aura
            1/CFL sauts de particule par maille.
            Ce nombre doit etre compris entre 1.E-8 et 1.(Defaut 0.05)
            (utilise uniquement dans le calcul 'CONVECTION_EXPLICITE').

     TAB2.'DELTAT_SAUVE' contient un reel : le pas de temps avec lequel
            on conserve les resultats pour un post-traitement. Si
            cette valeur est nulle, tous les temps de calcul seront
            sauvegardes. 

     TAB2.'TEMPS_SAUVES' contient un 'LISTREEL' : liste des temps
            auxquels les resultats sont sauvegardes pour un
            post-traitement.
            Cet indice 'ET' le precedent, 'DELTAT_SAUVE', ne peuvent
            pas co-exister. Si aucun des deux n'existe, tous les temps
            de calcul seront sauvegardes. 
            Disponible uniquement dans le calcul 'CONVECTION_DIFFUSION'.


     TAB2.'IMPERMEABLE' : objet MAILLAGE de points faces de type 'POI1'
            determinant les faces infranchissables par les particules.
            Disponible uniquement dans un calcul 'CONVECTION_DIFFUSION'.

     TAB2.'TEMPS_LACHER' : LISTREEL contenant les temps initiaux des
            differents lachers. Pour chacun des reels de cette liste,
            il doit y avoir un maillage de points de lacher (voir
            plus bas).

     TAB2.1 : maillage (elements 'POI1') representant les positions des
            particules qui seront lachees au premier temps de la liste
            TAB2.'TEMPS_LACHER' 
       "           "
       "           "
       "           "
     TAB2.i : maillage de POI1 representant les positions des
            particules qui seront lachees au ieme temps de la liste
            TAB2.'TEMPS_LACHER' 

     MODL4  est un modele, concatenation d'un modele elementaire par
            trajectoire. La formulation est 'TRAJECTOIRE'.
            Le maillage associe est un maillage de SEG2 appuyes sur
            les points des trajectoires (voir 'EXTR' 'ZONE' et 'COOR'
            pour en obtenir les coordonnees).

     MCH4   est un MCHAML base sur le modele MOD4, a valeurs aux
            noeuds des SEG2 des trajectoires. La composante
            'TMPS' contient le temps en chacun des points de chaque
            trajectoire. 

            
     Exemple:   (voir valitraj.dgibi et trajec.dgibi)
     --------

     si CARRE est un maillage et VN un champ de vitesse aux noeuds,
     et MODCAR un modele de type NAVIER_STOKES, on pourra coder la
     sequence suivante. 

*    DEFINITION DE LA TABLE DE LACHER
     LACHER = TABLE ;
     LACHER.TEMPS_LACHER = PROG 0.  2. ;
     LACHER.TEMPS_LIMITE =  180.;
     LACHER.CFL =  0.05 ;
     LACHER.DELTAT_SAUVE =  0. ;
*    points du premier lacher a t=0. :      
     LACHER.1 = (0.0 0.0) 'ET' (0.0 1.) 'ET' (0.0 2.0) ;
*    points du second lacher a t= 2. :      
     LACHER.2 = 'MANU' 'POI1' (0.0 1.)  ;
    
*    CALCUL DES TRAJECTOIRES
     MODL4 MCH4 = 'TRAJ' MODCAR VN LACHER ;
*    TRACE DES TRAJECTOIRES
     CROB = 'EXTR' MODL4 'MAIL' ;
     'TRACER' CROB  ;
*    TRACE DES ISOVALEURS DU TEMPS
     'TRACER' MCH4 MODL4 ;
      
*    IMPRESSION DU TEMPS DE SORTIE DE LA 1ERE PARTICULE
     TABZON = 'EXTR' MODL4 'ZONE' ;
     NBB    = 'NBEL' TABZON . 2 ;
     TT1    = 'EXTR' MCH4 'TMPS' 1 NBB 2 ;
     'LISTE' TT1 ;

        
     Remarques :
     -----------
     
    1. Les elements de maillages actuellement utilisables sont,
       en 2D : TRI3  QUA4
       en 3D : PRI6  CUB8  TET4
    2. Le nombre de Courant est utilise pour calculer le pas de
       temps dans chaque element du maillage par la formule
       dt = (L * C) / Vm
       ou L est 2 fois la plus petite distance entre le point centre
              et le milieu de faces de l'element.
          Vm est la vitesse moyenne sur l'element.
          dt est le pas de temps du calcul.
          C  est le nombre de Courant.
       C doit etre suffisament petit pour permettre de faire plusieurs
       pas par element. 
       Nous avons fait nos tests de validation avec C=0.05 .
     3. Lorsque TAB2.DELTAT_SAUVE est donne different de 0.,
        on conserve les positions des particules avec le pas de temps
        donne. Toutefois si la particule sort du domaine entre deux pas
        de temps on conserve la position et le temps du point de sortie.
     4. Dans le cas d'un calcul transitoire, la vitesse ou le flux sont
        interpoles lineairement en fonction du temps.