Télécharger asparam.procedur

Retour à la liste

Numérotation des lignes :

  1. * ASPARAM PROCEDUR MAGN 08/12/09 21:15:00 6189
  2. 'DEBPROC' asparam rxt*TABLE ;
  3.  
  4. TBT = rxt.'TBT' ;
  5. GEO = rxt.'GEO' ;
  6. *
  7. idim = vale 'DIME';
  8. DIM3D=FAUX ;
  9. Si(EGA idim 3) ;
  10. DIM3D=VRAI ;
  11. Finsi ;
  12. * MESS 'si VRAI dimension 3D' ;
  13. * LIST DIM3D ;
  14. *
  15. un = rxt.'TIC'.'UN' ;
  16. vn = rxt.'TIC'.'VN' ;
  17. xd = rxt.'TIC'.'XD' ;
  18. xdinj = 'IPOL' 0.0 t (rxt.'scenasp'.'xdinj') ;
  19. xdmin = xdinj * 0.01 ;
  20. ixd=(masq xd 'SUPERIEUR' xdmin ) ;
  21. tf = rxt.'TIC'.'TF' ;
  22. td = rxt.'TIC'.'TD' ;
  23. dd = rxt.'TIC'.'DD' ;
  24. t = rxt.'scenasp'.'t' ;
  25. ddinj = 'IPOL' 0.0 t (rxt.'scenasp'.'ddinj') ;
  26. ddmin = ddinj * 0.01 ;
  27. idd=(masq dd 'SUPERIEUR' ddmin ) ;
  28. dd = rxt.'TIC'.'DD' *idd ;
  29. rhogas = rxt.'TIC'.'RHO' ;
  30. rhovap = rxt.'TIC'.'RVP' ;
  31. mugas = rxt.'TIC'.'Mu' ;
  32. Cpgas = extr (rxt.'TIC'.'Cpm') (dime (rxt.'TIC'.'Cpm')) ;
  33. lamgas = TBT.'lamair' ;
  34. rhodis = rxt.'rod' ;
  35. Cpdis = rxt.'Cpd' ;
  36. $vtf = rxt.'GEO'.'$vtf' ;
  37. Lv = TBT.'Lv' ;
  38. dt = rxt.'DT0' ;
  39. *
  40. *corrections appel procedur rvst2.procedur
  41. * rhosat rhosp = rvsat td ;
  42. rgva = TBT.'Rgvap' ;
  43. Pvap = rgva * rhovap * (tf+273.0) + 0.000005 ;
  44. ind200= td masque 'EGINF' 200. ;
  45. ind20sup = td masque 'EGSUP' 20. ;
  46. indch= ind20sup * ind200 ;
  47. *
  48. **
  49. td1 = td * indch ;
  50. rhosat1 rhosp1 = rvsat td1 ;
  51. rhosp1 = rhosp1 * indch ;
  52. rhosat1 = rhosat1 * indch ;
  53. ind20inf = td masque 'EGINF' 19.999999999999999 ;
  54. ind0sup = td masque 'EGSUP' 0. ;
  55. indfr = ind20inf * ind0sup ;
  56. td2 = td * indfr ;
  57. rhosat2 rhosp2 = rvst2 td2 ;
  58. rhosat2 = rhosat2 * indfr ;
  59. rhosp2 = rhosp2 * indfr ;
  60. rhosat = rhosat1 + rhosat2 ;
  61. rhosp = rhosp1 + rhosp2 ;
  62. *
  63. *fin correction appel procedur
  64. *
  65. Pt = extr (dime rxt.'TIC'.'PT') rxt.'TIC'.'PT' ;
  66. Dvap = (((tf + 273.0)/273.0) ** 1.75) * 2.2 / Pt ;
  67. *
  68. drho = noel $vtf ((rhodis - rhogas)/rhodis) ;
  69. 'SI' DIM3D ;
  70. gdrx = nomc 'UX' (drho * 0.00) ;
  71. gdry = nomc 'UY' (drho * 0.00) ;
  72. gdrz = nomc 'UZ' (drho * 9.81) ;
  73. rxt.'TIC'.'GRD' = kcht $vtf vect centre (gdrx et gdry et gdrz) ;
  74. 'SINON' ;
  75. gdrx = nomc 'UX' (drho * 0.00) ;
  76. gdrz = nomc 'UY' (drho * 9.81) ;
  77. rxt.'TIC'.'GRD' = kcht $vtf vect centre (gdrx et gdrz) ;
  78. FINSI;
  79. *
  80. cmps = extr un comp ;
  81. ur = (psca (vn - un) (vn - un) cmps cmps) ** 0.5 ;
  82. Re = dd * ur * rhogas * (mugas** -1.) ;
  83. Pr = mugas * Cpgas * (lamgas** -1.) ;
  84. Sc = mugas * (rhogas ** -1.0) * (Dvap ** -1.0) ;
  85. * modif 1 2004
  86. Nu = 2.00 + (0.599 * ((Pr ** 0.33333) * (Re ** 0.5))) ;
  87. * Nu = (2.00 + 0.599 * ((Pr ** 0.33333) * (Re ** 0.5))) ;
  88. * modif 2 2004
  89. Sh = 2.00 + (0.599 * ((Sc ** 0.33333) * (Re ** 0.5))) ;
  90. * Sh = (2.00 + 0.599 * ((Sc ** 0.33333) * (Re ** 0.5))) ;
  91. *
  92. F_D = (1.00 + (0.15 * (Re ** 0.687))) ;
  93. *
  94. fk = 18.0 * mugas * F_D * (dd ** -2.0) ;
  95. fh = 6.0 * lamgas * Nu * (dd ** -2.0) ;
  96. fm = 6.0 * Dvap * Sh * (dd ** -2.0) * (rhovap - rhosat) ;
  97.  
  98.  
  99.  
  100. *fmax = (tf - td) * ((rhodis * Cpdis / dt) - (fh / Cpgas)) / Lv;
  101. *fmin = (-1.0*td) * ((rhodis * Cpdis / dt) - (fh / Cpgas)) / Lv;
  102. *
  103. num = ((rhovap - rhosat)/dt) * ((xd + 1.0e-10) ** -1.0) ;
  104. den = (rhosp * ((xd + 1.0e-10) ** -1.0)) ;
  105. den = (den * Lv / (Cpdis * rhodis)) + 1.0 ;
  106. fmax = abs (num * (den ** -1.0)) ;
  107. fmin = -1.0 * fmax ;
  108. *
  109. f1 = fm masque superieur fmax ;
  110. f2 = fm masque inferieur fmin ;
  111. f3 = (1.0 - f1) * (1.0 - f2) ;
  112. fm = (fmax * f1) + (fmin * f2) + (fm * f3) ;
  113. * fm = fm * 0.0 ;
  114.  
  115. kgas = fk * xd * (rhogas ** -1.0) * ixd ;
  116. kdis = fk / rhodis ;
  117. hgas = fh * xd / Cpgas * (rhogas ** -1.0) * ixd ;
  118. hdis = fh / Cpdis / rhodis ;
  119. sdia = fm * dd / 3.0 / rhodis ;
  120. sxd = fm * xd / rhodis ;
  121. srvp = fm * xd * (-1.0) ;
  122. *srho = fm * xd * (-1.0) ;
  123. std = fm * (Lv) / Cpdis / rhodis ;
  124.  
  125.  
  126. * rxt.'TIC'.'GRD' = kcht $vtf vect centre (gdrx et gdry et gdrz) ;
  127. rxt.'TIC'.'KGAS' = kcht $vtf scal sommet kgas ;
  128. rxt.'TIC'.'KDIS' = kcht $vtf scal sommet kdis ;
  129. rxt.'TIC'.'HGAS' = kcht $vtf scal sommet hgas ;
  130. rxt.'TIC'.'HDIS' = kcht $vtf scal sommet hdis ;
  131. rxt.'TIC'.'SDIA' = noel $vtf (kcht $vtf scal sommet sdia) ;
  132. rxt.'TIC'.'SXD' = noel $vtf (kcht $vtf scal sommet sxd ) ;
  133. *rxt.'TIC'.'SRHO' = noel $vtf (kcht $vtf scal sommet srho) ;
  134. rxt.'TIC'.'SRVP' = noel $vtf (kcht $vtf scal sommet srvp) ;
  135. rxt.'TIC'.'STD' = noel $vtf (kcht $vtf scal sommet std ) ;
  136. rxt.'TIC'.'fm' = kcht $vtf scal sommet fm ;
  137. rxt.'TIC'.'fmax' = kcht $vtf scal sommet fmax ;
  138.  
  139.  
  140.  
  141. * mess 'Mini/Max kdis ' (mini kdis) (maxi kdis) ;
  142. * mess 'Mini/Max hdis ' (mini hdis) (maxi hdis) ;
  143. * mess 'Mini/Max kgas ' (mini kgas) (maxi kgas) ;
  144. * mess 'Mini/Max hgas ' (mini hgas) (maxi hgas) ;
  145. * mess 'Mini/Max fk ' (mini fk ) (maxi fk ) ;
  146. * mess 'Mini/Max fh ' (mini fh ) (maxi fh ) ;
  147. * mess 'Mini/Max fm ' (mini fm ) (maxi fm ) ;
  148. * mess 'Mini/Max srvp ' (mini srvp) (maxi srvp) ;
  149. * mess 'Mini/Max xd ' (mini xd ) (maxi xd ) ;
  150. * mess 'Mini/Max rhogas' (mini rhogas) (maxi rhogas) ;
  151. * mess 'Cpgas' Cpgas 'Lv=' lv ;
  152. * mess 'Mini/Max mugas' (mini mugas) (maxi mugas) ;
  153. * mess 'Mini/Max lamgas' lamgas ;
  154. * mess 'Mini/Max Dvap ' (mini Dvap ) (maxi Dvap ) ;
  155.  
  156.  
  157. 'FINPROC' ;
  158.  
  159.  

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