Procedures | Cast3M

   1 : * FILTREKE  PROCEDUR  MAGN      00/03/14    21:15:02     3811           
   2 : *
   3 : * VERSION DU 16-12-97
   4 : *
   5 : 
   6 :  DEBP FILTREKE ;
   7 :  ARGU RX*TABLE ;
   8 : * Filtre sur K et Epsilon
   9 : * - Echelle de vitesse (K**0.5) inférieure à une fraction (alfk)
  10 : *   de Uref (vitesse caractéristique) (alfk=1 pour l'instant)
  11 : *   Uref=max(UN,U0) 
  12 : *  - K > K0 
  13 : *  - Epsilon tel que l'echelle de longueur reste inférieure
  14 : *   à (L0/a) où L0 = diamètre enceinte et a=f(Re)
  15 : *   => Nut < Uref*L0/a
  16 : 
  17 :  rv=rx.'EQEX' ;
  18 :  iarg=rx.'IARG' ;
  19 :  NASTOK = rv.'NAVISTOK' ;
  20 : 
  21 :  'SI' (('EGA' NASTOK 0) ou ('EGA' IARG 2)) ;
  22 : 
  23 : * Ancien Filtre sur K et Epsilon
  24 : *- Echelle de vitesse (K**0.5) inférieure à une fraction (alfk)
  25 : * de U0 (vitesse caractéristique) (alfk=1 pour l'instant)
  26 : * K > 0  Epsilon tel que l'echelle de longueur reste inférieure
  27 : * à L0 diamètre enceinte
  28 : *   => Nut < U0*L0
  29 : 
  30 :  si( non ( ega iarg 2)) ;
  31 :  mess 'Procedure FILTREKE : nombre d arguments incorrect ' iarg ;
  32 :  quitter FILTREKE  ;
  33 :  finsi ;
  34 :  si ( ega ('TYPE' rx.'ARG1') 'MOT     ') ;
  35 :  U1=rv.'INCO'.(rx.'ARG1') ;
  36 :  sinon ;
  37 :  si ( ega ('TYPE' (rx.'ARG1')) 'FLOTTANT') ;
  38 :  U1=rx.'ARG1' ;
  39 :  sinon ;
  40 :  mess 'Procedure FILTREKE : type argument 1 invalide ' ;
  41 :  quitter FILTREKE ;
  42 :  finsi ;
  43 :  finsi ;
  44 : 
  45 :  si ( ega ('TYPE' rx.'ARG2') 'MOT     ') ;
  46 :  L0=rv.'INCO'.(rx.'ARG2') ;
  47 :  sinon ;
  48 :  si ( ega ('TYPE' (rx.'ARG2')) 'FLOTTANT') ;
  49 :  L0=rx.'ARG2' ;
  50 :  sinon ;
  51 :  mess 'Procedure FILTREKE : type argument 2 invalide ' ;
  52 :  quitter FILTREKE ;
  53 :  finsi ;
  54 :  finsi ;
  55 : 
  56 : 
  57 :  nic=dime (rx.'LISTINCO') ;
  58 :  si( non ( ega nic 2)) ;
  59 :  mess 'Procedure FILTREKE : nombre d inconnues incorrect ' nic ;
  60 :  quitter FILTREKE ;
  61 :  finsi ;
  62 :  nomi1=extr 1 (rx.'LISTINCO');
  63 :  nomi2=extr 2 (rx.'LISTINCO');
  64 :  nom1= mot (text (chai nomi1)) ;
  65 :  nom2= mot (text (chai nomi2)) ;
  66 : 
  67 :  en=rv.'INCO'.nom2 ;
  68 :  kn=rv.'INCO'.nom1 ;
  69 : *Rec=100.;
  70 :  k0 = 1.e-10 ;
  71 : 
  72 : *mdu2=kops u1 'PSCA' u1 ;
  73 : *mdu=mdu2 ** 0.5 ;
  74 : *Re=kops (kops (kops mdu '*' L0) '/' nu) '+' (Rec / 10.) ;
  75 : *a= exp (kops Rec '/' Re ) ;
  76 : 
  77 :  mdu2=u1*u1 ;
  78 :  a=0.09 ;
  79 : 
  80 :  kn=kops kn '|<' k0 ;
  81 :  kn=kops kn '>|' mdu2 ;
  82 : 
  83 :  E0= kops (kops kn '**' 1.5) '*' (a / L0) ;
  84 :  en=kops en '|<' E0 ;
  85 : 
  86 :  rv.'INCO'.'EN'=en ;
  87 :  rv.'INCO'.'KN'=kn ;
  88 : 
  89 :  si( non ( ega NASTOK 0 ) ) ;
  90 :  as2 ama1  = 'KOPS' 'MATRIK' ;
  91 :  RESPRO as2 ama1 ;
  92 :  finsi ;
  93 : 
  94 : 'QUITTER' FILTREKE ;
  95 : 'FINSI'  ;
  96 : 
  97 :  si( non ( ega iarg 4)) ;
  98 :  mess 'Procedure FILTREKE : nombre d arguments incorrect ' iarg ;
  99 :  quitter FILTREKE  ;
 100 :  finsi ;
 101 :  si ( ega ('TYPE' rx.'ARG1') 'MOT     ') ;
 102 :  U1=rv.'INCO'.(rx.'ARG1') ;
 103 :  sinon ;
 104 :  si ( ega ('TYPE' (rx.'ARG1')) 'FLOTTANT') ;
 105 :  U1=rx.'ARG1' ;
 106 :  sinon ;
 107 :  mess 'Procedure FILTREKE : type argument 1 invalide ' ;
 108 :  quitter FILTREKE ;
 109 :  finsi ;
 110 :  finsi ;
 111 : 
 112 :  si ( ega ('TYPE' rx.'ARG2') 'MOT     ') ;
 113 :  L0=rv.'INCO'.(rx.'ARG2') ;
 114 :  sinon ;
 115 :  si ( ega ('TYPE' (rx.'ARG2')) 'FLOTTANT') ;
 116 :  L0=rx.'ARG2' ;
 117 :  sinon ;
 118 :  mess 'Procedure FILTREKE : type argument 2 invalide ' ;
 119 :  quitter FILTREKE ;
 120 :  finsi ;
 121 :  finsi ;
 122 : 
 123 :  si ( ega ('TYPE' rx.'ARG3') 'MOT     ') ;
 124 :  NU=rv.'INCO'.(rx.'ARG3') ;
 125 :  sinon ;
 126 :  si ( ega ('TYPE' (rx.'ARG3')) 'FLOTTANT') ;
 127 :  NU=rx.'ARG3' ;
 128 :  sinon ;
 129 :  mess 'Procedure FILTREKE : type argument 3 invalide ' ;
 130 :  quitter FILTREKE ;
 131 :  finsi ;
 132 :  finsi ;
 133 : 
 134 :  si ( ega ('TYPE' rx.'ARG4') 'MOT     ') ;
 135 :  UN=rv.'INCO'.(rx.'ARG4') ;
 136 :  sinon ;
 137 :  si ( ega ('TYPE' (rx.'ARG4')) 'CHPOINT') ;
 138 :  UN=rx.'ARG4' ;
 139 :  sinon ;
 140 :  mess 'Procedure FILTREKE : type argument 4 invalide ' ;
 141 :  quitter FILTREKE ;
 142 :  finsi ;
 143 :  finsi ;
 144 : 
 145 : 
 146 :  nic=dime (rx.'LISTINCO') ;
 147 :  si( non ( ega nic 2)) ;
 148 :  mess 'Procedure FILTREKE : nombre d inconnues incorrect ' nic ;
 149 :  quitter FILTREKE ;
 150 :  finsi ;
 151 : 
 152 :  nomi1=extr 1 (rx.'LISTINCO');
 153 :  nomi2=extr 2 (rx.'LISTINCO');
 154 :  nom1= mot (text (chai nomi1)) ;
 155 :  nom2= mot (text (chai nomi2)) ;
 156 : 
 157 :  en=rv.'INCO'.nom2 ;
 158 :  kn=rv.'INCO'.nom1 ;
 159 : 
 160 :  Rec=100.;
 161 :  k0 = 1.e-10 ;
 162 :  cnu=0.09;
 163 : 
 164 :  mdu=kops un 'PSCA' un ;
 165 :  mdu=mdu ** 0.5 ;
 166 : 
 167 :  Re=kops (kops (kops mdu '*' L0) '/' nu) '+' (Rec / 10.) ;
 168 :  a= exp (kops Rec '/' Re ) ;
 169 : 
 170 :  mdu = kops mdu '|<' u1 ;
 171 :  mdu2= kops mdu '*' mdu ;
 172 : 
 173 :  kn=kops kn '|<' k0 ;
 174 :  kn=kops kn '>|' mdu2 ;
 175 : 
 176 :  E0= kops (kops kn '**' 1.5) '*' (a / L0) ;
 177 :  en=kops en '|<' E0 ;
 178 : 
 179 :  rv.'INCO'.nom2=en ;
 180 :  rv.'INCO'.nom1=kn ;
 181 : 
 182 :  as2 ama1  = 'KOPS' 'MATRIK' ;
 183 :  RESPRO as2 ama1 ;
 184 : 
 185 :  FINPROC ;
 186 :