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.'vtf' = MT ;
  143.  
  144. RXT.'axe' = AXE ;
  145. *RXT.'diru1' = (0 1) ;
  146. RXT.'pi' = ( MT 'POIN' 'PROC' (0.0 2.0)) ;
  147. RXT.'IMPR' = 2 ;
  148.  
  149. RXT.'THERMP' = FAUX ;
  150.  
  151. *RXT.'TIMP1' = 'TABLE' ;
  152. *RXT.'TIMP1'.'MAILLAGE' = PCON ;
  153. *RXT.'TIMP1'.'t' = 'PROG' 0.0 1000 ;
  154. *RXT.'TIMP1'.'TIMP' = 'PROG' 120.0 120.0 ;
  155. *RXT.'TIMP1'.'ECHAN' = 25.0 ;
  156.  
  157. RXT.'VAPEUR' = VRAI ;
  158.  
  159. *RXT.'AIR' = VRAI ;
  160. RXT.'TF0' = 24.3 ;
  161. RXT.'PT0' = 1.01E5 ;
  162. *RXT.'Yvap0' = 0.4372 ;
  163. RXT.'Yvap0' = 0.0044 ;
  164.  
  165. RXT.'MODTURB' = 'NUTURB' ;
  166. RXT.'NUT' = 1.8E-2 ;
  167. *RXT.'MODTURB' = 'LMEL' ;
  168. *RXT.'LMEL' = 0.2 ;
  169. *RXT.'LMEL' = 0.04 ;
  170.  
  171. RXT.'epsi' = 1.E-5 ;
  172. RXT.'DISCR' = 'LINE' ;
  173. RXT.'KPRE' = 'MSOMMET' ;
  174. *RXT.'DISCR' = 'MACRO' ;
  175. *RXT.'KPRE' = 'CENTREP1' ;
  176. RXT.'GRAPH' = FAUX ;
  177.  
  178. RXT.'scenasp' = 'TABLE' ;
  179. RXT.'scenasp'.'t' = 'PROG' 0.0 25.5 ;
  180. RXT.'ASPER' = VRAI ;
  181.  
  182. ASPER = rxt.'ASPER' ;
  183. si ASPER ;
  184. **************description essai************************
  185. rop = 1.000e+3 ;
  186. C = 9.0e-2 ;
  187. Vi = 2.0 ;
  188. Dp = 6.0e-4 ;
  189. S = 3.14*((0.1)**2) ;
  190. C_s = C/S ;
  191. alpha = C_s/rop/Vi ;
  192. * alpha = 1.4e-04
  193. * list alpha ;
  194. * opti donn 5 ;
  195. rxt.'aspinj' = jd ;
  196. rxt.'toitf' = domf ;
  197. rxt.'rod' = rop ;
  198. rxt.'Cpd' = 4.180e+3 ;
  199. rxt.'scenasp'.'vzinj' = prog (-1.*Vi) (-1.*Vi) ;
  200. rxt.'scenasp'.'xdinj' = prog alpha alpha ;
  201. rxt.'scenasp'.'tdinj' = prog 4.000e+1 4.000e+1 ;
  202. rxt.'scenasp'.'ddinj' = prog Dp Dp ;
  203. finsi ;
  204.  
  205. *************************************
  206. ** EXECUTION **
  207. *************************************
  208. * temps caracteristiques a recalculer pour estimer le temps
  209. RXT.'DT0' = 0.1 ;
  210. NSAV = 0 ;
  211. NPAS = 5 ;
  212. execrxt NPAS RXT ;
  213. *opti donn 5;
  214. NSAV = NSAV '+' NPAS ;
  215.  
  216.  
  217.  
  218. *'OPTI' 'SAUVE' 'boite-asp.save' ;
  219. *'SAUVE' ;
  220.  
  221.  
  222. Maxe = chan axe QUAF ;
  223. Mdomf = chan domf QUAF ;
  224. $vtf= rxt.'GEO'.'$vtf' ;
  225. vtf = doma $vtf maillage ;
  226. * elim (Maxe et vtf) epsi ;
  227. * elim (Mdomf et vtf) epsi ;
  228. vtf = doma $vtf maillage ;
  229.  
  230.  
  231. un=rxt.'TIC'.'UN' ;
  232. vn=rxt.'TIC'.'VN' ;
  233. rho=rxt.'TIC'.'RHO';
  234. rvp=rxt.'TIC'.'RVP';
  235. tf=rxt.'TIC'.'TF' ;
  236. rvp=rxt.'TIC'.'RVP';
  237. td=rxt.'TIC'.'TD' ;
  238. xd=rxt.'TIC'.'XD' ;
  239. dd=rxt.'TIC'.'DD' ;
  240.  
  241. an= psca vn vn (mots ux uy uz)(mots ux uy uz) ;
  242. an= an **0.5 ;
  243. mess ' Norme de VN ' (mini an) (maxi an) ;
  244.  
  245. list rxt.'TIC'.'Tdmi' ;
  246. list rxt.'TIC'.'Qaspe';
  247. list rxt.'TIC'.'Mrest';
  248. list rxt.'TIC'.'LMAXU';
  249. *
  250.  
  251.  
  252.  
  253. lTdmi = prog
  254. 40.000 40.000 34.438 29.848 26.280 23.420 ;
  255.  
  256. lqaspe=prog
  257. 0.0000 0.0000 1.03552E-03 1.66385E-03 2.20096E-03 2.65541E-03 ;
  258.  
  259. lmrest=prog
  260. 0.35316 0.35316 0.35326 0.35343 0.35365 0.35392 ;
  261.  
  262.  
  263. * lmaxu=prog
  264. * 0.0000 0.0000 9.89819E-02 0.31752 0.51396 0.72466 ;
  265. *0.0000 0.0000 0.10029 0.33355 0.53667 0.75167 ;
  266. * Après correction Gounand 2012/12/20
  267. lmaxu=prog
  268. 0.0000 0.0000 8.34311E-02 0.26067 0.41732 0.68018;
  269.  
  270.  
  271. tic=rxt.'TIC' ;
  272. ****************
  273.  
  274.  
  275. ERti=SOMM( abs (ltdmi - tic.'Tdmi') )/ 80. ;
  276. ERQA=SOMM( abs (lqaspe - tic.'Qaspe' ) ) /3. ;
  277. ERMR=SOMM( abs (lmrest - tic.'Mrest' ) ) /6. ;
  278. ERum=SOMM( abs (Lmaxu - tic.'LMAXU' ) ) ;
  279.  
  280. Mess ' ERti=' ERti ' ERQA=' ERQA ' ERMR=' ERMR 'ERum=' ERum ;
  281.  
  282. Si (ERti '>' 1.e-2) ; erreur 5 ; Finsi ;
  283. Si (ERQA '>' 1.e-2) ; erreur 5 ; Finsi ;
  284. Si (ERMR '>' 1.e-2) ; erreur 5 ; Finsi ;
  285. Si (ERum '>' 1.e-2) ; erreur 5 ; Finsi ;
  286.  
  287. Finsi;
  288.  
  289. SI GRAPH ;
  290. uy1 = 'EXCO' 'UY' un ;
  291. evuy= evol 'CHPO' uy1 axe ;
  292. dess evuy 'TITRE' ' UY axe ';
  293.  
  294. evvy= evol 'CHPO' (exco vn uy) axe ;
  295. dess evvy 'TITRE' ' VY axe ';
  296.  
  297. evr= evol 'CHPO' rho axe ;
  298. dess evr 'TITRE' ' Rho axe ';
  299.  
  300. evrv= evol 'CHPO' rvp axe ;
  301. dess evrv 'TITRE' ' Rho vap axe ';
  302.  
  303. evtf= evol 'CHPO' tf axe ;
  304. dess evtf 'TITRE' ' Tf axe ';
  305. list evtf ;
  306.  
  307. evtd= evol 'CHPO' td axe ;
  308. dess evtd 'TITRE' ' Td axe ';
  309. list evtd ;
  310.  
  311. evxd= evol 'CHPO' xd axe ;
  312. dess evxd 'TITRE' ' Xd axe ';
  313.  
  314. evdd= evol 'CHPO' dd axe ;
  315. dess evdd 'TITRE' ' DD axe ';
  316.  
  317.  
  318. ixd=masq xd 'SUPERIEUR' 1.e-5 ;
  319. vn = ixd * vn ;
  320.  
  321. ung= vect un 1. ux uy jaune ;
  322. vng= vect vn 0.1 ux uy vert ;
  323. trace ung vtf (cont vtf);
  324. trace vng vtf ;
  325. *trace tp $vtp ;
  326. trace (ung et vng) vtf ;
  327.  
  328. *********
  329. *debit sur toutes surface
  330. Vparoi = vn * rxt.'rod' ;
  331. flimp = KOPS Xd '*' Vparoi ;
  332. nmfU = DOMA rxt.'GEO'.'$menvfU' 'NORMALEV' ;
  333. vnfU =vect nmfU 0.05 'UX' 'UY' 'ROUG' ;
  334. trace vnfU rxt.'GEO'.'menvfU' ;
  335. DepufU = 'DBIT' flimp rxt.'GEO'.'$menvfU' ;
  336. * list nmfu ;
  337. MESS 'debit sortie par les parois a t' ;
  338. list DepufU ;
  339. * masseU = DepufU * rxt.'TIC'.'LTPS' ;
  340. * MESS 'masse sortie par les parois a t' ;
  341. MESS 'temps final' ;
  342. list rxt.'TIC'.'LTPS' ;
  343.  
  344. denseau = 'KOPS' Xd '*' rop ;
  345. Diagvtf = 'DOMA' $vtf 'XXDIAGSI' ;
  346. masseeau = 'SOMT' (Diagvtf '*' denseau) ;
  347. MESS 'masse dans l enceinte a t' ;
  348. list masseeau ;
  349. trace denseau vtf TITR 'denseau' ;
  350. ev =evol chpo denseau domf ;
  351. dess ev ;
  352. ev =evol chpo denseau axe ;
  353. dess ev ;
  354. *
  355. *opti donn 5;
  356. un = rxt.'TIC'.'UN' ;
  357. uy1 = 'EXCO' 'UY' un ;
  358. vn = rxt.'TIC'.'VN' ;
  359. vy1 = 'EXCO' 'UY' vn ;
  360.  
  361.  
  362.  
  363. trace rvp vtf (cont vtf) TITR ' RVP ' ;
  364. trace tf vtf (cont vtf) TITR ' TF ' ;
  365. trace td vtf (cont vtf) TITR ' TD ' ;
  366. trace xd vtf (cont vtf) TITR ' XD ' ;
  367. *opti donn 5;
  368.  
  369.  
  370. 'FINSI' ;
  371.  
  372. 'FIN' ;
  373.  
  374.  
  375.  
  376.  
  377.  
  378.  
  379.  
  380.  
  381.  
  382.  
  383.  
  384.  
  385.  
  386.  
  387.  
  388.  
  389.  

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