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