Procedures | Cast3M

   1 : * CALCDISP  PROCEDUR  GBM       04/12/22    21:15:01     4973           
   2 : *---------------------------------------------------------------------
   3 : *---------PROCEDURE DE CALCUL DE LA DISPERSIVITE----------------------
   4 : *---------------------------------------------------------------------
   5 : 'DEBP'  CALCDISP QELEM*'CHPOINT' DISPL*'CHPOINT' DISPT*'CHPOINT'     ;
   6 : *  |-----------------------------------------------------------------|
   7 : *  | Généralités : CALCDISP calcule le tenseur de dispersion         |
   8 : *  |               du problème de transport convection-diffusion.    |
   9 : *  |-----------------------------------------------------------------|
  10 : *  |                                                                 |
  11 : *  |-----------------------------------------------------------------|
  12 : *  |                         ENTREES                                 |
  13 : *  |-----------------------------------------------------------------|
  14 : *  |                                                                 |
  15 : *  | ALPHAL   : coefficient longitudinal de dispersivité             |
  16 : *  |                                                                 |
  17 : *  | ALPHAT   : coefficient transverse de dispersivité               |
  18 : *  |                                                                 |
  19 : *  | QELEM    : vitesse au centre de chaque élément                  |
  20 : *  |-----------------------------------------------------------------|
  21 : *  |                         SORTIES                                 |
  22 : *  |-----------------------------------------------------------------|
  23 : *  |                                                                 |
  24 : *  |                                                                 |
  25 : *  | dif_disp : C'est le tenseur de dispersion                       |
  26 : *  |                                                                 |
  27 : *  |                                                                 |
  28 : **********************************************************************
  29 : 
  30 : 
  31 : 
  32 : * Seulement si présence de convection et présence de
  33 : * dispersivité et disp_l et disp_t
  34 : 
  35 : 'SI' (DISPERSI )                                                       ;
  36 : * vpetit = 1.D-20 sur le meme nombre de zone que DISPL
  37 :    vpetit =  'MANUEL' 'CHPO' ('EXTRAIRE' DISPL MAILLAGE) 'SCAL' 1.D-30 ;
  38 :    vpetit = 'NOMC' 'SCAL' vpetit 'NATU' DISCRET                        ;
  39 :    DISPL = 'NOMC' 'SCAL' DISPL                                         ;
  40 :    DISPT = 'NOMC' 'SCAL' DISPT                                         ;
  41 : 
  42 :  
  43 : *   on calcul la dispersivité
  44 : *   dimension 2
  45 :    'SI' ('EGA' ('VALEUR' 'DIME') 2)                                    ;
  46 :       V1 = 'EXCO' Qelem 'VX' 'SCAL'                                    ;
  47 :       V2 = EXCO Qelem 'VY' 'SCAL'                                      ;
  48 :       VC = (V1 * V1) '+' (V2 * V2)                                     ;
  49 :       VC = (VC**0.5)                                                   ;
  50 :       VC = VC + vpetit                                                 ;
  51 :       VCM1 = VC '**' -1                                                ;
  52 :       D11P = ((DISPL - DISPT) * V1*V1*VCM1) + (DISPT * VC)             ;
  53 :       D22P = ((DISPL - DISPT) * V2*V2*VCM1) + (DISPT * VC)             ; 
  54 :       D21P = (DISPL - DISPT) * V1*V2*VCM1                              ;   
  55 : 
  56 :       D11P = 'NOMC' K11 D11P                                           ;
  57 :       D21P = 'NOMC' K21 D21P                                           ; 
  58 :       D22P = 'NOMC' K22 D22P                                           ;
  59 :       dif_disp = D11P '+' D22P '+' D21P                                ;
  60 : 
  61 :    'FINSI'                                                             ;
  62 : *  dimension 3
  63 :    'SI' ('EGA' ('VALEUR' 'DIME') 3)                                    ;
  64 :       V1 = 'EXCO' Qelem 'VX' 'SCAL'                                    ;
  65 :       V2 = EXCO Qelem 'VY' 'SCAL'                                      ;
  66 :       V3 = EXCO Qelem 'VZ' 'SCAL'                                      ;
  67 :       VC = (V1 * V1) '+' (V2 * V2) '+' (V3 *  V3)                      ;
  68 :       VC = VC**0.5                                                     ;
  69 :       VC = VC + vpetit                                                 ;
  70 :       VCM1 =  VC '**' -1                                               ;
  71 :       D11P = ((DISPL - DISPT) * V1*V1*VCM1) + (DISPT * VC)             ;
  72 :       D22P = ((DISPL - DISPT) * V2*V2*VCM1) + (DISPT * VC)             ;
  73 :       D33P = ((DISPL - DISPT) * V3*V3*VCM1) + (DISPT * VC)             ;     
  74 :       D21P = (DISPL - DISPT) * V1*V2*VCM1                              ;
  75 :       D31P = (DISPL - DISPT) * V1*V3*VCM1                              ;
  76 :       D32P = (DISPL - DISPT) * V2*V3*VCM1                              ;
  77 : 
  78 :       D11P = 'NOMC' K11 D11P                                           ;
  79 :       D21P = 'NOMC' K21 D21P                                           ; 
  80 :       D22P = 'NOMC' K22 D22P                                           ;
  81 :       D31P = 'NOMC' K31 D31P                                           ;
  82 :       D32P = 'NOMC' K32 D32P                                           ; 
  83 :       D33P = 'NOMC' K33 D33P                                           ;
  84 : 
  85 :       dif_disp = D11P '+' D22P '+' D33P '+' D21P '+' D31P '+' D32P     ;        
  86 :    'FINSI'                                                             ;
  87 : 'FINSI'                                                                ;
  88 : 
  89 : 'FINPROC' dif_disp                                                     ;
  90 : 
  91 : 
  92 :