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 cyli pz sphe cd phx ;
  313. domm = basm proj cyli pz plan ph0 phx phy ;
  314.  
  315. l44s= l44 proj cyli pz sphe cd phx ;
  316. l55s= l55 proj cyli pz sphe cd phx ;
  317.  
  318. p44 = regle n4 l44 l44s ;
  319. p55 = regle n4 l55 l55s ;
  320.  
  321. p0s= p0 proj cyli 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.'vtf' = aria ;
  359. rxt.'epsi' = epsi ;
  360. rxt.'pi' = (0.0 0.0 1.1) ;
  361. rxt.'GRAPH' = GRAPH ;
  362. rxt.'FRPREC' = 3 ;
  363. *
  364. rxt.'VAPEUR' = VRAI ;
  365. rxt.'ASPER' = VRAI ;
  366.  
  367. *-- Conditions initiales dans l'enceinte de test
  368. rxt.'TF0' = 120.0 ;
  369. rxt.'PT0' = 2.50e5 ;
  370. rxt.'Yvap0' = 0.01 ;
  371.  
  372. rxt.'scenasp' = table ;
  373. rxt.'scenasp'.'t' = PROG 0.000 1000.000 ;
  374.  
  375. ASPER = rxt.'ASPER' ;
  376. si ASPER ;
  377. rxt.'aspinj' = jd ;
  378. rxt.'toitf' = domf ;
  379. rxt.'rod' = 1.000e+3 ;
  380. rxt.'Cpd' = 4.180e+3 ;
  381. rxt.'scenasp'.'vzinj' = prog 0. 0. ;
  382. rxt.'scenasp'.'xdinj' = prog 0. 0. ;
  383. rxt.'scenasp'.'tdinj' = prog 9.50e+1 9.50e+1 ;
  384. rxt.'scenasp'.'ddinj' = prog 1.400e-3 1.400e-3 ;
  385. finsi ;
  386. *
  387.  
  388. rxt.'DT0' = DT0 ;
  389.  
  390. *-- On impose la viscosite turbulente (m2/s)
  391. rxt.'MODTURB'='NUTURB' ;
  392. * nut/nuf=330 ;
  393. rxt.'NUT' = 5.e-3 ;
  394.  
  395. rxt.'IMPR' =1 ;
  396.  
  397.  
  398. **********************************************
  399. *-- Definition des capteurs de mesures dans l"enceinte
  400. R00 = 0. 0. 1.0 ;
  401. R01 = 0. 0. 1.5 ;
  402. R02 = 0. 0. 2.0 ;
  403.  
  404. PERSO=TABLE ;
  405. PERSO.CAPT0 = VRAI ;
  406. PERSO.TITRE = 'Tableperso' ;
  407. PERSO.TEMPS = 'PROG' ;
  408. PERSO.FREQ0 = 1 ;
  409. PERSO.LPOINTS= R00 et R01 et R02 ;
  410. PERSO.NPOINTS= 3 ;
  411. PERSO.CAPT='MOTS' 'CAP1'
  412. 'CAP2' 'CAP3' ;
  413. PERSO.COMP= MOTS
  414. 'TD' 'UY' 'UX' 'XVAP' 'VY' 'VX' 'DD' 'XD' 'TF' 'STD2' ;
  415.  
  416.  
  417.  
  418. *************************************
  419. ** Procedure PERSO pour **
  420. ** et suivi en ligne des capteurs **
  421. *************************************
  422. rxt.'PERSO' = VRAI ;
  423. rxt.'PRCPERSO'= 'VISUP' ;
  424. rxt.'TABPERSO'= PERSO ;
  425.  
  426.  
  427. EXECRXT 0 rxt ;
  428.  
  429. xdd = 'KCHT' (rxt.'GEO'.'$vtf') 'SCAL' 'SOMMET' 10.e-4 ;
  430. rxt.'TIC'.'XD' = xdd ;
  431. rxt.'TIC'.'XDM' = xdd ;
  432. dd = 'KCHT' (rxt.'GEO'.'$vtf') 'SCAL' 'SOMMET' 1.4e-3 ;
  433. rxt.'TIC'.'DD' = dd ;
  434. rxt.'TIC'.'DDM' = xdd ;
  435. td = 'KCHT' (rxt.'GEO'.'$vtf') 'SCAL' 'SOMMET' 95.00 ;
  436. rxt.'TIC'.'TD' = td ;
  437. rxt.'TIC'.'TDM' = td ;
  438.  
  439. EXECRXT nbit rxt ;
  440.  
  441.  
  442. *opti sauv 'asptest1.save';
  443. *sauv ;
  444.  
  445. Maxe = chan axe QUAF;
  446. $vtf= rxt.'GEO'.'$vtf';
  447. vtf = doma $vtf maillage ;
  448. elim (Maxe et vtf) epsi ;
  449. discr=rxt.'TBT'.'discr' ;
  450. $axe=mode Maxe 'NAVIER_STOKES' discr ;
  451. axe = doma $axe maillage ;
  452.  
  453. un=rxt.'TIC'.'UN';
  454. vn=rxt.'TIC'.'VN';
  455. rho=rxt.'TIC'.'RHO';
  456. rvp=rxt.'TIC'.'RVP';
  457. tf=rxt.'TIC'.'TF' ;
  458. rvp=rxt.'TIC'.'RVP';
  459. td=rxt.'TIC'.'TD' ;
  460. dd=rxt.'TIC'.'DD' ;
  461. xd=rxt.'TIC'.'XD' ;
  462. *std=rxt.'TIC'.'STD' ;
  463.  
  464. Si (NON COMPLET);
  465.  
  466. an= psca vn vn (mots ux uy uz)(mots ux uy uz) ;
  467. an= an **0.5 ;
  468. mess ' Norme de VN ' (mini an) (maxi an) ;
  469. test = abs ((maxi an) - 2.888 ) /4.;
  470. Si (test '>' 1.e-2) ; erreur 5 ; Finsi ;
  471.  
  472. list rxt.'TIC'.'Tfm' ;
  473. list rxt.'TIC'.'Tdmi' ;
  474. list rxt.'TIC'.'Tdma' ;
  475. list rxt.'TIC'.'PT' ;
  476. list rxt.'TIC'.'Qaspe';
  477. list rxt.'TIC'.'Mrest';
  478. list rxt.'TIC'.'LMAXU';
  479. *
  480.  
  481.  
  482. ltfm = prog
  483. 120.00 120.00 120.00 121.96
  484. 121.71 121.00 119.89 118.56
  485. 117.12 115.60 114.08 ;
  486.  
  487. lTdmi = prog
  488. 95.000 95.000 88.092 81.369
  489. 75.991 71.976 68.993 66.787
  490. 65.170 63.987 63.122 ;
  491.  
  492. lPT =prog
  493. 2.50000E+05 2.50000E+05 2.49708E+05 2.50937E+05
  494. 2.55002E+05 2.56878E+05 2.57203E+05 2.57406E+05
  495. 2.57546E+05 2.57374E+05 2.56982E+05;
  496.  
  497. lqaspe=prog
  498. 0.0000 0.0000 1.6970 1.6368
  499. 1.3683 1.0905 0.87348 0.70758
  500. 0.57940 0.48047 0.40380 ;
  501.  
  502. lmrest=prog
  503. 0.15186 0.15186 0.23671 0.31855
  504. 0.38697 0.44149 0.48516 0.52054
  505. 0.54951 0.57354 0.59373 ;
  506.  
  507. lmaxu=prog
  508. 0.0000 0.0000 1.60946E-03 4.07883E-03
  509. 6.43512E-03 8.23917E-03 9.37628E-03 9.89489E-03
  510. 9.92291E-03 9.60303E-03 9.05802E-03 ;
  511.  
  512.  
  513. tic=rxt.'TIC' ;
  514. ****************
  515.  
  516.  
  517. ERtf=somm( abs (ltfm - tic.'Tfm') )/ 80. ;
  518. ERPT=somm( abs (lPT - tic.'PT' ) ) /1.e5 ;
  519. ERQA=somm( abs (lqaspe - tic.'Qaspe' ) ) /3. ;
  520. ERMR=somm( abs (lmrest - tic.'Mrest' ) ) /6. ;
  521. ERum=somm( abs (Lmaxu - tic.'LMAXU' ) ) ;
  522.  
  523. Mess ' ERtf=' ERtf ' ERPT=' ERPT ' ERQasp=' ERQA 'ERum=' ERum ;
  524. Mess ' ERMR=' ERmr ;
  525. Si (ERtf '>' 1.e-2) ; erreur 5 ; Finsi ;
  526. Si (ERPT '>' 1.e-2) ; erreur 5 ; Finsi ;
  527. Si (ERQA '>' 1.e-2) ; erreur 5 ; Finsi ;
  528. Si (ERMR '>' 1.e-2) ; erreur 5 ; Finsi ;
  529. Si (ERum '>' 1.e-2) ; erreur 5 ; Finsi ;
  530.  
  531. finsi ;
  532. ************************************************
  533.  
  534. SI GRAPH ;
  535.  
  536.  
  537.  
  538. evuz= evol 'CHPO' (exco un uz) axe ;
  539. dess evuz 'TITRE' ' UZ axe ';
  540.  
  541. MESS 'COUCOU';
  542. opti donn 5;
  543.  
  544. evvz= evol 'CHPO' (exco vn uz) axe ;
  545. dess evvz 'TITRE' ' VZ axe ';
  546.  
  547. evr= evol 'CHPO' rho axe ;
  548. dess evr 'TITRE' ' Rho axe ';
  549.  
  550. evrv= evol 'CHPO' rvp axe ;
  551. dess evrv 'TITRE' ' Rho vap axe ';
  552.  
  553. evtf= evol 'CHPO' tf axe ;
  554. dess evtf 'TITRE' ' Tf axe ';
  555.  
  556. evtd= evol 'CHPO' td axe ;
  557. dess evtd 'TITRE' ' Td axe ';
  558.  
  559. evxd= evol 'CHPO' xd axe ;
  560. dess evxd 'TITRE' ' Xd axe ';
  561.  
  562. evdd= evol 'CHPO' dd axe ;
  563. dess evdd 'TITRE' ' DD axe ';
  564.  
  565. Mpco = chan (p44 et p55) QUAF ;
  566. $pco = mode Mpco 'NAVIER_STOKES' MACRO ;
  567. pco = doma $pco maillage ;
  568.  
  569.  
  570. elim (vtf et pco) epsi ;
  571.  
  572. ixd=masq xd 'SUPERIEUR' 1.e-7 ;
  573. vn = ixd * vn ;
  574. un1=redu un pco ;
  575. vn1=redu vn pco ;
  576. ung= vect un1 0.1 ux uy uz jaune ;
  577. vng= vect vn1 0.05 ux uy uz vert ;
  578. *trace ung pco (cont pco) ;
  579. *trace vng pco (cont pco) ;
  580. *trace (ung et vng) pco (cont pco) ;
  581.  
  582. trace rvp pco TITR ' RVP ' ;
  583. trace tf pco TITR ' TF ' ;
  584. trace td pco TITR ' TD ' ;
  585. trace xd pco TITR ' XD ' ;
  586.  
  587.  
  588.  
  589.  
  590.  
  591. *******************************************************************
  592. *---------------------------------------------------------------*
  593. * Trace de l'evolution de la temperature TD sur 1 capteurs *
  594. *---------------------------------------------------------------*
  595. TPS=rxt.TBT.TABPERSO.TEMPS;
  596. TD1=rxt.TBT.TABPERSO.CAP1.'TD';
  597.  
  598.  
  599. EV1TD=EVOL MANU TPS TD1;
  600.  
  601. *list ev1TD ;
  602.  
  603.  
  604. TAB1=TABLE;
  605. TAB1.'TITRE'=TABLE;
  606. TAB1.'TITRE'. 1= MOT 'Td z=1.0m';
  607.  
  608.  
  609. DESS (EV1TD) LEGE TAB1;
  610.  
  611. *******************************************************************
  612. *---------------------------------------------------------------*
  613. * Trace de l'evolution de la temperature du gaz TF sur les 4 capteurs *
  614. *---------------------------------------------------------------*
  615.  
  616. TF1=rxt.TBT.TABPERSO.CAP1.'TF';
  617.  
  618.  
  619. EV1TF=EVOL MANU TPS TF1;
  620.  
  621. *list ev1TF ;
  622.  
  623. TAB1=TABLE;
  624. TAB1.'TITRE'=TABLE;
  625. TAB1.'TITRE'. 1= MOT 'Tf z=1.0m';
  626.  
  627.  
  628.  
  629. DESS (EV1TF ) LEGE TAB1;
  630.  
  631. *******************************************************************
  632. *---------------------------------------------------------------*
  633. * Trace de l'evolution de la fraction volumique XD sur les 4 capteurs *
  634. *---------------------------------------------------------------*
  635.  
  636. XD1=rxt.TBT.TABPERSO.CAP1.'XD';
  637.  
  638.  
  639. EV1XD=EVOL MANU TPS XD1;
  640.  
  641. *list ev1xd ;
  642.  
  643.  
  644. TAB1=TABLE;
  645. TAB1.'TITRE'=TABLE;
  646. TAB1.'TITRE'. 1= MOT 'Xd z=1.0m';
  647.  
  648.  
  649. DESS (EV1XD ) LEGE TAB1;
  650.  
  651.  
  652. *******************************************************************
  653. *---------------------------------------------------------------*
  654. * Trace de l'evolution de la fraction volumique DD sur les 4 capteurs *
  655. *---------------------------------------------------------------*
  656.  
  657. DD1=rxt.TBT.TABPERSO.CAP1.'DD';
  658.  
  659. EV1DD=EVOL MANU TPS DD1;
  660.  
  661. *list ev1dd ;
  662.  
  663.  
  664. TAB1=TABLE;
  665. TAB1.'TITRE'=TABLE;
  666. TAB1.'TITRE'. 1= MOT 'Dd z=1.0m';
  667.  
  668. DESS (EV1DD) LEGE TAB1;
  669.  
  670. *******************************************************************
  671. *---------------------------------------------------------------*
  672. * Trace de l'evolution de STD sur les 4 capteurs *
  673. *---------------------------------------------------------------*
  674.  
  675. STD1 = rxt.TBT.TABPERSO.CAP1.'STD2';
  676.  
  677. EV1STD = EVOL MANU TPS STD1 ;
  678.  
  679. *list ev1std ;
  680.  
  681.  
  682. TAB1=TABLE;
  683. TAB1.'TITRE'=TABLE;
  684. TAB1.'TITRE'. 1= MOT 'STD z=1.0m';
  685.  
  686. DESS (EV1STD) LEGE TAB1;
  687. *OPTI DONN 5 ;
  688.  
  689. FINSI ;
  690.  
  691. 'FIN' ;
  692.  
  693.  
  694.  
  695.  
  696.  
  697.  
  698.  
  699.  
  700.  
  701.  
  702.  
  703.  
  704.  

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