Télécharger la plaquette Cast3M


fr en

Télécharger cormasse.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : cormasse.dgibi
  2. *
  3. *
  4. * Tests de la procédure cormasse qui corrige un chpo afin que ses valeurs
  5. * soient positives et inférieures à Maxrho, et que que son intégrale
  6. * soit égale à une valeur cible.
  7. * 
  8. 'OPTI' 'DIME' 2 'ELEM' 'QUA4' ;
  9. *
  10. Lref = 3. ;
  11. p1   = 0.   0.   ;
  12. p2   = Lref 0.   ;
  13. p3   = Lref Lref ;
  14. p4   = 0.   Lref ;
  15. vtf  = 'MANU' 'QUA4' P1 P2 P3 P4 ;
  16. Mvtf = 'CHAN' vtf 'QUAF' ;
  17. 'ELIM' Mvtf 1.D-5 ;
  18. $vtf = 'MODE' Mvtf 'NAVIER_STOKES' 'QUAF' ;
  19. vtf  = 'DOMA' $vtf 'MAILLAGE' ;
  20. *
  21. x y    = 'COOR' vtf ;
  22. MaxRho = 10000. ;
  23. *
  24. *
  25. * 1. Champ uniforme 2 : Mcfd=18
  26. *
  27. *
  28. RHOi = 'KCHT' $vtf 'SCAL' 'SOMMET' 2. ;
  29. *
  30. * 1.1/ Correction nulle
  31. mcible1 = 18. ;
  32. dRHO1 RHOk1 MiCFD = CORMASSE $vtf RHOi Mcible1 MaxRho ;
  33. *
  34. * 1.2/ Correction uniforme +9 / V = +1
  35. mcible2 = 27. ;
  36. dRHO2 RHOk2 MiCFD = CORMASSE $vtf RHOi Mcible2 MaxRho ;
  37. *
  38. * 1.3/Correction uniforme -9 / V = -1
  39. mcible3 = 9. ;
  40. dRHO3 RHOk3 MiCFD = CORMASSE $vtf RHOi Mcible3 MaxRho ;
  41. *
  42. *
  43. * 2. Champ non uniforme linéaire (y+1) : Mcfd=22.5
  44. *
  45. *
  46. RHOi = 'KCHT' $vtf 'SCAL' 'SOMMET' (y '+' 1.) ;
  47. *
  48. * 2.4/Correction uniforme -9/V = -1 en restant positif
  49. mcible4 = 13.500000001 ;
  50. dRHO4 RHOk4 MiCFD = CORMASSE $vtf RHOi Mcible4 MaxRho ;
  51. *
  52. * 2.5/Idem légèrement négatif
  53. mcible5 = 13.499999999 ;
  54. dRHO5 RHOk5 MiCFD = CORMASSE $vtf RHOi Mcible5 MaxRho ;
  55. *
  56. * 2.6/Correction non uniforme (redistribution du négatif)
  57. mcible6 = 9. ;
  58. dRHO6 RHOk6 MiCFD = CORMASSE $vtf RHOi Mcible6 MaxRho ;
  59. *
  60. * 2.7/Correction non uniforme (redistribution du négatif)
  61. mcible7 = 1. ;
  62. dRHO7 RHOk7 MiCFD = CORMASSE $vtf RHOi Mcible7 MaxRho ;
  63. *
  64. *
  65. * Tests de non-régression
  66. *
  67. *
  68. 'LIST' (rhok2 '-' 1.);
  69. 'LIST' (rhok2 '-' 3.);
  70. 'LIST' (rhok3 '-' 1.);
  71. 'LIST' (rhok4 '-' y);
  72. 'LIST' (rhok5 '-' y);
  73. 'LIST' (rhok6 '-' 2.4) ;
  74. 'LIST' (rhok7 '-' (2./3.));
  75. *
  76. erho1 = rhok1 '-' 2. ;
  77. val1 = 'MAXI' erho1 'ABS' ;
  78. erho2 = rhok2 '-' 3. ;
  79. val2 = 'MAXI' erho2 'ABS' ;
  80. erho3 = rhok3 '-' 1. ;
  81. val3 = 'MAXI' erho3 'ABS' ;
  82. erho4 = rhok4 '-' y ;
  83. val4 = 'MAXI' erho4 'ABS' ;
  84. erho5 = rhok5 '-' y ;
  85. val5 = 'MAXI' erho5 'ABS' ;
  86. erho6 = ('REDU' rhok6 (p3 'ET' p4)) '-' 2.4 ;
  87. val6 = 'MAXI' erho6 'ABS' ;
  88. erho7 = ('REDU' rhok7 (p3 'ET' p4)) '-' (2./3.) ;
  89. val7 = 'MAXI' erho7 'ABS' ;
  90. 'MESS' val1 val2 val3 val4 val5 val6 val7 ;
  91. *
  92. graph = faux ;
  93. 'SI' graph ;
  94.    'TRAC' vtf RHOk1 'TITR' 'Rho_i=2 , Rho_f=2' ;
  95.    'TRAC' vtf RHOk2 'TITR' 'Rho_i=2 , Rho_f=3' ;
  96.    'TRAC' vtf RHOk3 'TITR' 'Rho_i=2 , Rho_f=1' ;
  97.    'TRAC' vtf RHOk4 'TITR' 'Rho_i=y+1 , Rho_f=y^+' ;
  98.    'TRAC' vtf RHOk5 'TITR' 'Rho_i=y+1 , Rho_f=y^-' ;
  99.    'TRAC' vtf RHOk6
  100.    'TITR' 'Rho_i=y+1 , Rho_f=(0,0.9,2.4) en (bas,milieu,haut)' ;
  101.    'TRAC' vtf RHOk7
  102.    'TITR' 'Rho_i=y+1 , Rho_f=(0,0,2/3) en (bas,milieu,haut)' ;
  103. 'FINS' ;
  104. *
  105. vmax = 'MAXI' ('PROG' val1 val2 val3 val4 val5 val6 val7) ;
  106. 'SI' (vmax '>' 1.D-8) ;
  107.    'ERRE' 5 ;
  108. 'FINS' ;
  109. *
  110. 'FIN' ;
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  

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