Télécharger la plaquette Cast3M


fr en

Télécharger diffusion_advection_01.dgibi

Retour à la liste

Numérotation des lignes : 

  1.  
  2. OPTI DIME 2 'ELEM' 'TRI3';
  3. OPTI TRAC 'X';
  4.  
  5.  
  6. * PARAMETRAGE *
  7. Lx         = 1.  ;
  8. Ly         = 0.5 ;
  9. Nx         = 40  ;
  10. Ny         = 20  ;
  11.  
  12. Vx         = 10. ;
  13. Vy         = 2.  ;
  14. Val_Gauche = 0.  ;
  15. Val_Flux   = 1.  ;
  16. Alpha      = 0.1 ;
  17.  
  18.  
  19. * MAILLAGE *
  20. P1       = POIN 0. 0. ;
  21. P2       = P1 'PLUS' (Lx 0.);
  22.  
  23. L1       ='DROI' Nx P1 P2 ;
  24. L2       = L1 'PLUS' (0. Ly);
  25.  
  26. S1       ='REGL' Ny L1 L2 ;
  27. L_GAUCHE = POIN S1 DROI P1 (P1 PLUS (0. 1.)) 1.D-6 ;
  28. P_Mil    = POIN S1 PROC (P1 PLUS ((Lx/2.) (Ly/2.)));
  29.  
  30.  
  31. * Creation d'un MODELE 'THERMIQUE' 'ADVECTION' comme calcul de ref.
  32. MOD1 = MODE S1 'THERMIQUE' 'ADVECTION';
  33. MAT1 = MATE MOD1 'K' Alpha 'RHO' 1. 'C' 1. 'UX' Vx 'UY' Vy ;
  34. ADV1 = ADVE MOD1 MAT1 ;
  35. CON1 = COND MOD1 MAT1 ;
  36. BLO1 = BLOQ 'T' L_GAUCHE ;
  37. IMP1 = DEPI BLO1 Val_Gauche ;
  38. CHA1 = MANU 'CHPO' P_Mil 1 'Q' Val_Flux 'NATURE' 'DISCRET';
  39. RIG1 = ADV1 ET CON1 ET BLO1 ;
  40. SEC1 = CHA1 ET IMP1 ;
  41.  
  42.  
  43. * Creation d'un MODELE 'DIFFUSION' 'ADVECTION' (nouveau modele en 2022)
  44. MOD2 = MODE S1 'DIFFUSION' 'ADVECTION';
  45. MAT2 = MATE MOD2 'KD' Alpha 'CDIF' 1. 'UX' Vx 'UY' Vy ;
  46. ADV2 = ADVE MOD2 MAT2 ;
  47. CON2 = COND MOD2 MAT2 ;
  48. BLO2 = BLOQ 'CO' L_GAUCHE;
  49. IMP2 = DEPI BLO2 Val_Gauche ;
  50. CHA2 = MANU 'CHPO' P_Mil 1 'QCO' Val_Flux 'NATURE' 'DISCRET';
  51. RIG2 = ADV2 ET CON2 ET BLO2 ;
  52. SEC2 = CHA2 ET IMP2 ;
  53.  
  54.  
  55. * Creation d'un MODELE sans ADVECTION (pour tester ADVE avec un MODELE COMPLEXE)
  56. MOD0   = MODE  S1 'THERMIQUE' ;
  57. MAT0   = MATE  MOD0 'K' 1. 'RHO' 1. 'C' 1. ;
  58. MODTOT = MOD0 et MOD1 et MOD2 ;
  59. MATTOT = MAT0 et MAT1 et MAT2 ;
  60. RIGTOT = ADVE MODTOT MATTOT   ;
  61.  
  62.  
  63. * RESOLUTION *
  64. SOLa = RESO (RIG1 ET RIG2) (SEC1 ET SEC2) ;
  65. SolT ='EXCO' SOLa 'T'  ;
  66. SolC ='EXCO' SOLa 'CO' ;
  67.  
  68. TRAC S1 ('BORN' SOLa 'T' 'MAXI' (('MAXI' SolT) / 2.)) (CONT S1);
  69.  
  70.  
  71. * Comparaison des solutions *
  72. DIF1 ='MAXI' 'ABS' (SolT - SolC) / ('MAXI' 'ABS' SolT) ;
  73.  
  74. 'MESS' 'Ecart relatif MAXI entre les solutions THERMIQUE et DIFFUSION :' DIF1 ;
  75. 'SI' (DIF1 '>'('VALE' 'PREC'));
  76.    CHAI1 = CHAI 'Erreur dans le MODELE ADVECTION entre la THERMIQUE et la DIFFUSION';
  77.    ERRE CHAI1 ;
  78. 'FINS';
  79.  
  80. FIN;
  81.  
  82.  
  83.  

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