Télécharger tbsrc1.dgibi

Retour à la liste

Numérotation des lignes :

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

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