Télécharger la plaquette Cast3M


fr en

Télécharger broche.procedur

Retour à la liste

Numérotation des lignes : 

  1. * BROCHE    PROCEDUR  SERRE     21/06/28    21:15:01     11045          
  2. 'DEBPROC' broche KAS*'MOT' nab*'ENTIER' rxt*'TABLE';
  3. *mess '=================== DEBUT BROCHE ===================';
  4.  
  5.  GEO = rxt.'GEO'                                        ;
  6.  nub = nab                                              ;
  7.  ibr=index (rxt.'Sorties')                              ;
  8.  Briches =rxt.'Sorties'.(ibr.nub)                       ;
  9.  rxtsce  =Briches.'scenario'                            ;
  10.  GEO     =Briches                                       ;
  11.  CTinj   =chai 'Tinj'                                   ;
  12.  CKsi    =chai 'Ksi'                                    ;
  13.  CQlj    =chai 'Qlj'                                    ;
  14.  CHlj    =chai 'Hlj'                                    ;
  15.  
  16.   TTOUT=FAUX                                             ;
  17.   TTINJ=FAUX                                             ;
  18.   THINJ=FAUX                                             ;
  19.   TTINJ=exist (rxtsce ) 'tinj'                           ;
  20.   TTPEX=exist (rxtsce ) 'Pext'                           ;
  21.   TTQEX=exist (rxtsce ) 'Qext'                           ;
  22.   TTOUT=TTPEX ou TTQEX                                   ;
  23.  
  24.   Briches.'TTOUT'=TTOUT                                  ;
  25.  
  26. *mess ' TTOUT TTPEX TTQEX ******************************';
  27. *list ttout; list ttpex ; list ttqex;
  28. *mess ' KAS=================' KAS;
  29. *mess ' TTOUT TTPEX TTQEX ******************************';
  30.  
  31.   Si ((NON TTOUT) et (NON TTINJ))                        ;
  32.  Mess '*********************************************************' ;
  33.  Mess ' ERREUR   ERREUR   ERREUR   ERREUR   ERREUR   ERREUR     ' ;
  34.  Mess '                                                         ' ;
  35.  Mess 'il faut imposer au moins tinj ou pression !              ' ;
  36.  Mess '*********************************************************' ;
  37.  ERREUR 21                                                        ;
  38.   Finsi                                                           ;
  39.  
  40.   Si (TTPEX et TTQEX)                                             ;
  41.  Mess '*********************************************************' ;
  42.  Mess ' ERREUR   ERREUR   ERREUR   ERREUR   ERREUR   ERREUR     ' ;
  43.  Mess '                                                         ' ;
  44.  Mess 'il faut choisir entre Clapet de décharge et extraction   ' ;
  45.  Mess '*********************************************************' ;
  46.  ERREUR 21                                                    ;
  47.   Finsi                                                       ;
  48.  
  49. * Modifications ANTS du 06/01/2009
  50. * Evite d'avoir un debit d'injection nul
  51.  
  52. 'SI' (('NON' ('EXIS' rxtsce Qmodif)) et (TTINJ ou THINJ));
  53. 'MESS' '=============================================================' ;
  54. 'MESS' '=                       WARNING !!                          =' ;
  55. 'MESS' '=    LES DEBITS D INJECTION NULS ONT ETE IMPOSES A 1E-10    =' ;
  56. 'MESS' '=============================================================' ;
  57.   lstqxxx = 'MOTS' 'qeau' 'qliq' 'qair' 'qhe'
  58.                    'qh2'  'qo2'  'qn2'  'qco' 'qco2'    ;
  59.   lmots1 = 'MOTS'                                       ;
  60.   'REPE' bqxxx ('DIME' lstqxxx)                         ;
  61.     lmots = 'EXTR' lstqxxx &bqxxx                       ;
  62.     'SI' ('EXIS' rxtsce lmots)                          ;
  63.       lmots1 = lmots1 'ET' ('MOTS' lmots)               ;
  64.     'FINSI'                                             ;
  65.   'FIN' bqxxx                                           ;
  66.   mots1 = 'EXTR' lmots1 1                               ;
  67.   'REPE' bvxxx ('DIME' rxtsce.mots1)                    ;
  68.     indic = FAUX                                        ;
  69.     'REPETER' bq1xxx ('DIME' lmots1)                    ;
  70.       motsxx = 'EXTR' lmots1 &bq1xxx                    ;
  71.       valxxx = 'EXTR' rxtsce.motsxx &bvxxx              ;
  72.       'SI' ('>' valxxx 1.E-10)                          ;
  73.         indic = VRAI                                    ;
  74.         'QUIT' bq1xxx                                   ;
  75.       'FINSI'                                           ;
  76.     'FIN' bq1xxx                                        ;
  77.     'SI' ('NON' indic)                                  ;
  78.       'REMP' rxtsce.mots1 &bvxxx 1E-10                  ;
  79.     'FINSI'                                             ;
  80.   'FIN' bvxxx                                           ;
  81.   rxtsce.Qmodif = VRAI                                  ;
  82. 'FINSI' ;
  83. * Fin modifications ANTS du 06/01/2009
  84.  
  85. ************************************************************************
  86. ********* DEBUT KAS1 ***************************************************
  87. ************************************************************************
  88.  
  89. 'SI  KAS1' (EGA KAS 'KAS1');
  90.  
  91. * Initialisation des debits
  92.  aex   =1. ;
  93.  Ksi   =1. ;
  94.  qliq  =0. ; hliq  =0. ; Mliqj = 0.0 ;
  95.  qeau  =0. ; hvap  =0. ; Yvapj = 0.0 ;
  96.  qair  =0. ; hair  =0. ; Yairj = 0.0 ;
  97.  qhe   =0. ; hhe   =0. ; Yhej  = 0.0 ;
  98.  qh2   =0. ; hh2   =0. ; Yh2j  = 0.0 ;
  99.  qo2   =0. ; ho2   =0. ; Yo2j  = 0.0 ;
  100.  qn2   =0. ; hn2   =0. ; Yn2j  = 0.0 ;
  101.  qco   =0. ; hco   =0. ; Ycoj  = 0.0 ;
  102.  qco2  =0. ; hco2  =0. ; Yco2j = 0.0 ;
  103.  
  104. * Les donnees physiques
  105.  Cpvap  = TBT.'Cpvap' ;
  106.  Lv     = TBT.'Lv'    ;
  107.  Rabs   = TBT.'Rabs'  ;
  108.  Rgh2   = TBT.'Rgh2'  ; Mkh2   = TBT.'Mkh2'              ;
  109.  Rghe   = TBT.'Rghe'  ; Mkhe   = TBT.'Mkhe'              ;
  110.  Rgo2   = TBT.'Rgo2'  ; Mko2   = TBT.'Mko2'              ;
  111.  Rgn2   = TBT.'Rgn2'  ; Mkn2   = TBT.'Mkn2'              ;
  112.  Rgco2  = TBT.'Rgco2' ; Mkco2  = TBT.'Mkco2'             ;
  113.  Rgco   = TBT.'Rgco'  ; Mkco   = TBT.'Mkco'              ;
  114.  Rgvap  = TBT.'Rgvap' ; Mkvap  = TBT.'Mkvap'             ;
  115.  Rgair  = TBT.'Rgair' ; Mkair  = TBT.'Mkair'             ;
  116.  
  117.  
  118.   t        = rxtsce .'t'                                 ;
  119.  
  120. *------------------ TTOUT ----------------------------------------------
  121. *-- Formulation en Extraction ou clapet de décharge --------------------
  122. *------------------ TTOUT ----------------------------------------------
  123.   'SI   TTOUT KAS1' TTOUT                                     ;
  124.  
  125.   $breche = Briches.'$sreche';
  126.   brechei = Briches.'srechei';
  127.   brechec = Briches.'srechec';
  128.   breche  = Briches.'sreche' ;
  129.  
  130.   Rout=redu Breche tic.'RHO';
  131.   diru    =     Briches.'diru'            ;
  132.  
  133.  Si(EGA (VALE DIME) 2);
  134.  gjc   = 'MANU' 'CHPO' brechec 2 'UX' 0.0 'UY' 0.0
  135.                                  'NATURE' 'DISCRET'          ;
  136.  Sinon ;
  137.  gjc   = 'MANU' 'CHPO' brechec 3 'UX' 0.0 'UY' 0.0 'UZ' 0.
  138.                                  'NATURE' 'DISCRET'          ;
  139.  Finsi;
  140.  
  141. *==============================================================
  142. * C'est la qu'on definit l'extraction
  143. *
  144. ******* OUT INI deb ******************************************
  145. Si (NON (exist rxt.'TIC' 'OUT'));
  146. mess ' on definit l extraction OUT ';
  147. rxt.'TIC'.'OUT'=table;
  148. rxt.'TIC'.'OUT'.'Tps' =prog;
  149. rxt.'TIC'.'OUT'.'Tfm' =prog;
  150. rxt.'TIC'.'OUT'.'PT'  =prog;
  151. rxt.'TIC'.'OUT'.'Rhomn'=prog;
  152. rxt.'TIC'.'OUT'.'Uout' =prog;
  153. rxt.'TIC'.'OUT'.'Uoutn'=prog;
  154. rxt.'TIC'.'OUT'.'Qout'=prog;
  155. rxt.'TIC'.'OUT'.'Rdp' =prog;
  156. rxt.'TIC'.'OUT'.'Rsbr'=prog;
  157. Finsi;
  158. ******* OUT INI fin ******************************************
  159.  
  160.      Pt  =extr tic.'PT'  ndl                                 ;
  161.      Tfm =extr tic.'Tfm' ndl                                 ;
  162.      Rhomn=extr tic.'Rhomn' ndl                              ;
  163.  
  164.   Sout      = rxtsce.'Section'                               ;
  165.   Sbreche   = mesu breche                                    ;
  166.   Rsbr      = Sout/Sbreche                                   ;
  167.  
  168.  Si TTQEX                                                    ;
  170. mess ' +++++++++++++++++++++++++++++++++++++++++++++++++++++'
  171.      '+++++++++++++++++++++++++++++++'                       ;
  172. mess ' + Extraction d un débit imposé ++++++++++++++++++++++'
  173.      '+++++++++++++++++++++++++++++++'                       ;
  174.   Rdp       = 1.                                             ;
  175.   Pext      = 0.                                             ;
  176.   Kext      = 0.                                             ;
  177.   Bext      = 0.                                             ;
  178.  
  179.   Si (Sout < 1.e-10)                                         ;
  180.   Qout      = 0.                                             ;
  181.   Qvrai     = 0.                                             ;
  182.   Qout      = 0.                                             ;
  183.   Uout      = 0.                                             ;
  184.   uus = 'KCHT' $breche 'VECT' 'SOMMET' (0.*diru)  gjc        ;
  185.  mess 'ATTENTION: la section est < 1.e-10, on n extrait rien';
  186. mess ' +++++++++++++++++++++++++++++++++++++++++++++++++++++'
  187.      '+++++++++++++++++++++++++++++++'                       ;
  188.   Sinon                                                      ;
  189.   Qout      = 'IPOL' tps t (rxtsce .'Qext')                  ;
  190.   Uout      = Qout/Rhomn/Sout*Rsbr                           ;
  191.  
  192.   uus = 'KCHT' $breche 'VECT' 'SOMMET' (Rhomn*uout*diru)  gjc;
  193.   Qvrai     = DBIT $breche uus                               ;
  194.   alfa=qvrai/qout                                            ;
  195.   uoutn=uout/alfa                                            ;
  196.   uus = 'KCHT' $breche 'VECT' 'SOMMET' (Rhomn*uoutn*diru) gjc;
  197.   Briches.'gj' = uus                                         ;
  198.   Qvrai     = DBIT $breche uus                               ;
  199.  mess ' VERIFICATION (Qimpo) Qvrai=' Qvrai ' Qout=' Qout;
  200.  
  201. * Tinj      = Tout                                           ;
  202. *     rxtsce .'Tj'   = tinj                       ;
  203. mess ' Penceinte =' Pt ' Qout imposé=' Qout ' Uout=' Uoutn   ;
  204. mess ' Rhomn=' Rhomn '  Tfm=' Tfm ' Sout=' Sout              ;
  205.  
  206. mess ' +++++++++++++++++++++++++++++++++++++++++++++++++++++'
  207.      '+++++++++++++++++++++++++++++++'                       ;
  208.   Finsi                                                      ;
  209.  Finsi                                                       ;
  210.  
  211. mess 'µµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµ' ;
  212. mess 'µµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµ' ;
  213.  
  214.  
  215. *µµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµ
  216. *µµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµµ
  217.  
  218.  
  219. Si TTPEX                                                     ;
  221. mess ' +++++++++++++++++++++++++++++++++++++++++++++++++++++'
  222.      '+++++++++++++++++++++++++++++++'                       ;
  223. mess ' + Modèle de clapet de décharge ++++++++++++++++++++++'
  224.      '+++++++++++++++++++++++++++++++'                       ;
  225.  
  226.   Pext      = 'IPOL' tps t (rxtsce .'Pext')                  ;
  227.   Kext      = 'IPOL' tps t (rxtsce .'Kext')                  ;
  228.   Bext      = 'IPOL' tps t (rxtsce .'Bext')                  ;
  229.   Uout      = ((Pt - Pext)/Kext)**(1./Bext)                  ;
  230.   Qout      = Rhomn * Uout * Sout                            ;
  231.   Rdp       = Pt/Pext                                        ;
  232.  
  233. ******* OUT **************************************************
  234. rxt.'TIC'.'OUT'.'Tps' =rxt.'TIC'.'OUT'.'Tps' et (prog tps)   ;
  235. rxt.'TIC'.'OUT'.'Tfm' =rxt.'TIC'.'OUT'.'Tfm' et (prog Tfm)   ;
  236. rxt.'TIC'.'OUT'.'PT'  =rxt.'TIC'.'OUT'.'PT'  et (prog PT)    ;
  237. rxt.'TIC'.'OUT'.'Rhomn'=rxt.'TIC'.'OUT'.'Rhomn'et(prog Rhomn);
  238. rxt.'TIC'.'OUT'.'Rdp' =rxt.'TIC'.'OUT'.'Rdp'  et (prog Rdp ) ;
  239. rxt.'TIC'.'OUT'.'Rsbr'=rxt.'TIC'.'OUT'.'Rsbr' et (prog Rsbr) ;
  240.  
  241. ******* OUT **************************************************
  242.  
  243.   Si (Sout < 1.e-10)                                         ;
  244.   Qout      = 0.                                             ;
  245.   Qvrai     = 0.                                             ;
  246.   Uout      = 0.                                             ;
  247.   Uoutn     = 0.                                             ;
  248.   uus = 'KCHT' $breche 'VECT' 'SOMMET' (0.*diru)  gjc        ;
  249.  mess                                                        ;
  250.  mess 'ATTENTION: la section est < 1.e-10, on n extrait rien';
  251.  mess                                                        ;
  252.   Sinon                                                      ;
  253.  
  254.      gamma=1.29;
  255.      Pcr    = Pt * ( ( (2./(gamma + 1.) )**(gamma/(gamma - 1.))));
  256.      Co   = (gamma * Pcr / rhom)**0.5                        ;
  257.      Qcr  = Rhomn * Sout * Co                                ;
  258.  
  259.   Si ((Pext < Pcr) et (Sout > 1.e-10))                       ;
  260.      Uout = Co * Rsbr                                        ;
  261.      Qout      = Qcr                                         ;
  262.  
  263.   uus = 'KCHT' $breche 'VECT' 'SOMMET' (Rhomn*uout*diru) gjc ;
  264.   Qvrai     = DBIT $breche uus                               ;
  265.   alfa=qvrai/qout                                            ;
  266.   uoutn=uout/alfa                                            ;
  267.   uus = 'KCHT' $breche 'VECT' 'SOMMET' (Rhomn*uoutn*diru) gjc;
  268.   Briches.'gj' = uus                                         ;
  269.   Qvrai     = DBIT $breche uus                               ;
  270.  mess ' VERIFICATION (Clapet) Qvrai=' Qvrai ' Qout=' Qout;
  271.  
  272.   mess ' Il y a un col sonique : Débit critique = ' Qcr ' Uout=' Uout ;
  273.   Sinon                                                      ;
  274.     Si (Uout > Co)                                           ;
  275.      Uout = Co * Rsbr                                        ;
  276.      Qout      = Qcr                                         ;
  277.  
  278.   uus = 'KCHT' $breche 'VECT' 'SOMMET' (Rhomn*uout*diru) gjc ;
  279.   Qvrai     = DBIT $breche uus                               ;
  280.   alfa=qvrai/qout                                            ;
  281.   uoutn=uout/alfa                                            ;
  282.   uus = 'KCHT' $breche 'VECT' 'SOMMET' (Rhomn*uoutn*diru) gjc;
  283.   Briches.'gj' = uus                                         ;
  284.   Qvrai     = DBIT $breche uus                               ;
  285.  mess ' VERIFICATION (Clapet) Qvrai=' Qvrai ' Qout=' Qout;
  286.  
  287.   mess ' Il y a qd même un col sonique : Débit critique = '
  288.    Qcr ' Uout=' Uout ;
  289.     Sinon                                                    ;
  290.      Qout = Uout * Rsbr                                      ;
  291.   uus = 'KCHT' $breche 'VECT' 'SOMMET' (Rhomn*uout*diru) gjc ;
  292.   Qvrai     = DBIT $breche uus                               ;
  293.   alfa=qvrai/qout                                            ;
  294.   uoutn=uout/alfa                                            ;
  295.   uus = 'KCHT' $breche 'VECT' 'SOMMET' (Rhomn*uoutn*diru) gjc;
  296.   Briches.'gj' = uus                                         ;
  297.   Qvrai     = DBIT $breche uus                               ;
  298.  mess ' VERIFICATION (Clapet) Qvrai=' Qvrai ' Qout=' Qout;
  299.  
  300.   mess ' Il n y a pas de col sonique : Débit de sortie ' Qout
  301.      ' Uout=' Uout                                           ;
  302.     Finsi                                                    ;
  303.   Finsi                                                      ;
  304.  
  305. rxt.'TIC'.'OUT'.'Uout' =rxt.'TIC'.'OUT'.'Uout'  et (prog Uout)  ;
  306. ********rxt.'TIC'.'OUT'.'Uoutn'=rxt.'TIC'.'OUT'.'Uoutn' et (prog Uoutn) ;
  307. rxt.'TIC'.'OUT'.'Qout' =rxt.'TIC'.'OUT'.'Qout'  et (prog Qout) ;
  308. mess ' Pext            Penceinte       Rapport de pression   '
  309. 'Pression critique ';
  310. mess  Pext  Pt '  '  Rdp '      ' Pcr;
  311.  
  312. mess ' +++++++++++++++++++++++++++++++++++++++++++++++++++++'
  313.      '+++++++++++++++++++++++++++++++'                       ;
  315. *
  316. *==============================================================
  317.  Finsi                                                       ;
  318. FINSI                                                        ;
  319.  
  320.  
  321.  mess ' VERIFICATION Qvrai=' Qvrai ' Qout=' Qout;
  322.   Tout = redu tic.'TF' breche ;
  323.  
  324.  Si(Qout > 0.)                              ;
  325.   Tout = ('DBIT' $breche (Tout*uus)) / Qout ;
  326.   Rout = ('DBIT' $breche (Rout*uus)) / Qout ;
  327.   Roj = Rout                                ;
  328.   Tinj= Tout                                ;
  329.   Qj  = Qout                                ;
  330.  mess ' Tout=' Tout ' Rout=' Rout ' Qj=' Qj ;
  331.  Sinon                                      ;
  332.   Rout=1.;
  333.   Tout=0.;
  334.   Roj = 1.                                  ;
  335.   Tinj= 0.                                  ;
  336.   Qj  = 0.                                  ;
  337.  Finsi                                      ;
  338.  
  339.   src   = 'DBIT' $breche uus 'IMPR'                    ;
  340.  
  341.   Si (src < 0.);
  343.   mess 'L Extraction est dans le mauvais sens ! ';
  345.   erre 2;
  346.   quitter broche;
  347.   Finsi;
  348.  
  349. * $menvf  = rxt.'GEO'.'$menvf'        ;
  350. * ung=vect uus 0.5 ux uy uz jaune;
  351. * trace ung (doma $menvf maillage);
  352.  
  353.   RHO =tic.'RHO' redu Breche                     ;
  354.   YAIR=tic.'YAIR' redu Breche                    ;
  355.   Y   = YAIR                                     ;
  356.   qair      = 'DBIT' $breche (uus*RHO*YAIR)      ;
  357.   nqj=1                                          ;
  358.  
  359.    'SI' TBT.'THE' ;
  360.   nqj=nqj + 1                                          ;
  361.   YHE =tic.'YHE' redu Breche                           ;
  362.   qhe = 'DBIT' $breche (uus*YHE*RHO) ; Y=Y+YHE ; Finsi ;
  363.  
  364.    'SI' TBT.'TH2' ;
  365.   nqj=nqj + 1                                          ;
  366.   YH2 =tic.'YH2' redu Breche                           ;
  367.   qh2 = 'DBIT' $breche (uus*YH2*RHO) ; Y=Y+YH2 ; Finsi ;
  368.  
  369.    'SI' TBT.'TO2' ;
  370.   nqj=nqj + 1                                          ;
  371.   YO2 =tic.'YO2' redu Breche                           ;
  372.   qo2 = 'DBIT' $breche (uus*YO2*RHO) ; Y=Y+YO2 ; Finsi ;
  373.  
  374.    'SI' TBT.'TN2' ;
  375.   nqj=nqj + 1                                          ;
  376.   YN2 =tic.'YN2' redu Breche                           ;
  377.   qn2 = 'DBIT' $breche (uus*YN2*RHO) ; Y=Y+YN2 ; Finsi ;
  378.  
  379.    'SI' TBT.'TCO';
  380.   nqj=nqj + 1                                          ;
  381.   YCO =tic.'YCO' redu Breche                           ;
  382.   qco = 'DBIT' $breche (uus*YCO*RHO) ; Y=Y+YCO ; Finsi ;
  383.  
  384.    'SI' TBT.'TCO2' ;
  385.   nqj=nqj + 1                                          ;
  386.   YCO2 =tic.'YCO2' redu Breche                         ;
  387.   qco2 = 'DBIT' $breche (uus*YCO2*RHO) ; Y=Y+YCO2;Finsi;
  388.  
  389.    'SI' TBT.'VAPEUR'                                    ;
  390.   nqj=nqj + 1                                           ;
  391.   RVP =tic.'RVP' redu Breche                            ;
  392.   qeau      = 'DBIT' $breche (uus*RVP)                  ;
  393.   YVP= RVP * (inve RHO)                                 ;
  394.   Y = Y + YVP                                           ;
  395.    'FINSI'                                              ;
  396.  
  397.    Qj= qair + qhe + qh2 + qo2 + qn2 + qco + qco2 + qeau ;
  398.  
  399.   Si (Sout < 1.e-10)                                             ;
  400.   deltaq = 0.                                                    ;
  401.   Sinon                                                          ;
  402.    W1=0.                                                         ;
  403. *  deltaq=(w1*(Qj-Qout)/nqj) + ((1.-w1)*((Qout/Qj) - 1.)        );
  404.    deltaq=Qout/Qj     ;
  405.   Finsi                                                          ;
  406.  
  407.     qair = qair*deltaq                                ;
  408.    'SI' TBT.'THE' ; qhe = qhe*deltaq  ;Finsi ;
  409.    'SI' TBT.'TH2' ; qh2 = qh2*deltaq  ;Finsi ;
  410.    'SI' TBT.'TO2' ; qo2 = qo2*deltaq  ;Finsi ;
  411.    'SI' TBT.'TN2' ; qn2 = qn2*deltaq  ;Finsi ;
  412.    'SI' TBT.'TCO' ; qco = qco*deltaq  ;Finsi ;
  413.    'SI' TBT.'TCO2'; qco2= qco2*deltaq ;Finsi ;
  414.    'SI' TBT.'VAPEUR'; qeau = qeau*deltaq ;Finsi ;
  415.  
  416.    Qj= qair + qhe + qh2 + qo2 + qn2 + qco + qco2 + qeau ;
  417.  
  418.   mess '========> Qj=' Qj '=======> Qout=' Qout ' Deltaq=' deltaq ;
  419.  
  420.  'SI'('NON'('EGA' Qj 0.)) ; iQj=1./Qj ; 'SINON' ; iQj=0. ; 'FINSI' ;
  421.   Yvapj = qeau '*' iQj ; Yairj = qair '*' iQj ; Yhej  = qhe  '*' iQj ;
  422.   Yh2j  = qh2  '*' iQj ; Yo2j  = qo2  '*' iQj ; Yn2j  = qn2  '*' iQj ;
  423.   Yco2j = qco2 '*' iQj ; Ycoj  = qco  '*' iQj ;
  424.  
  425. * On ne change pas le signe de Qj mais celui de qeau et autres
  426. aex=-1.;
  427. qeau = aex * qeau ;
  428. qliq = aex * qliq ;
  429. qair = aex * qair ;
  430. qhe  = aex * qhe  ;
  431. qh2  = aex * qh2  ;
  432. qo2  = aex * qo2  ;
  433. qn2  = aex * qn2  ;
  434. qco  = aex * qco  ;
  435. qco2 = aex * qco2 ;
  436. mess 'qeau=' qeau ' qair=' qair ' qo2=' qo2 ' qn2=' qn2;
  437. *   mess ' MINI MAXI YSortie=' (mini Y) (maxi Y);
  438.  
  439. TTINJ=FAUX;
  440. Tinj=Tout;
  441.   'FINSI TTOUT KAS1'                                                   ;
  442. *-------------- FIN TTOUT ----------------------------------------------
  443. *-- Formulation en Extraction ou Clapet de décharge --------------------
  444.  
  445.  
  446.  
  447. *------------------ TTINJ ----------------------------------------------
  448.   'SI  TTINJ KAS1' TTINJ                                               ;
  449.  
  450.   $breche = Briches.'$sreche';
  451.   brechei = Briches.'srechei';
  452.   brechec = Briches.'srechec';
  453.   breche  = Briches.'sreche' ;
  454.  
  455.   Roj=redu Breche tic.'RHO';
  456.   diru    =Briches.'diru'            ;
  457.  
  458.  Si(EGA (VALE DIME) 2);
  459.  gjc   = 'MANU' 'CHPO' brechec 2 'UX' 0.0 'UY' 0.0
  460.                                  'NATURE' 'DISCRET'          ;
  461.  Sinon ;
  462.  gjc   = 'MANU' 'CHPO' brechec 3 'UX' 0.0 'UY' 0.0 'UZ' 0.
  463.                                  'NATURE' 'DISCRET'          ;
  464.  Finsi;
  465.  
  466.  $menvf  = rxt.'GEO'.'$menvf'        ;
  467.  
  468.   tinj      = 'IPOL' tps t (rxtsce .'tinj')                   ;
  469.   qair      = 'IPOL' tps t (rxtsce .'qair')                   ;
  470.    'SI' TBT.'THE' ;
  471.      qhe= 'IPOL' tps t (rxtsce .'qhe')   ; Finsi ;
  472.    'SI' TBT.'TH2' ;
  473.      qh2= 'IPOL' tps t (rxtsce .'qh2')   ; Finsi ;
  474.    'SI' TBT.'TO2' ;
  475.      qo2= 'IPOL' tps t (rxtsce .'qo2')   ; Finsi ;
  476.    'SI' TBT.'TN2' ;
  477.      qn2= 'IPOL' tps t (rxtsce .'qn2')   ; Finsi ;
  478.    'SI' TBT.'TCO';
  479.      qco= 'IPOL' tps t (rxtsce .'qco')   ; Finsi ;
  480.    'SI' TBT.'TCO2' ;
  481.      qco2='IPOL' tps t (rxtsce .'qco2')  ; Finsi ;
  482.  
  483.      Qj= qair + qhe + qh2 + qo2 + qn2 + qco + qco2      ;
  484.  
  485.    'SI' TBT.'VAPEUR'                                    ;
  486.      qeau = 'IPOL' tps t (rxtsce .'qeau')        ;
  487.      Qj = Qj + qeau                                     ;
  488.      Pt =extr tic.'PT'  ndl                             ;
  489.      Tfm=extr tic.'Tfm' ndl                             ;
  490.  
  491.    'SI'('NON'('EGA' Qj 0.)) ; iQj=1./Qj ; 'SINON' ; iQj=0. ; 'FINSI' ;
  492.     Yvapj = qeau '*' iQj ; Yairj = qair '*' iQj ; Yhej  = qhe  '*' iQj ;
  493.     Yh2j  = qh2  '*' iQj ; Yo2j  = qo2  '*' iQj ; Yn2j  = qn2  '*' iQj ;
  494.     Yco2j = qco2 '*' iQj ; Ycoj  = qco  '*' iQj ;
  495.  
  496.  
  497.   Mj=((Yvapj/Mkvap)+(Yairj/Mkair)+(Yhej/Mkhe)+(Yh2j/Mkh2)+(Yo2j/Mko2)+
  498.      (Yn2j*Mkn2)+(Yco2j/Mkco2)+(Ycoj/Mkco))**(-1.)      ;
  499.        Xvj=Yvapj*Mj/Mkvap                               ;
  500.        Pvj  = Pt * Xvj                                  ;
  501.  
  502. *      'SI' (TBT.'CONDMAS')                              ;
  503. *       Pt =extr tic.'PT'  ndl                           ;
  504. *       Tfm=extr tic.'Tfm' ndl                           ;
  505. *     Qj=qeau + qair + qhe + qh2 + qo2 + qn2 + qco + qco2;
  506. *a peu près Xvj = Yvj on fera mieux après
  507. *       Xvj  = qeau/Qt                                   ;
  508. *       Pvj  = Pt * Xvj                                  ;
  509. *       Hinj = VARI HVS Pvj (Tinj + 273.15)              ;
  510. *       Tsat = VARI TSATP  Pt                            ;
  511. *       Hliq = VARI HLS Pvj Tsat                         ;
  512. *       Hvap = VARI HVS Pvj Tsat                         ;
  513. *       mess ' Hinj=' Hinj ' Hliq=' Hliq ' Hvap=' Hvap   ;
  514. *       mess 'Pvj=' pvj ' Tsat = ' Tsat ' Tinj=' (Tinj + 273.15);
  515. *       Ksi=(Hinj - Hliq)/(Hvap - Hliq)                  ;
  516. *       Ksi=MINI (PROG 1. Ksi)                           ;
  517. *       mess ' Ksi  ' Ksi ' qeau= ' qeau ' qvap=' (qeau*Ksi);
  518. **?       Mliqbr= Mliqbr + (qeau*(1. - Ksi))               ;
  519. *       qeau = qeau * Ksi                                ;
  520. *      'FINSI'                                           ;
  521.  
  522.    'SINON'                                              ;
  523.  
  524.    'SI'('NON'('EGA' Qj 0.)) ; iQj=1./Qj ; 'SINON' ; iQj=0. ; 'FINSI' ;
  525.     Yvapj = qeau '*' iQj ; Yairj = qair '*' iQj ; Yhej  = qhe  '*' iQj ;
  526.     Yh2j  = qh2  '*' iQj ; Yo2j  = qo2  '*' iQj ; Yn2j  = qn2  '*' iQj ;
  527.     Yco2j = qco2 '*' iQj ; Ycoj  = qco  '*' iQj ;
  528.  
  529.  
  530.    'FINSI'                                              ;
  531.  
  532.    Si(EGA Qj 0.)                                        ;
  533.       Roj   = 0.                                        ;
  534.    Sinon                                                ;
  535.   Pt    = 'EXTR' (tic.'PT')    ndl                      ;
  536.   Rgj   = (TBT.'Rgvap' * Yvapj) + (TBT.'Rgair' * Yairj) +
  537.           (TBT.'Rghe'  * Yhej)  + (TBT.'Rgh2'  * Yh2j)  +
  538.           (TBT.'Rgo2'  * Yo2j)  + (TBT.'Rgn2'  * Yn2j)  +
  539.           (TBT.'Rgco'  * Ycoj)  + (TBT.'Rgco2' * Yco2j) ;
  540.  
  541.       Roj   = Pt '/' (Tinj '+' 273.15) '/' Rgj          ;
  542.    Finsi                                                ;
  543.  
  544.   'FINSI TTINJ KAS1'                                    ;
  545. *-------------- FIN TTINJ ----------------------------------------------
  546.  
  547.      rxtsce .'Yvapj'= Yvapj                      ;
  548.      rxtsce .'Yairj'= Yairj                      ;
  549.      rxtsce .'Yhej' = Yhej                       ;
  550.      rxtsce .'Yh2j' = Yh2j                       ;
  551.      rxtsce .'Yo2j' = Yo2j                       ;
  552.      rxtsce .'Yn2j' = Yn2j                       ;
  553.      rxtsce .'Yco2j'= Yco2j                      ;
  554.      rxtsce .'Ycoj' = Ycoj                       ;
  555.  
  556.      rxtsce .'aex'  = aex                        ;
  557.      rxtsce .'Qj'   = Qj                         ;
  558.      rxtsce .'Qlj'  = Qlj                        ;
  559.      rxtsce .'Tj'   = tinj                       ;
  560.      rxtsce .'Roj'  = Roj                        ;
  561.  
  562.   Briches.CTinj= Briches.CTinj 'ET' ('PROG' Tinj)        ;
  563.   Briches.CKsi = Briches.CKsi  'ET' ('PROG' Ksi)         ;
  564.   Briches.CQlj = Briches.CQlj  'ET' ('PROG' qliq)        ;
  565.   Briches.CHlj = Briches.CHlj  'ET' ('PROG' hliq)        ;
  566. *-------------------------------------------------------------
  567. *Caractéristiques de l'injection au pas de temps N
  568. * c'est a dire Dtn = tn - t(n-1) a l'instant n-1
  569.  
  570.  'SI' ( TBT.'IMPR' >EG 2)                                    ;
  571.  'MESS' '==> Détail des débits extraits (kg/s) Sortie ' nub  ;
  572.  'MESS' 'Qeau=' qeau 'Qliq=' qliq 'Qair=' qair 'Qhe =' qhe 'Qh2 =' qh2 ;
  573.  'MESS' 'Qo2 =' qo2  'Qn2 =' qn2  'Qco =' qco  'Qco2=' qco2  ;
  574.  'FINSI'                                                     ;
  575.  
  576. *mess '=================== FIN KAS1 =======================';
  577.   resproc qeau qliq qair qhe qh2 qo2 qn2 qco qco2         ;
  578.  
  579.    QUITTER BROCHE                                         ;
  580. 'FINSI KAS1'                                              ;
  581. ************************************************************************
  582. ********* FIN KAS1 *****************************************************
  583. ************************************************************************
  584.  
  585. ************************************************************************
  586. ********* DEBUT KAS2 ***************************************************
  587. ************************************************************************
  588.  
  589. 'SI  KAS2' (EGA KAS 'KAS2')          ;
  590.  
  591.   CSbreche =chai 'Ssreche'  ;
  592.   C$breche =chai '$sreche'  ;
  593.   C$brechc =chai '$srechec' ;
  594.   Cfacgeo  =chai 'facgeo'   ;
  595.   Cbrechei =chai 'srechei'  ;
  596.   Cbrechec =chai 'srechec'  ;
  597.   Cbreche  =chai 'sreche'   ;
  598.  
  599.  
  600.   diru    =Briches.'diru'            ;
  601.  
  602.  diru    = diru/(norm diru)          ;
  603.  Sbreche = GEO.CSbreche              ;
  604.  $breche = GEO.C$breche              ;
  605.  facgeo  = GEO.Cfacgeo               ;
  606.  brechei = GEO.Cbrechei              ;
  607.  brechec = GEO.Cbrechec              ;
  608.  breche  = GEO.Cbreche               ;
  609.  
  610.  $menvf  = rxt.'GEO'.'$menvf'        ;
  611.  $vtf    = rxt.'GEO'.'$vtf'          ;
  612.  
  613.      Roj = rxtsce .'Roj'                        ;
  614.      Qj  = rxtsce .'Qj'                         ;
  615.      Qlj = rxtsce .'Qlj'                        ;
  616.      tinj= rxtsce .'Tj'                         ;
  617.  
  618. **************************************************************
  619. **** Evolution de la source en fct du temps ******************
  620.  
  621.   t         = rxtsce .'t'                             ;
  622.  
  623. * Récupération des fractions massiques
  624.   Yvapj= rxtsce .'Yvapj'                              ;
  625.   Yairj= rxtsce .'Yairj'                              ;
  626.   Yhej = rxtsce .'Yhej'                               ;
  627.   Yh2j = rxtsce .'Yh2j'                               ;
  628.   Yo2j = rxtsce .'Yo2j'                               ;
  629.   Yn2j = rxtsce .'Yn2j'                               ;
  630.   Yco2j= rxtsce .'Yco2j'                              ;
  631.   Ycoj = rxtsce .'Ycoj'                               ;
  632.  
  633.  Qjf   = Qj '*' facgeo                                       ;
  634.  
  635.  'SI' ('EGA' Qj 0.0)                                         ;
  636.       Roj   = 0.0                                            ;
  637.       guj   = 0.0                                            ;
  638.  'SINON'                                                     ;
  639.       Roj   = rxtsce .'Roj'                                  ;
  640.       guj   = (Qjf '/' Sbreche '/' Roj )                     ;
  641.  'FINSI'                                                     ;
  642.  
  643.   Rvapj = Roj   '*' Yvapj                             ;
  644.   Rairj = Roj   '*' Yairj                             ;
  645.   Rhej  = Roj   '*' Yhej                              ;
  646.   Rh2j  = Roj   '*' Yh2j                              ;
  647.   Ro2j  = Roj   '*' Yo2j                              ;
  648.   Rn2j  = Roj   '*' Yn2j                              ;
  649.   Rco2j = Roj   '*' Yco2j                             ;
  650.   Rcoj  = Roj   '*' Ycoj                              ;
  651.  
  652. * Les donnees physiques
  653.  Cpvap  = TBT.'Cpvap' ;
  654.  Lv     = TBT.'Lv'    ;
  655.  Rabs   = TBT.'Rabs'  ;
  656.  Rgh2   = TBT.'Rgh2'  ; Mkh2   = TBT.'Mkh2'              ;
  657.  Rghe   = TBT.'Rghe'  ; Mkhe   = TBT.'Mkhe'              ;
  658.  Rgo2   = TBT.'Rgo2'  ; Mko2   = TBT.'Mko2'              ;
  659.  Rgn2   = TBT.'Rgn2'  ; Mkn2   = TBT.'Mkn2'              ;
  660.  Rgco2  = TBT.'Rgco2' ; Mkco2  = TBT.'Mkco2'             ;
  661.  Rgco   = TBT.'Rgco'  ; Mkco   = TBT.'Mkco'              ;
  662.  Rgvap  = TBT.'Rgvap' ; Mkvap  = TBT.'Mkvap'             ;
  663.  Rgair  = TBT.'Rgair' ; Mkair  = TBT.'Mkair'             ;
  664.  
  665. * Calcul Cp (Tj)   CALCP
  666.   Tj    = tinj                                               ;
  667.  
  668. Cph2j Cphej Cpo2j Cpn2j Cpco2j Cpcoj Cpairj = CALCP Tj       ;
  669.  
  670. * Calcul des Cvs pour l'injection
  671.   Cpvapj      = TBT.'Cpvap'                                  ;
  672.   Cvvapj      = Cpvapj '-'  Rgvap                            ;
  673.   Cvhej       = Cphej  '-'  Rghe                             ;
  674.   Cvh2j       = Cph2j  '-'  Rgh2                             ;
  675.   Cvo2j       = Cpo2j  '-'  Rgo2                             ;
  676.   Cvn2j       = Cpn2j  '-'  Rgn2                             ;
  677.   Cvco2j      = Cpco2j '-'  Rgco2                            ;
  678.   Cvcoj       = Cpcoj  '-'  Rgco                             ;
  679.   Cvairj      = Cpairj '-'  Rgair                            ;
  680.  
  681.   Cvj   = (Cvvapj * Yvapj)+(Cvairj * Yairj)+(Cvhej * Yhej)+
  682.           (Cvh2j  * Yh2j)+(Cvo2j * Yo2j)+(Cvn2j * Yn2j)+
  683.           (Cvcoj  * Ycoj)+(Cvco2j * Yco2j)                   ;
  684.   Cpj   = (Cpvapj * Yvapj)+(Cpairj * Yairj)+(Cphej * Yhej)+
  685.           (Cph2j  * Yh2j)+(Cpo2j * Yo2j)+(Cpn2j * Yn2j)+
  686.           (Cpcoj  * Ycoj)+(Cpco2j * Yco2j)                   ;
  687.   Rgj   = (Rgvap * Yvapj)+(Rgair * Yairj)+(Rghe * Yhej)+
  688.           (Rgh2  * Yh2j)+(Rgo2 * Yo2j)+(Rgn2 * Yn2j)+
  689.           (Rgco  * Ycoj)+(Rgco2 * Yco2j)                     ;
  690.  
  691.   Hj    = Cpj   '*' (Tj '+' 273.15)                          ;
  692.   Ej    = Cvj   '*' (Tj '+' 273.15)                          ;
  693.  
  694.   Lbreche=Sbreche**(1./IDIM)                                 ;
  695.   Knj=0.05 * guj                                             ;
  696.   Enj=(Knj**1.5)/Lbreche                                     ;
  697.  
  698.  'SI' ( TBT.'IMPR' >EG 2)                                    ;
  699.  'MESS' '==> Caracteristiques de l extraction ' nub          ;
  700.  'MESS' ' Débit total extrait Qj=' Qj 'kg/s'                 ;
  701.  'MESS' ' Vitesse extraction guj=' guj 'm/s'                 ;
  702.  'MESS' ' Densité extraction Roj=' Roj 'kg/m3'               ;
  703.  'MESS' ' R/M     extraction Rgj=' Rgj 'J/kg/K'              ;
  704.  'FINSI'                                                     ;
  705.  
  706.   src=0.;
  707.  
  708. 'SI' TBT.'FPAROI'   ;
  709.   brchp=breche      ;
  710. 'SINON'             ;
  711.   brchp=brechei     ;
  712. 'FINSI'             ;
  713.  
  714. 'SI' DIM3D                                                   ;
  715.  $brechec= GEO.C$brechc                                      ;
  716.   gj    = 'KCHT' $breche 'VECT' 'SOMMET' 'COMP'
  717.            ('MOTS' 1UN 2UN 3UN) (diru * guj)                 ;
  718.   gj    = 'REDU' gj brchp                                    ;
  719.   gjc   = 'KCHT' $brechec 'VECT' 'SOMMET' (0.0 0.0 0.0)      ;
  720.   gjx   = 'KCHT' $breche 'VECT' 'SOMMET' (diru *  guj)   gjc ;
  721. 'SINON'                                                      ;
  722.   gj    = 'KCHT' $breche 'VECT' 'SOMMET' 'COMP'
  723.           ('MOTS' 1UN 2UN) (diru * guj)                      ;
  724.   gj    = 'REDU' gj brchp                                    ;
  725.  gjc   = 'MANU' 'CHPO' brechec 2 'UX' 0.0 'UY' 0.0
  726.                                  'NATURE' 'DISCRET'          ;
  727.   gjx   = 'KCHT' $breche 'VECT' 'SOMMET' (diru * guj) gjc    ;
  728. 'FINSI'                                                      ;
  729.  
  730.   unvf0 = 'KCHT' $menvf  'VECT' 'SOMMET'  gjx                ;
  731.   src   = 'DBIT' unvf0 $menvf                                ;
  732.  
  733. *=============================================================
  734.  
  735.   Briches.'Rvapj'=Rvapj                                      ;
  736.   Briches.'Rairj'=Rairj                                      ;
  737.   Briches.'Rhej' =Rhej                                       ;
  738.   Briches.'Rh2j' =Rh2j                                       ;
  739.   Briches.'Ro2j' =Ro2j                                       ;
  740.   Briches.'Rn2j' =Rn2j                                       ;
  741.   Briches.'Rco2j'=Rco2j                                      ;
  742.   Briches.'Rcoj' =Rcoj                                       ;
  743.   Briches.'Rgj'  =Rgj                                        ;
  744. *  Briches.'Hj'   =Hj                                         ;
  745. *  Briches.'Ej'   =Ej                                         ;
  746.   Briches.'Tj'   =Tj                                         ;
  747. *  Briches.'Qj'   =Qj                                         ;
  748.   Briches.'Roj'  =Roj                                        ;
  749.   Briches.'gj'   =gj                                         ;
  750. *  Briches.'guj'  =guj                                        ;
  751.   Briches.'unvf0'=unvf0                                      ;
  752.   Briches.'src'  =src                                        ;
  753.      aex = rxtsce .'aex'                                     ;
  754. *mess '=================== FIN KAS2 =======================';
  755.  resproc
  756.  Rvapj Rairj Rhej Rh2j Ro2j Rn2j Rco2j Rcoj Rgj Hj Ej
  757.  Tj aex Qj Roj gj guj unvf0 src                              ;
  758.  
  759.    QUITTER BROCHE                                            ;
  760.  
  761. 'FINSI KAS2'                                                 ;
  762. ************************************************************************
  763. ********* FIN KAS2 *****************************************************
  764. ************************************************************************
  765.  
  766.  
  767. FINPROC ;
  768.  
  769.  
  770.  

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