Télécharger la plaquette Cast3M


fr en

Télécharger tubesrc2.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  tubesrc2.dgibi
  2. **  Exemple TUBESRC2
  3. **  ________________
  4. **
  5. **  --- 2 JUIN 1998 ---
  6. **
  7. **   Tube cylindrique Rayon R0=0.25  Longueur L0=16*R0
  8. **   test cas isotherme NS,FIMP KBBT en Implicite
  9. **
  10. **   On teste aussi un puits de masse uniformément réparti.
  11. **
  12.  
  13. GRAPH = 'N' ;
  14.  
  15.  
  16. DEBPROC TUBESRC ;
  17. ARGU TYPELT*MOT NH*ENTIER NV*ENTIER  GRAPH*MOT KPRESS*MOT
  18. MACRO*MOT SRC*FLOTTANT ;
  19.  
  20. OPTION DIME 2 ELEM TYPELT ;
  21.  option mode axis ;
  22.  
  23. r0=0.25 ; L0=16*R0 ;
  24. ae=0. ;
  25. P1=R0 0 ; p2=ae 0 ; p3 = ae (3.*L0/4.)  ;
  26. P4 = p3 plus (0 (L0/4));
  27. r1=R0/2. ;
  28. P7= P4 plus (R0 0) ;
  29. *P8=R0 (3.*L0/4.) ;
  30. entree= p1 d nh p2 ;
  31. axe=p2 d nv p4 ;
  32. sortie=p4 d nh p7 ;
  33. paroi=p7 d nv p1 ;
  34.  
  35. mt= entree axe sortie paroi daller ;
  36. *trace mt ;
  37. ent = chan entree poi1 ;
  38. ent=elem ent (lect 2 pas 1 (nbel ent)) ;
  39.  
  40.  
  41. mtq=chan mt quaf ;
  42. entree= chan entree quaf ;
  43. sortie= chan sortie quaf ;
  44. elim (mt et entree et sortie )1.e-5 ;
  45.  
  46. $mt= mode mtq 'NAVIER_STOKES' MACRO ; doma $mt 'IMPR' ;
  47.  
  48. NU=1.5E-2;
  49. uE=1. ;
  50. KPRESS='CENTREP1' ;
  51.  
  52. RU=eqex 'OMEGA' 0.9 'NITER' 5
  53.   'OPTI' 'EF' 'IMPL' 'SUPG' KPRESS 'ALE'
  54.   ZONE  $mt       OPER    KBBT 1.          INCO 'UN' 'PRES'
  55.   ZONE  $mt       OPER    NS 1. 'UN' NU  'WN'   INCO 'UN'
  56.   'OPTI'  INCOD KPRESS
  57.   ZONE  $mt       OPER    FIMP src       INCO 'PRES'
  58.   CLIM
  59.  UN UIMP entree 0.  UN VIMP ent    UE
  60.  UN UIMP axe 0. UN UIMP paroi 0.
  61.  UN VIMP paroi 0.  ;
  62.  
  63. ru.inco=table 'INCO' ;
  64. ru.'INCO'.'UN'=kcht $mt vect sommet (0. 1. ) ;
  65. ru.'INCO'.'WN'=kcht $mt vect sommet (0. 0. ) ;
  66. ru.'INCO'.'PRES' =  kcht $mt  scal KPRESS 0.     ;
  67.  
  68.  exec ru ;
  69.  
  70.  $entree=mode entree 'NAVIER_STOKES' MACRO  ;
  71.  $sortie=mode sortie 'NAVIER_STOKES' MACRO  ;
  72.  
  73.  un=ru.'INCO'.'UN';
  74.  qe=dbit un $entree ;
  75.  qs=dbit un $sortie ;
  76.  
  77.  si ('EGA' graph 'O' );
  78.  
  79.  mt=doma $mt 'MAILLAGE' ;
  80.  ung1=vect un 5.e-2 ux uy jaune ;
  81.  trace ung1 mt ;
  82.  
  83.  pn=elno $mt (ru.'INCO'.'PRES') KPRESS ;
  84.  trace pn mt ;
  85.  
  86.  srti=doma $sortie 'MAILLAGE' ;
  87.  evolV = EVOL 'CHPO' (ru.'INCO'.'UN') UY (srti  ) ;
  88.  
  89.  TAB1=TABLE;
  90.  TAB1.'TITRE'=TABLE ;
  91.  TAB1 . 1 ='MARQ      REGU '     ;
  92.  TAB1.'TITRE' . 1 = mot 'Composante_UX '     ;
  93.  DESS evolV 'TITX' 'R (m)' 'TITY' 'V (m/s)' LEGE  TAB1 ;
  94.  
  95.  FINSI ;
  96.  
  97.  V=doma $mt 'VOLUME' ;
  98.  VT=somt V ;
  99.  
  100. FINPROC RV qe qs vt ;
  101.  
  102.  
  103. **************   MACRO  CENTREP1 SRCE=-0.1 ******************
  104.  nv=15 ; nh= 3  ;
  105.  MACRO= 'MACRO';
  106.  kpress='CENTREP1';
  107.  typelt=qua8 ;
  108.  src= 0.1 ;
  109.  
  110. RV qe qs vt= TUBESRC typelt nh nv graph kpress macro src ;
  111.  
  112.  dq= (abs qe ) - qs  ;
  113.  mess ' qe ' qe ' qs ' qs ' dq ' dq ;
  114.  mess 'VT= ' vt ' src= ' src  ' VT*src= ' (vt*src) ;
  115.  
  116.  erq= abs (dq-(vt*src));
  117.  err1=8.e-7 ;
  118.  mess ' Erreur sur le debit : dq-vt*src ' (dq-(vt*src)) ;
  119. * Attention sur les bilans effet du parametre de relaxation 0.9
  120. * si ( erq > err1 ) ; erreur 5 ; finsi ;
  121.  
  122.  err2 =1.e-7 ;
  123.  dq1=dq - 7.85390E-02 ;
  124.  dq1=abs dq1 ;
  125.  mess 'DQ1=' DQ1 ;
  126.  si ( dq1  > err2  ) ; erreur 5 ; finsi ;
  127.  
  128.  
  129.  FIN ;
  130.  
  131.  
  132.  
  133.  
  134.  
  135.  
  136.  
  137.  
  138.  
  139.  
  140.  

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