Télécharger la plaquette Cast3M


fr en

Télécharger test-aspH.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  test-aspH.dgibi
  2. * KTEST: changement de phase sur gouttes vers un régime permanent 
  3. *        - pas de gravité
  4. *        - pas d'injection
  5. *        - option COMPLET régime permanent 
  6. *        - (rapport de Stage Antoine BOUQUET (voir L. Blumenfeld)) 
  7. *        - repartition homogène au depart
  8. *$$$$ VISUP
  9. ** PROCEDURE DE VISUALISATION DES RESULTATS 3D
  10. DEBPROC VISUP tps*'FLOTTANT' rxt*'TABLE'           ;
  11.  
  12. NUPADT = rxt.'TIC'.'NUPADT';
  13. LTPS   = rxt.'TIC'.'LTPS';
  14.  
  15. LOG1=FAUX;
  16. 'SI' (EXIST rxt.'TBT'.'TABPERSO' 'CAPT0');
  17.   LOG1=rxt.'TBT'.'TABPERSO'.'CAPT0';
  18. 'FINSI';
  19. 'SI' LOG1;
  20.   TEMPS = rxt.'TBT'.'TABPERSO'.'TEMPS';
  21.   NTPS1= DIME TEMPS;
  22.   FREQ0  = rxt.'TBT'.'TABPERSO'.'FREQ0' ;
  23.   NULOC0 = (NTPS1+1)*FREQ0  ;
  24.   'SI' (NUPADT EGA NULOC0);
  25.     MESS 'Execution du captage de donnees ';
  26.     CAPTP rxt;
  27.     TI1 = EXTR LTPS (NULOC0+1) ;
  28.     rxt.'TBT'.'TABPERSO'.'TEMPS' = 
  29.            rxt.'TBT'.'TABPERSO'.'TEMPS' et (PROG TI1);
  30.   'FINSI';
  31. 'FINSI';
  32.  
  33. LOG2=FAUX;
  34. 'SI' (EXIST rxt.'TBT'.'TABPERSO' 'CH3D0');
  35.   LOG2=rxt.'TBT'.'TABPERSO'.'CH3D0';
  36. 'FINSI';
  37. 'SI' LOG2;
  38.   TPS3D = rxt.'TBT'.'TABPERSO'.'TPS3D';
  39.   NTPS2= DIME TPS3D;
  40.   FREQ3D = rxt.'TBT'.'TABPERSO'.'FREQ3D';
  41.   NULOC3D= (NTPS2+1)*FREQ3D ;
  42.   'SI' (NUPADT EGA NULOC3D);
  43.     MESS 'Execution du trace de champ 3D suivant une coupe ';
  44.     CH3D rxt;
  45.     TI2 = EXTR LTPS (NULOC3D+1) ;
  46.     rxt.'TBT'.'TABPERSO'.'TPS3D' = 
  47.            rxt.'TBT'.'TABPERSO'.'TPS3D' et (PROG TI2);
  48.   'FINSI';
  49. 'FINSI';
  50.  
  51.  
  52. FINPROC       ;
  53.  
  54.  
  55. *$$$$ CAPTP
  56. ** Procedure pour capteur de Temp en des points donnes au cours du temps
  57. 'DEBPROC' CAPTP rxt*'TABLE';
  58.  
  59. L1= rxt.'TBT'.'TABPERSO'.'COMP';
  60. N1= DIME L1;
  61.  
  62. PTS =rxt.'TBT'.'TABPERSO'.'LPOINTS';
  63. M1  =rxt.'TBT'.'TABPERSO'.'NPOINTS';
  64. CAPT=rxt.'TBT'.'TABPERSO'.'CAPT'   ;
  65.  
  66. vtf=rxt.'GEO'.'vtf';
  67.  
  68. TEMPS = rxt.'TBT'.'TABPERSO'.'TEMPS';
  69. LTEMPS= DIME TEMPS;
  70.  
  71.  
  72. *   Initialisation des tables
  73. 'SI' (LTEMPS EGA 0 );
  74.   'REPETER' BOU1 M1;
  75.     PT1=PTS 'POIN' &BOU1;
  76.     CAPI= EXTR CAPT &BOU1;
  77.     rxt.'TBT'.'TABPERSO'.CAPI=TABLE;
  78.        'REPETER' BOU2 N1;
  79.           MOI = EXTR L1 &BOU2 ;
  80.           rxt.'TBT'.'TABPERSO'.CAPI.MOI=PROG;
  81.        'FIN' BOU2;
  82.   'FIN' BOU1;
  83. 'FINSI';
  84.  
  85.  
  86.  
  87. 'REPETER' BOU1 M1;
  88.  
  89.   PT1=PTS 'POIN' &BOU1;
  90.   CAPI= EXTR CAPT &BOU1;
  91.  
  92.   'REPETER' BOU2 N1;
  93.  
  94.     MOI = EXTR L1 &BOU2 ;
  95.  
  96.     'SI' ('EGA' MOI 'TF');
  97.         TF=rxt.'TIC'.'TF';
  98.         PT2=vtf 'POINT' 'PROCHE' PT1;
  99.         TF1=REDU TF PT2;
  100.         X1 = 'EXTR' TF1 'SCAL' PT2 ;
  101.  
  102.         rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  103.                  rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  104.     'FINSI';
  105.  
  106.   'SI' ('EGA' MOI 'XVAP');
  107.         YVAP = rxt.'TIC'.'YVAP'                    ;
  108.          l29 = yvap*0.029                          ;
  109.          l18 = yvap*0.018                          ;
  110.          lc = 0.018                                ;
  111.          xvap = (l29/(lc+l29-l18))                 ;
  112.  
  113.         PT2=vtf 'POINT' 'PROCHE' PT1                ;
  114.         Xvap1=REDU Xvap PT2                         ;
  115.         X1 = 'EXTR' Xvap1 'SCAL' PT2                ;
  116.  
  117.         rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  118.                  rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  119.     'FINSI';
  120.  
  121.     'SI' ('EGA' MOI 'XD');
  122.         XD=rxt.'TIC'.'XD';
  123.         PT2=vtf 'POINT' 'PROCHE' PT1;
  124.         XD1=REDU XD PT2;
  125.         X1 = 'EXTR' XD1 'SCAL' PT2 ;
  126.  
  127.         rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  128.                  rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  129.      'FINSI';
  130.  
  131.  
  132.  
  133.     'SI' ('EGA' MOI 'TD');
  134.         TD=rxt.'TIC'.'TD';
  135.         PT2=vtf 'POINT' 'PROCHE' PT1;
  136.         TD1=REDU TD PT2;
  137.         X1 = 'EXTR' TD1 'SCAL' PT2 ;
  138.  
  139.         rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  140.                  rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  141.      'FINSI';
  142.  
  143.  
  144.      'SI' ('EGA' MOI 'STD2');
  145.         STD2= elno $vtf rxt.'TIC'.'STD';
  146.         PT2=vtf 'POINT' 'PROCHE' PT1;
  147.         STD1=REDU STD2 PT2;
  148.         X1 = 'EXTR' STD1 'SCAL' PT2 ;
  149.  
  150.         rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  151.                  rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  152.      'FINSI';
  153.  
  154.  
  155.  
  156.  
  157.  
  158.       'SI' ('EGA' MOI 'DD');
  159.         DD=rxt.'TIC'.'DD';
  160.         PT2=vtf 'POINT' 'PROCHE' PT1;
  161.         DD1=REDU DD PT2;
  162.         X1 = 'EXTR' DD1 'SCAL' PT2 ;
  163.  
  164.         rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  165.                  rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  166.      'FINSI';
  167.  
  168.     'SI' ('EGA' MOI 'UX');
  169.       UN=rxt.'TIC'.'UN';
  170.       UX=EXCO 'UX' UN;
  171.       PT2=vtf 'POINT' 'PROCHE' PT1;
  172.       UX1=REDU UX PT2;
  173.       X1 = 'EXTR' UX1 'SCAL' PT2 ;
  174.  
  175.       rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  176.                  rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  177.     'FINSI';
  178.  
  179.     'SI' ('EGA' MOI 'UY');
  180.       UN=rxt.'TIC'.'UN';
  181.       UY=EXCO 'UY' UN;
  182.       PT2=vtf 'POINT' 'PROCHE' PT1;
  183.       UY1=REDU UY PT2;
  184.       X1 = 'EXTR' UY1 'SCAL' PT2 ;
  185.  
  186.       rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  187.                  rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  188.     'FINSI';
  189.  
  190.     IDIM = vale 'DIME';
  191.     'SI' ('EGA' IDIM 3);
  192.       'SI' ('EGA' MOI 'UZ');
  193.         UN=rxt.'TIC'.'UN';
  194.         UZ=EXCO 'UZ' UN;
  195.         PT2=vtf 'POINT' 'PROCHE' PT1;
  196.         UZ1=REDU UZ PT2;
  197.         X1 = 'EXTR' UZ1 'SCAL' PT2 ;
  198.  
  199.         rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  200.                  rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  201.       'FINSI';
  202.     'FINSI';
  203.  
  204. 'SI' ('EGA' MOI 'VX');
  205.       VN=rxt.'TIC'.'VN';
  206.       VX=EXCO 'UX' VN;
  207.       PT2=vtf 'POINT' 'PROCHE' PT1;
  208.       VX1=REDU VX PT2;
  209.       X1 = 'EXTR' VX1 'SCAL' PT2 ;
  210.  
  211.       rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  212.                  rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  213.     'FINSI';
  214.  
  215.     'SI' ('EGA' MOI 'VY');
  216.       VN=rxt.'TIC'.'VN';
  217.       VY=EXCO 'UY' VN;
  218.       PT2=vtf 'POINT' 'PROCHE' PT1;
  219.       VY1=REDU VY PT2;
  220.       X1 = 'EXTR' VY1 'SCAL' PT2 ;
  221.  
  222.       rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  223.                  rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  224.     'FINSI';
  225.  
  226.     IDIM = vale 'DIME';
  227.     'SI' ('EGA' IDIM 3);
  228.       'SI' ('EGA' MOI 'VZ');
  229.         VN=rxt.'TIC'.'VN';
  230.         VZ=EXCO 'UZ' VN;
  231.         PT2=vtf 'POINT' 'PROCHE' PT1;
  232.         VZ1=REDU VZ PT2;
  233.         X1 = 'EXTR' VZ1 'SCAL' PT2 ;
  234.  
  235.         rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  236.                  rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  237.       'FINSI';
  238.     'FINSI';
  239.  
  240.   'FIN' BOU2;
  241.  
  242. 'FIN' BOU1;
  243.  
  244.  
  245. 'FINPROC' ;
  246.  
  247.  opti dime 3 elem cu20                                        ;
  248.  
  249.  COMPLET=VRAI                                                 ;
  250.  COMPLET=FAUX                                                 ;
  251.  GRAPH=VRAI                                                   ;
  252.  GRAPH=FAUX                                                   ;
  253.  opti TRACE 'PSC'                                             ;
  254.  
  255.  Si COMPLET                                                   ;
  256.  DT0          =   2.e-2                                       ;
  257. * nbit=20                                                     ;
  258.  nbit=200                                                     ;
  259.  n1   =  2          ; n2   =  2          ; n3   =  1          ;
  260. * n4   =  -1 * 10          ;  
  261.   n4   =  8                ;                                  
  262.  Sinon                                                        ;
  263. * nbit=100                                                    ;
  264.  DT0          =   5.e-2                                       ;
  265. *DT0          =   1.e-1                                       ;
  266.  nbit = 10                                                    ;
  267.  n1   =  1          ; n2   =  2          ; n3   =  1          ;
  268.  n4   =  4          ;                                         ;
  269.  Finsi                                                        ;
  270.  
  271.  ri   = 0.75      ; sp   =  1.000      ; hc   = 4             ;
  272.  epsi =  1.000e-6   ;                                         ;
  273.  
  274.  p0   =     0.000             0.000                0.000      ;
  275.  px   = -1000.000             0.000                0.000      ;
  276.  py   =     0.000         -1000.000                0.000      ;
  277.  pz   =     0.000             0.000             1000.000      ;
  278.  cd   =     0.000             0.000              -20.000      ;
  279.  ph0  =     0.000             0.000                hc         ;
  280.  phx  =     ri                0.000                hc         ;
  281.  phy  =     0.000             ri                   hc         ;
  282.  
  283.  p1   =  (ri*0.6)             0.000                0.000      ;
  284.  p2   =  (ri*0.6)             (ri*0.6)             0.000      ;
  285.  p3   =  0.000                (ri*0.6)             0.000      ;
  286.  p4   =   ri                  0.000                0.000      ;
  287.  p5   =  0.000                 ri                  0.000      ;
  288.  p6   =  (ri+sp)              0.000                0.000      ;
  289.  p7   =  0.000                (ri+sp)              0.000      ;
  290.  
  291.  l1   =  droi p0 p1 n1                                        ;
  292.  l2   =  droi p1 p2 n1                                        ;
  293.  l3   =  droi p2 p3 n1                                        ;
  294.  l4   =  droi p3 p0 n1                                        ;
  295.  l5   =  cerc p4 p0 p5 (2*n1) 'DINI' 0.4 'DFIN'  0.01         ;
  296.  l6   =  cerc p6 p0 p7 (2*n1)                                 ;
  297.  
  298.  basf =  dall l1 l2 l3 l4                                     ;
  299.  basf =  basf et (regl (l2 et l3) l5 n2) ; elim basf epsi     ;
  300.  basf =  basf et (syme  basf droi p0 p3) ; elim basf epsi     ;
  301.  basf =  basf et (syme  basf droi p0 p1) ; elim basf epsi     ;
  302.  
  303.  l44  =  (l4 et (p3 d n2 p5) )                                ;
  304.  l44  =  l44 et (syme l44 droi p0 p1) ; elim l44 epsi         ;
  305.  l44  =  l44 coul rouge                                       ;
  306.  l55  =  l44 tour 90. p0 (p0 plus (0 0 100)) ;
  307.  
  308.  basm =  regl l5 l6 n3                                        ;
  309.  basm =  basm et (syme  basm droi p0 p3) ; elim basm epsi     ;
  310.  basm =  basm et (syme  basm droi p0 p1) ; elim basm epsi     ;
  311.  
  312.  domf =  basf proj dire pz sphe cd phx                        ;
  313.  domm =  basm proj dire pz plan  ph0 phx phy                  ;
  314.  
  315.  l44s= l44 proj dire pz sphe cd phx ;
  316.  l55s= l55 proj dire pz sphe cd phx ;
  317.  
  318.  p44 = regle n4 l44 l44s ;
  319.  p55 = regle n4 l55 l55s ;
  320.  
  321.  p0s= p0 proj dire pz sphe cd phx ;
  322.  axe = droit n4 p0  p0s 'DINI' 0.4 'DFIN' 0.05 ;
  323.  elim (axe et p55) epsi ;
  324. *opti donn 5;
  325.  
  326.  
  327.  aria =  volu basf domf n4 'DINI' 0.4 'DFIN' 0.05             ;
  328.  muro = (volu basm domm n4 'DINI' 0.4 'DFIN' 0.05) coul roug  ;
  329.  elim (aria et muro) epsi                                     ;
  330.  
  331.  elim (aria et p44 et p55 ) epsi ;
  332.  pjg  =  poin basf proc (0.000    0.000    0.000)             ;
  333. *pjg  =  poin basf proc ((ri/2.0) 0.000    0.000)             ;
  334.  jg   = (elem basf 'APPUIE' 'LARGEMENT' pjg) coul vert        ;
  335.  
  336.  pjd  =  poin domf proc (0.000 ri    hc)                      ;
  337. *pjd  =  poin domf proc ((ri/2.0) 0.000    hc)                ;
  338.  jd   =  (elem domf 'APPUIE' 'LARGEMENT'  pjd) coul vert      ;
  339.  
  340.  mesjd = MESU jd ;
  341.  'MESS' mesjd ; 
  342.  
  343.  
  344.  
  345.  enva = enve aria                                             ;
  346.  enva = orie enva point (0.0 0.0 (hc/2.0))                    ;
  347.  envm = enve muro                                             ;
  348. * pare = inte enva envm                                        ;
  349.  
  350.  mesaria = MESU aria ;
  351.  'MESS' mesaria ; 
  352. * trac aria;
  353. * opti donn 5;
  354.  
  355.  
  356. * DT0          =   1e-1                                       ;
  357.  rxt          =   table                                       ;
  358.  rxt.'VERSION'= 'V0'                                          ;
  359.  rxt.'vtf'    =   aria                                        ;
  360.  rxt.'epsi'   =   epsi                                        ;
  361.  rxt.'pi'     =  (0.0 0.0 1.1)                                ;
  362.  rxt.'GRAPH'     =  GRAPH                                     ;
  363.  rxt.'FRPREC'    = 3                                          ;
  364. *
  365.  rxt.'VAPEUR'    =  VRAI                                      ;
  366.  rxt.'ASPER'     =  VRAI                                      ;
  367.  
  368. *-- Conditions initiales dans l'enceinte de test
  369.  rxt.'TF0'    =        120.0                                  ;
  370.  rxt.'PT0'    =        2.50e5                                 ;
  371.  rxt.'Yvap0'  =        0.01                                   ;
  372.  
  373.  rxt.'scenasp'          = table                               ;
  374.  rxt.'scenasp'.'t'      =  PROG   0.000  1000.000             ;
  375.  
  376.  ASPER  = rxt.'ASPER'                                         ;
  377.  si ASPER                                                     ;
  378.   rxt.'aspinj' =   jd                                         ;
  379.   rxt.'toitf'  =   domf                                       ;
  380.   rxt.'rod'    =   1.000e+3                                   ;
  381.   rxt.'Cpd'    =   4.180e+3                                   ;
  382.   rxt.'scenasp'.'vzinj'  = prog   0.     0.                   ;
  383.   rxt.'scenasp'.'xdinj'  = prog   0.     0.                   ;
  384.   rxt.'scenasp'.'tdinj'  = prog   9.50e+1     9.50e+1         ;
  385.   rxt.'scenasp'.'ddinj'  = prog   1.400e-3    1.400e-3        ;
  386.  finsi                                                        ;
  387. *
  388.  
  389.  rxt.'DT0'    =  DT0                                          ;
  390.  
  391. *-- On impose la viscosite turbulente (m2/s)
  392.  rxt.'MODTURB'='NUTURB'                                       ;
  393. * nut/nuf=330                                                 ; 
  394.  rxt.'NUT'    =  5.e-3                                        ;
  395.  
  396.  rxt.'IMPR' =1                                                ;
  397.  
  398.  
  399. **********************************************
  400. *-- Definition des capteurs de mesures dans l"enceinte                  
  401. R00 =   0.  0.  1.0          ;                                                     
  402. R01 =   0.  0.  1.5          ;                                                    
  403. R02 =   0.  0.  2.0          ;
  404.  
  405. PERSO=TABLE                ;                                                            
  406. PERSO.CAPT0 = VRAI         ;                                                    
  407. PERSO.TITRE = 'Tableperso' ;                                            
  408. PERSO.TEMPS = 'PROG'       ;                                          
  409. PERSO.FREQ0 = 1            ;                    
  410. PERSO.LPOINTS= R00 et R01 et R02 ;
  411. PERSO.NPOINTS= 3                 ;
  412. PERSO.CAPT='MOTS' 'CAP1' 
  413.  'CAP2' 'CAP3'  ;
  414. PERSO.COMP= MOTS  
  415. 'TD' 'UY' 'UX' 'XVAP' 'VY' 'VX' 'DD' 'XD' 'TF'  'STD2' ;
  416.  
  417.  
  418.  
  419. *************************************                                   
  420. ** Procedure PERSO pour            **
  421. ** et suivi en ligne des capteurs  **                                   
  422. ************************************* 
  423. rxt.'PERSO'   = VRAI              ; 
  424. rxt.'PRCPERSO'= 'VISUP'           ;     
  425. rxt.'TABPERSO'= PERSO             ;
  426.  
  427.  
  428.  EXECRXT 0  rxt                                              ;
  429.  
  430.  xdd = 'KCHT' (rxt.'GEO'.'$vtf') 'SCAL' 'SOMMET' 10.e-4       ;
  431.  rxt.'TIC'.'XD' =  xdd                                        ;  
  432.  rxt.'TIC'.'XDM' =  xdd                                       ; 
  433.  dd = 'KCHT' (rxt.'GEO'.'$vtf') 'SCAL' 'SOMMET' 1.4e-3        ;
  434.  rxt.'TIC'.'DD' =  dd                                         ;  
  435.  rxt.'TIC'.'DDM' =  xdd                                       ; 
  436.  td = 'KCHT' (rxt.'GEO'.'$vtf') 'SCAL' 'SOMMET' 95.00         ;
  437.  rxt.'TIC'.'TD' =  td                                         ;  
  438.  rxt.'TIC'.'TDM' =  td                                        ; 
  439.  
  440.  EXECRXT nbit  rxt              ;
  441.  
  442.  
  443. *opti sauv  'asptest1.save';
  444. *sauv                    ;
  445.  
  446.  Maxe = chan axe QUAF;
  447.  $vtf= rxt.'GEO'.'$vtf';
  448.  vtf = doma $vtf maillage ;
  449.  elim (Maxe et vtf) epsi ;
  450.  discr=rxt.'TBT'.'discr' ;
  451.  $axe=mode Maxe 'NAVIER_STOKES' discr ;
  452.  axe = doma $axe maillage ;
  453.  
  454. un=rxt.'TIC'.'UN';
  455. vn=rxt.'TIC'.'VN';
  456. rho=rxt.'TIC'.'RHO';
  457. rvp=rxt.'TIC'.'RVP';
  458. tf=rxt.'TIC'.'TF' ;
  459. rvp=rxt.'TIC'.'RVP';
  460. td=rxt.'TIC'.'TD' ;
  461. dd=rxt.'TIC'.'DD' ;
  462. xd=rxt.'TIC'.'XD' ;
  463. *std=rxt.'TIC'.'STD' ;
  464.  
  465. Si (NON COMPLET);
  466.  
  467. an= psca vn vn (mots ux uy uz)(mots ux uy uz) ;
  468. an= an **0.5 ;
  469. mess ' Norme de VN ' (mini an) (maxi an) ;
  470. test = abs ((maxi an) - 2.888  ) /4.;
  471. Si (test '>' 1.e-2) ; erreur 5 ; Finsi ;
  472.  
  473. list rxt.'TIC'.'Tfm' ;
  474. list rxt.'TIC'.'Tdmi' ;
  475. list rxt.'TIC'.'Tdma' ;
  476. list rxt.'TIC'.'PT'  ;
  477. list rxt.'TIC'.'Qaspe';
  478. list rxt.'TIC'.'Mrest';
  479. list rxt.'TIC'.'LMAXU';
  480. *
  481.  
  482.  
  483. ltfm = prog 
  484.    120.00       120.00       120.00       121.96 
  485.    121.71       121.00       119.89       118.56
  486.    117.12       115.60       114.08 ;
  487.  
  488. lTdmi = prog
  489.    95.000       95.000       88.092       81.369
  490.    75.991       71.976       68.993       66.787
  491.    65.170       63.987       63.122 ;
  492.  
  493. lPT =prog
  494.   2.50000E+05  2.50000E+05  2.49708E+05  2.50937E+05
  495.   2.55002E+05  2.56878E+05  2.57203E+05  2.57406E+05
  496.   2.57546E+05  2.57374E+05  2.56982E+05;
  497.  
  498. lqaspe=prog
  499.    0.0000       0.0000      1.6970       1.6368
  500.    1.3683       1.0905      0.87348      0.70758
  501.    0.57940      0.48047     0.40380 ;
  502.  
  503. lmrest=prog
  504.   0.15186      0.15186      0.23671      0.31855
  505.   0.38697      0.44149      0.48516      0.52054
  506.   0.54951      0.57354      0.59373 ;
  507.  
  508. lmaxu=prog
  509.   0.0000       0.0000       1.60946E-03  4.07883E-03
  510.   6.43512E-03  8.23917E-03  9.37628E-03  9.89489E-03
  511.   9.92291E-03  9.60303E-03  9.05802E-03 ;
  512.  
  513.  
  514. tic=rxt.'TIC' ;
  515. ****************
  516.  
  517.  
  518.  ERtf=somm( abs (ltfm - tic.'Tfm') )/ 80.     ;
  519.  ERPT=somm( abs (lPT  - tic.'PT' ) ) /1.e5    ;
  520.  ERQA=somm( abs (lqaspe - tic.'Qaspe' ) ) /3. ;
  521.  ERMR=somm( abs (lmrest - tic.'Mrest' ) ) /6. ;
  522.  ERum=somm( abs (Lmaxu  - tic.'LMAXU' ) )     ;
  523.  
  524.  Mess ' ERtf=' ERtf ' ERPT=' ERPT ' ERQasp=' ERQA 'ERum=' ERum ;
  525.  Mess ' ERMR=' ERmr                                            ;
  526.  Si (ERtf '>' 1.e-2) ; erreur 5 ; Finsi ;
  527.  Si (ERPT '>' 1.e-2) ; erreur 5 ; Finsi ;
  528.  Si (ERQA '>' 1.e-2) ; erreur 5 ; Finsi ;
  529.  Si (ERMR '>' 1.e-2) ; erreur 5 ; Finsi ;
  530.  Si (ERum '>' 1.e-2) ; erreur 5 ; Finsi ;
  531.  
  532.  finsi ;
  533. ************************************************
  534.  
  535.   SI GRAPH ;
  536.  
  537.  
  538.  
  539. evuz= evol 'CHPO' (exco un uz) axe ;
  540. dess evuz 'TITRE' ' UZ axe ';
  541.  
  542. MESS 'COUCOU';
  543. opti donn 5;
  544.  
  545. evvz= evol 'CHPO' (exco vn uz) axe ;
  546. dess evvz 'TITRE' ' VZ axe ';
  547.  
  548. evr= evol 'CHPO' rho axe ;
  549. dess evr 'TITRE' ' Rho axe ';
  550.  
  551. evrv= evol 'CHPO' rvp axe ;
  552. dess evrv 'TITRE' ' Rho vap axe ';
  553.  
  554. evtf= evol 'CHPO' tf axe ;
  555. dess evtf 'TITRE' ' Tf axe ';
  556.  
  557. evtd= evol 'CHPO' td axe ;
  558. dess evtd 'TITRE' ' Td  axe ';
  559.  
  560. evxd= evol 'CHPO' xd axe ;
  561. dess evxd 'TITRE' ' Xd axe ';
  562.  
  563. evdd= evol 'CHPO' dd axe ;
  564. dess evdd 'TITRE' ' DD axe ';
  565.  
  566.  Mpco = chan (p44 et p55) QUAF ;
  567.  $pco = mode Mpco 'NAVIER_STOKES' MACRO ;
  568.  pco = doma $pco maillage ;
  569.  
  570.  
  571.  elim (vtf et pco) epsi ;
  572.  
  573. ixd=masq xd 'SUPERIEUR' 1.e-7 ;
  574. vn = ixd *  vn ;
  575. un1=redu un pco ;
  576. vn1=redu vn pco ;
  577. ung= vect un1 0.1 ux uy uz jaune ;
  578. vng= vect vn1 0.05 ux uy uz vert  ;
  579. *trace ung pco (cont pco) ;
  580. *trace vng pco (cont pco) ;
  581. *trace (ung et vng) pco (cont pco) ;
  582.  
  583. trace rvp pco TITR ' RVP ' ;
  584. trace tf pco TITR ' TF ' ;
  585. trace td pco TITR ' TD ' ;
  586. trace xd pco TITR ' XD ' ;
  587.  
  588.  
  589.  
  590.  
  591.  
  592. *******************************************************************
  593. *---------------------------------------------------------------*
  594. * Trace de l'evolution de la temperature TD  sur 1 capteurs *
  595. *---------------------------------------------------------------*
  596. TPS=rxt.TBT.TABPERSO.TEMPS;
  597. TD1=rxt.TBT.TABPERSO.CAP1.'TD';
  598.  
  599.  
  600. EV1TD=EVOL MANU TPS TD1;
  601.  
  602. *list ev1TD ;
  603.  
  604.  
  605. TAB1=TABLE;
  606. TAB1.'TITRE'=TABLE;
  607. TAB1.'TITRE'. 1= MOT 'Td z=1.0m';
  608.  
  609.  
  610. DESS (EV1TD) LEGE TAB1;
  611.  
  612. *******************************************************************
  613. *---------------------------------------------------------------*
  614. * Trace de l'evolution de la temperature du gaz TF sur les 4 capteurs *
  615. *---------------------------------------------------------------*
  616.  
  617. TF1=rxt.TBT.TABPERSO.CAP1.'TF';
  618.  
  619.  
  620. EV1TF=EVOL MANU TPS TF1;
  621.  
  622. *list ev1TF ;
  623.  
  624. TAB1=TABLE;
  625. TAB1.'TITRE'=TABLE;
  626. TAB1.'TITRE'. 1= MOT 'Tf z=1.0m';
  627.  
  628.  
  629.  
  630. DESS (EV1TF ) LEGE TAB1;
  631.  
  632. *******************************************************************
  633. *---------------------------------------------------------------*
  634. * Trace de l'evolution de la fraction volumique XD  sur les 4 capteurs *
  635. *---------------------------------------------------------------*
  636.  
  637. XD1=rxt.TBT.TABPERSO.CAP1.'XD';
  638.  
  639.  
  640. EV1XD=EVOL MANU TPS XD1;
  641.  
  642. *list ev1xd ;
  643.  
  644.  
  645. TAB1=TABLE;
  646. TAB1.'TITRE'=TABLE;
  647. TAB1.'TITRE'. 1= MOT 'Xd z=1.0m';
  648.  
  649.  
  650. DESS (EV1XD ) LEGE TAB1;
  651.  
  652.  
  653. *******************************************************************
  654. *---------------------------------------------------------------*
  655. * Trace de l'evolution de la fraction volumique DD  sur les 4 capteurs *
  656. *---------------------------------------------------------------*
  657.  
  658. DD1=rxt.TBT.TABPERSO.CAP1.'DD';
  659.  
  660. EV1DD=EVOL MANU TPS DD1;
  661.  
  662. *list ev1dd ;
  663.  
  664.  
  665. TAB1=TABLE;
  666. TAB1.'TITRE'=TABLE;
  667. TAB1.'TITRE'. 1= MOT 'Dd z=1.0m';
  668.  
  669. DESS (EV1DD) LEGE TAB1;
  670.  
  671. *******************************************************************
  672. *---------------------------------------------------------------*
  673. * Trace de l'evolution de STD  sur les 4 capteurs *
  674. *---------------------------------------------------------------*
  675.  
  676. STD1 = rxt.TBT.TABPERSO.CAP1.'STD2';
  677.  
  678. EV1STD = EVOL MANU TPS STD1       ;
  679.  
  680. *list ev1std ;
  681.  
  682.  
  683. TAB1=TABLE;
  684. TAB1.'TITRE'=TABLE;
  685. TAB1.'TITRE'. 1= MOT 'STD z=1.0m';
  686.  
  687. DESS (EV1STD) LEGE TAB1;
  688. *OPTI DONN 5 ;
  689.  
  690.  FINSI ;
  691.  
  692.   'FIN' ;
  693.  
  694.  
  695.  
  696.  
  697.  
  698.  
  699.  
  700.  
  701.  
  702.  
  703.  
  704.  
  705.  
  706.  
  707.  
  708.  
  709.  

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