Télécharger la plaquette Cast3M


fr en

Télécharger adve_06.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : adve_06.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Cas-test de l'operateur ADVEction dans la formulation DIFFUSION      *
  5. *                                                                      *
  6. * On resout : v.gradC = 4,                                             *
  7. *                                                                      *
  8. * Avec : vx=1, vy=1                                                    *
  9. * on a : dC/dx + dC/dy = 4                                             *
  10. *                                                                      *
  11. * De + : C(1,0)=0, C(0,1)=0                                            *
  12. * Sur 1 element QUA4 de geometrie carree et de cote 1.                 *
  13. *                                                                      *
  14. * La solution analytique est : C(x,y) = 2(x+y-1)                       *
  15. * Soit : C(0,0)=C(P1)=-2 et C(1,1)=C(P3)=2                             *
  16. *                                                                      *
  17. ***************************************************
  18. ** TESTING FILE FOR THE OPERATOR ADVE BUILDING A  *
  19. ** NON-SYMMETRIC STIFFNESS MATRIX                 *
  20. ***************************************************
  21.  
  22. OPTI DIME 2 ELEM qua4 MODE PLAN CONT ;
  23. p1     = 0. 0. ;
  24. p2     = 1. 0. ;
  25. p3     = 1. 1. ;
  26. p4     = 0. 1. ;
  27. LI1    = droite 1 p1 p2 ;
  28. LI2    = droite 1 p2 p3 ;
  29. LI3    = droite 1 p3 p4 ;
  30. LI4    = droite 1 p1 p4 ;
  31.  
  32. SU     = DALL LI1 LI2 LI3 LI4 ;
  33. MOD1   = MODE SU DIFFUSION ISOTROPE ADVECTION ;
  34. jj     = MATE MOD1 'VITX' 1. 'VITY' 1. ;
  35. matno1 = ADVE mod1 JJ ;
  36. list matno1;
  37.  
  38. * Matrice elementaire numero   : 1 ( ligne1,ligne2,ligne3...)
  39. * -0.33333E+00  0.83333E-01  0.16667E+00  0.83333E-01 ;
  40. * -0.25000E+00  0.00000E+00  0.25000E+00  0.00000E+00 ;
  41. * -0.16667E+00 -0.83333E-01  0.33333E+00 -0.83333E-01 ;
  42. * -0.25000E+00  0.00000E+00  0.25000E+00  0.00000E+00 ;
  43.  
  44. * Le flux surfacique :
  45. vec1 ='SOUR' MOD1 SU 4. ;
  46.  
  47. * Les concentrations imposees :
  48. BL2  = BLOQ 'CO' P2 ;
  49. BL4  = BLOQ 'CO' P4 ;
  50.  
  51. *A linear system is solved with a non-symmetric matrix
  52. SOL7 = RESO (matno1 et (bl2 et bl4)) vec1 ;
  53.  
  54. * Verification :
  55. CP1  = -2. ;
  56. CP3  = +2. ;
  57. ERR1 = ABS((EXTRAIRE SOL7 'CO' P1) - CP1) ;
  58. ERR2 = ABS((EXTRAIRE SOL7 'CO' P3) - CP3) ;
  59. SI ((ERR1+ERR2) > 1.E-10);
  60.   ERRE 5 ;
  61. * MESS '###### TEST NOT PASSED ######';
  62. SINON ;
  63.   ERRE 0 ;
  64. * MESS '###### TEST PASSED ######';
  65. FINSI ;
  66.  
  67. *****************************************
  68. matno1 = ADVE mod1 JJ 'SYMM';
  69. list matno1;
  70.  
  71. * Matrice elementaire numero   : 1 ( ligne1,ligne2,ligne3...)
  72. * -0.33333E+00 -0.83333E-01 0.00000E+00 -0.83333E-01
  73. * -0.83333E-01  0.00000E+00 0.83333E-01  0.00000E+00
  74. *  0.00000E+00  0.83333E-01 0.33333E+00  0.83333E-01
  75. * -0.83333E-01  0.00000E+00 0.83333E-01  0.00000E+00
  76.  
  77. FIN;
  78.  
  79.  
  80.  

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