Télécharger la plaquette Cast3M


fr en

Télécharger timp_echanp.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  timp_echanp.dgibi
  2. 'OPTI' 'DIME' 2 'ELEM' 'QUA4' ;
  3.  
  4. DISCR = 'LINE' ;
  5. KPRES = 'MSOMMET' ;
  6.  
  7. EPSI  = 1.0E-5 ;
  8. NPAS  = 20 ;
  9.  
  10. *** cavité soumise à une différence de température fluide/paroi
  11. *** avec une stratification thermique initiale
  12.  
  13. L1 = 1.0 ;
  14. H1 = 2.0 ;
  15.  
  16. P1 = 0.0 0.0 ;
  17. P2 =  L1 0.0 ;
  18. P3 =  L1  H1 ;
  19. P4 = 0.0  H1 ;
  20.  
  21. NX1 =  5 ;
  22. NY1 = 10 ;                                     
  23.  
  24. D1 = 'DROIT' P1 P2 NX1 ;
  25. D2 = 'DROIT' P2 P3 NY1 ;
  26. D3 = 'DROIT' P3 P4 NX1 ;
  27. D4 = 'DROIT' P4 P1 NY1 ;
  28.  
  29. S1 = 'DALLER' D1 D2 D3 D4 ;
  30.  
  31. *********************
  32. *** partie CALCUL ***
  33. *********************
  34.  
  35. *** calcul de la façon standard de type TIMP1
  36. RXT1 = 'TABLE' ;
  37. rxt1.'VERSION'= 'V0'                                          ;
  38.  
  39. RXT1.'vtf' = S1 ;
  40. RXT1.'pi' = S1 'POIN' 'PROC' ((L1/2.0) (H1/2.0)) ;
  41. RXT1.'IMPR' = 2 ;
  42.  
  43. TF0 = 100.0 ;
  44.  
  45. RXT1.'TF0' = TF0 ;
  46. RXT1.'PT0' = 1.0E5 ;
  47.  
  48. RXT1.'MODTURB' = 'NUTURB' ;
  49. RXT1.'NUT' = 200.0*2.0E-5 ;
  50.  
  51. RXT1.'epsi'  = EPSI  ;
  52. RXT1.'DISCR' = DISCR ;
  53. RXT1.'KPRE'  = KPRES ;
  54. RXT1.'GRAPH' = FAUX  ;
  55.  
  56. RXT1 . 'TIMP' = 'TABLE' ;
  57. RXT1 . 'TIMP' .'TIMP1' = 'TABLE' ;
  58. RXT1 . 'TIMP' .'TIMP1'.'MAILLAGE' = D2 'ET' D4 ;
  59. RXT1 . 'TIMP' .'TIMP1'.'ECHAN'    = 10.0 ;
  60. RXT1 . 'TIMP' .'TIMP1'.'t'        = 'PROG'    0.0   1000.0 ;
  61. RXT1 . 'TIMP' .'TIMP1'.'TIMP'     = 'PROG'   20.0     20.0 ;
  62.  
  63. RXT1.'DT0' = 1.0 ;
  64.  
  65. *** calcul de la façon ECHANP
  66. RXT2 = 'TABLE' ;
  67. rxt2.'VERSION'= 'V0'                                          ;
  68.  
  69. RXT2.'vtf' = S1 ;
  70. RXT2.'pi' = S1 'POIN' 'PROC' ((L1/2.0) (H1/2.0)) ;
  71. RXT2.'IMPR' = 2 ;
  72.  
  73. RXT2.'TF0' = TF0 ;
  74. RXT2.'PT0' = 1.0E5 ;
  75.  
  76. RXT2.'MODTURB' = 'NUTURB' ;
  77. RXT2.'NUT' = 200.0*2.0E-5 ;
  78.  
  79. RXT2.'epsi'  = EPSI  ;
  80. RXT2.'DISCR' = DISCR ;
  81. RXT2.'KPRE'  = KPRES ;
  82. RXT2.'GRAPH' = FAUX  ;
  83.  
  84. RXT2.'ECHANP' = 'TABLE' ;
  85. RXT2.'ECHANP'.'MAILLAGE' = D2 'ET' D4 ;
  86. RXT2.'ECHANP'.'ECHAN' = 10.0 ;
  87. RXT2.'ECHANP'.'TMUR'  = 20.0 ;
  88.  
  89. RXT2.'DT0' = 1.0 ;
  90.  
  91. EXECRXT  RXT1 0 ;
  92. EXECRXT  RXT2 0 ;
  93.  
  94. EPS1 = 30.0 ;
  95. *** définition d'une stratification thermique instable
  96. ZZ1 = 'COORD' 2 RXT1.'GEO'.'vtf' ;
  97. TF_Z1 = (-2.0 * EPS1 / H1 * ZZ1) + (TF0 + EPS1) ;
  98. RXT1.'TIC'.'TF' = 'KCHT' RXT1.'GEO'.'$vtf' 'SCAL' 'SOMMET' TF_Z1 ;
  99.  
  100. ZZ2 = 'COORD' 2 RXT2.'GEO'.'vtf' ;
  101. TF_Z2 = (TF0 + EPS1) - ((2.0 * EPS1 / H1) * ZZ2) ;
  102. RXT2.'TIC'.'TF' = 'KCHT' RXT2.'GEO'.'$vtf' 'SCAL' 'SOMMET' TF_Z2 ;
  103.  
  104. *** itération en temps
  105.  
  106. EXECRXT  RXT1 NPAS ;
  107. EXECRXT  RXT2 NPAS ;
  108.  
  109. *** tests de non-régression
  110. *** la température fluide calculée par les deux péthodes doit être
  111. *** strictement la même
  112.  
  113. vtf1 = RXT1.'GEO'.'vtf' ;
  114. vtf2 = RXT2.'GEO'.'vtf' ;
  115. 'ELIM' (vtf1 'ET' vtf2) EPSI ;
  116.  
  117. TF1  = RXT1.'TIC'.'TF'  ;
  118. TF2  = RXT2.'TIC'.'TF'  ;
  119.  
  120. DIFF1 = 'ABS' (TF2 - TF1) ;
  121.  
  122. MAX1  = 'MAXI' DIFF1 ;
  123. ERROR = 0 ;
  124.  
  125. 'SI' (MAX1 '>' 1.E-10) ;
  126. ERROR = ERROR + 1 ;
  127. 'MESS' 'Erreur dans le calcul de la température fluide' ;
  128. 'MESS' 'en utilisant la méthode TIMP1 ou ECHANP' ;
  129. 'MESS' '|TF1 - TF2| = ' MAX1 ' > 1.E-10 ' ;
  130. 'FINSI' ;
  131.  
  132. 'SI' (ERROR > 0) ;
  133. 'MESS' 'Il y a des erreurs !' ;
  134. 'ERRE' 5 ;
  135. 'FINSI' ;
  136.  
  137. 'FIN' ;
  138.  
  139.  
  140.  
  141.  
  142.  
  143.  
  144.  
  145.  
  146.  

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