Télécharger injairA.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : injairAC.dgibi
  2. *
  3. * Enceinte 2D Axisymetrique
  4. * Injection d'air en adiabatique
  5. * Comparaison Castem 2D axi / Cathare / Solution analytique
  6. 'OPTI' 'DIME' 2 'ELEM' 'QUA4' ;
  7. 'OPTI' 'MODE' 'AXIS' ;
  8. opti trace 'X';
  9. opti EPTR 10 ;
  10.  
  11. COMPLET= VRAI ;
  12. COMPLET= FAUX ;
  13. GRAPH = VRAI ;
  14. GRAPH = FAUX ;
  15.  
  16. Tmax=500.;
  17.  
  18. Si COMPLET;
  19. DTC=0.;
  20. n1=5 ;
  21. n2=10;
  22. n3=20;
  23. Sinon;
  24. DTC=2.;
  25. n1=3;
  26. n2=3;
  27. n3=3;
  28. Finsi;
  29.  
  30. *=======================================================================
  31. * Feuille inox de protection
  32. e4=(2.36-2.3588)/2.;
  33. lb4=17.;
  34. h4i=e4/lb4;
  35. * Isolant
  36. e1=(2.46-2.36)/2.;
  37. lb1=0.025;
  38. h1i=e1/lb1;
  39. * paroi inox Riva
  40. e3=0.02;
  41. lb3=18.;
  42. h3i=e3/lb3;
  43. * Echange externe
  44. hext=10.;
  45. Text=20.;
  46. * H equivalent
  47. Heq=(h4i + h1i + h3i + (1./hext))**(-1.);
  48. mess 'Heq=' Heq;
  49.  
  50. *=======================================================================
  51. *=======================================================================
  52. * Solution analytique
  53. DEBPROC SOREF Tmax*'FLOTTANT' DT0*'FLOTTANT' QM*'FLOTTANT' Vr*'FLOTTANT'
  54. Sr*'FLOTTANT' Pnm*'FLOTTANT' tetai*'FLOTTANT' INCOND*'MOT'
  55. HEXT*'FLOTTANT' TEXT*'FLOTTANT' ;
  56. * Les températures son en °C
  57.  
  58. Tps=0.;
  59. nbit = 50000;
  60.  
  61. Si('EGA' INCOND 'AIR');
  62. a1 a2 a3 a4 a5 a6 Cp=calcp tetai ;
  63. Rinc=287.1;
  64. SINON;
  65. Si('EGA' INCOND 'N2');
  66. a1 a2 a3 Cp a5 a6 a7=calcp tetai ;
  67. Rinc=296.9;
  68. SINON;
  69. mess ' Incondensable inconnu';
  70. QUITTER SOREF;
  71. FINSI;
  72. FINSI;
  73. * Rinc=292.; Valeur pour faire mieux fiter la solution analytique avec
  74. * cathare en air.
  75.  
  76. mess ' Incondensable R=' Rinc;
  77. Rhonm=Pnm/Rinc/(tetai + 273.15);
  78. a=Cp/Rinc-1.;
  79.  
  80. Lrho= prog Rhonm;
  81. Ltps= prog Tps ;
  82. Lp = prog Pnm ;
  83. LT = prog tetai;
  84.  
  85. repeter BLOC nbit;
  86. tps=tps + dt0;
  87. Si (tps > Tmax); quitter Bloc ; Finsi ;
  88. Rhon=Rhonm + (DT0*Qm/Vr);
  89. Fd=(Qm**3.)/(Rhon**2.)/(Vr**(10./3.));
  90.  
  91. A = Cp/Rinc - 1. ;
  92. Alfa = Sr*hext/Rhon/Rinc;
  93. Beta = Sr*hext*(Text + 273.15) ;
  94. AA = DT0/(Alfa*DT0+A) ;
  95.  
  96. Pn = (A*DT0/(A/DT0 + Alfa)*Pnm) + (AA/Vr*Qm*Cp*(tetai + 273.15))
  97. + (AA*Beta);
  98.  
  99. **Pn =Pnm + (DT0*Qm/Vr*(Rinc*Cp)/(Cp-Rinc)*(tetai + 273.15));
  100.  
  101. Tn =(Pn/Rhon/Rinc) - 273.15;
  102.  
  103. Si('EGA' INCOND 'AIR');
  104. a1 a2 a3 a4 a5 a6 Cp=calcp tn;
  105. Finsi;
  106. Si('EGA' INCOND 'N2');
  107. a1 a2 a3 Cp a5 a6 a7=calcp tn;
  108. FINSI;
  109.  
  110. Lrho= Lrho et (prog Rhon);
  111. Ltps= Ltps et (prog Tps) ;
  112. Lp = Lp et (prog Pn) ;
  113. LT = LT et (prog Tn) ;
  114. Rhonm = Rhon;
  115. Pnm = pn;
  116.  
  117. Fin BLOC;
  118.  
  119. evrho1 = (evol manu ltps lrho) coul rouge;
  120. evp1 = (evol manu ltps lp) coul rouge ;
  121. evt1 = (evol manu ltps lt) coul rouge ;
  122. evmair1= evrho1* Vr ;
  123. FINPROC evp1 evt1 evrho1 evmair1 ;
  124.  
  125. DT0=1. ;
  126. Qm=0.2 ;
  127. Vr=17.4;
  128. Sr=35.914;
  129. Pnm=1.e5;
  130. tetai=300. - 273.15;
  131. hext=0.; Text = 0.;
  132. *hext=Heq; Text = 20.;
  133. evp1 evt1 evrho1 evmair1 = SOREF Tmax DT0 QM Vr Sr Pnm tetai 'AIR'
  134. HEXT TEXT ;
  135.  
  136.  
  137. *=======================================================================
  138. Si GRAPH;
  139. TAB1 =TABLE ;
  140. TAB1.'TITRE'=TABLE ;
  141. tab1 . 1 = 'MOT' 'MARQ CROI REGU ';
  142. tab1 . 'TITRE' . 1 = 'Sol Ref';
  143.  
  144. tit= chai 'Pression enceinte' ;
  145. dess evp1 titr tit
  146. 'GRIL' titx 's' tity 'Pa' lege tab1;
  147.  
  148. tit= chai 'Température enceinte' ;
  149. dess evt1 titr tit
  150. 'GRIL' titx 's' tity '°C' lege tab1;
  151.  
  152. tit = chai 'Densité moyenne enceinte';
  153. dess evrho1 titr tit
  154. 'GRIL' titx 's' tity 'kg/m3' lege tab1;
  155.  
  156. tit = chai 'Masse d air';
  157. dess evmair1 titr tit
  158. 'GRIL' titx 's' tity 'kg' lege tab1;
  159. Finsi;
  160.  
  161.  
  162. epsi=1.e-5;
  163.  
  164. episo=0.044;
  165. db=0.1 ;
  166. R1=1.2 ;
  167. R1=1.2 - episo;
  168. h1=1.605;
  169. h2=2.342 - episo;
  170. hb=0.835;
  171. PA0=R1 0.;
  172. PAb=PA0 moins (0. hb);
  173. PAbm=PAb moins ((R1/2.) 0.);
  174. PA1=R1 h1;
  175. PB1=R1 ((-1.)*h1);
  176. PA1b=PA1 moins (0. hb);
  177. CA1=0. 0.275 ;
  178. CA2=(0.33*R1) H1 ;
  179. PAX2=0. h2;
  180. PA20= PAX2 tour -5. CA1;
  181. PA2= PA20 tour -15. CA1;
  182. p0=0 0 ;
  183. paxb= p0 moins (0. hb);
  184. p2=db 0. ;
  185. p2b=p2 moins (0. hb);
  186. breche = paxb 'DROI' n1 p2b ;
  187. bas= p2b 'DROI' dini (db/n1) dfin (R1/10.) PABM
  188. 'DROI' dini (R1/10.) dfin (db/n1) PAB;
  189. n2= nbel bas;
  190. plaf0 = PAX2 C n1 CA1 PA20;
  191. plaf1 = PA20 C (n2/2) CA1 PA2;
  192. plaf2 = PA1 C (n2/2) CA2 PA2;
  193.  
  194. plaf =plaf0 et plaf1 et plaf2;
  195. axeh = PAX2 'DROI' n3 paxb ;
  196. PAXM=(PAX2 plus PAXB)*0.5;
  197.  
  198. axeh = PAX2 DROI dini (h2/n3/4.) dfin (h1/n3) PAXM
  199. DROI dini (h1/n3) dfin (h2/n3/4.) PAXB;
  200. n3n= nbel axeh ;
  201. paroih=PAb d n3n PA1;
  202.  
  203. pla=plaf et paroih ;
  204. pla = pla et (syme pla droit p0 pax2);
  205. Si GRAPH;
  206. trace (pla et ca1 et ca2 et PA1);
  207. trace (plaf et paroih et axeh et breche et bas
  208. et ca1);
  209. Finsi;
  210. mth= 'DALLER' (breche 'ET' bas) paroih plaf axeh;
  211.  
  212. entb = breche plus (0. 0.);
  213. DEPL entb tour 15. P2b;
  214. pbx2= pax2 syme 'DROITE' P0 PA0;
  215. p0b= poin 1 entb ;
  216. axeb = PBX2 'DROI' n3 p0b;
  217. paroib =pab d n3 pb1;
  218. fond = plaf SYME 'DROITE' P0 PA0;
  219. elim (entb 'ET' bas et paroib et fond et axeb) epsi;
  220. mtb= 'DALLER' (entb 'ET' bas) paroib fond axeb;
  221. axe=axeb et axeh;
  222. mt = mth et mtb ;
  223. elim (mtb et bas et entb et bas) epsi;
  224. Parois = paroib et paroih et plaf et fond et entb ;
  225. parext = paroib et paroih et plaf et fond ;
  226. Si Graph;
  227. trace (elim (mt et (syme mt 'DROIT' P0 pax2)) epsi);
  228. Finsi;
  229. ************************************************************************
  230. * Petite disgression pour calculer la surface de l'enceinte
  231. Si FAUX;
  232.  
  233. la= prog;
  234. ln=prog;
  235. nn=25;
  236. repeter BLOC 40;
  237. a= rota parext nn 90. pax2 pbx2;
  238. sa=((mesu a 'SURF' ) * 4.);
  239. la = la et (prog sa);
  240. ln = ln et (prog nn);
  241. nn= nn + 25;
  242. fin bloc;
  243. evmes= evol manu ln la ;
  244. dess evmes GRIL TITR ' Surface interne de l enceinte RIVA';
  245. *opti donn 5;
  246. Finsi ;
  247. ************************************************************************
  248.  
  249. rxt = 'TABLE' ;
  250. rxt.'vtf'=mt ;
  251. rxt.'axe' = axe ;
  252. rxt.'breche'=breche ;
  253. rxt.'diru1'=0. 1. ;
  254. rxt.'epsi'= 1.e-4 ;
  255. rxt.'pi' = 0.5 0.5 ;
  256. rxt.'DT0' = (maxi (prog DT0 DTC));
  257. rxt.'DETMAT'= VRAI ;
  258. rxt.'RENU' = 'RIEN' ;
  259. rxt.'DISCR'= 'LINE';
  260. rxt.'KPRE' = 'MSOMMET';
  261. rxt.'MODTURB' = 'LMEL' ;
  262. rxt.'LMEL' = 0.01 ;
  263. rxt.'GRAPH' = GRAPH ;
  264.  
  265.  
  266. TF0 = Tetai ;
  267. rxt.'TF0' = TF0 ;
  268. rxt.'PT0' = 1.0e5 ;
  269.  
  270. Qinc = Qm ;
  271. rxt.'scenario' = table ;
  272. rxt.'scenario'.'t' = prog 0.0 1000. ;
  273. rxt.'scenario'.'qair' = prog Qinc Qinc ;
  274. rxt.'scenario'.'tinj' = prog TF0 TF0 ;
  275.  
  276. rxt.'GRAPH' = FAUX ;
  277. execrxt 0 rxt ;
  278.  
  279. Si COMPLET;
  280. rxt.'TBT'.'RTF'.'1TSCA'.'KOPT'.'CMD'=0.;
  281. Finsi;
  282.  
  283. Diag=rxt.'GEO'.'Diag';
  284. mess ' Volume de RIVA ' (somt diag);
  285. rxt.'GRAPH' = GRAPH ;
  286.  
  287. nbit = ENTI (Tmax / rxt.'DT0') ;
  288. execrxt nbit rxt ;
  289.  
  290. TAB1 =TABLE ;
  291. TAB1.'TITRE'=TABLE ;
  292. tab1 . 1 = 'MOT' 'MARQ CROI REGU ';
  293. tab1 . 2 = 'MOT' 'MARQ CARR REGU TIRC';
  294. tab1 . 3 = 'MOT' 'MARQ LOSA REGU TIRC';
  295. tab1 . 4 = 'MOT' 'MARQ TRIU REGU TIRC';
  296. tab1 . 5 = 'MOT' 'MARQ TRID REGU TIRC';
  297. tab1 . 'TITRE' . 1 = 'Sol Ref';
  298. tab1 . 'TITRE' . 2 = 'Nautilus CFD';
  299. *.................................................................
  300. evrho2= evol manu (rxt.'TIC'.'LTPS') (rxt.'TIC'. 'Rhom');
  301. evp2= evol manu (rxt.'TIC'.'LTPS') (rxt.'TIC'. 'PT');
  302. evt2= evol manu (rxt.'TIC'.'LTPS') (rxt.'TIC'. 'Tfm');
  303. evmair2=evrho2 * Vr;
  304.  
  305. Si GRAPH;
  306. tit= chai 'Pression enceinte' ;
  307. dess (evp1 et evp2) titr tit
  308. 'GRIL' titx 's' tity 'Pa' lege tab1;
  309.  
  310. tit= chai 'Température moyenne enceinte';
  311. dess (evt1 et evt2) titr tit
  312. 'GRIL' titx 's' tity '°C' lege tab1;
  313.  
  314. tit= chai 'Densité moyenne enceinte ' ;
  315. dess (evrho1 et evrho2) titr tit
  316. 'GRIL' titx 's' tity 'kg/m3' lege tab1;
  317.  
  318. tit= chai 'Masse d air ' ;
  319. dess (evmair1 et evmair2) titr tit
  320. 'GRIL' titx 's' tity 'kg/m3' lege tab1;
  321. Finsi;
  322.  
  323. devp=evp1 - evp2 abs;
  324. erp= maxi (extr devp 'ORDO');
  325. mess ' ERP =' ERP;
  326. Si (ERP '>' 5000. ) ; erreur 5 ; Finsi ;
  327.  
  328. devt=evt1 - evt2 abs;
  329. ert= maxi (extr devt 'ORDO');
  330. mess ' ERT =' ERT;
  331. Si (ERT '>' 1.2 ) ; erreur 5 ; Finsi ;
  332.  
  333.  
  334. FIN ;
  335.  
  336.  
  337.  
  338.  
  339.  

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