Télécharger la plaquette Cast3M


fr en

Télécharger ccar3d.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  ccar3d.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. **
  5. **  --- 12 OCTOBRE 1998 ---
  6. **
  7. **   TEST CAVITE CUBIQUE
  8. **
  9. **   teste KCCT NS en 3D + le Bi CG
  10.  
  11. GRAPH = 'N'                   ;
  12. err1=5.e-3;
  13.  
  14. option dime 2 elem tri6 ;
  15. opti isov suli ;
  16. p1= 0 0 ; p12=0.5 0. ; p2= 1 0 ;
  17.  
  18. ds1=0.051  ; ds2=0.12 ;
  19. *ds1=0.01 ; ds2=0.1 ;
  20.  ds1=2.   ; ds2=2.  ;
  21.  ds1=0.25 ; ds2=0.2 ;
  22.  ds1=0.3  ; ds2=0.3 ;
  23. ab=p1 d dini ds1 dfin ds2 p12 d dini ds2 dfin ds1 p2 ;
  24. ab12= p12 d dini ds1 dfin ds2 (0.5 0.5) d dini ds2
  25. dfin ds1 (0.5 1.) ;
  26. mt= ab trans dini ds1 dfin ds2 (0 0.5) trans dini ds2
  27. dfin ds1 (0 0.5) ;
  28.  bc=cote 2 mt ;
  29.  cd=cote 3 mt ;
  30.  da=cote 4 mt ;
  31. ct=ab et bc et cd et da ;
  32. elim ct 1.e-3 ;
  33.  mt=ab bc cd da daller ;
  34.  
  35. mt=orie mt ;
  36.  
  37. opti dime 3 elem cu20;
  38. mth=mt plus (0 0 1) ;
  39. ab= ab12 plus (0 0 1) ;
  40. oeil = 10 10 100 ;
  41. nph=2 ;
  42. cav=mt volu  nph mth ;
  43. *trace oeil cav;
  44. macro= '   ' ;
  45. f1=face 1 cav ;
  46. f2=face 2 cav ;
  47. f3=face 3 cav ;
  48. cav= chan cav quaf ;
  49. f1=chan f1 quaf ;
  50. f2=chan f2 quaf ;
  51. f3=chan f3 quaf ;
  52. ab=chan ab quaf ;
  53. psup=cd trans nph (0 0 1) ;
  54. psup= chan psup quaf ;
  55. elim (cav et f1 et f2 et f3 et psup et ab ) 1.e-5 ;
  56. elim (f3 et psup) 1.e-5 ;
  57. paroi=psup diff f3 ;
  58.  
  59.  $mt=mode cav 'NAVIER_STOKES' QUAF;
  60.  $AB=mode AB 'NAVIER_STOKES'  QUAF;
  61.  
  62. MU=1. ;
  63. RO= 400. ;
  64. kpress='CENTREP1' ;
  65. prep1=doma $mt kpress;
  66. bcp=elem  prep1 POI1 (lect 1) ;
  67.  
  68.   rv= eqex 'OMEGA' 0.3  'NITER' 3
  69.   'OPTI' 'EF' 'IMPL' KPRESS 'SUPG' 'DIV2'
  70.   ZONE  $mt       OPER    KBBT (1.)         INCO 'UN' 'PRES'
  71.   ZONE  $mt       OPER    NS 1. 'UN'  (MU/RO)       INCO 'UN'        ;
  72.  
  73.   rv=eqex rv
  74.   CLIM
  75.   UN VIMP (F3 ET F1) 0. UN WIMP (F1 et F2 et F3)  0.
  76.   UN UIMP PSUP 1.  UN UIMP (PAROI ET F1) 0. PRES TIMP bcp 0.  ;
  77.  
  78.  
  79. rv.inco= table inco ;
  80. rv.inco.tn =  kcht $mt  scal sommet 5. ;
  81. rv.inco.un =  kcht $mt  vect sommet (1.e-5 1.e-5 1.e-5) ;
  82.  
  83. rv.inco.pres =  kcht $mt  scal KPRESS 1.e-5  ;
  84.  
  85.  rv.'METHINV'.TYPINV=3 ;
  86.  rv.'METHINV'.IMPINV=1 ;
  87.  rv.'METHINV'.NITMAX=1000;
  88.  rv.'METHINV'.PRECOND=3 ;
  89.  rv.'METHINV'.RESID  =1.e-10 ;
  90.  
  91.  exec rv ;
  92.  
  93.  srti=doma $AB 'MAILLAGE'        ;
  94.  evolV = EVOL 'CHPO' (rv.'INCO'.'UN') UX (srti  ) ;
  95.  evx=extr evolV 'ORDO' ;
  96.  list evx ;
  97.  evy=extr evolV 'ABSC' ;
  98.  evolV= evol 'MANU' 'Vitesse' evx 'Hauteur' evy ;
  99.  
  100.  si ('EGA' graph 'O' );
  101.  TAB1=TABLE;
  102.  TAB1.'TITRE'=TABLE ;
  103.  TAB1 . 1 ='MARQ      REGU '     ;
  104.  
  105.  TAB1.'TITRE' . 1 = mot 'Composante_UX '     ;
  106.  DESS evolV 'TITX' 'R (m)' 'TITY' 'V (m/s)' LEGE  TAB1 ;
  107.  
  108.  c1=vect (rv.inco.'UN') 0.3 ux uy uz jaune ;
  109.  trace c1 mt ;
  110.  pn=elno $mt (rv.inco.'PRES') kpress;
  111.  trace pn mt ;
  112.  finsi ;
  113.  evx=abs evx ;
  114.  lrr='PROG'  3.43000E-06 +3.74144E-02 +6.40402E-02 +6.91359E-02
  115.  +6.26670E-02 +3.46496E-02  5.43639E-02   .17102       .65700  ;
  116.  ER=SOMM( (((abs (evx - lrr))*0.125)**2.) / 1. ) ;
  117.  ER=SOMM( abs (evx - lrr) ) ;
  118.  mess ' Ecart sur CENTREP1 TRI7 QUADR ' er ;
  119.  si ( er > err1 ) ; erreur 5 ; finsi ;
  120.  FIN ;
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  
  131.  
  132.  
  133.  
  134.  
  135.  
  136.  
  137.  
  138.  
  139.  
  140.  
  141.  
  142.  
  143.  
  144.  

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