Télécharger la plaquette Cast3M


fr en

Télécharger adve_02.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : adve_02.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. * Cas-test de l'operateur ADVEction dans la formulation THERMIQUE      *
  5. *                                                                      *
  6. * Ce cas-test verifie que le produit de la rigidite d'advection avec un*
  7. * champ de temperature donne : v.gradT est egal a la solution attendue.*
  8. *                                                                      *
  9. * Attention, il s'agit du champ solution integre sur l'element fini    *
  10. * et non pas la solution analytique.                                   *
  11. *                                                                      *
  12. *----------------------------------------------------------------------*
  13. 'OPTI' 'DIME' 2 'ELEM' 'QUA4' ;
  14. *                                                                      *
  15. * Commentez cette ligne pour voir les traces :                         *
  16. 'OPTI' 'TRAC' 'PSC' ;
  17. *                                                                      *
  18. O1       = 0 0 ;
  19. X1       = 1 0 ;
  20. Y1       = 0 1 ;
  21. *                                                                      *
  22. S1       =(O1 'DROI' 1 X1) 'TRAN' 1 Y1 ;
  23. L1       = S1 'COTE' 1 'COUL' 'ROUG' ;
  24. L3       = S1 'COTE' 3 ;
  25. L4       = S1 'COTE' 4 'COUL' 'ROUG' ;
  26. 'TITR' ' Maillage : carre de cote 1...' ;
  27. 'TRAC' 'QUAL' S1 ;
  28. *                                                                      *
  29. mo1      = 'MODE' S1 thermique advection ;
  30. ma1      = 'MATE' mo1 'K' 1. 'RHO' 1.5 'C' 1. 'VITX' 0.22 'VITY' 1.0 ;
  31. *                                                                      *
  32. cht1     = (s1 'COOR' 1) + (S1 'COOR' 2) 'NOMC' T ;
  33. *                                                                      *
  34. KA1      = 'ADVE' mo1 ma1 ;
  35. chq1     = ka1 * cht1 ;
  36. *                                                                      *
  37. chqref   = 'MANU' 'CHPO' s1 1 'Q' 0.305 'NATURE' 'DISCRET' ;
  38. *                                                                      *
  39. err1     =  'MAXI' 'ABS' (chq1 - chqref / chqref) ;
  40. *                                                                      *
  41. 'OPTI' 'ECHO' 0 ;
  42. 'MESS' ;
  44. ' > Ecart relatif entre champs calcule et reference =' err1 ;
  45. 'MESS' ;
  46. 'OPTI' 'ECHO' 1 ;
  47. *                                                                      *
  48. 'SI' (err1 > 1.e-5) ;
  49.   'ERRE' 5 ;
  50. 'SINO' ;
  51.   'OPTI' 'ECHO' 0 ;
  52.   'MESS' ;
  53.   'MESS' ' > Test reussi ! ' ;
  54.   'MESS' ;
  55. 'OPTI' 'ECHO' 1 ;
  56.  
  57. 'FINS' ;
  58. *                                                                      *
  59. 'FIN' ;
  60.  
  61.  
  62.  

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