Télécharger la plaquette Cast3M


fr en

Télécharger htc_wtr.procedur

Retour à la liste

Numérotation des lignes : 

  1. * HTC_WTR   PROCEDUR  AM        11/09/13    21:15:04     7126           
  2. *-----------------------------------------------------------------------
  3. *23456789012345678901234567890123456789012345678901234567890123456789012
  4. *-----------------------------------------------------------------------
  5. *-------------------------------------------------
  6. * HTC_WTR Properties of water 
  7. * Pressione di saturazione, densita,
  8. * calore latente di evaporazione.
  9. * Date: 16-03-98
  10. *---------------------------------------------------------------
  11. 'DEBP' HTC_WTR TAV1*'TABLE' M1*'MOT' TA*'CHPOINT' PA/'CHPOINT' ;
  12. *---------------------------------------------------------------
  13. *
  14.  FLG1= TAV1.'FLG1';
  15. 'SI' FLG1; 'MESS' 'PROCEDURE-HTC_WTR'; 'FINSI'; 
  16. 'SI'( 'EGA' M1 'SATU');
  17. *  RESO1 = TA 'MASQUE' 'EGINFE' 374.15;
  18. *  RESO2 = 1. - RESO1 ;
  19. *
  20. * Saturation pressure in Mpa
  21. *
  22. 'SI'('NON'('EXIS' TAV1 'EVPS')); 
  23. TPR1= 'PROG' 0. 4. 5. 6. 8. 10. 'PAS' 1. 36.
  24.             38. 40. 'PAS' 5. 100. 'PAS' 10. 300.
  25.             320. 340. 360. 374.14 2000.;
  26. PS21='PROG' .000611 .000813 .000872 .000935 .001072 .001228 .001312
  27.     .001402 .001497 .001598 .001705 .001818 .001938 .002064 .002198
  28.     .002339 .002487 .002645 .002810 .002985 .003169 .003363 .003567;
  29. PS22='PROG' .003782 .004008 .004246 .004496 .004759 .005034 .005324
  30.     .005628 .005947 .006632 .007384 .009593 .01235  .01576  .01994    
  31.     .02503  .03119  .03858  .04739  .05783  .07014  .08455  .10133 
  32.     .1433   .1985   .2701   .3613   .4758   .6178   .7917  1.002    
  33.    1.254   1.554   1.906   2.318   2.795   3.344   3.973   4.688   
  34.    5.499   6.412   7.436   8.581  11.27   14.59   18.65   22.09  22.09;
  35. PS2=PS21 'ET' PS22;
  36. TAV1.'EVPS'= 'EVOL' 'MANU' TPR1 'SCAL' PS2 'SCAL' ;
  37. 'FINSI';
  38. OUT= 'NOMC' 'SCAL' ('VARI' (TAV1.'EVPS') TA);           
  39. *
  40. *  AAA=33661.805;
  41. *  BBB=4744.058;
  42. *  OUT1= RESO1 *(AAA*('EXP'(((TA + 273.15)**-1)*(-1.*BBB))));
  43. *  OUT2= RESO2 *(AAA*('EXP'(((647.3)**-1)*(-1.*BBB)))) ;
  44. *  OUT = OUT1 + OUT2; 
  45. * 
  46. 'SINON'; 
  47.  'SI'( 'EGA' M1 'RHOW');
  48. * Density as fun. of T and P
  49.    A0= 4.8863E-7 ;A1=-1.6528E-9 ;A2= 1.8621E-12;
  50.    A3= 2.4266E-13;A4=-1.5996E-15;A5= 3.3703E-18;
  51.    B0= 1.0213E3  ;B1=-7.7377E-1 ;B2= 8.7696E-3 ;
  52.    B3=-9.2118E-5 ;B4= 3.3534E-7 ;B5=-4.4034E-10;
  53.    PR= 2.E7;
  54.    B0=B0 - (PR*A0);B1=B1 - (PR*A1);B2=B2 - (PR*A2);
  55.    B3=B3 - (PR*A3);B4=B4 - (PR*A4);B5=B5 - (PR*A5);
  56.    RO1=TA*(TA*(TA*(TA*(TA*B5+B4)+B3)+B2)+B1)+B0;
  57.    OUT=PA*(TA*(TA*(TA*(TA*(TA*A5+A4)+A3)+A2)+A1)+A0)+RO1;
  58. * mess 'density' (maxi out);
  59.  'SINON';
  60.   'SI'( 'EGA' M1 'HFG');
  61. * Heat of vaporization J/Kg
  62. *
  63.     RESO1 = TA 'MASQUE' 'EGINFE' 374.15;           
  64. *
  65.     OUT   = RESO1 * 350000.*(('ABS'(374.15 - TA))**(1./3.));
  66.   'SINON';
  67.    'SI'( 'EGA' M1 'CW');
  68. * Specific heat J/Kg            
  69.      OUT   = 4100. ;
  70.    'FINSI';
  71.   'FINSI'; 
  72.  'FINSI';
  73. 'FINSI'; 
  74. *
  75. 'FINPROC' OUT;
  76.  
  77.  

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