Télécharger test-aspH.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : test-aspH.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * KTEST: changement de phase sur gouttes vers un régime permanent
  5. * - pas de gravité
  6. * - pas d'injection
  7. * - option COMPLET régime permanent
  8. * - (rapport de Stage Antoine BOUQUET (voir L. Blumenfeld))
  9. * - repartition homogène au depart
  10. *$$$$ VISUP
  11. ** PROCEDURE DE VISUALISATION DES RESULTATS 3D
  12. DEBPROC VISUP tps*'FLOTTANT' rxt*'TABLE' ;
  13.  
  14. NUPADT = rxt.'TIC'.'NUPADT';
  15. LTPS = rxt.'TIC'.'LTPS';
  16.  
  17. LOG1=FAUX;
  18. 'SI' (EXIST rxt.'TBT'.'TABPERSO' 'CAPT0');
  19. LOG1=rxt.'TBT'.'TABPERSO'.'CAPT0';
  20. 'FINSI';
  21. 'SI' LOG1;
  22. TEMPS = rxt.'TBT'.'TABPERSO'.'TEMPS';
  23. NTPS1= DIME TEMPS;
  24. FREQ0 = rxt.'TBT'.'TABPERSO'.'FREQ0' ;
  25. NULOC0 = (NTPS1+1)*FREQ0 ;
  26. 'SI' (NUPADT EGA NULOC0);
  27. MESS 'Execution du captage de donnees ';
  28. CAPTP rxt;
  29. TI1 = EXTR LTPS (NULOC0+1) ;
  30. rxt.'TBT'.'TABPERSO'.'TEMPS' =
  31. rxt.'TBT'.'TABPERSO'.'TEMPS' et (PROG TI1);
  32. 'FINSI';
  33. 'FINSI';
  34.  
  35. LOG2=FAUX;
  36. 'SI' (EXIST rxt.'TBT'.'TABPERSO' 'CH3D0');
  37. LOG2=rxt.'TBT'.'TABPERSO'.'CH3D0';
  38. 'FINSI';
  39. 'SI' LOG2;
  40. TPS3D = rxt.'TBT'.'TABPERSO'.'TPS3D';
  41. NTPS2= DIME TPS3D;
  42. FREQ3D = rxt.'TBT'.'TABPERSO'.'FREQ3D';
  43. NULOC3D= (NTPS2+1)*FREQ3D ;
  44. 'SI' (NUPADT EGA NULOC3D);
  45. MESS 'Execution du trace de champ 3D suivant une coupe ';
  46. CH3D rxt;
  47. TI2 = EXTR LTPS (NULOC3D+1) ;
  48. rxt.'TBT'.'TABPERSO'.'TPS3D' =
  49. rxt.'TBT'.'TABPERSO'.'TPS3D' et (PROG TI2);
  50. 'FINSI';
  51. 'FINSI';
  52.  
  53.  
  54. FINPROC ;
  55.  
  56.  
  57. *$$$$ CAPTP
  58. ** Procedure pour capteur de Temp en des points donnes au cours du temps
  59. 'DEBPROC' CAPTP rxt*'TABLE';
  60.  
  61. L1= rxt.'TBT'.'TABPERSO'.'COMP';
  62. N1= DIME L1;
  63.  
  64. PTS =rxt.'TBT'.'TABPERSO'.'LPOINTS';
  65. M1 =rxt.'TBT'.'TABPERSO'.'NPOINTS';
  66. CAPT=rxt.'TBT'.'TABPERSO'.'CAPT' ;
  67.  
  68. vtf=rxt.'GEO'.'vtf';
  69.  
  70. TEMPS = rxt.'TBT'.'TABPERSO'.'TEMPS';
  71. LTEMPS= DIME TEMPS;
  72.  
  73.  
  74. * Initialisation des tables
  75. 'SI' (LTEMPS EGA 0 );
  76. 'REPETER' BOU1 M1;
  77. PT1=PTS 'POIN' &BOU1;
  78. CAPI= EXTR CAPT &BOU1;
  79. rxt.'TBT'.'TABPERSO'.CAPI=TABLE;
  80. 'REPETER' BOU2 N1;
  81. MOI = EXTR L1 &BOU2 ;
  82. rxt.'TBT'.'TABPERSO'.CAPI.MOI=PROG;
  83. 'FIN' BOU2;
  84. 'FIN' BOU1;
  85. 'FINSI';
  86.  
  87.  
  88.  
  89. 'REPETER' BOU1 M1;
  90.  
  91. PT1=PTS 'POIN' &BOU1;
  92. CAPI= EXTR CAPT &BOU1;
  93.  
  94. 'REPETER' BOU2 N1;
  95.  
  96. MOI = EXTR L1 &BOU2 ;
  97.  
  98. 'SI' ('EGA' MOI 'TF');
  99. TF=rxt.'TIC'.'TF';
  100. PT2=vtf 'POINT' 'PROCHE' PT1;
  101. TF1=REDU TF PT2;
  102. X1 = 'EXTR' TF1 'SCAL' PT2 ;
  103.  
  104. rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  105. rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  106. 'FINSI';
  107.  
  108. 'SI' ('EGA' MOI 'XVAP');
  109. YVAP = rxt.'TIC'.'YVAP' ;
  110. l29 = yvap*0.029 ;
  111. l18 = yvap*0.018 ;
  112. lc = 0.018 ;
  113. xvap = (l29/(lc+l29-l18)) ;
  114.  
  115. PT2=vtf 'POINT' 'PROCHE' PT1 ;
  116. Xvap1=REDU Xvap PT2 ;
  117. X1 = 'EXTR' Xvap1 'SCAL' PT2 ;
  118.  
  119. rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  120. rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  121. 'FINSI';
  122.  
  123. 'SI' ('EGA' MOI 'XD');
  124. XD=rxt.'TIC'.'XD';
  125. PT2=vtf 'POINT' 'PROCHE' PT1;
  126. XD1=REDU XD PT2;
  127. X1 = 'EXTR' XD1 'SCAL' PT2 ;
  128.  
  129. rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  130. rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  131. 'FINSI';
  132.  
  133.  
  134.  
  135. 'SI' ('EGA' MOI 'TD');
  136. TD=rxt.'TIC'.'TD';
  137. PT2=vtf 'POINT' 'PROCHE' PT1;
  138. TD1=REDU TD PT2;
  139. X1 = 'EXTR' TD1 'SCAL' PT2 ;
  140.  
  141. rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  142. rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  143. 'FINSI';
  144.  
  145.  
  146. 'SI' ('EGA' MOI 'STD2');
  147. STD2= elno $vtf rxt.'TIC'.'STD';
  148. PT2=vtf 'POINT' 'PROCHE' PT1;
  149. STD1=REDU STD2 PT2;
  150. X1 = 'EXTR' STD1 'SCAL' PT2 ;
  151.  
  152. rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  153. rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  154. 'FINSI';
  155.  
  156.  
  157.  
  158.  
  159.  
  160. 'SI' ('EGA' MOI 'DD');
  161. DD=rxt.'TIC'.'DD';
  162. PT2=vtf 'POINT' 'PROCHE' PT1;
  163. DD1=REDU DD PT2;
  164. X1 = 'EXTR' DD1 'SCAL' PT2 ;
  165.  
  166. rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  167. rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  168. 'FINSI';
  169.  
  170. 'SI' ('EGA' MOI 'UX');
  171. UN=rxt.'TIC'.'UN';
  172. UX=EXCO 'UX' UN;
  173. PT2=vtf 'POINT' 'PROCHE' PT1;
  174. UX1=REDU UX PT2;
  175. X1 = 'EXTR' UX1 'SCAL' PT2 ;
  176.  
  177. rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  178. rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  179. 'FINSI';
  180.  
  181. 'SI' ('EGA' MOI 'UY');
  182. UN=rxt.'TIC'.'UN';
  183. UY=EXCO 'UY' UN;
  184. PT2=vtf 'POINT' 'PROCHE' PT1;
  185. UY1=REDU UY PT2;
  186. X1 = 'EXTR' UY1 'SCAL' PT2 ;
  187.  
  188. rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  189. rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  190. 'FINSI';
  191.  
  192. IDIM = vale 'DIME';
  193. 'SI' ('EGA' IDIM 3);
  194. 'SI' ('EGA' MOI 'UZ');
  195. UN=rxt.'TIC'.'UN';
  196. UZ=EXCO 'UZ' UN;
  197. PT2=vtf 'POINT' 'PROCHE' PT1;
  198. UZ1=REDU UZ PT2;
  199. X1 = 'EXTR' UZ1 'SCAL' PT2 ;
  200.  
  201. rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  202. rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  203. 'FINSI';
  204. 'FINSI';
  205.  
  206. 'SI' ('EGA' MOI 'VX');
  207. VN=rxt.'TIC'.'VN';
  208. VX=EXCO 'UX' VN;
  209. PT2=vtf 'POINT' 'PROCHE' PT1;
  210. VX1=REDU VX PT2;
  211. X1 = 'EXTR' VX1 'SCAL' PT2 ;
  212.  
  213. rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  214. rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  215. 'FINSI';
  216.  
  217. 'SI' ('EGA' MOI 'VY');
  218. VN=rxt.'TIC'.'VN';
  219. VY=EXCO 'UY' VN;
  220. PT2=vtf 'POINT' 'PROCHE' PT1;
  221. VY1=REDU VY PT2;
  222. X1 = 'EXTR' VY1 'SCAL' PT2 ;
  223.  
  224. rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  225. rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  226. 'FINSI';
  227.  
  228. IDIM = vale 'DIME';
  229. 'SI' ('EGA' IDIM 3);
  230. 'SI' ('EGA' MOI 'VZ');
  231. VN=rxt.'TIC'.'VN';
  232. VZ=EXCO 'UZ' VN;
  233. PT2=vtf 'POINT' 'PROCHE' PT1;
  234. VZ1=REDU VZ PT2;
  235. X1 = 'EXTR' VZ1 'SCAL' PT2 ;
  236.  
  237. rxt.'TBT'.'TABPERSO'.CAPI.MOI =
  238. rxt.'TBT'.'TABPERSO'.CAPI.MOI et ('PROG' X1);
  239. 'FINSI';
  240. 'FINSI';
  241.  
  242. 'FIN' BOU2;
  243.  
  244. 'FIN' BOU1;
  245.  
  246.  
  247. 'FINPROC' ;
  248.  
  249. opti dime 3 elem cu20 ;
  250.  
  251. COMPLET=VRAI ;
  252. COMPLET=FAUX ;
  253. GRAPH=VRAI ;
  254. GRAPH=FAUX ;
  255. opti TRACE 'PSC' ;
  256.  
  257. Si COMPLET ;
  258. DT0 = 2.e-2 ;
  259. * nbit=20 ;
  260. nbit=200 ;
  261. n1 = 2 ; n2 = 2 ; n3 = 1 ;
  262. * n4 = -1 * 10 ;
  263. n4 = 8 ;
  264. Sinon ;
  265. * nbit=100 ;
  266. DT0 = 5.e-2 ;
  267. *DT0 = 1.e-1 ;
  268. nbit = 10 ;
  269. n1 = 1 ; n2 = 2 ; n3 = 1 ;
  270. n4 = 4 ; ;
  271. Finsi ;
  272.  
  273. ri = 0.75 ; sp = 1.000 ; hc = 4 ;
  274. epsi = 1.000e-6 ; ;
  275.  
  276. p0 = 0.000 0.000 0.000 ;
  277. px = -1000.000 0.000 0.000 ;
  278. py = 0.000 -1000.000 0.000 ;
  279. pz = 0.000 0.000 1000.000 ;
  280. cd = 0.000 0.000 -20.000 ;
  281. ph0 = 0.000 0.000 hc ;
  282. phx = ri 0.000 hc ;
  283. phy = 0.000 ri hc ;
  284.  
  285. p1 = (ri*0.6) 0.000 0.000 ;
  286. p2 = (ri*0.6) (ri*0.6) 0.000 ;
  287. p3 = 0.000 (ri*0.6) 0.000 ;
  288. p4 = ri 0.000 0.000 ;
  289. p5 = 0.000 ri 0.000 ;
  290. p6 = (ri+sp) 0.000 0.000 ;
  291. p7 = 0.000 (ri+sp) 0.000 ;
  292.  
  293. l1 = droi p0 p1 n1 ;
  294. l2 = droi p1 p2 n1 ;
  295. l3 = droi p2 p3 n1 ;
  296. l4 = droi p3 p0 n1 ;
  297. l5 = cerc p4 p0 p5 (2*n1) 'DINI' 0.4 'DFIN' 0.01 ;
  298. l6 = cerc p6 p0 p7 (2*n1) ;
  299.  
  300. basf = dall l1 l2 l3 l4 ;
  301. basf = basf et (regl (l2 et l3) l5 n2) ; elim basf epsi ;
  302. basf = basf et (syme basf droi p0 p3) ; elim basf epsi ;
  303. basf = basf et (syme basf droi p0 p1) ; elim basf epsi ;
  304.  
  305. l44 = (l4 et (p3 d n2 p5) ) ;
  306. l44 = l44 et (syme l44 droi p0 p1) ; elim l44 epsi ;
  307. l44 = l44 coul rouge ;
  308. l55 = l44 tour 90. p0 (p0 plus (0 0 100)) ;
  309.  
  310. basm = regl l5 l6 n3 ;
  311. basm = basm et (syme basm droi p0 p3) ; elim basm epsi ;
  312. basm = basm et (syme basm droi p0 p1) ; elim basm epsi ;
  313.  
  314. domf = basf proj cyli pz sphe cd phx ;
  315. domm = basm proj cyli pz plan ph0 phx phy ;
  316.  
  317. l44s= l44 proj cyli pz sphe cd phx ;
  318. l55s= l55 proj cyli pz sphe cd phx ;
  319.  
  320. p44 = regle n4 l44 l44s ;
  321. p55 = regle n4 l55 l55s ;
  322.  
  323. p0s= p0 proj cyli pz sphe cd phx ;
  324. axe = droit n4 p0 p0s 'DINI' 0.4 'DFIN' 0.05 ;
  325. elim (axe et p55) epsi ;
  326. *opti donn 5;
  327.  
  328.  
  329. aria = volu basf domf n4 'DINI' 0.4 'DFIN' 0.05 ;
  330. muro = (volu basm domm n4 'DINI' 0.4 'DFIN' 0.05) coul roug ;
  331. elim (aria et muro) epsi ;
  332.  
  333. elim (aria et p44 et p55 ) epsi ;
  334. pjg = poin basf proc (0.000 0.000 0.000) ;
  335. *pjg = poin basf proc ((ri/2.0) 0.000 0.000) ;
  336. jg = (elem basf 'APPUIE' 'LARGEMENT' pjg) coul vert ;
  337.  
  338. pjd = poin domf proc (0.000 ri hc) ;
  339. *pjd = poin domf proc ((ri/2.0) 0.000 hc) ;
  340. jd = (elem domf 'APPUIE' 'LARGEMENT' pjd) coul vert ;
  341.  
  342. mesjd = MESU jd ;
  343. 'MESS' mesjd ;
  344.  
  345.  
  346.  
  347. enva = enve aria ;
  348. enva = orie enva point (0.0 0.0 (hc/2.0)) ;
  349. envm = enve muro ;
  350. * pare = inte enva envm ;
  351.  
  352. mesaria = MESU aria ;
  353. 'MESS' mesaria ;
  354. * trac aria;
  355. * opti donn 5;
  356.  
  357.  
  358. * DT0 = 1e-1 ;
  359. rxt = table ;
  360. rxt.'vtf' = aria ;
  361. rxt.'epsi' = epsi ;
  362. rxt.'pi' = (0.0 0.0 1.1) ;
  363. rxt.'GRAPH' = GRAPH ;
  364. rxt.'FRPREC' = 3 ;
  365. *
  366. rxt.'VAPEUR' = VRAI ;
  367. rxt.'ASPER' = VRAI ;
  368.  
  369. *-- Conditions initiales dans l'enceinte de test
  370. rxt.'TF0' = 120.0 ;
  371. rxt.'PT0' = 2.50e5 ;
  372. rxt.'Yvap0' = 0.01 ;
  373.  
  374. rxt.'scenasp' = table ;
  375. rxt.'scenasp'.'t' = PROG 0.000 1000.000 ;
  376.  
  377. ASPER = rxt.'ASPER' ;
  378. si ASPER ;
  379. rxt.'aspinj' = jd ;
  380. rxt.'toitf' = domf ;
  381. rxt.'rod' = 1.000e+3 ;
  382. rxt.'Cpd' = 4.180e+3 ;
  383. rxt.'scenasp'.'vzinj' = prog 0. 0. ;
  384. rxt.'scenasp'.'xdinj' = prog 0. 0. ;
  385. rxt.'scenasp'.'tdinj' = prog 9.50e+1 9.50e+1 ;
  386. rxt.'scenasp'.'ddinj' = prog 1.400e-3 1.400e-3 ;
  387. finsi ;
  388. *
  389.  
  390. rxt.'DT0' = DT0 ;
  391.  
  392. *-- On impose la viscosite turbulente (m2/s)
  393. rxt.'MODTURB'='NUTURB' ;
  394. * nut/nuf=330 ;
  395. rxt.'NUT' = 5.e-3 ;
  396.  
  397. rxt.'IMPR' =1 ;
  398.  
  399.  
  400. **********************************************
  401. *-- Definition des capteurs de mesures dans l"enceinte
  402. R00 = 0. 0. 1.0 ;
  403. R01 = 0. 0. 1.5 ;
  404. R02 = 0. 0. 2.0 ;
  405.  
  406. PERSO=TABLE ;
  407. PERSO.CAPT0 = VRAI ;
  408. PERSO.TITRE = 'Tableperso' ;
  409. PERSO.TEMPS = 'PROG' ;
  410. PERSO.FREQ0 = 1 ;
  411. PERSO.LPOINTS= R00 et R01 et R02 ;
  412. PERSO.NPOINTS= 3 ;
  413. PERSO.CAPT='MOTS' 'CAP1'
  414. 'CAP2' 'CAP3' ;
  415. PERSO.COMP= MOTS
  416. 'TD' 'UY' 'UX' 'XVAP' 'VY' 'VX' 'DD' 'XD' 'TF' 'STD2' ;
  417.  
  418.  
  419.  
  420. *************************************
  421. ** Procedure PERSO pour **
  422. ** et suivi en ligne des capteurs **
  423. *************************************
  424. rxt.'PERSO' = VRAI ;
  425. rxt.'PRCPERSO'= 'VISUP' ;
  426. rxt.'TABPERSO'= PERSO ;
  427.  
  428.  
  429. EXECRXT 0 rxt ;
  430.  
  431. xdd = 'KCHT' (rxt.'GEO'.'$vtf') 'SCAL' 'SOMMET' 10.e-4 ;
  432. rxt.'TIC'.'XD' = xdd ;
  433. rxt.'TIC'.'XDM' = xdd ;
  434. dd = 'KCHT' (rxt.'GEO'.'$vtf') 'SCAL' 'SOMMET' 1.4e-3 ;
  435. rxt.'TIC'.'DD' = dd ;
  436. rxt.'TIC'.'DDM' = xdd ;
  437. td = 'KCHT' (rxt.'GEO'.'$vtf') 'SCAL' 'SOMMET' 95.00 ;
  438. rxt.'TIC'.'TD' = td ;
  439. rxt.'TIC'.'TDM' = td ;
  440.  
  441. EXECRXT nbit rxt ;
  442.  
  443.  
  444. *opti sauv 'asptest1.save';
  445. *sauv ;
  446.  
  447. Maxe = chan axe QUAF;
  448. $vtf= rxt.'GEO'.'$vtf';
  449. vtf = doma $vtf maillage ;
  450. elim (Maxe et vtf) epsi ;
  451. discr=rxt.'TBT'.'discr' ;
  452. $axe=mode Maxe 'NAVIER_STOKES' discr ;
  453. axe = doma $axe maillage ;
  454.  
  455. un=rxt.'TIC'.'UN';
  456. vn=rxt.'TIC'.'VN';
  457. rho=rxt.'TIC'.'RHO';
  458. rvp=rxt.'TIC'.'RVP';
  459. tf=rxt.'TIC'.'TF' ;
  460. rvp=rxt.'TIC'.'RVP';
  461. td=rxt.'TIC'.'TD' ;
  462. dd=rxt.'TIC'.'DD' ;
  463. xd=rxt.'TIC'.'XD' ;
  464. *std=rxt.'TIC'.'STD' ;
  465.  
  466. Si (NON COMPLET);
  467.  
  468. an= psca vn vn (mots ux uy uz)(mots ux uy uz) ;
  469. an= an **0.5 ;
  470. mess ' Norme de VN ' (mini an) (maxi an) ;
  471. test = abs ((maxi an) - 2.888 ) /4.;
  472. Si (test '>' 1.e-2) ; erreur 5 ; Finsi ;
  473.  
  474. list rxt.'TIC'.'Tfm' ;
  475. list rxt.'TIC'.'Tdmi' ;
  476. list rxt.'TIC'.'Tdma' ;
  477. list rxt.'TIC'.'PT' ;
  478. list rxt.'TIC'.'Qaspe';
  479. list rxt.'TIC'.'Mrest';
  480. list rxt.'TIC'.'LMAXU';
  481. *
  482.  
  483. ltfm = prog
  484. 120.00 120.00 120.00 121.96
  485. 121.97 121.55 120.52 119.20
  486. 117.74 116.21 114.65 ;
  487.  
  488. lTdmi = prog
  489. 95.000 95.000 88.092 81.369
  490. 75.991 71.978 69.003 66.807
  491. 65.197 64.019 63.160 ;
  492.  
  493. lPT =prog
  494. 2.50000E+05 2.50000E+05 2.49708E+05 2.50937E+05
  495. 2.55002E+05 2.57084E+05 2.57655E+05 2.57874E+05
  496. 2.57950E+05 2.57771E+05 2.57383E+05 ;
  497.  
  498. lqaspe=prog
  499. 0.00000E+00 0.00000E+00 1.6970 1.6368
  500. 1.3683 1.0913 .87450 .70836
  501. .58036 .48159 .40494 ;
  502.  
  503. lmrest=prog
  504. .15186 .15186 .23671 .31855
  505. .38697 .44153 .48526 .52068
  506. .54969 .57377 .59402 ;
  507.  
  508. lmaxu=prog
  509. 0.00000E+00 0.00000E+00 1.60946E-03 4.07883E-03
  510. 6.43512E-03 8.24107E-03 9.38082E-03 9.89954E-03
  511. 9.92635E-03 9.60522E-03 9.05903E-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.  

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