Télécharger la plaquette Cast3M


fr en

Télécharger testkfpt.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  testkfpt.dgibi
  2. opti dime 2 elem QUA4;
  3.  
  4. nx = 1000;ny = 5;
  5.  
  6. p0=0 0;
  7. p1 = 30. 0. ;
  8. co1= p0 d dini 0.001 dfin 1. p1;
  9.  
  10. DISCR= QUAF;
  11. *DISCR= MACRO;
  12.  
  13. COMPLET = FAUX ;
  14. Si (NON COMPLET);
  15.  DISCR= LINE;
  16. p1 = 1. 0. ;
  17. co1= p0 d 5 p1;
  18. FINSI;
  19.  
  20.  
  21.  
  22. Mco1 = chan QUAF co1;
  23.  
  24. $co1 = mode Mco1  'NAVIER_STOKES' DISCR;
  25.  
  26. co1 = doma $co1  maillage ;
  27.  
  28. YP  = 0.001;
  29.  
  30.  
  31.  uet = ((coor 1 co1)*0.5+1.e-5);
  32.  teta= (coor 1 co1);
  33.  RO =  1.  ;
  34.  CP =  1.  ;
  35.  mu = 1.e-6 ;
  36.  lb = 1.e-6;
  37.  alpha = lb / ro*cp ;
  38.  h =  kfpt $co1 RO mu  CP lb uet yp ;
  39.  yplus = ro * yp * uet / mu;
  40.  tetap = ro * cp * uet *(inve h);
  41.  y1 = extr  (evol chpo yplus co1) 'ORDO';
  42.  t1 = extr  (evol chpo tetap co1) 'ORDO';
  43.  ev1 = evol 'MANU' 'Y+' y1    'TETA+' t1    ;
  44.  pr = mu/alpha;
  45.  Si COMPLET;
  46.  dess ev1  TITR (chai 'Prandtl ' pr)'XBOR' 1. 1.e3 'LOGX';
  47.  Finsi ;
  48.  
  49.  
  50.  
  51.  mu = 1.e-7;
  52.   h =  kfpt $co1 RO mu  CP lb uet yp ;
  53.  yplus = ro * yp * uet / mu;
  54.  tetap = ro * cp * uet *(inve h);
  55.  y1 = extr  (evol chpo yplus co1) 'ORDO';
  56.  t1 = extr  (evol chpo tetap co1) 'ORDO';
  57.  ev01 = evol 'MANU' 'Y+' y1    'TETA+' t1    ;
  58.  pr = mu/alpha;
  59.  
  60.  
  61.  mu = 1.e-5;
  62.   h =  kfpt $co1 RO mu  CP lb uet yp ;
  63.  yplus = ro * yp * uet / mu;
  64.  tetap = ro * cp * uet *(inve h);
  65.  y1 = extr  (evol chpo yplus co1) 'ORDO';
  66.  t1 = extr  (evol chpo tetap co1) 'ORDO';
  67.  ev10 = evol 'MANU' 'Y+' y1    'TETA+' t1    ;
  68.  pr = mu/alpha;
  69.  
  70.  
  71.  
  72.  Si COMPLET;
  73.  dess (ev1 et ev01 et ev10) 'XBOR' 1.  1.e3 'LOGX';
  74.  Finsi ;
  75.  
  76.  
  77. FIN ;
  78.  
  79.  
  80.  
  81.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales