Télécharger la plaquette Cast3M


fr en

Télécharger enc2D-therco.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  enc2D-therco.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. *  Enceinte 2D Axisymetrique
  6. *  Refroidissement d'un gaz chaud (AIR)
  7. *  Interaction GAz/paroi/temperature exterieure
  8. *  Teste thermique paroi et coefficient d'echange
  9.  
  10. 'OPTI' 'DIME' 2 'ELEM' 'CU20' ;
  11.  'OPTI' 'MODE' 'AXIS' ;
  12.  
  13.  COMPLET= VRAI                                          ;
  14.  COMPLET= FAUX                                          ;
  15.  GRAPH  = VRAI                                          ;
  16.  GRAPH  = FAUX                                          ;
  17.  opti trac 'PSC';
  18.  
  19.  'SI' COMPLET ;
  20.   nbit=200    ;
  21.   DT0 = 2.    ;
  22.   n1=3 ;
  23.   n2=5;
  24.   n3=10;
  25.   ds1= 0.02 ;
  26.   ds2= 0.2  ;
  27.  'SINON'      ;
  28.   nbit=20     ;
  29.   DT0 = 2.    ;
  30.   n1=3 ;
  31.   n2=5;
  32.   n3=10;
  33.   ds1= 0.02 ;
  34.   ds2= 0.2  ;
  35.  'FINSI'      ;
  36.  
  37.  
  38. p1=5.e-4 0 ;
  39. p2=0.5 0. ;
  40. p3=1 0 ;
  41. p4=1 2;
  42. p41=0.5 2. ;
  43. p5=5.e-4 2 ;
  44.  
  45.  
  46. bas= p1 'DROI' DINI ds1 DFIN ds2 p2 d DINI ds2 DFIN ds1 p3 ;
  47. paroid=p3 'DROI' n3 p4 ;
  48. plaf=inve (bas plus (0. 2.) );
  49. axe = p5 'DROI' n3 p1 ;
  50. elim (bas et plaf et paroid et axe) 1.e-5 ;
  51. vtp = (paroid trans 2 (0.1  0.)) coul bleu;
  52. parext = (cote 3 vtp) coul rouge ;
  53.  
  54.  
  55. mt= 'DALLER'  bas paroid plaf axe ;
  56. Si GRAPH;
  57. trace (mt et vtp);
  58. Finsi ;
  59. rxt = 'TABLE' ;
  60. rxt.'VERSION'= 'V0'                                           ;
  61. rxt.'vtf'=mt ;
  62. rxt.'axe'   = axe ;
  63.  
  64.  rxt.'TF0'    =       100.                        ;
  65.  rxt.'PT0'    =       1.0e5                       ;
  66.  
  67.  rxt.'THERMP'=VRAI ;
  68.  rxt.'THERCO'=VRAI ;
  69.  rxt.'vtp'=vtp;
  70.  rxt.'ROCP'=1.e4 ;
  71.  rxt.'LAMBDA'=40.  ;
  72.  rxt.'Tp0'     = 30. ;
  73.  rxt.'ECHAN'=50. ;
  74.  
  75.  rxt.'ECHEXT'  = VRAI ;
  76.  rxt.'HEXT'=50. ;
  77.  rxt.'TPEXT'=20. ;
  78.  rxt.'parext'=parext ;
  79.  
  80. rxt.'epsi'= 1.e-4 ;
  81. rxt.'pi'  = 0.5 0.5 ;
  82. rxt.'DT0'  = DT0     ;
  83. rxt.'RENU'  = 'RIEN' ;
  84. rxt.'DISCR'= 'QUAF';
  85. rxt.'KPRE' = 'MSOMMET';
  86. rxt.'MODTURB'  = 'NUTURB'  ;
  87. rxt.'NUT'  = 0.005  ;
  88. rxt.'GRAPH'  = GRAPH ;
  89.  
  90. EXECRXT 0    rxt ;
  91. $vtp=rxt.'GEO'.'$vtp' ;
  92. $vtf=rxt.'GEO'.'$vtf' ;
  93. *rxt.'LAMBDA'=kcht $vtf scal sommet 40. ;
  94.  
  95. EXECRXT nbit rxt ;
  96.  
  97. un=rxt.'TIC'.'UN';
  98. tf =rxt.'TIC'.'TF' ;
  99.  
  100. $vtf=rxt.'GEO'.'$vtf' ;
  101. $vtp=rxt.'GEO'.'$vtp' ;
  102. vtf=doma $vtf maillage;
  103. vtp=doma $vtp maillage;
  104.  
  105. rho=rxt.'TIC'.'RHO' ;
  106. un =rxt.'TIC'.'UN';
  107. tp =rxt.'TIC'.'TP';
  108.  
  109. ixe=inve axe ;
  110. evauz = evol chpo (exco un 'UY') ixe ;
  111. evarh = evol chpo rho ixe ;
  112. evatf = evol chpo tf  ixe ;
  113.  
  114.  
  115. Si (NON COMPLET) ;
  116. list rxt.TIC.'Tfm' ;
  117. list rxt.TIC.'PT'  ;
  118. list rxt.TIC.'LMAXU';
  119.  
  120. ltfm=Prog
  121.    100.00       79.731       76.588       71.335       67.287
  122.    63.939       60.724       57.698       54.980       52.548    
  123.    50.370       48.421       46.657       45.037       43.536
  124.    42.140       40.845       39.644       38.530       37.493    
  125.    36.528 ;   
  126.  
  127. lPT =Prog
  128.  1.00000E+05    98275.       92267.       91693.       91538.
  129.    90412.       89390.       88653.       87949.       87272.    
  130.    86679.       86162.       85693.       85262.       84864.
  131.    84492.       84145.       83823.       83523.       83246.    
  132.    82987. ;   
  133.  
  134. Lmaxu=Prog
  135.    0.0000       0.0000      0.26029      0.44149      0.51786
  136.    0.48432      0.42517     0.35634      0.29129      0.23483    
  137.    0.19800      0.17967     0.16889      0.16210      0.15651
  138.    0.15082      0.14500     0.13929      0.13407      0.12947    
  139.    0.12546 ;
  140.  
  141. tic=rxt.'TIC' ;
  142. ERtf=SOMM( abs (ltfm - tic.'Tfm') )/ 80. ;
  143. ERPT=SOMM( abs (lPT  - tic.'PT' ) ) /1.e5 ;
  144. ERum=SOMM( abs (Lmaxu  - tic.'LMAXU' ) )  ;
  145.  
  146. Mess ' ERtf=' ERtf ' ERPT=' ERPT ' ERum=' ERum ;
  147.  
  148. Si (ERtf '>' 2.e-4) ; erreur 5 ; Finsi ;
  149. Si (ERPT '>' 2.e-4) ; erreur 5 ; Finsi ;
  150. Si (ERum '>' 2.e-3) ; erreur 5 ; Finsi ;
  151.  
  152. evuz= extr evauz ordo ;
  153. list evuz ;
  154. lrr= prog
  155.   8.68197E-39  1.81018E-02  5.70173E-02  8.89906E-02  0.10936
  156.   0.12289      0.12546      0.12244      0.10940      9.93363E-02
  157.   8.06895E-02  7.14766E-02  5.53056E-02  5.39819E-02  4.18936E-02
  158.   4.67526E-02  3.46966E-02  3.43026E-02  2.18271E-02  3.34970E-03
  159.  -9.75217E-38 ;
  160.  
  161. ER=SOMM( (evuz - lrr)*(evuz - lrr) ) *0.05;
  162.  mess ' Ecart  ' er ;
  163. Si ( er > 2.e-4) ; erreur 5 ; finsi ;
  164.  
  165. Finsi ;
  166.  
  167.  
  168. Si GRAPH ;
  169.  
  170.  dess evauz 'TITRE' 'Vitesse sur l axe';
  171.  dess evarh 'TITRE' 'Densite melange sur l axe';
  172.  dess evatf 'TITRE' 'Temperature sur l axe';
  173.  
  174. ung= vect un 0.1 ux uy jaune;
  175. trace ung vtf;
  176.  trace tp  vtp 'TITRE' ' Tp ' ;
  177.  opti isov suli ;
  178. trace tf  vtf 'TITRE' ' Tf ' ;
  179. trace rho vtf 'TITRE' ' Rho' ;
  180. Finsi ;
  181.  
  182.  FIN ;
  183.  
  184.  
  185.  
  186.  
  187.  
  188.  
  189.  
  190.  
  191.  
  192.  
  193.  
  194.  
  195.  
  196.  
  197.  
  198.  
  199.  

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