Télécharger la plaquette Cast3M


fr en

Télécharger dzvdt.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  dzvdt.dgibi
  2. * -------------- VARI option DZVDT -------------------------------
  3. * Test de l'opérateur VARI DZVDT(P,T)
  4. * dervivee partielle du facteur de compressibilite de la vapeur
  5. * par rapport a la temperature
  6. * Les données sont un FLOTTANT, un LISTREEL ou un CHPO
  7. * notation : DT**** derivee partielle
  8. * notation : DC**** derivee calculee par accroissement
  9.  
  10. 'OPTI' 'DIME' 2 'ELEM' 'QUA4' ;
  11. 'OPTI' 'ECHO' 0 ;
  12. DELTAT=0.1D0 ;
  13. DELTAP=1000D0;
  14. *---------------- FLOTTANT----------------------------------------
  15. P1 = 1.0133D5 ;
  16. T1=373.15 ; 
  17. P2 = 5.021D5 ;
  18. T2=425.15; 
  19. P3 = 10.027D5;
  20. T3=453.15 ; 
  21.  
  22. DTPOINT1 = VARI DZVDT P1 T1 ;
  23. DTPOINT2 = VARI DZVDT P2 T2 ;
  24. DTPOINT3 = VARI DZVDT P3 T3 ;
  25.  
  26. R2=(VARI ROVAP (P1) (T1+DELTAT));
  27. R1=(VARI ROVAP (P1) T1);
  28. R4=(VARI ROVAP (P2) (T2+DELTAT));
  29. R3=(VARI ROVAP (P2) T2);
  30. R6=(VARI ROVAP (P3) (T3+DELTAT));
  31. R5=(VARI ROVAP (P3) T3);
  32.  
  33. DCPOINT4=((VARI ZVAP R2 (T1+DELTAT))-(VARI ZVAP R1 T1))/DELTAT;
  34. DCPOINT5=((VARI ZVAP R4 (T2+DELTAT))-(VARI ZVAP R3 T2))/DELTAT;
  35. DCPOINT6=((VARI ZVAP R6 (T3+DELTAT))-(VARI ZVAP R5 T3))/DELTAT;
  36.  
  37.  
  38.  
  39.  
  40. *----------------- LISTREEL----------------------------------------
  41.  
  42. T1 = 'PROG' 373.15D0 425.15D0 453.15D0 ;
  43. P1 = 'PROG' 1.0133D5 5.021D5 10.027D5 ;
  44. DELTAP2 = 'PROG' 1000D0 1000D0 1000D0 ;
  45. DELTAT2 = 'PROG' 0.1D0 0.1D0 0.1D0 ;
  46. R6= VARI ROVAP P1 T1;
  47. R7 = VARI ROVAP (P1) (T1+DELTAT2) ;
  48. DTLIST1 = VARI DZVDT P1 T1 ;
  49.  
  50. DCLIST1= ((VARI ZVAP (R7) (T1+DELTAT2))-(VARI ZVAP R6 T1))/DELTAT2;
  51. *-----------------CHAMPOINT-----------------------------------------
  52. P11 = 1.0133D5 0. ;
  53. P22 = 5.021D5 0. ;
  54. P33 = 10.027D5 0. ;
  55. P11P33 = P11 'DROI' 1 P22 'DROI' 1 P33 ;
  56. P11P33_1 = CHAN 'POI1' P11P33 ;
  57. CT1= MANU 'CHPO' P11P33_1 'REAL' T1 ;
  58. CP1= MANU 'CHPO' P11P33_1 'REAL' P1;
  59. CT2 = MANU 'CHPO' P11P33_1 'REAL' (T1+DELTAT2);
  60. R8= VARI ROVAP CP1 CT1 ;
  61. R9 = VARI ROVAP CP1 CT2 ;
  62. DCCHPO1= ((VARI ZVAP R9 CT2)-(VARI ZVAP R8 CT1))/DELTAT;
  63. DTCHPO1 = VARI DZVDT CP1 CT1 ;
  64. *--------------------------------------------------------------------
  65. *----RESULTATS ET STOCKAGES DES DONNEES------------------
  66. *----------extraction des listreal-----------------------
  67. LPOINT1=  'EXTR' DTLIST1 1 ;
  68. LPOINT2 = 'EXTR' DTLIST1 2 ;
  69. LPOINT3=  'EXTR' DTLIST1 3 ;
  70. LPOINT4 = 'EXTR' DCLIST1 1 ;
  71. LPOINT5 = 'EXTR' DCLIST1 2 ;
  72. LPOINT6=  'EXTR' DCLIST1 3 ;
  73. *----------extraction des champoints---------------------
  74. CPOINT4 = 'EXTR' DCCHPO1 'REAL' P11 ;
  75. CPOINT5 = 'EXTR' DCCHPO1 'REAL' P22 ;
  76. CPOINT6 = 'EXTR' DCCHPO1 'REAL' P33 ;
  77. CPOINT1 = 'EXTR' DTCHPO1 'REAL' P11 ;
  78. CPOINT2 = 'EXTR' DTCHPO1 'REAL' P22 ;
  79. CPOINT3=  'EXTR' DTCHPO1 'REAL' P33 ;
  80. *-----------------------------------------------------------------------
  81. *----CALCULS VALEURS REFERENCES XREFi CASTEM----------------------------
  82. *-----------------------------------------------------------------------
  83. XREF1=(DCPOINT4+LPOINT4+CPOINT4)/3;
  84. XREF2=(DCPOINT5+LPOINT5+CPOINT5)/3;
  85. XREF3=(DCPOINT6+LPOINT6+CPOINT6)/3;
  86. *MESS DCPOINT4 LPOINT4 CPOINT4;
  87.  
  88.  
  89. *------------CALCULS DES ERREURS---------------------------------------
  90. CTRL1 = DTPOINT1+  LPOINT1 + CPOINT1 / 3.d0 - DTPOINT1 / DTPOINT1 ;
  91. CTRL2 = DTPOINT2 + LPOINT2 + CPOINT2 / 3.d0 - DTPOINT2 / DTPOINT2 ;
  92. CTRL3 = DTPOINT3 + LPOINT3 + CPOINT3 / 3.d0 - DTPOINT3 / DTPOINT2 ;
  93. CTRL4 = CTRL1 + CTRL2 + CTRL3 ;
  94. CTRL5 = DTPOINT1 + LPOINT1 + CPOINT1 / 3.d0 - XREF1 / XREF1 ;
  95. CTRL6 = DTPOINT2 + LPOINT2 + CPOINT2 / 3.d0 - XREF2 / XREF2 ;
  96. CTRL7 = DTPOINT3 + LPOINT3 + CPOINT3 / 3.d0 - XREF3 / XREF3 ;
  97. CTRL8 = CTRL5 + CTRL6 + CTRL7 ;
  98. *-------------------------------------------> Affichage
  99. 'MESS' '     ' ;
  100. 'MESS' '     ' ;
  101. 'MESS' 'Pt 1  VARI DZVDT ----->' DTPOINT1 LPOINT1 CPOINT1 XREF1 ;
  102. 'MESS' 'Pt 2  VARI DZVDT ----->' DTPOINT2 LPOINT2 CPOINT2 XREF2 ;
  103. 'MESS' 'Pt 3  VARI DZVDT ----->' DTPOINT3 LPOINT3 CPOINT3 XREF3 ;
  104. 'MESS' '     ' ;
  105. 'MESS' 'Comparaison calculs rel.----->' CTRL4 ;
  106. 'MESS' 'Erreur rel. VDI         ----->' CTRL8 ;
  107. 'MESS' '     ' ;
  108. 'MESS' '     ' ;
  109. 'MESS' '     ' ;
  110. 'MESS' '     ' ;
  111. 'MESS' '     ' ;
  112. *------------------------------------------> Compte-rendu et sortie
  113. EPS4  = 1.E-14 ;
  114. EPS8  = 3.E-3  ;
  115. CTRL4 = ABS CTRL4 ;
  116. CTRL8 = ABS CTRL8 ;
  117. LOG4 = CTRL4 > EPS4 ;
  118. LOG8 = CTRL8 > EPS8 ;
  119. L0   = LOG4 'OU' LOG8 ;
  120. 'SI' L0        ;
  121.    'ERREUR' 5  ;
  122. 'SINON'        ;
  123.    'ERREUR' 0  ;
  124. 'FINSI'        ;
  125. 'FIN' ;
  126.  
  127.  
  128.  
  129.  
  130.  

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