Télécharger la plaquette Cast3M


fr en

Télécharger precipite4EFMH.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  precipite4EFMH.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *
  5. ********************* CAS TEST : precipite4.dgibi ********************
  6. *
  7. GRAPH = 'N'   ;
  8. 'OPTI' 'ECHO' 1 ;
  9. CRIT1 = 1.D-6 ;
  10. 'SAUT' 'PAGE' ;
  11. *
  12. *---------------------------------------------------------------------
  13. * TEST PRECIPITE4
  14. * precipité uniforme initial de 1 (par unité de volume solide) partout.
  15. * vitesse vx = 1 vy = 0 uniforme
  16. * maillage 2D pseudo 1D.
  17. * concini = 0.
  18. * concentration à gauche nulle, flux diff nul à droite.
  19. * Le précipité fond, le frond se déplace de gauche à droite selon
  20. * l'expression VX * limsol * tps '/' F00 '/' (1.D0 '-' POROS)
  21. * ou poros est la porosité, F00 la concentration de precipité, tps
  22. * le temps écoulé et limsol la limite de solubilité
  23. *
  24. * utilise TRANSGEN
  25. *---------------------------------------------------------------------
  26. *
  27. *------------------
  28. * Options generales
  29. *------------------
  30. *
  31. 'OPTI' 'DIME' 2 'ELEM' 'QUA4' ;
  32. 'OPTI' 'ISOV' 'SURF'          ;
  33. *
  34. *
  35. *=========
  36. * MAILLAGE
  37. *=========
  38. *
  39. *
  40. *- Création des points supports du contour du domaine, et des droites
  41. *- passant par les centres et les faces pour le post-traitement.
  42. *
  43. L     = 100.D0          ;
  44. H     = 1.D0            ;
  45. HS2   = H/2.D0          ;
  46. X0    = 0.D0            ;
  47. X1    = X0 + L          ;
  48. Y0    = 0.D0            ;
  49. Y1    = Y0 + H          ;
  50. INUMX = 200             ;
  51. INUMY = 1               ;
  52. INUM1 = INUMX - 1       ;
  53.  
  54. DX    = X1 - X0 / INUMX ;
  55.  
  56. A1    = X0 Y0   ;
  57. A3    = X1 Y0   ;
  58. D1    = X0 Y1   ;
  59. D3    = X1 Y1   ;
  60.  
  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 = DX / (5. * INUMX) ;
  69. *
  70. *- Creation maillage GEOMETRIQUE
  71. *
  72. PTOT1  = 'DALL' DRBAS DRGAU DRHAU DRDRO ;
  73. PTOT1  = 'ORIE' PTOT1                   ;
  74. *
  75. *- Creation maillage HYBRIDE y compris sous-objets (cond. limites)
  76. *
  77. QFTOT= 'CHAN' PTOT1  QUAF ;
  78. QFGAU= 'CHAN' DRGAU  QUAF ;
  79. QFDRO= 'CHAN' DRDRO  QUAF ;
  80.  ELIM PELIM (QFTOT ET QFGAU ET QFDRO) ;
  81. *
  82. *================
  83. * INITIALISATIONS
  84. *================
  85. *
  86. *                                                    ----------------
  87. *                                                    = MODELISATION =
  88. *                                                    ----------------
  89. MODHYB = MODE QFTOT 'DARCY' 'ANISOTROPE' ;
  90. MODDRO = MODE QFDRO 'DARCY' 'ISOTROPE' ;
  91. MODGAU = MODE QFGAU 'DARCY' 'ISOTROPE' ;
  92.  
  93. *
  94. *                                               ---------------------
  95. *                                               = Donnees physiques =
  96. *                                               ---------------------
  97. *
  98.  
  99. *     ---------------------
  100. *          DIFFUSIVITE
  101. *     ---------------------
  102.  
  103. VK    = 1.D-6                ;
  104. MATI2 = MANU 'CHPO' ('DOMA' MODHYB 'CENTRE')  'K' VK ;
  105.  
  106.  
  107. *     ---------------------
  108. *          VITESSE
  109. *     ---------------------
  110.  
  111. VX1   = 5.D0 ; VY1 = 0.D0    ;
  112. *
  113. SPEED = 'MANU' 'CHPO' ('DOMA' MODHYB 'FACE')  2
  114.         'VX' VX1 'VY' VY1 'NATURE' 'DISCRET' ;
  115. SPEEDC = 'MANU' 'CHPO' ('DOMA' MODHYB 'CENTRE')  2
  116.         'VX' VX1 'VY' VY1 'NATURE' 'DISCRET' ;
  117. *
  118. *    -----------------------
  119. *      Donnees transitoire 
  120. *    -----------------------
  121. *
  122. TETA   = 1.00D0       ;
  123. TMIN   = 0.D0         ;
  124. TMAX   = 0.25D2      ;
  125. TSUP   = 1.2D0 * TMAX ;
  126. DELTAT = 0.250D0        ;
  127. *
  128. LICALC = 'PROG' TMIN 'PAS' DELTAT TMAX ;
  129. LISAUV = LICALC                        ;
  130.  
  131.  
  132. *
  133. *      ------------------------
  134. *       Conditions initiales 
  135. *      ------------------------
  136.  
  137. T0 = 0.D0 ;
  138. CINI =  'MANU' 'CHPO' ('DOMA' MODHYB 'CENTRE') 1 'I129' T0
  139.                                   'NATURE' 'DISCRET' ;
  140.  
  141.  
  142. *         --------------
  143. *          T imposée  
  144. *         --------------
  145.  
  146. TGAUCHE = 0.D0  ;
  147. T1 = TGAUCHE;
  148. TGAUCHE = 'KCHT' MODGAU SCAL CENTRE (1.D0 * TGAUCHE);
  149. TGAUCHE = NOMC I129 TGAUCHE;
  150. TGAUCHE = 'CHAR' TGAUCHE
  151.         ('EVOL' 'MANU' ('PROG' 0. TSUP) ('PROG' 1. 1.));
  152.  
  153. *
  154. *                                          ---------------------------
  155. *                                          = Table DARCY_TRANSITOIRE =
  156. *                                          ---------------------------
  157. *
  158. *
  159. *
  160.  
  161. *-- Table de transport :
  162. Transp                         = 'TABLE';
  163. Transp . 'MODELE'              = MODHYB ;
  164. Transp.'TEMPS'                 = 'TABLE';
  165. Transp.'CONCENTRATION'         = 'TABLE';
  166. Transp.'FLUXDIFF'              = 'TABLE';
  167. Transp.'FLUXCONV'              = 'TABLE';
  168. Transp.'CARACTERISTIQUES'      = MATI2;
  169. Transp.'POROSITE'              = 0.2D0;
  170. Transp.'COEF_RETARD'           = 10.D0;
  171. Transp.'CONVECTION'            = 1.D0 * SPEED ;
  172. Transp.'VITELEM'               = 1.D0 * SPEEDC ;
  173. *Transp . 'ALPHAL'              = MANU 'CHPO' (doma MODHYB centre)
  174. *                                 'SCAL' 1.D0;
  175. *Transp . 'ALPHAT'              = MANU 'CHPO' (doma MODHYB centre)
  176.  
  177. *Transp.'DECROISSANCE'          = ('LOG' 2.D0) '/' 5.D0;
  178.  
  179. *                                 'SCAL' 0.D0;
  180.  
  181. * Conditions initiales :
  182. Transp.'TEMPS'. 0              = TMIN ;
  183. Transp.'CONCENTRATION'. 0      = CINI ;
  184. Transp.'FLUXDIFF'. 0           = 'MANUEL' 'CHPO'
  185.                         (doma MODHYB face)  'I129' 0.D0 ;
  186. Transp.'FLUXCONV'. 0           = 'MANUEL' 'CHPO'
  187.                         (doma MODHYB face)  'I129' 0.D0 ;
  188.  
  189. * Conditions aux limites DIRICHLET à gauche:
  190. Transp . 'TRACE_IMPOSE' =  TGAUCHE       ;
  191.  
  192. * CONDITION FLUX MIXTE à gauche
  193. *Transp . 'FLUXTOT_IMP' =  TGAUCHE       ;
  194.  
  195. * Fluxdiffusif a gauche
  196. *Transp . 'FLUX_IMPOSE' =  TGAUCHE        ;
  197.  
  198.  
  199. * Paramètres numériques :
  200. Transp.'THETA_DIFF'            = 1.D0;
  201. Transp.'THETA_CONV'            = 1.D0;
  202. Transp.'THETA_DEC'             = 0.5D0;
  203. Transp.'LUMP'                  = FAUX;
  204. Transp.'TYPDISCRETISATION'     = 'EFMH';
  205. Transp.'DECENTR'          = VRAI;
  206. Transp.'NUM_PECLET'       = 2.D0;
  207. Transp.'SEUILCALC'          = 1.D-14;
  208.  
  209. TABRES = table METHINV;
  210. TABRES . 'TYPINV' = 1;
  211. TABRES . 'PRECOND' = 3;
  212.  
  213. Transp . 'METHINV' = TABRES;
  214.  
  215.  
  216. Transp.'TEMPS_CALCULES'        = LiCalc;
  217. Transp.'TEMPS_SAUVES'          = LiSauv;
  218.  
  219.  
  220. Transp . INTCONC = TABLE;
  221. Transp . INTCONC . 0 = 0.D0 * CINI;
  222.  
  223.  
  224. * Precipité
  225.  
  226. * CE SONT DES CHAMPS . A CORRIGER
  227. Transp . 'LIMITE_SOLUBILITE'   =  MANU 'CHPO' (doma MODHYB centre)
  228.                                  'I129' 0.25D0;
  229.  
  230. Transp . 'COEF_DISSOLUTION' =  MANU 'CHPO' (doma MODHYB centre)
  231.                                  'I129' 2.5D8 ;
  232.  
  233. Transp . 'PRECIPITE' = 'TABLE' ;
  234. Transp . 'PRECIPITE' . 0 =  MANU 'CHPO' (doma MODHYB centre)
  235.                                  'I129' 1.D0;
  236.  
  237. Transp . INTPREC = TABLE;
  238. Transp . INTPREC . 0 = 0.D0 * CINI;                                 
  239.  
  240. * IL FAUT UN INDICE DISSOLUTION
  241.  
  242. Transp . 'DISSOLUTION' = 'TABLE'  ;
  243. Transp . 'DISSOLUTION' . 0 = 0.D0 * CINI;
  244.  
  245.  
  246.  
  247. *     TESTER UN TERME SOURCE; TESTER UNE SEULE MAILLE,
  248. *     TESTER CONVECTION PURE ...... TRES IMPORTATN
  249. *     TESTER DECROISSANCE '-' filiation
  250. *     TESTER FLUX imposés, dissolution, precipitation
  251. *     retard poros, diffusion seule. 
  252. *     changer poros ne devrait rien change 'SAUF' niveau source
  253. *                                                            ==========
  254. *                                                            | CALCUL |
  255. *                                                            ==========
  256.  
  257.  
  258. *=======================
  259. * Resolution transitoire
  260. *=======================
  261. *
  262. TRANSGEN TRANSP ;
  263. *
  264. *
  265. *=================
  266. * POST-TRAITEMENT
  267. *=================
  268.  
  269. TRANS2 = TABLE TRANSP;
  270.  
  271. *
  272. *--------------------------------------------------------------------
  273. *      Dans chaque cas on trace
  274. *         La trace de concentration le long de DRMID
  275. *         La concentration le long de DRMIC
  276. *--------------------------------------------------------------------
  277. *      Tests de NON-REGRESSION :
  278. *         Principe du maximum
  279. *         Position du centre de gravité du champ de concentration
  280. *--------------------------------------------------------------------
  281. *
  282. *                                                Critères numériques
  283. CFL = VX1 * DELTAT / DX                   ;
  284. PEK = VX1 * DX / (2. * VK)                ;
  285. FOU = 2. * VK * DELTAT / (DX * DX)        ;
  286. 'SAUT' 1 'LIGNE'                          ;
  287. 'MESS' '   Critères numériques             ' ;
  288. 'MESS' '    CFL     ' CFL ;
  289. 'MESS' '    PECLET  ' PEK ;
  290. 'MESS' '    FOURIER ' FOU ;
  291. 'MESS' '                                    ' ;
  292. *
  293. XC YC = 'COOR' (DOMA MODHYB 'CENTRE')   ;
  294. ISOR1 = INDEX ( TRANS2 . 'TEMPS') ;
  295. NTSOR = DIME ISOR1                ;
  296.  
  297.  
  298. NTSO1 = NTSOR - 1                 ;
  299. IOK   = FAUX                      ;
  300. IRESU = 1                         ;
  301. *
  302. *-----------------------
  303. REPETER VISURESU NTSO1 ;
  304. *-----------------------
  305. *
  306. IRESU = IRESU + 1                        ;
  307. INDI1 = ISOR1.IRESU                      ;
  308. TTRA  = TRANS2 . 'TEMPS' . INDI1         ;
  309. EPR1  = TRANS2 . 'CONCENTRATION' . INDI1        ;
  310.  
  311.  
  312. poros = TRANSP . 'POROSITE';
  313. limsol = 'MAXIMUM' TRANSP . 'LIMITE_SOLUBILITE';
  314. F00    = maximum TRANSP . PRECIPITE . 0 ;
  315.  
  316. solanal = VX1 * limsol * TTRA '/' F00 '/' (1.D0 '-' POROS);
  317.  
  318. 'MESSAGE' solanal;
  319.  
  320. *
  321. 'SI' ('NEG' GRAPH 'N' )                        ;
  322.    LTI2 = 'CHAINE' 'Front 1D-h temps ' TTRA    ;
  323.    'TITR' LTI2                                 ;
  324.    drmil = 'QUELCONQUE' 'SEG2' ('EXTRAIRE' epr1 maillage) ;
  325.    drmil = 'INVERSE' drmil                     ;
  326.    AV2  = 'EVOL' 'ROUG' 'CHPO' EPR1 'I129'  drmil ;
  327.    'DESS' AV2 'MIMA'  'NCLK'                         ;
  328.    toto = 'KCHA' modhyb epr1 'CHAM';
  329. *   'TRACER' modhyb toto 'NCLK';
  330. 'FINS'                                         ;
  331.  
  332.  
  333.  
  334.  
  335. * concentration au temps t
  336. *
  337. *-------------
  338. FIN VISURESU ;
  339. *-------------
  340.  
  341.  
  342. IRESU = 1                         ;
  343. *
  344. *-----------------------
  345. REPETER VISURES NTSO1 ;
  346. *-----------------------
  347. *
  348. IRESU = IRESU + 1                        ;
  349. INDI1 = ISOR1.IRESU                      ;
  350. TTRA  = TRANS2 . 'TEMPS' . INDI1         ;
  351. EPR1  = TRANS2 . 'PRECIPITE' . INDI1        ;
  352.  
  353. *
  354. 'SI' ('NEG' GRAPH 'N' )                        ;
  355.    LTI2 = 'CHAINE' 'Front 1D-h temps ' TTRA    ;
  356.    'TITR' LTI2                                 ;
  357.    drmil = 'QUELCONQUE' 'SEG2' ('EXTRAIRE' epr1 maillage) ;
  358.    drmil = 'INVERSE' drmil                     ;
  359.    AV2  = 'EVOL' 'ROUG' 'CHPO' EPR1 'I129'  drmil ;
  360.    'DESS' AV2 'MIMA'  'NCLK'                         ;
  361.    toto = 'KCHA' modhyb epr1 'CHAM';
  362. *   'TRACER' modhyb toto 'NCLK';
  363. 'FINS'                                         ;
  364.  
  365. * valeur non nulle de concentration
  366. roro = 'MASQUE' epr1 SUPERIEUR 0.1D0;
  367. * coord x du maillage
  368. zozo = 'COORDONNEE' 1 roro;
  369. * multiplie pour trouver valeur du front
  370. tutu = roro * zozo;
  371. tutu = 'KCHA' tutu modhyb 'CHAM';
  372. titi = 'ELEM' tutu SUPERIEUR 1.D-10;
  373. * on extrait la coordonnée x du front
  374. roro = 'MINIMUM' ('COORDONNEE' 1 titi);
  375.  
  376. 'LISTE' roro;
  377.  
  378. * concentration au temps t
  379. *
  380. *-------------
  381. FIN VISURES ;
  382. *-------------
  383.  
  384. 'SI' (roro 'EGA' solanal (2.D0 * DX)) ;
  385.    'ERREUR' 0 ;
  386. 'SINON'                      ;
  387.    'ERREUR' 5 ;
  388. 'FINSI' ;
  389.  
  390.  
  391. *
  392. 'FIN' ;
  393.  
  394.  
  395.  
  396.  
  397.  
  398.  
  399.  

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