Télécharger la plaquette Cast3M


fr en

Télécharger adve_04.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : adve_04.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Cas-test de l'operateur ADVEction dans la formulation DIFFUSION      *
  5. *                                                                      *
  6. * Comparaison a une solution analytique.                               *
  7. *                                                                      *
  8. * On calcule la concentration d'un fluide qui s'ecoule dans un tuyau   *
  9. * chauffe sur toute sa longueur.                                       *
  10. * On verifie la concentration du fluide en sortie du tuyau.            *
  11. *                                                                      *
  12. * Le flux lineique est de         : 500 mol/m-3/m-1                    *
  13. * La concentration en entree est de : 20mol.m-3                        *
  14. * La vitesse d'ecoulement est de  : 25 cm/s                            *
  15. * La section du tuyau est de      : 5.e-5 m2, diametre ~ 8 mm          *
  16. * La longueur du tuyau est de     : 5 m                                *
  17. *                                                                      *
  18. * La concentration attendue en sortie est de : 70mol.m-3               *
  19. *                                                                      *
  20. * On utilise les elements lineaires et quadratiques (TUY2,TUY3)        *
  21. *                                                                      *
  22. * On met une conduction legere pour eviter les problemes de pivot nul  *
  23. * dans la matrice                                                      *
  24. *                                                                      *
  25. *----------------------------------------------------------------------*
  26. 'OPTI' 'DIME' 3 'ELEM' 'SEG2' ;
  27. *                                                                      *
  28. * Commentez cette ligne pour voir les traces :                         *
  29. 'OPTI' 'TRAC' 'PSC' ;
  30. *                                                                      *
  31. O1       = 0 0 0 ;
  32. X1       = 1 0 0 ;
  33. P1       = 2.5 * X1 ;
  34. P2       = 5.0 * X1 ;
  35. Y1       = 0 1 0 ;
  36. Z1       = 0 0 1 ;
  37. *                                                                      *
  38. L1       = O1 'DROI' 5 P1 'COUL' 'VERT' ;
  39. L2       = P1 'DROI' 10 P2 ;
  40. L2       = 'CHAN' 'QUAD' L2 'COUL' 'BLEU' ;
  41. 'TITR' ' Maillage du tuyau, longueur = 5 m (Vert=SEG2,Bleu=SEG3) ' ;
  42. 'TRAC' 'QUAL' (L1 'ET' L2) ;
  43. *                                                                      *
  44. mo1      ='MODE' L1 diffusion advection 'TUY2' ;
  45. mo2      ='MODE' L2 diffusion advection 'TUY3' ;
  46. mo1      = mo1 'ET' mo2 ;
  47. *                                                                      *
  48. L0       = L1 'ET' L2 ;
  49. *
  50. * Le flux advectee dans un tuyau est J = Cdif.Sect.V :                 *
  51. ma1      = 'MATE' mo1 'CDIF' 4.e6 'VITE' 0.25 'SECT' 5.e-5 'K' 1d-6 ;
  52. *                                                                      *
  53. sq1      = 'SOUR' mo1 ma1 L0 (5.e2 / 5.e-5) ;
  54. list (MAXI (resu sq1));
  55. 'TITR' ' Terme source le long du tuyau ' ;
  56. 'DESS' ('EVOL' 'VERT' 'CHPO' sq1 L0 'QCO') 'YBOR' 0. 300. ;
  57. *                                                                      *
  58. KA1      = 'ADVE' mo1 ma1 ;
  59. KC1      = 'COND' mo1 ma1 ;
  60. *                                                                      *
  61. cl1      = 'BLOQ' 'CO' O1 ;
  62. dcl1     = 'DEPI' cl1 20. ;
  63. *                                                                      *
  64. chC2     = 'RESO' (ka1 'ET' KC1 'ET' CL1) (sq1 'ET' dcl1);
  65. *                                                                      *
  66. 'TITR' ' Evolution de la concentration le long du tuyau ' ;
  67. 'DESS' ('EVOL' 'VERT' 'CHPO' chC2 L0 'CO') ;
  68. *                                                                      *
  69. * Test :
  70. cref1    = 70. ;
  71. Cp1      = 'EXTR' chC2 'CO' P2 ;
  72. err1     = 'ABS' (cREF1 - Cp1) / Cp1 ;
  73. *                                                                      *
  74. 'OPTI' 'ECHO' 0 ;
  75. 'MESS' ;
  76. mot1     = 'CHAI' ' > concentration calculee en  sortie = '
  77.   'FORMAT' '(F3.0)' Cp1 'mol.m-3 pour 70mol.m-3 attendu ' ;
  78. 'MESS' mot1 ;
  79. 'MESS' ' > Erreur relative = ' err1 ;
  80. 'MESS' ;
  81. 'OPTI' 'ECHO' 1 ;
  82. *                                                                      *
  83. 'SI' (err1 > 1.e-5) ;
  84.   'ERRE' 5 ;
  85. 'SINO' ;
  86.   'OPTI' 'ECHO' 0 ;
  87.   'MESS' ;
  88.   'MESS' ' > Test reussi ' ;
  89.   'MESS' ;
  90. 'FINS' ;
  91. *                                                                      *
  92. 'FIN' ;
  93.  
  94.  
  95.  
  96.  
  97.  

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