Télécharger la plaquette Cast3M


fr en

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