Télécharger la plaquette Cast3M


fr en

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