Télécharger la plaquette Cast3M


fr en

Télécharger adve_07.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : adve_07.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. ************************************************************************
  6. *           Test adve_07.dgibi: Jeux de données                        *
  7. *           -----------------------------------                        *
  8. ************************************************************************
  9. *__________________________________________________________________
  10. *   TEST ADVE_07                                                   *
  11. *                                                                  *
  12. *   COMPARAISON NUMERIQUE ENTRE LES MODELES                        *
  13. *   THERMIQUE ADVECTION ET DIFFUSION ADVECTION                     *
  14. *                                                                  *
  15. *     - Maillage 2D rectangulaire                                  *
  16. *     - Les inconnues ('T' et 'CO') sont bloquees a gauche         *
  17. *     - Un flux generalise est injecte au milieu du domaine        *
  18. *     - Le milieu est animé d'une vitesse de transport {Vx;Vy}     *
  19. *     - Le calcul est stationnaire                                 *
  20. *                                                                  *
  21. *   On s'interesse a la solution donnee par RESO pour la THERMIQUE *
  22. *   et pour la DIFFUSION                                           *
  23. *                                                                  *
  24. *       L'ecart entre les valeurs donnes par les deux solutions    *
  25. *   doit etre de l'ordre de 'VALE' 'PREC' (XZPREC).                *
  26. *                                                                  *
  27. *__________________________________________________________________*
  28. *
  29. OPTI DIME 2 'ELEM' 'TRI3';
  30. OPTI TRAC 'PSC';
  31.  
  32.  
  33. * PARAMETRAGE *
  34. Lx         = 1.  ;
  35. Ly         = 0.5 ;
  36. Nx         = 40  ;
  37. Ny         = 20  ;
  38.  
  39. Vx         = 10. ;
  40. Vy         = 2.  ;
  41. Val_Gauche = 0.  ;
  42. Val_Flux   = 1.  ;
  43. Alpha      = 0.1 ;
  44.  
  45.  
  46. * MAILLAGE *
  47. P1       = POIN 0. 0. ;
  48. P2       = P1 'PLUS' (Lx 0.);
  49.  
  50. L1       ='DROI' Nx P1 P2 ;
  51. L2       = L1 'PLUS' (0. Ly);
  52.  
  53. S1       ='REGL' Ny L1 L2 ;
  54. L_GAUCHE = POIN S1 DROI P1 (P1 PLUS (0. 1.)) 1.D-6 ;
  55. P_Mil    = POIN S1 PROC (P1 PLUS ((Lx/2.) (Ly/2.)));
  56.  
  57.  
  58. * Creation d'un MODELE 'THERMIQUE' 'ADVECTION' comme calcul de ref.
  59. MOD1 = MODE S1 'THERMIQUE' 'ADVECTION';
  60. MAT1 = MATE MOD1 'K' Alpha 'RHO' 1. 'C' 1. 'VITX' Vx 'VITY' Vy ;
  61. ADV1 = ADVE MOD1 MAT1 ;
  62. CON1 = COND MOD1 MAT1 ;
  63. BLO1 = BLOQ 'T' L_GAUCHE ;
  64. IMP1 = DEPI BLO1 Val_Gauche ;
  65. CHA1 = MANU 'CHPO' P_Mil 1 'Q' Val_Flux 'NATURE' 'DISCRET';
  66. RIG1 = ADV1 ET CON1 ET BLO1 ;
  67. SEC1 = CHA1 ET IMP1 ;
  68.  
  69.  
  70. * Creation d'un MODELE 'DIFFUSION' 'ADVECTION' (nouveau modele en 2022)
  71. MOD2 = MODE S1 'DIFFUSION' 'ADVECTION';
  72. MAT2 = MATE MOD2 'KD' Alpha 'CDIF' 1. 'VITX' Vx 'VITY' Vy ;
  73. ADV2 = ADVE MOD2 MAT2 ;
  74. CON2 = COND MOD2 MAT2 ;
  75. BLO2 = BLOQ 'CO' L_GAUCHE;
  76. IMP2 = DEPI BLO2 Val_Gauche ;
  77. CHA2 = MANU 'CHPO' P_Mil 1 'QCO' Val_Flux 'NATURE' 'DISCRET';
  78. RIG2 = ADV2 ET CON2 ET BLO2 ;
  79. SEC2 = CHA2 ET IMP2 ;
  80.  
  81.  
  82. * Creation d'un MODELE sans ADVECTION (pour tester ADVE avec un MODELE COMPLEXE)
  83. MOD0   = MODE  S1 'THERMIQUE' ;
  84. MAT0   = MATE  MOD0 'K' 1. 'RHO' 1. 'C' 1. ;
  85. MODTOT = MOD0 et MOD1 et MOD2 ;
  86. MATTOT = MAT0 et MAT1 et MAT2 ;
  87. RIGTOT = ADVE MODTOT MATTOT   ;
  88.  
  89.  
  90. * RESOLUTION *
  91. SOLa = RESO (RIG1 ET RIG2) (SEC1 ET SEC2) ;
  92. SolT ='EXCO' SOLa 'T'  ;
  93. SolC ='EXCO' SOLa 'CO' ;
  94.  
  95. TRAC S1 ('BORN' SOLa 'T' 'MAXI' (('MAXI' SolT) / 2.)) (CONT S1);
  96.  
  97.  
  98. * Comparaison des solutions *
  99. DIF1 ='MAXI' 'ABS' (SolT - SolC) / ('MAXI' 'ABS' SolT) ;
  100.  
  101. 'MESS' 'Ecart relatif MAXI entre les solutions THERMIQUE et DIFFUSION :' DIF1 ;
  102. 'SI' (DIF1 '>'('VALE' 'PREC'));
  103.    CHAI1 = CHAI 'Erreur dans le MODELE ADVECTION entre la THERMIQUE et la DIFFUSION';
  104.    ERRE CHAI1 ;
  105. 'FINS';
  106.  
  107. FIN;
  108.  
  109.  
  110.  

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