Télécharger la plaquette Cast3M


fr en

Télécharger decroissanceVF.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  decroissanceVF.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. GRAPH = 'O'   ;
  6. 'OPTI' 'ECHO' 1 ;
  7. CRIT1 = 1.D-6 ;
  8. 'SAUT' 'PAGE' ;
  9. *
  10. *---------------------------------------------------------------------
  11. * On test en 0D la décroissance radioactive en partant d'une concentration
  12. * initiale T0 (unité). FLux nul aux bords du domaine. Plusieurs mailles
  13. * car cela ne coute rien d'essayer.
  14. * diffusion discrétisée (obligatoire en EFMH) mais inutile car C constante.
  15. * retard et porosité differentes de 1.
  16. *
  17. *---------------------------------------------------------------------
  18. *
  19. *------------------
  20. * Options generales
  21. *------------------
  22. *
  23. 'OPTI' 'DIME' 2 'ELEM' 'QUA4' ;
  24. 'OPTI' 'ISOV' 'SURF'          ;
  25. *
  26. *
  27. *=========
  28. * MAILLAGE
  29. *=========
  30. *
  31. *
  32. *- Création des points supports du contour du domaine, et des droites
  33. *- passant par les centres et les faces pour le post-traitement.
  34. *
  35. L     = 100.D0          ;
  36. LS2   = L / 2.D0        ;
  37. H     = 1.D0            ;
  38. HS2   = H / 2.D0        ;
  39. X0    = 0.D0            ;
  40. X1    = X0 + L          ;
  41. Y0    = 0.D0            ;
  42. Y1    = Y0 + H          ;
  43. INUMX = 100             ;
  44. INUMY = 1               ;
  45. INUM1 = INUMX - 1       ;
  46. Y01   = Y0 + Y1 * 0.5D0 ;
  47. DX    = X1 - X0 / INUMX ;
  48. DX1   = DX / 2.D0       ;
  49. XG    = X0 + DX1        ;
  50. XD    = X1 - DX1        ;
  51. *
  52. A1    = X0 Y0   ;
  53. A3    = X1 Y0   ;
  54. D1    = X0 Y1   ;
  55. D3    = X1 Y1   ;
  56. B1    = X0 Y01  ;
  57. B3    = X1 Y01  ;
  58. C1    = XG Y01  ;
  59. C3    = XD Y01  ;
  60. P6    = LS2 Y01 ;
  61. *
  62. *- Création des DROITES frontieres
  63. *
  64. DRBAS = A3 'DROI' INUMX A1 ;
  65. DRGAU = A1 'DROI' INUMY D1 ;
  66. DRHAU = D1 'DROI' INUMX D3 ;
  67. DRDRO = D3 'DROI' INUMY A3 ;
  68. PELIM = DX1 / (5. * INUMX) ;
  69. *
  70. *- Creation maillage GEOMETRIQUE
  71. *
  72. PTOT1  = 'DALL' DRBAS DRGAU DRHAU DRDRO ;
  73. PTOT2  = 'ORIE' PTOT1                   ;
  74. *
  75. *- Creation maillage HYBRIDE y compris sous-objets (cond. limites)
  76. *
  77. DRMID = B1 'DROI' INUMX B3              ;
  78. DRMIC = C1 'DROI' INUM1 C3              ;
  79. EXT1  = 'MANU' 'POI1' B1                ;
  80. *
  81. QFTOT= 'CHAN' PTOT2  QUAF ;
  82. QFGAU= 'CHAN' DRGAU  QUAF ;
  83. QFDRO= 'CHAN' DRDRO  QUAF ;
  84.  ELIM PELIM (QFTOT ET QFGAU ET QFDRO ET DRMID ET DRMIC ET EXT1) ;
  85. *
  86. *================
  87. * INITIALISATIONS
  88. *================
  89. *
  90. *                                                    ----------------
  91. *                                                    = MODELISATION =
  92. *                                                    ----------------
  93. MODHYB = MODE QFTOT 'DARCY' 'ANISOTROPE' ;
  94. MODDRO = MODE QFDRO 'DARCY' 'ANISOTROPE' ;
  95. MODGAU = MODE QFGAU 'DARCY' 'ANISOTROPE' ;
  96. CHYB1 = 'DOMA' MODHYB 'SURFACE'          ;
  97. CHYB2 = 'DOMA' MODHYB 'NORMALE'          ;
  98.  
  99.  
  100.  
  101. *
  102. *                                               ---------------------
  103. *                                               = Donnees physiques =
  104. *                                               ---------------------
  105. *
  106. T0    = 1.D0 ; T1  = 1.D0    ;
  107. VK    = 1.D-2                ;
  108. VK1 = manu CHPO (doma modhyb centre) 1 'K11' VK NATURE DISCRET;
  109. VK2 = manu CHPO (doma modhyb centre) 1 'K22' VK NATURE DISCRET;
  110. VK12 = manu CHPO (doma modhyb centre) 1 'K21' 0.D0 NATURE DISCRET;
  111. MATI2 = VK1 et VK2 et VK12;
  112. *
  113. *
  114. *                                             -----------------------
  115. *                                             = Donnees transitoire =
  116. *                                             -----------------------
  117. *  TETA   : Parametre de le theta-méthode
  118. *  TMAX   : Temps final
  119. *  TSUP   : Temps pour conditions aux limites
  120. *  DELTAT : Pas de temps
  121. *
  122. TETA   = 0.00D0       ;
  123. TMIN   = 0.D0         ;
  124. TMAX2  = 15.00D0      ;
  125. TMAX   = 30.00D0      ;
  126. TSUP   = 1.2D0 * TMAX ;
  127. DELTAT = 1.D0        ;
  128. *
  129. LICALC = 'PROG' TMIN 'PAS' (DELTAT/2.D0) TMAX2 ;
  130. LICALC = LICALC ET ('PROG' (TMAX2 + DELTAT) 'PAS' DELTAT TMAX) ;
  131. LISAUV = 'PROG' TMAX                   ;
  132. *
  133. *                                             ------------------------
  134. *                                             = Conditions initiales =
  135. *                                             ------------------------
  136. TCHYB = 'MANU' 'CHPO' ('DOMA' MODHYB 'CENTRE') 1 'H' T0 
  137.                                              'NATURE' 'DISCRET' ;
  138. *
  139. *                                                       --------------
  140. *                                                       = T imposée  =
  141. *                                                       --------------
  142. *
  143. *                                          ---------------------------
  144. *                                          = Table DARCY_TRANSITOIRE =
  145. *                                          ---------------------------
  146. *
  147. *
  148. *
  149.  
  150. *-- Table de transport :
  151. Transp                         = 'TABLE';
  152. Transp . 'MODELE'              = MODHYB ;
  153. Transp.'TEMPS'                 = 'TABLE';
  154. Transp.'CONCENTRATION'         = 'TABLE';
  155. Transp.'FLUXDIFF'              = 'TABLE';
  156. Transp.'FLUXCONV'              = 'TABLE';
  157. Transp.'CARACTERISTIQUES'      = MATI2;
  158.  
  159. Transp . 'COEF_RETARD' = 100.D0;
  160. Transp . 'POROSITE'    = 0.1D0;
  161.  
  162. * Conditions initiales :
  163. Transp.'TEMPS'. 0              = TMIN ;
  164. Transp.'CONCENTRATION'. 0      = TCHYB;
  165. Transp.'FLUXDIFF'. 0           = 0.D0 * TCHYB;
  166. Transp.'FLUXCONV'. 0           = 0.D0 * TCHYB;
  167.  
  168. * Conditions aux limites :
  169. * flux diffusif nul donc rien
  170.  
  171.  
  172. * Paramètres numériques :
  173. Transp.'THETA_DIFF'            = 0.5D0;
  174. Transp.'THETA_CONVECTION'      = 1.0D0;
  175. Transp.'LUMP'                  = FAUX;
  176. Transp.'TYPDISCRETISATION'     = 'VF';
  177. Transp.'DECENTR'          = FAUX;
  178.  
  179. TABRES = table METHINV;
  180. TABRES . 'TYPINV' = 1;
  181. TABRES . 'PRECOND' = 3;
  182.  
  183. Transp . 'METHINV' = TABRES;
  184.  
  185.  
  186. Transp.'TEMPS_CALCULES'        = LiCalc;
  187. Transp.'TEMPS_SAUVES'          = LiSauv;
  188.  
  189.  
  190. Transp . INTCONC = TABLE;
  191. Transp . INTCONC . 0 = 0.D0 * TCHYB;
  192. *                                                            ==========
  193. Transp . 'DECROISSANCE' = (log 2.D0) / 10.D0;
  194. *                                                            | CALCUL |
  195. *                                                            ==========
  196.  
  197.  
  198. *=======================
  199. * Resolution transitoire
  200. *=======================
  201. *
  202. TRANSGEN TRANSP ;
  203. *
  204. *
  205. *=================
  206. * POST-TRAITEMENT
  207. *=================
  208.  
  209.  
  210. ****** comparaison sol analytique ************
  211. **   cini * exp(-lambda t)
  212.  
  213. concEFMH =  Transp . CONCENTRATION . (('DIME' transp . TEMPS) '-' 1);
  214. solan =  Transp . CONCENTRATION . 0 * ('EXP'
  215.              (-1.D0 * TRANSP . 'DECROISSANCE' * TMAX));
  216.  
  217. toto = concEFMH '-' solan ;
  218. toto = 'KOPS' toto * toto;
  219. toto = ('DOMA' modhyb volume) * toto;
  220. toto = 'MAXIMUM' (('RESULT' toto));
  221. titi = ('DOMA' modhyb volume) * ('KOPS' solan * solan);
  222. titi = 'MAXIMUM' ('RESULT' titi);
  223. toto = (toto '/' titi)**0.5D0;
  224.  
  225. 'MESSAGE' 'erreur relative L2' (toto); 
  226.  
  227. *
  228.  
  229. 'SI' (toto > 1.D-3) ;
  230.    'ERREUR' 5 ;
  231. 'SINON' ;
  232.     'ERREUR' 0;
  233. 'FINSI' ;
  234.  
  235.  
  236. *
  237. 'FIN' ;
  238.  
  239.  
  240.  
  241.  
  242.  
  243.  
  244.  

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