Télécharger test-asp2D.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : test-asp2D.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. *****************************************************************
  6. ** Jeu de données - maillage 2D
  7. ** injection de gouttes à 40°C
  8. ** dans une enceinte remplie d'air à 24°C
  9. ** modèle à 7 équations (diphasique-1 phase dispersée)
  10. **
  11. ** Date : L. Blumenfeld 29/09/05
  12. *****************************************************************
  13. 'OPTI' 'DIME' 2 ;
  14. 'OPTI' 'TRAC' 'X' ;
  15. 'OPTI' 'MODE' 'AXIS' ;
  16. 'OPTI' 'ELEM' 'QUA4' ;
  17.  
  18. *****************************************************************
  19. ** MAILLAGE ***
  20. *****************************************************************
  21. GRAPH = VRAI ;
  22. GRAPH = FAUX ;
  23. HB = 7.294 ;
  24. HA = 1.294 ;
  25. RA = 1.9 ;
  26. RJ = 0.1 ;
  27. RM = RA/2. ;
  28. HF = (HB+HA)/2. ;
  29. *................................................
  30. * Nombre d'elements injection et cotés *
  31. * en MACRO, N réel = n x 2 *
  32. *...............................................*
  33. *___________________CAS N°1_______________________
  34. nj = -4 ;
  35. nr = -20 ;
  36. nm = nr/2 ;
  37. nhb = -50 ;
  38. nha = -10 ;
  39. * nhe = -1 ;
  40. nf = nhb/2 ;
  41.  
  42. *_______________coordonées des points_______________*
  43. A0 = 0. HA ;
  44. A1 = RA HA ;
  45. AM = RM HA ;
  46. AJ = RJ HA ;
  47. C0 = 0. 0. ;
  48. CJ = RJ 0. ;
  49. C1 = RA 0. ;
  50. CM = RM 0. ;
  51. B0 = 0. HB ;
  52. BJ = RJ HB ;
  53. B1 = RA HB ;
  54. BM = RM HB ;
  55. * E0 = 0. HE ;
  56. * EJ = RJ HE ;
  57. * E1 = RA HE ;
  58. * EM = RM HE ;
  59. F0 = 0. HF ;
  60. F1 = RA HF ;
  61. FJ = RJ HF ;
  62. *______________________densité________________________*
  63. dh0 = 1. ;
  64. *dhj = 1. ;
  65. dhm1 = 1. ;
  66. dhm2 = 1. ;
  67. dh1 = 1. ;
  68. *dh1=0.2 ;
  69. dvc = 1. ;
  70. dvb = 1. ;
  71. dva1 = 1. ;
  72. dva2 = 1. ;
  73. dve1 = 1. ;
  74. dve2 = 1. ;
  75. df1 = 1. ;
  76. df2 = 1. ;
  77. dhj = 1. ;
  78. * ______________Définition des segments_______________*
  79.  
  80. CJC0 = 'DROIT' nj BJ B0 'DINI' dhj 'DFIN' dh0 ;
  81. CJCM = 'DROIT' nr CJ CM 'DINI' dhj 'DFIN' dhm1 ;
  82. CMC1 = 'DROIT' nm CM C1 'DINI' dhm2 'DFIN' dh1 ;
  83. C1A1 = 'DROIT' nha C1 A1 'DINI' dvc 'DFIN' dva1 ;
  84. A1AM = 'DROIT' nm A1 AM 'DINI' dh1 'DFIN' dhm2 ;
  85. AMAJ = 'DROIT' nr AM AJ 'DINI' dhm1 'DFIN' dhj ;
  86. AJCJ = 'DROIT' nha AJ CJ 'DINI' dva2 'DFIN' dvc ;
  87.  
  88. A1F1 = 'DROIT' nf A1 F1 'DINI' dve2 'DFIN' df1 ;
  89. F1B1 = 'DROIT' nf F1 B1 'DINI' df2 'DFIN' dvb ;
  90. B1BM = 'DROIT' nm B1 BM 'DINI' dh1 'DFIN' dhm2 ;
  91. BMBJ = 'DROIT' nr BM BJ 'DINI' dhm1 'DFIN' dhj ;
  92. BJFJ = 'DROIT' nf BJ FJ 'DINI' dvb 'DFIN' df2 ;
  93. FJAJ = 'DROIT' nf FJ AJ 'DINI' df1 'DFIN' dve2 ;
  94.  
  95. BJB0 = 'DROIT' nj BJ B0 'DINI' dhj 'DFIN' dh0 ;
  96. B0F0 = 'DROIT' nf B0 F0 'DINI' dvb 'DFIN' df2 ;
  97. F0A0 = 'DROIT' nf F0 A0 'DINI' df1 'DFIN' dve2 ;
  98. A0AJ = 'DROIT' nj A0 AJ 'DINI' dh0 'DFIN' dhj ;
  99.  
  100. CJC1 = CJCM 'ET' CMC1 ;
  101. A1AJ = A1AM 'ET' AMAJ ;
  102. A1B1 = A1F1 'ET' F1B1 ;
  103. B1BJ = B1BM 'ET' BMBJ ;
  104. A1B1 = A1F1 'ET' F1B1 ;
  105. B0A0 = B0F0 'ET' F0A0 ;
  106. A1AJ = A1AM 'ET' AMAJ ;
  107. B1B0 = B1BJ 'ET' BJB0 ;
  108. BJAJ = BJFJ 'ET' FJAJ ;
  109.  
  110. ZONE_1 = 'DALLER' CJC1 C1A1 A1AJ AJCJ PLAN ;
  111. ZONE_2 = 'DALLER' A1B1 B1BJ BJAJ ('INVERSE' A1AJ) PLAN ;
  112. ZONE_3 = 'DALLER' BJB0 B0A0 A0AJ ('INVERSE' BJAJ) PLAN ;
  113.  
  114. 'ELIMINATION' 1e-5 ZONE_1 'ET' ZONE_2;
  115. 'ELIMINATION' 1e-5 ZONE_2 'ET' ZONE_3;
  116. MT = ZONE_2 'ET' ZONE_3;
  117. *'TRACE' MT ;
  118. * opti donn 5;
  119.  
  120.  
  121. * Lignes pour conditions aux limites et post-traitement
  122.  
  123. AXE = ('INVERSE' B0A0) ;
  124. * INJ = A0AJ ;
  125. jd = ('INVERSE' BJB0) coul roug ;
  126. domf = ('INVERSE' B1B0) coul vert;
  127. BAS = A1AJ ET ('INVERSE' A0AJ) ;
  128. * BAS = INJ 'ET' AJCJ 'ET' BAS2 ;
  129. * HAUT = ('INVERSE' B1B0) ;
  130. PCON = A1B1 ;
  131. * COTE = C1A1 ;
  132. PTOT = BAS 'ET' PCON ;
  133. aria = PTOT et jd et domf ;
  134. * elim aria PTOT 1.e-5 ;
  135. * elim jd domf 1.e-5 ;
  136. * trac aria ;
  137.  
  138. *opti donn 5;
  139. * Mise en place de la procedure
  140.  
  141. RXT = 'TABLE' ;
  142. rxt.'VERSION'= 'V0' ;
  143. RXT.'vtf' = MT ;
  144.  
  145. RXT.'axe' = AXE ;
  146. *RXT.'diru1' = (0 1) ;
  147. RXT.'pi' = ( MT 'POIN' 'PROC' (0.0 2.0)) ;
  148. RXT.'IMPR' = 2 ;
  149.  
  150. RXT.'THERMP' = FAUX ;
  151.  
  152. *RXT.'TIMP1' = 'TABLE' ;
  153. *RXT.'TIMP1'.'MAILLAGE' = PCON ;
  154. *RXT.'TIMP1'.'t' = 'PROG' 0.0 1000 ;
  155. *RXT.'TIMP1'.'TIMP' = 'PROG' 120.0 120.0 ;
  156. *RXT.'TIMP1'.'ECHAN' = 25.0 ;
  157.  
  158. RXT.'VAPEUR' = VRAI ;
  159.  
  160. *RXT.'AIR' = VRAI ;
  161. RXT.'TF0' = 24.3 ;
  162. RXT.'PT0' = 1.01E5 ;
  163. *RXT.'Yvap0' = 0.4372 ;
  164. RXT.'Yvap0' = 0.0044 ;
  165.  
  166. RXT.'MODTURB' = 'NUTURB' ;
  167. RXT.'NUT' = 1.8E-2 ;
  168. *RXT.'MODTURB' = 'LMEL' ;
  169. *RXT.'LMEL' = 0.2 ;
  170. *RXT.'LMEL' = 0.04 ;
  171.  
  172. RXT.'epsi' = 1.E-5 ;
  173. RXT.'DISCR' = 'LINE' ;
  174. RXT.'KPRE' = 'MSOMMET' ;
  175. *RXT.'DISCR' = 'MACRO' ;
  176. *RXT.'KPRE' = 'CENTREP1' ;
  177. RXT.'GRAPH' = FAUX ;
  178.  
  179. RXT.'scenasp' = 'TABLE' ;
  180. RXT.'scenasp'.'t' = 'PROG' 0.0 25.5 ;
  181. RXT.'ASPER' = VRAI ;
  182.  
  183. ASPER = rxt.'ASPER' ;
  184. si ASPER ;
  185. **************description essai************************
  186. rop = 1.000e+3 ;
  187. C = 9.0e-2 ;
  188. Vi = 2.0 ;
  189. Dp = 6.0e-4 ;
  190. S = 3.14*((0.1)**2) ;
  191. C_s = C/S ;
  192. alpha = C_s/rop/Vi ;
  193. * alpha = 1.4e-04
  194. * list alpha ;
  195. * opti donn 5 ;
  196. rxt.'aspinj' = jd ;
  197. rxt.'toitf' = domf ;
  198. rxt.'rod' = rop ;
  199. rxt.'Cpd' = 4.180e+3 ;
  200. rxt.'scenasp'.'vzinj' = prog (-1.*Vi) (-1.*Vi) ;
  201. rxt.'scenasp'.'xdinj' = prog alpha alpha ;
  202. rxt.'scenasp'.'tdinj' = prog 4.000e+1 4.000e+1 ;
  203. rxt.'scenasp'.'ddinj' = prog Dp Dp ;
  204. finsi ;
  205.  
  206. *************************************
  207. ** EXECUTION **
  208. *************************************
  209. * temps caracteristiques a recalculer pour estimer le temps
  210. RXT.'DT0' = 0.1 ;
  211. NSAV = 0 ;
  212. NPAS = 5 ;
  213. execrxt NPAS RXT ;
  214. *opti donn 5;
  215. NSAV = NSAV '+' NPAS ;
  216.  
  217.  
  218.  
  219. *'OPTI' 'SAUVE' 'boite-asp.save' ;
  220. *'SAUVE' ;
  221.  
  222.  
  223. Maxe = chan axe QUAF ;
  224. Mdomf = chan domf QUAF ;
  225. $vtf= rxt.'GEO'.'$vtf' ;
  226. vtf = doma $vtf maillage ;
  227. * elim (Maxe et vtf) epsi ;
  228. * elim (Mdomf et vtf) epsi ;
  229. vtf = doma $vtf maillage ;
  230.  
  231.  
  232. un=rxt.'TIC'.'UN' ;
  233. vn=rxt.'TIC'.'VN' ;
  234. rho=rxt.'TIC'.'RHO';
  235. rvp=rxt.'TIC'.'RVP';
  236. tf=rxt.'TIC'.'TF' ;
  237. rvp=rxt.'TIC'.'RVP';
  238. td=rxt.'TIC'.'TD' ;
  239. xd=rxt.'TIC'.'XD' ;
  240. dd=rxt.'TIC'.'DD' ;
  241.  
  242. an= psca vn vn (mots ux uy uz)(mots ux uy uz) ;
  243. an= an **0.5 ;
  244. mess ' Norme de VN ' (mini an) (maxi an) ;
  245.  
  246. list rxt.'TIC'.'Tdmi' ;
  247. list rxt.'TIC'.'Qaspe';
  248. list rxt.'TIC'.'Mrest';
  249. list rxt.'TIC'.'LMAXU';
  250. *
  251.  
  252.  
  253.  
  254. lTdmi = prog
  255. 40.000 40.000 34.438 29.848 26.280 23.420 ;
  256.  
  257. lqaspe=prog
  258. 0.0000 0.0000 1.03552E-03 1.66385E-03 2.20096E-03 2.65541E-03 ;
  259.  
  260. lmrest=prog
  261. 0.35316 0.35316 0.35326 0.35343 0.35365 0.35392 ;
  262.  
  263.  
  264. * lmaxu=prog
  265. * 0.0000 0.0000 9.89819E-02 0.31752 0.51396 0.72466 ;
  266. *0.0000 0.0000 0.10029 0.33355 0.53667 0.75167 ;
  267. * Après correction Gounand 2012/12/20
  268. lmaxu=prog
  269. 0.0000 0.0000 8.34311E-02 0.26067 0.41732 0.68018;
  270.  
  271.  
  272. tic=rxt.'TIC' ;
  273. ****************
  274.  
  275.  
  276. ERti=SOMM( abs (ltdmi - tic.'Tdmi') )/ 80. ;
  277. ERQA=SOMM( abs (lqaspe - tic.'Qaspe' ) ) /3. ;
  278. ERMR=SOMM( abs (lmrest - tic.'Mrest' ) ) /6. ;
  279. ERum=SOMM( abs (Lmaxu - tic.'LMAXU' ) ) ;
  280.  
  281. Mess ' ERti=' ERti ' ERQA=' ERQA ' ERMR=' ERMR 'ERum=' ERum ;
  282.  
  283. Si (ERti '>' 1.e-2) ; erreur 5 ; Finsi ;
  284. Si (ERQA '>' 1.e-2) ; erreur 5 ; Finsi ;
  285. Si (ERMR '>' 1.e-2) ; erreur 5 ; Finsi ;
  286. Si (ERum '>' 1.e-2) ; erreur 5 ; Finsi ;
  287.  
  288. Finsi;
  289.  
  290. SI GRAPH ;
  291. uy1 = 'EXCO' 'UY' un ;
  292. evuy= evol 'CHPO' uy1 axe ;
  293. dess evuy 'TITRE' ' UY axe ';
  294.  
  295. evvy= evol 'CHPO' (exco vn uy) axe ;
  296. dess evvy 'TITRE' ' VY axe ';
  297.  
  298. evr= evol 'CHPO' rho axe ;
  299. dess evr 'TITRE' ' Rho axe ';
  300.  
  301. evrv= evol 'CHPO' rvp axe ;
  302. dess evrv 'TITRE' ' Rho vap axe ';
  303.  
  304. evtf= evol 'CHPO' tf axe ;
  305. dess evtf 'TITRE' ' Tf axe ';
  306. list evtf ;
  307.  
  308. evtd= evol 'CHPO' td axe ;
  309. dess evtd 'TITRE' ' Td axe ';
  310. list evtd ;
  311.  
  312. evxd= evol 'CHPO' xd axe ;
  313. dess evxd 'TITRE' ' Xd axe ';
  314.  
  315. evdd= evol 'CHPO' dd axe ;
  316. dess evdd 'TITRE' ' DD axe ';
  317.  
  318.  
  319. ixd=masq xd 'SUPERIEUR' 1.e-5 ;
  320. vn = ixd * vn ;
  321.  
  322. ung= vect un 1. ux uy jaune ;
  323. vng= vect vn 0.1 ux uy vert ;
  324. trace ung vtf (cont vtf);
  325. trace vng vtf ;
  326. *trace tp $vtp ;
  327. trace (ung et vng) vtf ;
  328.  
  329. *********
  330. *debit sur toutes surface
  331. Vparoi = vn * rxt.'rod' ;
  332. flimp = KOPS Xd '*' Vparoi ;
  333. nmfU = DOMA rxt.'GEO'.'$menvfU' 'NORMALEV' ;
  334. vnfU =vect nmfU 0.05 'UX' 'UY' 'ROUG' ;
  335. trace vnfU rxt.'GEO'.'menvfU' ;
  336. DepufU = 'DBIT' flimp rxt.'GEO'.'$menvfU' ;
  337. * list nmfu ;
  338. MESS 'debit sortie par les parois a t' ;
  339. list DepufU ;
  340. * masseU = DepufU * rxt.'TIC'.'LTPS' ;
  341. * MESS 'masse sortie par les parois a t' ;
  342. MESS 'temps final' ;
  343. list rxt.'TIC'.'LTPS' ;
  344.  
  345. denseau = 'KOPS' Xd '*' rop ;
  346. Diagvtf = 'DOMA' $vtf 'XXDIAGSI' ;
  347. masseeau = 'SOMT' (Diagvtf '*' denseau) ;
  348. MESS 'masse dans l enceinte a t' ;
  349. list masseeau ;
  350. trace denseau vtf TITR 'denseau' ;
  351. ev =evol chpo denseau domf ;
  352. dess ev ;
  353. ev =evol chpo denseau axe ;
  354. dess ev ;
  355. *
  356. *opti donn 5;
  357. un = rxt.'TIC'.'UN' ;
  358. uy1 = 'EXCO' 'UY' un ;
  359. vn = rxt.'TIC'.'VN' ;
  360. vy1 = 'EXCO' 'UY' vn ;
  361.  
  362.  
  363.  
  364. trace rvp vtf (cont vtf) TITR ' RVP ' ;
  365. trace tf vtf (cont vtf) TITR ' TF ' ;
  366. trace td vtf (cont vtf) TITR ' TD ' ;
  367. trace xd vtf (cont vtf) TITR ' XD ' ;
  368. *opti donn 5;
  369.  
  370.  
  371. 'FINSI' ;
  372.  
  373. 'FIN' ;
  374.  
  375.  
  376.  
  377.  
  378.  
  379.  
  380.  
  381.  
  382.  
  383.  
  384.  
  385.  
  386.  
  387.  
  388.  
  389.  
  390.  
  391.  
  392.  

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