Télécharger precipite1EFMH.dgibi

Retour à la liste

Numérotation des lignes :

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

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