Télécharger prodt.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : prodt.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Teste la procedure PRODT production d'energie turbulente
  5. *
  6.  
  7. opti dime 3 elem cub8 ;
  8.  
  9. DISCR = 'MACRO' ;
  10.  
  11. a = 1.13 ;
  12. b = 2.16 ;
  13. n1= 2;
  14. n2= 3;
  15. n3= 5;
  16. P0 = 0 0 0;
  17. p1 = a 0 0;
  18. p2 = a b 0;
  19. p3 = 0 b 0;
  20. cmt1= p0 d n1 p1 d n2 p2 d n1 p3 d n2 p0;
  21. mt1 = surf 'PLAN' cmt1 ;
  22. mt= mt1 volu n3 'TRAN' (0 0 3.51);
  23. Mmt= chan mt QUAF ;
  24. $mt= mode Mmt 'NAVIER_STOKES' DISCR;
  25. mt=doma $mt maillage ;
  26. x=coor 1 mt ;
  27. y=coor 2 mt ;
  28. z=coor 3 mt ;
  29.  
  30. u1= x + (2.*y) + (3.*z);
  31. u2= (4.*x) + (5.*y) + (6.*z);
  32. u3= (7.*x) + (8.*y) + (9.*z);
  33. u=(nomc 'UX' u1) et (nomc 'UY' u2) et (nomc 'UZ' u3);
  34. t=(8*x) + (16*y) + (32*z) ;
  35.  
  36.  
  37. ** Cas ou dt/dz >0 est de meme signe que gb
  38. ** Stratification instable => G=0.
  39. gb=(1 1 1);
  40. P=prodt u $mt gb t ;
  41. mi= mini P ;
  42. ma= maxi P ;
  43. d=abs ( mi + ma - 1092.) ;
  44. mess mi ma 'Solution : 546 d=' d ;
  45. Si(d > 1.e-10);
  46. erreur 5 ;
  47. Finsi ;
  48.  
  49.  
  50. ** Cas ou dt/dz >0 est de signe contraire a gb (<0)
  51. ** Stratification stable => (P + G)(1 + C3*Rf)
  52. ** Rf=-G /(P+G)
  53. gb=(1 1 1)*(-1.);
  54. P=prodt u $mt gb t ;
  55. mi= mini P ;
  56. ma= maxi P ;
  57. d=abs ( mi + ma - 1060.) ;
  58. mess mi ma 'Solution : 530 d=' d ;
  59. Si(d > 1.e-10);
  60. erreur 5 ;
  61. Finsi ;
  62. FIN ;
  63.  
  64.  
  65.  
  66.  
  67.  

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