Télécharger prim_errord.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : prim_errord.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ************************************************************************
  5. **** APPROCHE VF "Cell-Centred Formulation" pour la solution des ****
  6. **** Equations d'Euler pour un gaz parfait ****
  7. **** OPERATEUR PRIM ****
  8. **** Test erreur ordre composantes vitesse-fractions massiques ****
  9. **** ****
  10. **** A. BECCANTINI DRN/DMT/SEMT/TTMF FEVRIER 2000 ****
  11. ************************************************************************
  12. *
  13. 'OPTION' 'DIME' 2 ;
  14. 'OPTION' 'ELEM' QUA4 ;
  15. 'OPTION' 'ECHO' 0 ;
  16. 'OPTION' 'TRAC' 'X';
  17. *
  18. *
  19. **** GRAPH
  20. *
  21. *
  22. GRAPH = FAUX ;
  23. * GRAPH = VRAI ;
  24. *
  25.  
  26.  
  27. *) Degre de polynome = 4
  28.  
  29. ***********************
  30. **** LA TABLE PGAZ ****
  31. ***********************
  32.  
  33. PGAZ = 'TABLE' ;
  34.  
  35. *
  36. **** Ordre des polynoms
  37. *
  38.  
  39. PGAZ . 'NORD' = 4 ;
  40.  
  41. *
  42. **** Especes qui sont dans les equations d'Euler
  43. *
  44.  
  45. PGAZ . 'ESPEULE' = 'MOTS' 'H2 ' 'N2 ' ;
  46.  
  47. *
  48. **** Espece qui n'y est pas
  49. *
  50.  
  51.  
  52. PGAZ . 'ESPNEULE' = 'O2 ' ;
  53.  
  54. *
  55.  
  56. PGAZ . 'H2 ' = 'TABLE' ;
  57. PGAZ . 'O2 ' = 'TABLE' ;
  58. PGAZ . 'N2 ' = 'TABLE' ;
  59.  
  60. *
  61. **** R (J/Kg/K)
  62. *
  63.  
  64. PGAZ . 'H2 ' . 'R' = 4130.0 ;
  65. PGAZ . 'O2 ' . 'R' = 259.8 ;
  66. PGAZ . 'N2 ' . 'R' = 296.8 ;
  67.  
  68.  
  69. *
  70. **** Regressions polynomials
  71. *
  72.  
  73. PGAZ . 'H2 ' . 'A' = 'PROG' 9834.91866 0.54273926 0.000862203836
  74. -2.37281455E-07 1.84701105E-11 ;
  75. PGAZ . 'O2 ' . 'A' = 'PROG' 575.012333 0.350522002 -0.000128294865
  76. 2.33636971E-08 -1.53304905E-12;
  77. PGAZ . 'N2 ' . 'A' = 'PROG' 652.940766 0.288239099 -7.80442298E-05
  78. 8.78233606E-09 -3.05514485E-13 ;
  79.  
  80. *
  81. **** "Enthalpies" (ou energies) de formations a OK (J/Kg)
  82. *
  83.  
  84. PGAZ . 'H2 ' . 'H0K' = -4.195D6 ;
  85. PGAZ . 'O2 ' . 'H0K' = -2.634D5 ;
  86. PGAZ . 'N2 ' . 'H0K' = -2.953D5 ;
  87.  
  88. *
  89. *** Fin PGAZ
  90. *
  91.  
  92.  
  93. ***************************
  94. ***** DOMAINE SPATIAL ****
  95. ***************************
  96.  
  97.  
  98. A1 = 0.0D0 0.0D0;
  99. A2 = 2.0D0 0.0D0;
  100. A3 = 3.0D0 0.0D0;
  101. A4 = 4.0D0 1.0D0;
  102.  
  103. L12 = A1 'DROIT' 21 A2;
  104. L23 = A2 'DROIT' 32 A3;
  105. L34 = A3 'DROIT' 43 A4;
  106. L41 = A4 'DROIT' 51 A1;
  107.  
  108. LDOM1 = L12 'ET' L23 'ET' L34 'ET' L41 ;
  109. DOM1 = 'SURFACE' LDOM1 'PLANE';
  110.  
  111. 'SI' GRAPH ;
  112. 'TRACER' DOM1 'TITRE' 'Domaine' ;
  113. 'FINSI' ;
  114.  
  115. $DOM1 = 'DOMA' DOM1 ;
  116.  
  117. MAIL = $DOM1 . 'CENTRE' ;
  118.  
  119. *
  120. **** Energie, densite, temperature, vitesse, fractions massiques
  121. *
  122.  
  123. RN = 'BRUI' 'BLAN' 'UNIF' 1.0 0.1 MAIL ;
  124. TN = 'BRUI' 'BLAN' 'UNIF' 5000 1000 MAIL ;
  125.  
  126. YN2 = 'BRUI' 'BLAN' 'UNIF' 0.3 0.1 MAIL ;
  127. YO2 = 'BRUI' 'BLAN' 'UNIF' 0.2 0.1 MAIL ;
  128. RYN2 = 'NOMC' (RN * YN2) 'N2 ' 'NATU' 'DISCRET' ;
  129. RYH2 = 'NOMC' (RN * YN2) 'H2 ' 'NATU' 'DISCRET' ;
  130. RYN = RYH2 'ET' RYN2 ;
  131.  
  132. VNX = 'BRUI' 'BLAN' 'UNIF' 1.0 0.1 MAIL ;
  133. VNY = 'BRUI' 'BLAN' 'UNIF' 1.0 0.1 MAIL ;
  134.  
  135. VN = ('NOMC' VNX 'UX' 'NATU' 'DISCRET' ) 'ET'
  136. ('NOMC' VNY 'UY' 'NATU' 'DISCRET' ) ;
  137.  
  138. GN = ('NOMC' (VNX * RN) 'UX' 'NATU' 'DISCRET' ) 'ET'
  139. ('NOMC' (VNY * RN) 'UY' 'NATU' 'DISCRET' ) ;
  140.  
  141. * Energie de l'ordre de (0.5 * 1^2) '+' (*10^6)
  142.  
  143. RET = 'BRUI' 'BLAN' 'UNIF' 1.0D6 1.0D5 MAIL ;
  144.  
  145. *
  146. **** PRIM
  147. *
  148.  
  149. LCOMRYN = 'EXTRAIRE' RYN 'COMP' ;
  150. LCOMRVN = 'EXTRAIRE' GN 'COMP' ;
  151.  
  152. VIT PRES TEMP YN GAMN
  153. = 'PRIM' 'PERFTEMP' PGAZ RN GN RET RYN ;
  154.  
  155. *
  156. **** On inverse l'ordre des composantes des fractions massiques
  157. *
  158.  
  159. RYN = RYN2 'ET' RYH2 ;
  160.  
  161. VIT1 PRES1 TEMP1 YN1 GAMN1
  162. = 'PRIM' 'PERFTEMP' PGAZ RN GN RET RYN ;
  163.  
  164.  
  165. ERRV = ('MAXIMUM' (VIT1 '-' VIT) 'ABS')
  166. '/' ('MAXIMUM' VN 'ABS') ;
  167.  
  168. ERRP = 'MAXIMUM' ((PRES1 '-' PRES) '/' PRES) 'ABS';
  169.  
  170. ERRT = 'MAXIMUM' ((TEMP1 '-' TEMP) '/' TEMP) 'ABS';
  171.  
  172. ERRG = 'MAXIMUM' ((GAMN1 '-' GAMN) '/' GAMN) 'ABS';
  173.  
  174. ERRY = ('MAXIMUM' (YN1 '-' YN) 'ABS')
  175. '/' ('MAXIMUM' YN 'ABS') ;
  176.  
  177. ERRO = 'MAXIMUM' ('PROG' ERRV ERRP ERRT ERRG ERRY );
  178.  
  179. 'SI' (ERRO > 1.0D-16) ;
  180. 'MESSAGE' ('CHAINE' 'Erreur maximum ' ERRO) ;
  181. 'ERREUR' 5;
  182. 'FINSI' ;
  183.  
  184.  
  185. * On inverse l'ordre des composante de vitesse
  186.  
  187. GN = ('NOMC' (VNY * RN) 'UY' 'NATU' 'DISCRET' ) 'ET'
  188. ('NOMC' (VNX * RN) 'UX' 'NATU' 'DISCRET' ) ;
  189.  
  190. VIT1 PRES1 TEMP1 YN1 GAMN1
  191. = 'PRIM' 'PERFTEMP' PGAZ RN GN RET RYN ;
  192.  
  193.  
  194.  
  195. ERRV = ('MAXIMUM' (VIT1 '-' VIT) 'ABS')
  196. '/' ('MAXIMUM' VN 'ABS') ;
  197.  
  198. ERRP = 'MAXIMUM' ((PRES1 '-' PRES) '/' PRES) 'ABS';
  199.  
  200. ERRT = 'MAXIMUM' ((TEMP1 '-' TEMP) '/' TEMP) 'ABS';
  201.  
  202. ERRG = 'MAXIMUM' ((GAMN1 '-' GAMN) '/' GAMN) 'ABS';
  203.  
  204. ERRY = ('MAXIMUM' (YN1 '-' YN) 'ABS')
  205. '/' ('MAXIMUM' YN 'ABS') ;
  206.  
  207. ERRO = 'MAXIMUM' ('PROG' ERRV ERRP ERRT ERRG ERRY );
  208.  
  209.  
  210. 'SI' (ERRO > 1.0D-16) ;
  211. 'MESSAGE' ('CHAINE' 'Erreur maximum ' ERRO) ;
  212. 'ERREUR' 5;
  213. 'FINSI' ;
  214.  
  215. 'FIN' ;
  216.  
  217.  
  218.  
  219.  
  220.  

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