Télécharger la plaquette Cast3M


fr en

Télécharger tubesrc1.dgibi

Retour à la liste

Numérotation des lignes : 

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

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