Télécharger la plaquette Cast3M


fr en

Télécharger mdrecomb.procedur

Retour à la liste

Numérotation des lignes : 

  1. * MDRECOMB  PROCEDUR  SERRE     21/06/28    21:15:04     11045          
  2. 'DEBPROC' MDRECOMB RXT*'TABLE   ' ;
  3. *
  4. *
  5. * Couplage du modèle 0D du PAR avec le modèle fluide CFD
  6. *
  7. * E/    : rxt    : Table contenant tous les data pour les fluides
  8. *  /S   : Qo_i   : Différence de débit entre la sortie et l'entrée des PARs
  9. *                  pour l'espèce i (en kg/s)
  10. *  /S   : drhpar : Différence de débit enthalpique entre la sortie et
  11. *                  l'entrée des PARs (en W)
  12. *
  13. *  /S : Sauvegarde pour post-traitement local dans tic . 'RECOMB' . &BCLREC
  14. *  /S : Sauvegarde pour C.L. de l'inconnue 'XXX' dans tic . 'RECOMB' . 'XXX'
  15. *
  16. *
  17. 'MESS' '$$$$ RECOMBINEUR $$$$' ;
  18. tbt = RXT . 'TBT' ;
  19. tic = RXT . 'TIC' ;
  20. geo = RXT . 'GEO' ;
  21. *
  22. idim = 'VALE' 'DIME' ;
  23. 'SI' ('EGA' idim 2) ;
  24.      LMOT1 = 'MOTS' '1UN' '2UN' ;
  25. 'SINO' ;
  26.      LMOT1 = 'MOTS' '1UN' '2UN' '3UN' ;
  27. 'FINS' ;
  28. Dt  = tic . 'DT' ;
  29. Tps = tic . 'Tps' ;
  30. Pe = 'EXTR' tic . 'PT' ('DIME' tic . 'PT') ;
  31. *
  32. * Initialisation pour les bilans de masse et d'energie 0D
  33. QIN_H2   = 0. ;
  34. QIN_H2O  = 0. ;
  35. QIN_O2   = 0. ;
  36. QIN_N2   = 0. ;
  37. QIN_HE   = 0. ;
  38. QIN_CO2  = 0. ;
  39. QIN_CO   = 0. ;
  40. *
  41. QOUT_H2  = 0. ;
  42. QOUT_H2O = 0. ;
  43. QOUT_O2  = 0. ;
  44. QOUT_N2  = 0. ;
  45. QOUT_HE  = 0. ;
  46. QOUT_CO2 = 0. ;
  47. QOUT_CO  = 0. ;
  48. *
  49. drhpar   = 0. ;
  50. *
  51. * Initialisation des conditions aux limites scalaire à la sortie des PARs
  52. cltf mat1 = 'KOPS' 'MATRIK' ;
  53. clrvp  = 'COPI' cltf ;
  54. clrh2  = 'COPI' cltf ;
  55. clro2  = 'COPI' cltf ;
  56. clrn2  = 'COPI' cltf ;
  57. clrhe  = 'COPI' cltf ;
  58. clrco  = 'COPI' cltf ;
  59. clrco2 = 'COPI' cltf ;
  60. clrair = 'COPI' cltf ;
  61. *
  62. * Initialisation des CL vectorielles à l'entrée et à la sortie des PARs
  63. clupar = 'COPI' cltf ;
  64. *
  65. *
  66. *- Boucle principale sur les PARs
  67. *
  68. *
  69. NBREC = 'DIME' rxt . 'RECOMB' ;
  70. 'REPE' BCLREC NBREC ;
  71.    RECI = rxt . 'RECOMB' . &BCLREC ;
  72.    RTRI = tic . 'RECOMB' . &BCLREC ;
  73.    TYPR     = RECI . 'TYPR' ;
  74.    NSECTION = RECI . 'NSECTION' ;
  75. *
  76. ** Modèle, surface, facteur géométrique et normale à l'entrée et à la sortie du PAR
  77.    NAME = 'CHAI' '$REC'  'ent' &BCLREC ; $ENTREE = GEO . NAME ;
  78.    NAME = 'CHAI' 'SREC'  'ent' &BCLREC ; SENTREE = GEO . NAME ;
  79.    NAME = 'CHAI' 'fgREC' 'ent' &BCLREC ; FG_ENT  = GEO . NAME ;
  80.    NAME = 'CHAI' 'DREC'  'ent' &BCLREC ; dire    = GEO . NAME ;
  81.    NAME = 'CHAI' 'REC'  'entI' &BCLREC ; ENTREI  = GEO . NAME ;
  82. *
  83.    NAME = 'CHAI' '$REC'  'sor' &BCLREC ; $SORTIE = GEO . NAME ;
  84.    NAME = 'CHAI' 'SREC'  'sor' &BCLREC ; SSORTIE = GEO . NAME ;
  85.    NAME = 'CHAI' 'fgREC' 'sor' &BCLREC ; FG_SOR  = GEO . NAME ;
  86.    NAME = 'CHAI' 'DREC'  'sor' &BCLREC ; dirs    = GEO . NAME ;
  87.    NAME = 'CHAI' 'REC'  'sorI' &BCLREC ; SORTII  = GEO . NAME ;
  88. *
  89. ** Surfaces des mailles duales
  90.    Sduae = 'DOMA' $ENTREE 'XXDIAGSI' ;
  91.    Sduas = 'DOMA' $SORTIE 'XXDIAGSI' ;
  92. *
  93. ** Fractions massiques moyennes à l'entrée du PAR
  94. ** (avec décomposition de l'air en N2 et O2)
  95.    Rgair = tbt . 'Rgair' ;
  96.    Rgo2  = tbt . 'Rgo2' ;
  97.    Rgn2  = tbt . 'Rgn2' ;
  98.    YO2e  = 0.21 '*' Rgair '/' Rgo2 '*' tic . 'YAIR' ;
  99.    YN2e  = 0.79 '*' Rgair '/' Rgn2 '*' tic . 'YAIR' ;
  100.    'SI' TBT . 'TO2' ;
  101.       YO2e  = 'KCHT' $ENTREE 'SCAL' 'SOMMET' (tic . 'YO2' '+' YO2e) ;
  102.       YO2e  = ('SOMT' (YO2e  '*' Sduae)) '/' SENTREE ;
  103.    'SINO' ;
  104.       YO2e  = 'KCHT' $ENTREE 'SCAL' 'SOMMET' YO2e ;
  105.       YO2e  = ('SOMT' (YO2e  '*' Sduae)) '/' SENTREE ;
  106.    'FINSI' ;
  107.    'SI' TBT . 'TN2' ;
  108.       YN2e  = 'KCHT' $ENTREE 'SCAL' 'SOMMET' (tic . 'YN2' '+' YN2e) ;
  109.       YN2e  = ('SOMT' (YN2e  '*' Sduae)) '/' SENTREE ;
  110.    'SINO' ;
  111.       YN2e  = 'KCHT' $ENTREE 'SCAL' 'SOMMET' YN2e ;
  112.       YN2e  = ('SOMT' (YN2e  '*' Sduae)) '/' SENTREE ;
  113.    'FINSI' ;
  114. *
  115.    'SI' TBT . 'TH2' ;
  116.       YH2e  = 'KCHT' $ENTREE 'SCAL' 'SOMMET' tic . 'YH2' ;
  117.       YH2e  = ('SOMT' (YH2e  '*' Sduae)) '/' SENTREE ;
  118.    'SINO' ;
  119.       YH2e  = 0. ;
  120.    'FINSI' ;
  121.    'SI' TBT . 'VAPEUR' ;
  122.       YH2Oe = 'KCHT' $ENTREE 'SCAL' 'SOMMET' tic . 'YVAP' ;
  123.       YH2Oe = ('SOMT' (YH2Oe '*' Sduae)) '/' SENTREE ;
  124.    'SINO' ;
  125.       YH2Oe = 0. ;
  126.    'FINSI' ;
  127.    'SI' TBT . 'THE' ;
  128.       YHEe  = 'KCHT' $ENTREE 'SCAL' 'SOMMET' tic . 'YHE' ;
  129.       YHEe  = ('SOMT' (YHEe  '*' Sduae)) '/' SENTREE ;
  130.    'SINO' ;
  131.       YHEe  = 0. ;
  132.    'FINSI' ;
  133.    'SI' TBT . 'TCO' ;
  134.       YCOe  = 'KCHT' $ENTREE 'SCAL' 'SOMMET' tic . 'YCO' ;
  135.       YCOe  = ('SOMT' (YCOe  '*' Sduae)) '/' SENTREE ;
  136.    'SINO' ;
  137.       YCOe  = 0. ;
  138.    'FINSI' ;
  139.    'SI' TBT . 'TCO2' ;
  140.       YCO2e = 'KCHT' $ENTREE 'SCAL' 'SOMMET' tic . 'YCO2' ;
  141.       YCO2e = ('SOMT' (YCO2e '*' Sduae)) '/' SENTREE ;
  142.    'SINO' ;
  143.       YCO2e = 0. ;
  144.    'FINSI' ;
  145. *
  146. ** Tgaz, Tplaque à l'entrée du PAR
  147.    Te = 'KCHT' $ENTREE 'SCAL' 'SOMMET' tic . 'TF' ;   
  148.    TIN = ('SOMT' (Te '*' Sduae)) '/' SENTREE ;
  149.    TINK = TIN '+' 273.15 ;
  150.    TP = 'EXTR' rtri . 'Tpla' ('DIME' rtri . 'Tpla') ;
  151.    TPK = TP '+' 273.15 ;
  152. *
  153. *----------------------------
  154. * Recombineur de type SIEMENS
  155. *----------------------------
  156. *
  157.    'SI' ('EGA' typr 'SIEMENS') ;
  158.      AA1    = RECI . 'A' ;
  159.      BB1    = RECI . 'B' ;
  160.      XH2ONI = RECI . 'XH2ON' ;
  161.      XH2OFI = RECI . 'XH2OFF' ;
  162.      BETAI  = RECI . 'BETA' ;
  163.      YH2s YO2s YN2s YH2Os YHEs YCO2s YCOs TsK TPK DEB EFF XH2ONF
  164.      = 'FLAM' 'RECOMBIN'
  165.      YH2e YO2e YN2e YH2Oe YHEe YCO2e YCOe
  166.      TINK TPK Pe Dt 0 typr AA1 BB1 XH2ONI XH2OFI BETAI ;
  167.      RECI . 'XH2ON' = XH2ONF ;
  168. *
  169. *---------------------------
  170. * Recombineur de type HEATER
  171. *---------------------------
  172. *
  173.    'SINO' ;
  174.      PUI0 = 'IPOL' Tps RECI . 'PUISSANCE' ;
  175.      PUI1 = 'IPOL' (Tps '-' Dt) RECI . 'PUISSANCE' ;
  176.      A1 = PUI0 '-' PUI1 '/' DT ;
  177.      FUITE = RECI . 'QHOUSING' ;
  178.      HRAPP = RECI . 'HRAPP' ;
  179.      CRAPP = RECI . 'CRAPP' ;
  180.      YH2s YO2s YN2s YH2Os YHEs YCO2s YCOs TsK TPK DEB
  181.      = 'FLAM' 'RECOMBIN'
  182.      YH2e YO2e YN2e YH2Oe YHEe YCO2e YCOe
  183.      TINK TPK Pe Dt 0 typr A1 fuite hrapp crapp ;
  184.      EFF = 0.D0 ;
  185.    'FINSI' ;
  186. *
  187.    TS = TSK '-' 273.15 ;
  188.    TP = TPK '-' 273.15 ;
  189.    NDEB = DEB '*' NSECTION ; 
  190. *
  191. * Vitesses débitantes à l'entrée et à la sortie du PAR (en m/s)
  192. * avec prise en compte du facteur géométrique et du nombre de sections
  193.    RG_IN  = (TBT . 'Rgvap' '*' YH2Oe)
  194.         '+' (TBT . 'Rgh2'  '*' YH2e) 
  195.         '+' (TBT . 'Rgo2'  '*' YO2e) 
  196.         '+' (TBT . 'Rgn2'  '*' YN2e) 
  197.         '+' (TBT . 'Rghe'  '*' YHEe) 
  198.         '+' (TBT . 'Rgco'  '*' YCOe) 
  199.         '+' (TBT . 'Rgco2' '*' YCO2e) ;
  200.    RG_IN  = ('SOMT' (RG_IN '*' Sduae)) '/' SENTREE ;
  201.    RO_IN  = Pe '/' (RG_IN '*' TINK) ;
  202.    UN_IN  = ndeb '*' FG_ENT '/' (SENTREE '*' RO_IN) ;
  203. *
  204.    RG_OUT = (TBT . 'Rgvap' '*' YH2Os) 
  205.         '+' (TBT . 'Rgh2'  '*' YH2s) 
  206.         '+' (TBT . 'Rgo2'  '*' YO2s) 
  207.         '+' (TBT . 'Rgn2'  '*' YN2s) 
  208.         '+' (TBT . 'Rghe'  '*' YHEs) 
  209.         '+' (TBT . 'Rgco'  '*' YCOs) 
  210.         '+' (TBT . 'Rgco2' '*' YCO2s) ; 
  211.    RG_OUT = ('SOMT' (RG_OUT '*' Sduas)) '/' SSORTIE ;
  212.    RO_OUT = Pe '/' (RG_OUT '*' TSK) ;
  213.    UN_OUT = ndeb '*' FG_SOR '/' (SSORTIE '*' RO_OUT) ;
  214. *
  215. ** Densité des constituants à l'entrée et à la sortie du PAR (en kg/m3)
  216.    RH2E  = RO_IN '*' YH2e  ;
  217.    RH2OE = RO_IN '*' YH2Oe ;
  218.    RO2E  = RO_IN '*' YO2e  ;
  219.    RN2E  = RO_IN '*' YN2e  ;
  220.    RHEE  = RO_IN '*' YHEe  ;
  221.    RCO2E = RO_IN '*' YCO2e ;
  222.    RCOE  = RO_IN '*' YCOe  ;
  223. *
  224.    RH2S  = RO_OUT '*' YH2s  ;
  225.    RH2OS = RO_OUT '*' YH2Os ;
  226.    RO2S  = RO_OUT '*' YO2s  ;
  227.    RN2S  = RO_OUT '*' YN2s  ;
  228.    RHES  = RO_OUT '*' YHEs  ;
  229.    RCO2S = RO_OUT '*' YCO2s ;
  230.    RCOS  = RO_OUT '*' YCOs  ;
  231. *
  232. ** Enthalpie massique à l'entrée et à la sortie du PAR (en J/kg)
  233.    Cph2_e Cphe_e Cpo2_e Cpn2_e Cpco2_e Cpco_e Cpair_e
  234.    = CALCP TIN ;
  235.    Cph2_s Cphe_s Cpo2_s Cpn2_s Cpco2_s Cpco_s Cpair_s
  236.    = CALCP TS ;
  237.    Cpvap = TBT . 'Cpvap' ;
  238.    CP_IN  = Cpvap '*' YH2Oe  
  239.      '+' (Cph2_e  '*' YH2e)  
  240.      '+' (Cpo2_e  '*' YO2e)  
  241.      '+' (Cpn2_e  '*' YN2e)  
  242.      '+' (Cphe_e  '*' YHEe)  
  243.      '+' (Cpco2_e '*' YCO2e) 
  244.      '+' (Cpco_e  '*' YCOe)  ;
  245.    H_IN  = CP_IN '*' TINK ;
  246. *
  247.    CP_OUT = Cpvap '*' YH2Os
  248.      '+' (Cph2_s  '*' YH2s)  
  249.      '+' (Cpo2_s  '*' YO2s)  
  250.      '+' (Cpn2_s  '*' YN2s)  
  251.      '+' (Cphe_s  '*' YHEs)  
  252.      '+' (Cpco2_s '*' YCO2s) 
  253.      '+' (Cpco_s  '*' YCOs)  ;
  254.    H_OUT = CP_OUT '*' TSK ;
  255. *
  256. ** Contribution du PAR aux bilans massiques 0D (Q en kg/s)
  257.    QIN_H2   = ndeb '*' YH2e  '+' QIN_H2   ;
  258.    QIN_H2O  = ndeb '*' YH2Oe '+' QIN_H2O  ;
  259.    QIN_O2   = ndeb '*' YO2e  '+' QIN_O2   ;
  260.    QIN_N2   = ndeb '*' YN2e  '+' QIN_N2   ;
  261.    QIN_HE   = ndeb '*' YHEe  '+' QIN_HE   ;
  262.    QIN_CO2  = ndeb '*' YCO2e '+' QIN_CO2  ;
  263.    QIN_CO   = ndeb '*' YCOe  '+' QIN_CO   ;
  264. *
  265.    QOUT_H2  = ndeb '*' YH2S  '+' QOUT_H2  ;
  266.    QOUT_H2O = ndeb '*' YH2OS '+' QOUT_H2O ;
  267.    QOUT_O2  = ndeb '*' YO2S  '+' QOUT_O2  ;
  268.    QOUT_N2  = ndeb '*' YN2S  '+' QOUT_N2  ;
  269.    QOUT_HE  = ndeb '*' YHES  '+' QOUT_HE  ;
  270.    QOUT_CO2 = ndeb '*' YCO2S '+' QOUT_CO2 ;
  271.    QOUT_CO  = ndeb '*' YCOS  '+' QOUT_CO  ;
  272. *
  273. ** Contribution des PARs au bilan d'energie 0D (drhpar en W)
  274.    drhpar = H_OUT '-' H_IN '*' ndeb '+' drhpar ;
  275. *
  276. ** Conditions aux limites en vitesse (en m/s)
  277.    toto1 = 'KCHT' $ENTREE 'VECT' 'SOMMET' 'COMP' LMOT1 (dire * UN_IN) ;
  278.    tutu1 = 'KCHT' $SORTIE 'VECT' 'SOMMET' 'COMP' LMOT1 (dirs * UN_OUT) ;
  279.    'SI' tbt . 'FPAROI' ;
  280.      NAME = 'CHAI' 'REC'  'ent' &BCLREC ; ENTREE  = GEO . NAME ;
  281.      NAME = 'CHAI' 'REC'  'sor' &BCLREC ; SORTIE  = GEO . NAME ;
  282.      toto1 = 'REDU' toto1 ENTREE ;
  283.      tutu1 = 'REDU' tutu1 SORTIE ;
  284.    'SINO' ;
  285.      toto1 = 'REDU' toto1 ENTREI ;
  286.      tutu1 = 'REDU' tutu1 SORTII ;
  287.    'FINS' ;
  288.    toto1 = 'KCHT' $ENTREE 'VECT' 'SOMMET' 'COMP' LMOT1 toto1 ;
  289.    tutu1 = 'KCHT' $SORTIE 'VECT' 'SOMMET' 'COMP' LMOT1 tutu1 ;
  290.    clupar = clupar 'ET' toto1 'ET' tutu1 ;
  291. *
  292. ** Conditions aux limites en température (en C)
  293.    cltf   = cltf 'ET'
  294.    ('MANU' 'CHPO' SORTII 1 'TF' TS 'NATURE' 'DISCRET') ;
  295. *
  296. ** Conditions aux limites en densités partielles (en kg/m3)
  297.    clrvp  = clrvp  'ET'
  298.    ('MANU' 'CHPO' SORTII 1 'RVP'  RH2Os 'NATURE' 'DISCRET') ;
  299.    clrh2  = clrh2  'ET'
  300.    ('MANU' 'CHPO' SORTII 1 'RH2'  RH2s  'NATURE' 'DISCRET') ;
  301.    clro2  = clro2  'ET'   
  302.    ('MANU' 'CHPO' SORTII 1 'RO2'  RO2s  'NATURE' 'DISCRET') ;
  303.    clrn2  = clrn2  'ET'  
  304.    ('MANU' 'CHPO' SORTII 1 'RN2'  RN2s  'NATURE' 'DISCRET') ;
  305.    clrhe  = clrhe  'ET'  
  306.    ('MANU' 'CHPO' SORTII 1 'RHE'  RHEs  'NATURE' 'DISCRET') ;
  307.    clrco  = clrco  'ET'  
  308.    ('MANU' 'CHPO' SORTII 1 'RCO'  RCOs  'NATURE' 'DISCRET') ;
  309.    clrco2 = clrco2 'ET'
  310.    ('MANU' 'CHPO' SORTII 1 'RCO2' RCO2s 'NATURE' 'DISCRET') ;
  311.    clrair = clrair 'ET'
  312.    ('MANU' 'CHPO' SORTII 1 'RAIR' 0.    'NATURE' 'DISCRET') ;
  313. *
  314. * Sauvegarde dans TIC des data pour post-traitement local
  315. * (T en C, R en kg/m3, DEB en kg/s, Efficiency et Y sans unité, U en m/s)
  316.    rtri . 'Tpla'  = rtri . 'Tpla'  'ET' ('PROG' TP)  ;
  317.    rtri . 'Tin'   = rtri . 'Tin'   'ET' ('PROG' TIN) ;
  318.    rtri . 'Tout'  = rtri . 'Tout'  'ET' ('PROG' TS)  ;
  319. *
  320.    rtri . 'RH2E'  = rtri . 'RH2E'  'ET' ('PROG' RH2E ) ;
  321.    rtri . 'RH2OE' = rtri . 'RH2OE' 'ET' ('PROG' RH2OE) ;
  322.    rtri . 'RO2E'  = rtri . 'RO2E'  'ET' ('PROG' RO2E ) ;
  323.    rtri . 'RN2E'  = rtri . 'RN2E'  'ET' ('PROG' RN2E ) ;
  324.    rtri . 'RHEE'  = rtri . 'RHEE'  'ET' ('PROG' RHEE ) ;
  325.    rtri . 'RCO2E' = rtri . 'RCO2E' 'ET' ('PROG' RCO2E) ;
  326.    rtri . 'RCOE'  = rtri . 'RCOE'  'ET' ('PROG' RCOE ) ;
  327. *
  328.    rtri . 'RH2S'  = rtri . 'RH2S'  'ET' ('PROG' RH2S ) ;
  329.    rtri . 'RH2OS' = rtri . 'RH2OS' 'ET' ('PROG' RH2OS) ;
  330.    rtri . 'RO2S'  = rtri . 'RO2S'  'ET' ('PROG' RO2S ) ;
  331.    rtri . 'RN2S'  = rtri . 'RN2S'  'ET' ('PROG' RN2S ) ;
  332.    rtri . 'RHES'  = rtri . 'RHES'  'ET' ('PROG' RHES ) ;
  333.    rtri . 'RCO2S' = rtri . 'RCO2S' 'ET' ('PROG' RCO2S) ;
  334.    rtri . 'RCOS'  = rtri . 'RCOS'  'ET' ('PROG' RCOS ) ;
  335. *   
  336.    rtri . 'DEB'   = rtri . 'DEB'   'ET' ('PROG' ndeb)  ;
  337.    rtri . 'EFF'   = rtri . 'EFF'   'ET' ('PROG' EFF)   ;
  338. *
  339.    rtri . 'YH2E'  = rtri . 'YH2E'  'ET' ('PROG' YH2e)  ;
  340.    rtri . 'YH2OE' = rtri . 'YH2OE' 'ET' ('PROG' YH2Oe) ;
  341.    rtri . 'YO2E'  = rtri . 'YO2E'  'ET' ('PROG' YO2e)  ;
  342.    rtri . 'YN2E'  = rtri . 'YN2E'  'ET' ('PROG' YN2e)  ;
  343.    rtri . 'YHEE'  = rtri . 'YHEE'  'ET' ('PROG' YHEe)  ;
  344.    rtri . 'YCO2E' = rtri . 'YCO2E' 'ET' ('PROG' YCO2e) ;   
  345.    rtri . 'YCOE'  = rtri . 'YCOE'  'ET' ('PROG' YCOe)  ;
  346. *   
  347.    rtri . 'YH2S'  = rtri . 'YH2S'  'ET' ('PROG' YH2s)  ;
  348.    rtri . 'YH2OS' = rtri . 'YH2OS' 'ET' ('PROG' YH2Os) ;
  349.    rtri . 'YO2S'  = rtri . 'YO2S'  'ET' ('PROG' YO2s)  ;
  350.    rtri . 'YN2S'  = rtri . 'YN2S'  'ET' ('PROG' YN2s)  ;
  351.    rtri . 'YHES'  = rtri . 'YHES'  'ET' ('PROG' YHEs)  ;
  352.    rtri . 'YCO2S' = rtri . 'YCO2S' 'ET' ('PROG' YCO2s) ;
  353.    rtri . 'YCOS'  = rtri . 'YCOS'  'ET' ('PROG' YCOs)  ;
  354. *
  355.    rtri . 'Uin'   = rtri . 'Uin'   'ET' ('PROG' UN_IN)  ;
  356.    rtri . 'Uout'  = rtri . 'Uout'  'ET' ('PROG' UN_OUT) ;
  357. *
  358.    rtri . 'Hin'   = rtri . 'Hin'   'ET' ('PROG' H_IN)  ;
  359.    rtri . 'Hout'  = rtri . 'Hout'  'ET' ('PROG' H_OUT) ;
  360. 'FIN' BCLREC ;
  361. *
  362. ** Sauvegarde des conditions aux limites
  363. tic . 'RECOMB' . 'UN'   = clupar ;
  364. tic . 'RECOMB' . 'TF'   = cltf   ;
  365. tic . 'RECOMB' . 'RVP'  = clrvp  ;
  366. tic . 'RECOMB' . 'RH2'  = clrh2  ;
  367. tic . 'RECOMB' . 'RO2'  = clro2  ;
  368. tic . 'RECOMB' . 'RN2'  = clrn2  ;
  369. tic . 'RECOMB' . 'RHE'  = clrhe  ;
  370. tic . 'RECOMB' . 'RCO'  = clrco  ;
  371. tic . 'RECOMB' . 'RCO2' = clrco2 ;
  372. tic . 'RECOMB' . 'RAIR' = clrair ;
  373. *
  374. ** Contribution des PARs aux bilans de masse 0D (Qo_beta en kg/s)
  375. QO_H2  = QOUT_H2  '-' QIN_H2  ;
  376. QO_O2  = QOUT_O2  '-' QIN_O2  ;
  377. QO_H2O = QOUT_H2O '-' QIN_H2O ;
  378. QO_N2  = QOUT_N2  '-' QIN_N2  ;
  379. *
  380. * Affichage pour vérification
  381. QO_OUT = QOUT_H2  '+' QOUT_O2 '+' QOUT_H2O '+' QOUT_N2 '+' QOUT_HE '+'
  382.          QOUT_CO2 '+' QOUT_CO ; 
  383. QO_IN  = QIN_H2   '+' QIN_O2  '+' QIN_H2O  '+' QIN_N2  '+' QIN_HE  '+'
  384.          QIN_CO2  '+' QIN_CO  ;
  385. 'MESS' 'DQ / Qin / Qout (kg/s)' ;
  386. 'MESS' '-------------------------------------------------------------' ;
  387. 'MESS' 'Balance :' ' ' (QO_OUT '-' QO_IN) ' ' QO_OUT ' ' QO_IN ;
  388. 'MESS' '-------------------------------------------------------------' ;
  389. 'MESS' 'H2      :' ' ' QO_H2  ' ' QIN_H2  ' ' QOUT_H2  ;
  390. 'MESS' 'H2O     :' ' ' QO_H2O ' ' QIN_H2O ' ' QOUT_H2O ;
  391. 'MESS' 'O2      :' ' ' QO_O2  ' ' QIN_O2  ' ' QOUT_O2  ;
  392. 'MESS' 'N2      :' ' ' QO_N2  ' ' QIN_N2  ' ' QOUT_N2  ;
  393. 'MESS' ' ' ;
  394. 'MESS' 'HE      :' ' ' (QIN_HE  '-' QOUT_HE)  ' ' QIN_HE  ' ' QOUT_HE  ;
  395. 'MESS' 'CO2     :' ' ' (QIN_CO2 '-' QOUT_CO2) ' ' QIN_CO2 ' ' QOUT_CO2 ;
  396. 'MESS' 'CO      :' ' ' (QIN_CO  '-' QOUT_CO)  ' ' QIN_CO  ' ' QOUT_CO  ;
  397. 'MESS' '-------------------------------------------------------------' ;
  398. *
  399. 'RESP' QO_H2 QO_O2 QO_H2O QO_N2 drhpar ;
  400. 'FINP' ;
  401.  

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