Télécharger la plaquette Cast3M


fr en

Télécharger pq1-lref.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  pq1-lref.dgibi
  2. *---------------------------------------------------------
  3. * Cas test de non régression basé sur pq1
  4. *---------------------------------------------------------
  5. *     Pressurisation d'une enceinte type Phébus
  6. *
  7. * On ajoute LREF pour prise en compte convection forcée
  8. * Pour lref = 0.  on doit retrouver pq1
  9. * Pour lref = 0.1 cela doit bouger
  10. *
  11. *                                 revisite Mars     2001
  12. *--------------------------------------------------------
  13. COMPLET = FAUX ;
  14. *COMPLET = VRAI ;
  15.  GRAPH =  VRAI                                          ;
  16.  GRAPH =  FAUX                                          ;
  17.  'OPTI' 'TRAC' 'X'                                ;
  18. *'OPTI' 'TRAC' 'PSC'                              ;
  19.  
  20.  'SI' COMPLET ;
  21.  
  22.  nbit=100;
  23.  DT0 = 1. ;
  24.  n1   =  1          ; n2   =  4          ; n3   =  4          ;
  25.  n4   =  8          ; nn   =  2          ;
  26.  
  27.  'SINON' ;
  28.  
  29.  nbit= 5 ;
  30.  DT0 = 10. ;
  31.  n1   =  1          ; n2   =  2          ; n3   =  4          ;
  32.  n4   =  4          ; nn   =  1          ;
  33.  
  34.  'FINSI' ;
  35.  
  36. *--------------------------------------------------------
  37. * Definition du maillage de l'enceinte cylindrique
  38. *
  39.  'OPTI' 'DIME' 3 'ELEM' 'CU20'                                ;
  40.  
  41.  ri   = 1.052      ; sp   =  0.10     ; hc   = 4.163      ;
  42.  
  43.  
  44.  
  45.  epsi =  1.000e-2   ;                                         ;
  46.  
  47.  p0   =     0.000             0.000                0.000      ;
  48.  px   = -1000.000             0.000                0.000      ;
  49.  py   =     0.000         -1000.000                0.000      ;
  50.  pz   =     0.000             0.000             1000.000      ;
  51.  cd   =     0.000             0.000              -20.000      ;
  52.  ph0  =     0.000             0.000                hc         ;
  53.  phx  =     ri                0.000                hc         ;
  54.  phy  =     0.000             ri                   hc         ;
  55.  
  56.  fg1 = 0.25  ;
  57.  fg2 = fg1 * (2.0 ** 0.5) / 2. ;
  58.  
  59.  p1   =  (ri*fg1)             0.000                0.000      ;
  60.  p2   =  (ri*fg2)             (ri*fg2)             0.000      ;
  61.  p3   =  0.000                (ri*fg1)             0.000      ;
  62.  p4   =   ri                  0.000                0.000      ;
  63.  p5   =  0.000                 ri                  0.000      ;
  64.  p6   =  (ri+sp)              0.000                0.000      ;
  65.  p7   =  0.000                (ri+sp)              0.000      ;
  66.  
  67. * Hauteur de l'enceinte
  68.   h1 = 4.163                                                    ;
  69. * Vecteur de translation
  70.   v1 = 0. 0. h1                                               ;
  71.  
  72.  l1   =  'DROI' p0 p1 n1                                        ;
  73.  l2   =  'DROI' p1 p2 n1                                        ;
  74.  l3   =  'DROI' p2 p3 n1                                        ;
  75.  l4   =  'DROI' p3 p0 n1                                        ;
  76.  l5   =  'CERC' p4 p0 p5 (2*n1)                                 ;
  77.  l6   =  'CERC' p6 p0 p7 (2*n1)                                 ;
  78.  
  79.  basf0=  'DALL' l1 l2 l3 l4 'PLAN'                              ;
  80.  basf1=('REGL' (l2 'ET' l3) l5 n2) ;
  81.  l44= cote 2 basf1;
  82.  ax4= (inve l4) et l44 ;
  83.  l11= cote 4 basf1;
  84.  ax1= l11 et (inve l1) ;
  85.  
  86.  basf =  basf0 'ET' ('REGL' (l2 'ET' l3) l5 n2) ;
  87.  'ELIM' basf epsi ;
  88.  basf =  basf 'ET' ('SYME'  basf 'DROI' p0 p3) ;
  89.  ax11 =  ('SYME'  ax1  'DROI' p0 p3) 'ET' (inve ax1) ;
  90.  'ELIM' basf epsi ;
  91.  basf =  basf 'ET' ('SYME'  basf 'DROI' p0 p1) ;
  92.  ax44 =  (inve ax4)  'ET' ('SYME'  ax4  'DROI' p0 p4) ;
  93.  'ELIM' basf epsi ;
  94.  
  95.  basm =  'REGL' l5 l6 n3                                        ;
  96.  basm =  basm 'ET' ('SYME'  basm 'DROI' p0 p3) ;
  97.  'ELIM' basm epsi ;
  98.  basm =  basm 'ET' ('SYME'  basm 'DROI' p0 p1) ;
  99.  'ELIM' basm epsi ;
  100.  
  101. * Creation du volume
  102.  
  103.  dx = ri / 2. ;
  104.  nz1 = ('ENTIER' ( h1 / dx ))*nn ;
  105.  
  106.  bas = basf 'VOLU' nz1 'TRAN' v1 ;
  107.  mbas = basm 'VOLU' nz1 'TRAN' v1 ;
  108.  mbas = mbas 'COUL' 'ROUG' ;
  109.  plan1 = ax11 'TRAN' nz1 v1 ;
  110.  plan4 = ax44 'TRAN' nz1 v1 ;
  111.  'ELIM' (bas et plan1 et plan4) epsi ;
  112.  
  113.  mt = bas ;
  114.  wall = mbas ;
  115.  elim (mt et wall) epsi ;
  116.  
  117. * Localisation d'une brèche éventuelle au bas de l'enceinte
  118.  
  119.  pjg  =  'POIN' basf 'PROC' (0.000    0.000    0.000)             ;
  120.  jg   = ('ELEM' basf 'APPUIE' 'LARGEMENT' pjg) 'COUL' 'VERT'      ;
  121.  
  122. *--------------------------------------------------------------------
  123. * Fin de la définition du maillage
  124. *--------------------------------------------------------------------
  125.  
  126.  
  127. *--------------------------------------------------------------------
  128. * Début de l'initialisation de la procédure ENCEINTE : table RXT
  129. *--------------------------------------------------------------------
  130.  
  131.  rxt             = 'TABLE' ;
  132.  rxt . 'VERSION' = 'V0' ;
  133.  
  134. *-- Nom du volume fluide
  135.  rxt . 'vtf'     = mt ;
  136.  rxt . 'epsi'    = epsi ;
  137.  
  138. *-- Definition des murs de l'enceinte : ici un seul mur
  139. *-- en ACIER dont on traitera la thermique dans l'épaisseur
  140. *-- et que l'on initialise a 40oC
  141.  
  142. *-- On definit d'abord le matériau ACIER avec sa conductivite
  143. *-- thermique LAMBDA (W/m/K) et le produit ro*Cp (J/m3/K)
  144.  rxt . 'THERMP' = VRAI                                        ;
  145.  rxt . 'vtp'    = wall                                        ;
  146.  rxt . 'LAMBDA' = 15.                                         ;
  147.  rxt . 'ROCP'   = 3.9E6                                       ;
  148.  rxt . 'Tp0'    = 60.                                         ;
  149.  rxt . 'ECHAN'  = 10.                                         ;
  150.  
  151. *-- On definit un point interne au maillage pour imposer la valeur de
  152. *-- la pression
  153.  rxt . 'pi'   =  (0.0 0.0 0.5)         ;
  154.  
  155. *-- On indique que le calcul comporte de la vapeur d'eau
  156.  rxt . 'VAPEUR' = VRAI                                        ;
  157.  rxt . 'Yvap0'  = 0.0023                      ;
  158.  
  159. *-- Conditions initiales dans l'enceinte de test
  160.  rxt . 'TF0'    = 40.0                        ;
  161.  rxt . 'PT0'    = 1.00000e5                   ;
  162.  
  163. *-- Definition du scenario à la breche
  164.  
  165.  rxt . 'Breches' = 'TABLE' ;
  166.  rxt . 'Breches' . 'A' = 'TABLE' ;
  167.  rxt . 'Breches' . 'A' . 'Maillage' =  jg ;
  168.  rxt . 'Breches' . 'A' . 'diru'     =  (0. 0. 1.) ;
  169.  rxt . 'Breches' . 'A' . 'scenario' = 'TABLE' ;
  170.  rxt . 'Breches' . 'A' . 'scenario' . 't'    =  PROG 0.0  1000. ;
  171.  rxt . 'Breches' . 'A' . 'scenario' . 'qair' =  PROG 0.0  0.0   ;
  172.  rxt . 'Breches' . 'A' . 'scenario' . 'qeau' =  PROG 0.05 0.05  ;
  173.  rxt . 'Breches' . 'A' . 'scenario' . 'tinj' =  PROG 150. 150.  ;
  174.  
  175. *-- On impose la viscosite turbulente (m2/s)
  176.  rxt . 'MODTURB' = 'NUTURB' ;
  177.  rxt . 'NUT'     =  1.e-2   ;
  178.  
  179. *-- On impose le pas de temps (s)
  180.  rxt . 'DT0'    =  DT0                                          ;
  181.  
  182. *-- On active le recalcul automatique du préconditionnement
  183. *-- toutes les 5 itérations
  184.  rxt . 'FRPREC' = 5    ;
  185.  rxt . 'DETMAT' = VRAI ;
  186. *rxt.'TRTF'   = FAUX                                       ;
  187.  
  188.  
  189. Si VRAI;
  190. *-- On lance le calcul sur 20 itérations d'une seconde
  191.  rxt . 'GRAPH' = GRAPH ;
  192.  
  193.  EXECRXT  2    rxt                              ;
  194.  cont1 = 'ENVE' rxt . 'GEO' . 'vtf' ;
  195.  rxt . 'TIC' . 'LREF' = 'MANU' 'CHPO' cont1 0. 'SCAL' ;
  196.  EXECRXT  (nbit - 2) rxt                        ;
  197.  
  198. list rxt.'TIC'.'Tfm' ;
  199. list rxt.'TIC'.'PT'  ;
  200. list rxt.'TIC'.'Qc'  ;
  201. list rxt.'TIC'.'LMAXU';
  202.  
  203. ltfm = 'PROG' 
  204. 40.000   65.434   73.712   81.158   86.809   90.721  ;
  205. lPT = 'PROG' 
  206. 1.00000E+05 1.06009E+05 1.21199E+05 1.29592E+05 1.36461E+05 1.43987E+05 ;
  207. Lqc = 'PROG' 
  208. 0.0000   0.0000   0.0000  4.06133E-04  2.65247E-03  4.88852E-03 ;
  209. Lmaxu = 'PROG' 
  210. 0.0000  0.81320   2.0813   2.6901   2.2350   2.5277 ;
  211.  
  212.  
  213. tic=rxt.'TIC' ;
  214. ERtf=somm( abs (ltfm - tic.'Tfm') )/ 80. ;
  215. ERPT=somm( abs (lPT  - tic.'PT' ) ) /1.e5 ;
  216. ERQc=somm( abs (lqc  - tic.'Qc' ) )  ;
  217. ERum=somm( abs (Lmaxu  - tic.'LMAXU' ) )/2. ;
  218.  
  219. Mess 'ERtf=' ERtf 'ERPT=' ERPT 'ERQc=' ERQc 'ERum=' ERum ;
  220.  
  221. ierr = 0 ;
  222. Si (ERtf '>' 1.e-4) ; ierr = ierr '+' 1 ; Finsi ;
  223. Si (ERPT '>' 1.e-3) ; ierr = ierr '+' 1 ; Finsi ;
  224. Si (ERQc '>' 1.e-4) ; ierr = ierr '+' 1 ; Finsi ;
  225. Si (ERum '>' 1.e-2) ; ierr = ierr '+' 1 ; Finsi ;
  226.  
  227. ******************************************************************
  228.  
  229.  rxt = enlev rxt 'GEO';
  230.  rxt = enlev rxt 'TBT';
  231.  rxt = enlev rxt 'TIC';
  232.  
  233.  EXECRXT  2    rxt                              ;
  234.  cont1 = 'ENVE' rxt . 'GEO' . 'vtf' ;
  235.  rxt . 'TIC' . 'LREF' = 'MANU' 'CHPO' cont1 1. 'SCAL' ;
  236.  EXECRXT  (nbit - 2) rxt                        ;
  237.  
  238. list rxt.'TIC'.'Tfm' ;
  239. list rxt.'TIC'.'PT'  ;
  240. list rxt.'TIC'.'Qc'  ;
  241. list rxt.'TIC'.'LMAXU';
  242.  
  243. ltfm = 'PROG' 
  244. *40.000    65.434    73.712    81.158    86.809    90.721  ;
  245.  40.000    65.434    73.712    81.158    86.787    90.532  ;
  246. lPT = 'PROG' 
  247. *1.00000E+05 1.06009E+05 1.21199E+05 1.29592E+05 1.36461E+05 1.43987E+05;
  248.  1.00000E+05 1.06009E+05 1.21199E+05 1.29592E+05 1.36453E+05 1.43899E+05;
  249. Lqc = 'PROG' 
  250. *0.0000    0.0000    0.0000   4.06133E-04  2.65247E-03 4.88852E-03 ;
  251.  0.0000    0.0000    0.0000   4.98153E-04  3.47975E-03 5.32965E-03 ;
  252. Lmaxu = 'PROG' 
  253. *0.0000   0.81320    2.0813    2.6901    2.2350    2.5277 ;
  254.  0.0000   0.81320    2.0813    2.6901    2.2350    2.5277 ;
  255.  
  256.  
  257. tic=rxt.'TIC' ;
  258. ERtf=somm( abs (ltfm - tic.'Tfm') )/ 80. ;
  259. ERPT=somm( abs (lPT  - tic.'PT' ) ) /1.e5 ;
  260. ERQc=somm( abs (lqc  - tic.'Qc' ) )  ;
  261. ERum=somm( abs (Lmaxu  - tic.'LMAXU' ) )/2. ;
  262.  
  263. Mess 'ERtf=' ERtf 'ERPT=' ERPT 'ERQc=' ERQc 'ERum=' ERum ;
  264.  
  265. ierr = 0 ;
  266. Si (ERtf '>' 1.e-4) ; ierr = ierr '+' 1 ; Finsi ;
  267. Si (ERPT '>' 1.e-3) ; ierr = ierr '+' 1 ; Finsi ;
  268. Si (ERQc '>' 1.e-4) ; ierr = ierr '+' 1 ; Finsi ;
  269. Si (ERum '>' 1.e-2) ; ierr = ierr '+' 1 ; Finsi ;
  270.  
  271. ******************************************************************
  272.  
  273.  rxt.'vtf'    =   (chan 'MACRO' mt)                     ;
  274.  rxt.'vtp'    =   (chan 'MACRO' wall)                   ;
  275.  rxt . 'Breches' . 'A' . 'Maillage' = (chan 'MACRO' jg) ;
  276.  
  277.  rxt = enlev rxt 'GEO';
  278.  rxt = enlev rxt 'TBT';
  279.  rxt = enlev rxt 'TIC';
  280.  
  281.  EXECRXT  2    rxt                              ;
  282.  cont1 = 'ENVE' rxt . 'GEO' . 'vtf' ;
  283.  rxt . 'TIC' . 'LREF' = 'MANU' 'CHPO' cont1 1. 'SCAL' ;
  284.  EXECRXT  (nbit - 2) rxt                        ;
  285.  
  286. list rxt.'TIC'.'Tfm' ;
  287. list rxt.'TIC'.'PT'  ;
  288. list rxt.'TIC'.'Qc'  ;
  289. list rxt.'TIC'.'LMAXU';
  290.  
  291. ltfm  = 'PROG' 
  292. *40.000 65.434    73.716    81.163    86.907    90.942 ; 
  293.  40.000 65.434    73.716    81.163    86.890    90.769 ;
  294. lPT   = 'PROG' 
  295. *1.00000E+05 1.06006E+05 1.21166E+05 1.29578E+05 1.36479E+05 1.44098E+05;
  296.  1.00000E+05 1.06006E+05 1.21166E+05 1.29578E+05 1.36473E+05 1.44021E+05;
  297. Lqc   = 'PROG'
  298. *0.0000 0.0000    0.0000   3.11548E-04 2.41832E-03 4.44576E-03 ;
  299.  0.0000 0.0000    0.0000   3.81417E-04 3.18170E-03 4.84069E-03 ;
  300. Lmaxu = 'PROG'
  301. *0.0000 0.81320    2.0840    2.6866    2.2313    2.5160 ;
  302.  0.0000 0.81320    2.0840    2.6866    2.2313    2.5160 ;
  303.  
  304.  
  305. tic=rxt.'TIC' ;
  306. ERtf=somm( abs (ltfm - tic.'Tfm') )/ 80. ;
  307. ERPT=somm( abs (lPT  - tic.'PT' ) ) /1.e5 ;
  308. ERQc=somm( abs (lqc  - tic.'Qc' ) )  ;
  309. ERum=somm( abs (Lmaxu  - tic.'LMAXU' ) )/2. ;
  310.  
  311. Mess 'ERtf=' ERtf 'ERPT=' ERPT 'ERQc=' ERQc 'ERum=' ERum ;
  312.  
  313. Si (ERtf '>' 1.e-4) ; ierr = ierr '+' 1 ; Finsi ;
  314. Si (ERPT '>' 1.e-3) ; ierr = ierr '+' 1 ; Finsi ;
  315. Si (ERQc '>' 1.e-4) ; ierr = ierr '+' 1 ; Finsi ;
  316. Si (ERum '>' 1.e-2) ; ierr = ierr '+' 1 ; Finsi ;
  317.  
  318.  
  319. ******************************************************************
  320.  
  321.  rxt.'vtf'    =   (chan 'QUAF' mt)                         ;
  322.  rxt.'vtp'    =   (chan 'QUAF' wall)                        ;
  323.  rxt . 'Breches' . 'A' . 'Maillage' = (chan 'QUAF' jg) ;
  324.  rxt = enlev rxt 'GEO';
  325.  rxt = enlev rxt 'TBT';
  326.  rxt = enlev rxt 'TIC';
  327.  
  328.  EXECRXT  2    rxt                              ;
  329.  cont1 = 'ENVE' rxt . 'GEO' . 'vtf' ;
  330.  rxt . 'TIC' . 'LREF' = 'MANU' 'CHPO' cont1 1. 'SCAL' ;
  331.  EXECRXT  (nbit - 2) rxt                        ;
  332.  
  333. list rxt.'TIC'.'Tfm' ;
  334. list rxt.'TIC'.'PT'  ;
  335. list rxt.'TIC'.'Qc'  ;
  336. list rxt.'TIC'.'LMAXU';
  337.  
  338. ltfm  = 'PROG' 
  339. *40.000    65.434    73.712    81.158    86.809    90.721  ;
  340.  40.000    65.434    73.712    81.158    86.787    90.532  ;
  341. lPT   = 'PROG' 
  342. *1.00000E+05 1.06009E+05 1.21199E+05 1.29592E+05 1.36461E+05 1.43987E+05;
  343.  1.00000E+05 1.06009E+05 1.21199E+05 1.29592E+05 1.36453E+05 1.43899E+05;
  344. Lqc   = 'PROG'
  345. *0.0000    0.0000    0.0000   4.06133E-04 2.65247E-03 4.88852E-03;
  346.  0.0000    0.0000    0.0000   4.98153E-04 3.47975E-03 5.32965E-03;
  347. Lmaxu = 'PROG'
  348. *0.0000   0.81320    2.0813    2.6901    2.2350    2.5277  ;
  349.  0.0000   0.81320    2.0813    2.6901    2.2350    2.5277  ;
  350.  
  351. tic=rxt.'TIC' ;
  352. ERtf=somm( abs (ltfm - tic.'Tfm') )/ 80. ;
  353. ERPT=somm( abs (lPT  - tic.'PT' ) ) /1.e5 ;
  354. ERQc=somm( abs (lqc  - tic.'Qc' ) )  ;
  355. ERum=somm( abs (Lmaxu  - tic.'LMAXU' ) )/2. ;
  356.  
  357. Mess 'ERtf=' ERtf 'ERPT=' ERPT 'ERQc=' ERQc 'ERum=' ERum ;
  358.  
  359. Si (ERtf '>' 1.e-4) ; ierr = ierr '+' 1 ; Finsi ;
  360. Si (ERPT '>' 1.e-3) ; ierr = ierr '+' 1 ; Finsi ;
  361. Si (ERQc '>' 1.e-4) ; ierr = ierr '+' 1 ; Finsi ;
  362. Si (ERum '>' 1.e-2) ; ierr = ierr '+' 1 ; Finsi ;
  363.  
  364. finsi;
  365.  
  366. ******************************************************************
  367.  rxt.'vtf'    =   (chan 'LINE' mt)                         ;
  368.  rxt.'vtp'    =   (chan 'LINE' wall)                        ;
  369.  rxt . 'Breches' . 'A' . 'Maillage' = (chan 'LINE' jg) ;
  370.  rxt = enlev rxt 'GEO';
  371.  rxt = enlev rxt 'TBT';
  372.  rxt = enlev rxt 'TIC';
  373.  
  374.  EXECRXT  2    rxt                              ;
  375.  cont1 = 'ENVE' rxt . 'GEO' . 'vtf' ;
  376.  rxt . 'TIC' . 'LREF' = 'MANU' 'CHPO' cont1 1. 'SCAL' ;
  377.  EXECRXT  (nbit - 2) rxt                        ;
  378.  
  379. list rxt.'TIC'.'Tfm' ;
  380. list rxt.'TIC'.'PT'  ;
  381. list rxt.'TIC'.'Qc'  ;
  382. list rxt.'TIC'.'LMAXU';
  383.  
  384. ltfm  = 'PROG' 
  385. *40.000    67.363    75.733    83.257    88.699    92.375  ;
  386.  40.000    67.363    75.733    83.257    88.644    92.093  ;
  387. lPT   = 'PROG' 
  388. *1.00000E+05 1.06452E+05 1.22785E+05 1.31710E+05 1.39047E+05 1.46818E+05;
  389.  1.00000E+05 1.06452E+05 1.22785E+05 1.31710E+05 1.39028E+05 1.46663E+05;
  390. Lqc   = 'PROG'
  391. *0.0000    0.0000    0.0000   1.04818E-03 3.74210E-03 6.39247E-03 ;
  392.  0.0000    0.0000    0.0000   1.26761E-03 4.89033E-03 6.99583E-03 ;
  393. Lmaxu = 'PROG'
  394. *0.0000   0.81320    2.0604    2.5919    2.1446    2.4526  ;
  395.  0.0000   0.81320    2.0604    2.5919    2.1446    2.4526  ;
  396.  
  397. tic=rxt.'TIC' ;
  398. ERtf=somm( abs (ltfm - tic.'Tfm') )/ 80. ;
  399. ERPT=somm( abs (lPT  - tic.'PT' ) ) /1.e5 ;
  400. ERQc=somm( abs (lqc  - tic.'Qc' ) )  ;
  401. ERum=somm( abs (Lmaxu  - tic.'LMAXU' ) )/2. ;
  402.  
  403. Mess 'ERtf=' ERtf 'ERPT=' ERPT 'ERQc=' ERQc 'ERum=' ERum ;
  404.  
  405. Si (ERtf '>' 1.e-4) ; ierr = ierr '+' 1 ; Finsi ;
  406. Si (ERPT '>' 1.e-3) ; ierr = ierr '+' 1 ; Finsi ;
  407. Si (ERQc '>' 1.e-4) ; ierr = ierr '+' 1 ; Finsi ;
  408. Si (ERum '>' 1.e-2) ; ierr = ierr '+' 1 ; Finsi ;
  409.  
  410.  
  411.  
  412. Si GRAPH ;
  413.  
  414. $vtf=rxt.'GEO'.'$vtf' ;
  415. mt  =doma $vtf maillage ;
  416. Mpl1=chan 'QUAF' plan1 ;
  417. Mpl4=chan 'QUAF' plan4 ;
  418. ELIM (mt et Mpl1 et Mpl4) epsi ;
  419. $mpl1= mode Mpl1 'NAVIER_STOKES' MACRO ;
  420. $mpl4= mode Mpl4 'NAVIER_STOKES' MACRO ;
  421. plan1= doma $mpl1 maillage ;
  422. plan4= doma $mpl4 maillage ;
  423. plan=plan1 et plan4 ;
  424.  
  425. paroif = rxt.'GEO'.'paroif';
  426. rho=rxt.'TIC'.'RHO';
  427. rvp=rxt.'TIC'.'RVP';
  428. tf=rxt.'TIC'.'TF';
  429. un=rxt.'TIC'.'UN' ;
  430. un1=redu un plan ;
  431. ung= vect un1 1. ux uy uz jaune ;
  432. trace ung plan ;
  433. opti isov suli ;
  434. trace rho plan 'TITRE'' Rho' ;
  435. trace rvp plan 'TITRE'' Rvp' ;
  436. trace tf  plan 'TITRE'' Tf ' ;
  437.  
  438. trace rho paroif  'TITRE'' Rho' ;
  439. trace rvp paroif  'TITRE'' Rvp' ;
  440. trace Tf  paroif  'TITRE'' Tf ' ;
  441. Fcond = tic.'Fcondw';
  442. trace Fcond paroif  'TITRE' ' Flux de condensation Kg / m**2 ' ;
  443.  
  444.  axe  = p0 d nz1 (p0 plus v1) ;
  445.  axe = chan axe 'QUAF' ;
  446.  elim (axe et mt) epsi ;
  447.  
  448. evr= evol 'CHPO' rho axe ;
  449. dess evr 'TITRE' ' Rho axe ';
  450.  
  451. evt = evol 'CHPO' Tf  axe ;
  452. dess evt 'TITRE' ' Température axe ';
  453.  
  454. evvp= evol 'CHPO' rvp axe ;
  455. dess evvp 'TITRE' ' rvp axe ';
  456.  
  457.  
  458.  
  459. 'FINSI' ;
  460.  
  461. 'SI' (IERR '>' 0) ;
  462.   'MESS' 'Il y a des problemes !!!' ;
  463.   'ERRE' 5 ;
  464. 'SINO';
  465.   'MESS' 'Tout s est bien passe!' ;
  466. 'FINS' ;
  467.  
  468.  
  469.  
  470.  'FIN' ;
  471.  
  472.  
  473.  
  474.  
  475.  
  476.  
  477.  
  478.  
  479.  
  480.  
  481.  
  482.  
  483.  
  484.  
  485.  
  486.  
  487.  
  488.  
  489.  

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