Télécharger lbdaliq.dgibi

Retour à la liste

Numérotation des lignes :

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

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