Télécharger la plaquette Cast3M


fr en

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