Télécharger la plaquette Cast3M


fr en

Télécharger lbdaliq.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : lbdaliq.dgibi
  2. *
  3. *
  4. * Conductivité thermique de l'eau
  5. * Valeurs de référence : tables VDI
  6. * Comparaisons de différentes formules et tests de 'VARI' 'LBDALIQ'
  7. *
  8. *
  9. 'OPTI' 'DIME' 2 'ELEM' 'QUA4' 'TRAC' 'PSC' 'EPTR' 10 'DENS' 1. ;
  10. graph = faux ;
  11. *
  12. *- Calcul de la conductivité thermique de l'eau (formule Cathare)
  13. *
  14. 'DEBP' lambdak TW0*'LISTREEL' PW0*'LISTREEL' ;
  15. *
  16. UN  = TW0 / TW0 ;
  17. PS1 = 'VARI' 'PSATT' TW0 ;
  18. HL1 = 'VARI' 'HLS'   PW0 TW0 ;
  19. *
  20. a0 = 0.573738622   ;
  21. a1 = 0.2536103551  ;
  22. a2 = -0.145468269   ;
  23. a3 = 0.01387472485 ;
  24. *
  25. x  = 1.7196904D-6 * HL1 ;
  26. L1 = (a0 * un) + (a1 * x) + (a2 * x * x) + (a3 * x * x * x) ;
  27. *
  28. 'FINP' L1 ;
  29. *
  30. *- Calcul de la conductivité thermique de l'eau (formule VDI)
  31. *
  32. 'DEBP' lambdav TW0*'LISTREEL' PW0*'LISTREEL' ;
  33. *
  34. UN  = TW0 / TW0 ;
  35. HL1 = 'VARI' 'HLS'   PW0 TW0 ;
  36. *
  37. a0 = -0.92247   ;
  38. a1 =  6.728934102  ;
  39. a2 =-10.11230521 ;
  40. a3 =  6.996953832 ;
  41. a4 = -2.316062510 ;
  42. *
  43. b0 =-0.20954276 ;
  44. b1 = 1.320227345 ;
  45. b2 =-2.485904388 ;
  46. b3 = 1.517081933 ;
  47. *
  48. c0 = 0.08104183147 ;
  49. c1 =-0.4513858027 ;
  50. c2 = 0.8057261332 ;
  51. c3 =-0.4668315566 ;
  52. *
  53. PC0 = 221.2D5 * UN ;
  54. TC0 = (374.15 + 273.15) * UN ;
  55. BETA = PW0 / PC0 ;
  56. TETA = TW0 / TC0 ;
  57. PS1 = 'VARI' 'PSATT' TW0 ;
  58. BETAS = PS1 / PC0 ;
  59. *
  60. L1 = (a0*un) + (a1 * TETA) + (a2 * TETA * TETA)
  61.    + (a3 * TETA * TETA * TETA) + (a4 * TETA * TETA * TETA * TETA) ;
  62. L2 = (b0*un) + (b1 * TETA) + (b2 * TETA * TETA)
  63.    + (b3 * TETA * TETA * TETA) ;
  64. L3 = (c0*un) + (c1 * TETA) + (c2 * TETA * TETA)
  65.    + (c3 * TETA * TETA * TETA) ;
  66. L4 = BETA - BETAS ;
  67. L5 = L1 + (L4*L2) + (L4*L4*L3) ; 
  68. *
  69. 'FINP' L5 ;
  70. *
  71. * Température et pression 
  72. *
  73. TW1 = ('PROG' 20. 40. 60. 80. 90.) + 273.15 ;
  74. UN  = TW1 / TW1 ;
  75. PW1 = 1.D5 * UN ;
  76. LREF1 = 'PROG' 0.604 0.632 0.654 0.670 0.676 ;
  77. L10 = LAMBDAK TW1 PW1 ;
  78. L11 = LAMBDAV TW1 PW1 ;
  79. *
  80. TW2 = ('PROG' 10. 30. 50. 100. 130. 150. 170.) + 273.15 ;
  81. UN  = TW2 / TW2 ;
  82. PW2 = 10.D5 * UN ;
  83. LREF2 = 'PROG' 0.588 0.619 0.644 0.681 0.688 0.687 0.681 ;
  84. L20 = LAMBDAK TW2 PW2 ;
  85. L21 = LAMBDAV TW2 PW2 ;
  86. *
  87. TW3 = ('PROG' 10. 50. 100. 150. 200. 220.) + 273.15 ;
  88. UN  = TW3 / TW3 ;
  89. PW3 = 25.D5 * UN ;
  90. LREF3 = 'PROG' 0.589 0.645 0.682 0.688 0.665 0.649 ;
  91. L30 = LAMBDAK TW3 PW3 ;
  92. L31 = LAMBDAV TW3 PW3 ;
  93. *
  94. * Test de VARI avec T sous forme de FLOTTANT
  95. 'DEBP' testf p*'LISTREEL' t*'LISTREEL' ;
  96.    nb     = 'DIME' t ;
  97.    listlb = 'PROG' ;
  98.    'REPE' blo1 nb ;
  99.       valp   = 'EXTR' p &blo1 ;
  100.       valt   = 'EXTR' t &blo1 ;
  101.       vallb  = 'VARI' 'LBDALIQ' valp valt ; 
  102.       listlb = listlb 'ET' ('PROG' vallb) ;
  103.    'FIN' blo1 ;
  104. 'FINP' listlb ;
  105. test_f1 = testf pw1 tw1 ;
  106. test_f2 = testf pw2 tw2 ;
  107. test_f3 = testf pw3 tw3 ;
  108. *
  109. * Test de VARI avec T sous forme de LISTREEL
  110. test_l1 = vari 'LBDALIQ' pw1 tw1 ;
  111. test_l2 = vari 'LBDALIQ' pw2 tw2 ;
  112. test_l3 = vari 'LBDALIQ' pw3 tw3 ;
  113. *
  114. * Test de VARI avec T sous forme de CHPO
  115. 'DEBP' testc p*'CHPOINT' t*'CHPOINT' mail1*'MAILLAGE' ;
  116.    test1  = 'VARI' 'LBDALIQ' p t ;
  117.    listlb = 'EXTR' test1 'VALE' mail1 ;
  118. 'FINP' listlb ;
  119. P1 = 1. 0. ; P2 = 2. 0. ;
  120. GEO1 = P1 'DROI' 4 P2 ; PT1 = 'CHAN' 'POI1' GEO1 ;
  121. GEO2 = P1 'DROI' 6 P2 ; PT2 = 'CHAN' 'POI1' GEO2 ;
  122. GEO3 = P1 'DROI' 5 P2 ; PT3 = 'CHAN' 'POI1' GEO3 ;
  123. T1 = 'MANU' 'CHPO' PT1 'SCAL' TW1 ; 
  124. T2 = 'MANU' 'CHPO' PT2 'SCAL' TW2 ; 
  125. T3 = 'MANU' 'CHPO' PT3 'SCAL' TW3 ; 
  126. P1 = 'MANU' 'CHPO' PT1 'SCAL' PW1 ; 
  127. P2 = 'MANU' 'CHPO' PT2 'SCAL' PW2 ; 
  128. P3 = 'MANU' 'CHPO' PT3 'SCAL' PW3 ; 
  129. test_c1 = testc p1 t1 pt1 ;
  130. test_c2 = testc p2 t2 pt2 ;
  131. test_c3 = testc p3 t3 pt3 ;
  132. *
  133. * Comparaisons mise en oeuvre Esope/Gibi pour Cathare 
  134. d10 = 'EVOL'        'MANU' tw1 (l10 - test_L1 / l10 * 100.) ;
  135. d20 = 'EVOL' 'VERT' 'MANU' tw2 (l20 - test_l2 / l20 * 100.) ;
  136. d30 = 'EVOL' 'ROUG' 'MANU' tw3 (l30 - test_l3 / l30 * 100.) ;
  137. *
  138. *-------------------------------------
  139. *- Sorties Graphiques
  140. *-------------------------------------
  141. *
  142. 'SI' graph ;
  143. *
  144. ev10   = 'EVOL'        'MANU' tw1 l10 ;
  145. ev20   = 'EVOL' 'VERT' 'MANU' tw2 l20 ;
  146. ev30   = 'EVOL' 'ROUG' 'MANU' tw3 l30 ;
  147. ev11   = 'EVOL'        'MANU' tw1 l11 ;
  148. ev21   = 'EVOL' 'VERT' 'MANU' tw2 l21 ;
  149. ev31   = 'EVOL' 'ROUG' 'MANU' tw3 l31 ;
  150. evr1   = 'EVOL'        'MANU' tw1 lref1 ;
  151. evr2   = 'EVOL' 'VERT' 'MANU' tw2 lref2 ;
  152. evr3   = 'EVOL' 'ROUG' 'MANU' tw3 lref3 ;
  153. rela10 = 'EVOL'        'MANU' tw1 (l10 - lref1 / lref1 * 100.) ;
  154. rela20 = 'EVOL' 'VERT' 'MANU' tw2 (l20 - lref2 / lref2 * 100.) ;
  155. rela30 = 'EVOL' 'ROUG' 'MANU' tw3 (l30 - lref3 / lref3 * 100.) ;
  156. rela11 = 'EVOL'        'MANU' tw1 (l11 - lref1 / lref1 * 100.) ;
  157. rela21 = 'EVOL' 'VERT' 'MANU' tw2 (l21 - lref2 / lref2 * 100.) ;
  158. rela31 = 'EVOL' 'ROUG' 'MANU' tw3 (l31 - lref3 / lref3 * 100.) ;
  159. *
  160. TAB1 = TABLE ;
  161. TAB1 . 'TITRE' = TABLE ;
  162. TAB1 .  1 = ' MARQ PLEIN ROND NOLI' ;
  163. TAB1 .  2 = ' MARQ PLEIN CARR NOLI' ;
  164. TAB1 .  3 = ' MARQ PLEIN LOSA NOLI' ;
  165. TAB1 .  4 = ' MARQ ROND NOLI' ;
  166. TAB1 .  5 = ' MARQ CARR NOLI' ;
  167. TAB1 .  6 = ' MARQ LOSA NOLI' ;
  168. *
  169. TIT1 = 'Thermal conductivity of water (Cathare)' ;
  170. TIT2 = ' C' ;
  171. TIT3 = 'W/m2/K' ;
  172. TAB1 . 'TITRE' .  1 = 'MOT' ' 1 bar' ;
  173. TAB1 . 'TITRE' .  2 = 'MOT' '10 bar' ;
  174. TAB1 . 'TITRE' .  3 = 'MOT' '25 bar' ;
  175. TAB1 . 'TITRE' .  4 = 'MOT' 'ref  1' ;
  176. TAB1 . 'TITRE' .  5 = 'MOT' 'ref 10' ;
  177. TAB1 . 'TITRE' .  6 = 'MOT' 'ref 25' ;
  178. 'DESS' (ev10 'ET' ev20 'ET' ev30 'ET' evr1 'ET' evr2 'ET' evr3)
  179. 'MIMA' 'TITR' TIT1 'TITX' TIT2 'TITY' TIT3 'LEGE' TAB1 ;
  180. *
  181. TIT1 = 'Thermal conductivity of water (VDI)' ;
  182. TIT2 = ' C' ;
  183. TIT3 = 'W/m2/K' ;
  184. TAB1 . 'TITRE' .  1 = 'MOT' ' 1 bar' ;
  185. TAB1 . 'TITRE' .  2 = 'MOT' '10 bar' ;
  186. TAB1 . 'TITRE' .  3 = 'MOT' '25 bar' ;
  187. TAB1 . 'TITRE' .  4 = 'MOT' 'ref  1' ;
  188. TAB1 . 'TITRE' .  5 = 'MOT' 'ref 10' ;
  189. TAB1 . 'TITRE' .  6 = 'MOT' 'ref 25' ;
  190. 'DESS' (ev11 'ET' ev21 'ET' ev31 'ET' evr1 'ET' evr2 'ET' evr3)
  191. 'MIMA' 'TITR' TIT1 'TITX' TIT2 'TITY' TIT3 'LEGE' TAB1 ;
  192. *
  193. TIT1 = 'Thermal conductivity of water. Relative error (x100).' ;
  194. TIT2 = ' C' ;
  195. TIT3 = ' %' ;
  196. TAB1 . 'TITRE' .  1 = 'MOT' 'Cat.  (1 bar)' ;
  197. TAB1 . 'TITRE' .  2 = 'MOT' 'Cat. (10 bar)' ;
  198. TAB1 . 'TITRE' .  3 = 'MOT' 'Cat. (25 bar)' ;
  199. TAB1 . 'TITRE' .  4 = 'MOT' 'VDI.  (1 bar)' ;
  200. TAB1 . 'TITRE' .  5 = 'MOT' 'VDI. (10 bar)' ;
  201. TAB1 . 'TITRE' .  6 = 'MOT' 'VDI. (25 bar)' ;
  202. 'DESS' (rela10 'ET' rela20 'ET' rela30 'ET'
  203.         rela11 'ET' rela21 'ET' rela31)
  204. 'MIMA' 'TITR' TIT1 'TITX' TIT2 'TITY' TIT3 'LEGE' TAB1 ;
  205. *
  206. TIT1 = 'Thermal conductivity of water. Relative error (Cath., x100).' ;
  207. TIT2 = ' C' ;
  208. TIT3 = ' %' ;
  209. TAB1 . 'TITRE' .  1 = 'MOT' 'Cat.  (1 bar)' ;
  210. TAB1 . 'TITRE' .  2 = 'MOT' 'Cat. (10 bar)' ;
  211. TAB1 . 'TITRE' .  3 = 'MOT' 'Cat. (25 bar)' ;
  212. 'DESS' (rela10 'ET' rela20 'ET' rela30)
  213. 'MIMA' 'TITR' TIT1 'TITX' TIT2 'TITY' TIT3 'LEGE' TAB1 ;
  214. *
  215. TIT1 = 'Thermal conductivity of water. Relative error (VDI, x100).' ;
  216. TIT2 = ' C' ;
  217. TIT3 = ' %' ;
  218. TAB1 . 'TITRE' .  1 = 'MOT' 'VDI.  (1 bar)' ;
  219. TAB1 . 'TITRE' .  2 = 'MOT' 'VDI. (10 bar)' ;
  220. TAB1 . 'TITRE' .  3 = 'MOT' 'VDI. (25 bar)' ;
  221. 'DESS' (rela11 'ET' rela21 'ET' rela31)
  222. 'MIMA' 'TITR' TIT1 'TITX' TIT2 'TITY' TIT3 'LEGE' TAB1 ;
  223. *
  224. TIT1 =
  225. 'Esope/Gibi relative error for Cathare formula (x100)' ;
  226. TIT2 = ' C' ;
  227. TIT3 = ' %' ;
  228. TAB1 . 'TITRE' .  1 = 'MOT' 'Esope/Gibi (1 bar)' ;
  229. TAB1 . 'TITRE' .  2 = 'MOT' 'Esope/Gibi (10 bar)' ;
  230. TAB1 . 'TITRE' .  3 = 'MOT' 'Esope/Gibi (25 bar)' ;
  231. 'DESS' (d10 'ET' d20 'ET' d30)
  232. 'MIMA' 'TITR' TIT1 'TITX' TIT2 'TITY' TIT3 'LEGE' TAB1 ;
  233. 'FINSI' ;
  234. *
  235. *----------------------------------------------------------------------
  236. *- Tests de bon fonctionnement
  237. *- (Comparaisons des différents types et de la mise en oeuvre en gibi)
  238. *----------------------------------------------------------------------
  239. *
  240. diff1 = 'MAXI' (test_f1 - test_l1 / l10) 'ABS' ;
  241. diff2 = 'MAXI' (test_f2 - test_l2 / l20) 'ABS' ;
  242. diff3 = 'MAXI' (test_f3 - test_l3 / l30) 'ABS' ;
  243. diff4 = 'MAXI' (test_f1 - test_c1 / l10) 'ABS' ;
  244. diff5 = 'MAXI' (test_f2 - test_c2 / l20) 'ABS' ;
  245. diff6 = 'MAXI' (test_f3 - test_c3 / l30) 'ABS' ;
  246. diff7 = 'MAXI' (test_f1 - l10 / l10)     'ABS' ;
  247. diff8 = 'MAXI' (test_f2 - l20 / l20)     'ABS' ;
  248. diff9 = 'MAXI' (test_f3 - l30 / l30)     'ABS' ;
  249. didi  = diff1 + diff2 + diff3 + diff4 + diff5
  250.       + diff6 + diff7 + diff8 + diff9 ;
  251. 'MESS' 'Sommes du max des différences relatives ' didi ;
  252. 'SI' (didi > 1.D-14) ;
  253.    'ERRE' 5 ;
  254. 'FINSI' ;
  255. *
  256. 'FIN' ;
  257.  
  258.  
  259.  
  260.  
  261.  
  262.  
  263.  
  264.  
  265.  

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