Télécharger tubesrc.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : tubesrc.dgibi
  2. *$$$ TUBESRC
  3.  
  4.  
  5. ** Exemple TUBESRC
  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. ************** QUADRATIQUE CENTREP1 SRCE=-0.1 *************
  106. nv=15 ; nh= 3 ;
  107. MACRO= 'QUAF' ;
  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-5 ;
  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-5 ;
  125. dq1=dq - 7.85390E-02 ;
  126. dq1=abs dq1 ;
  127. mess 'DQ1=' DQ1 ;
  128. si ( dq1 > err2 ) ; erreur 5 ; finsi ;
  129.  
  130. ************** QUADRATIQUE CENTREP0 SRCE=-0.1 ******************
  131. nv=15 ; nh= 3 ;
  132. MACRO= 'QUAF';
  133. kpress='CENTREP0';
  134. typelt=qua8 ;
  135. src= 0.1 ;
  136.  
  137. RV qe qs vt= TUBESRC typelt nh nv graph kpress macro src ;
  138.  
  139. dq= (abs qe ) - qs ;
  140. mess ' qe ' qe ' qs ' qs ' dq ' dq ;
  141. mess 'VT= ' vt ' src= ' src ' VT*src= ' (vt*src) ;
  142.  
  143. erq= abs (dq-(vt*src));
  144. err1=8.e-5 ;
  145. mess ' Erreur sur le debit : dq-vt*src ' (dq-(vt*src)) ;
  146. * Attention sur les bilans effet du parametre de relaxation 0.9
  147. si ( erq > err1 ) ; erreur 5 ; finsi ;
  148.  
  149. err2 =1.e-5 ;
  150. dq1=dq - 7.85390E-02 ;
  151. dq1=abs dq1 ;
  152. mess 'DQ1=' DQ1 ;
  153. si ( dq1 > err2 ) ; erreur 5 ; finsi ;
  154.  
  155. ************** QUADRATIQUE CENTREP1 SRCE= 0. ******************
  156. nv=15 ; nh= 3 ;
  157. MACRO= 'QUAF';
  158. kpress='CENTREP1';
  159. typelt=qua8 ;
  160. src=0. ;
  161.  
  162. RV qe qs vt= TUBESRC typelt nh nv graph kpress macro src ;
  163.  
  164. dq= (abs qe ) - qs ;
  165. mess ' qe ' qe ' qs ' qs ' dq ' dq ;
  166. mess 'VT= ' vt ' src= ' src ' VT*src= ' (vt*src) ;
  167.  
  168. erq= abs (dq-(vt*src));
  169. err1=8.e-5 ;
  170. mess ' Erreur sur le debit : dq-vt*src ' (dq-(vt*src)) ;
  171. * Attention sur les bilans effet du parametre de relaxation 0.9
  172. si ( erq > err1 ) ; erreur 5 ; finsi ;
  173.  
  174. err2 =1.e-5 ;
  175. dq1=dq ;
  176. dq1=abs dq1 ;
  177. mess 'DQ1=' DQ1 ;
  178. si ( dq1 > err2 ) ; erreur 5 ; finsi ;
  179.  
  180. ************** QUADRATIQUE CENTREP1 SRCE= 0.1 TRI7 *************
  181. nv=15 ; nh= 3 ;
  182. MACRO= 'QUAF';
  183. kpress='CENTREP1';
  184. typelt=TRI6 ;
  185. src=-0.1;
  186.  
  187. RV qe qs vt= TUBESRC typelt nh nv graph kpress macro src ;
  188.  
  189. dq= (abs qe ) - qs ;
  190. mess ' qe ' qe ' qs ' qs ' dq ' dq ;
  191. mess 'VT= ' vt ' src= ' src ' VT*src= ' (vt*src) ;
  192.  
  193. erq= abs (dq-(vt*src));
  194. err1=8.e-5 ;
  195. mess ' Erreur sur le debit : dq-vt*src ' (dq-(vt*src)) ;
  196. * Attention sur les bilans effet du parametre de relaxation 0.9
  197. si ( erq > err1 ) ; erreur 5 ; finsi ;
  198.  
  199. err2 =1.e-5 ;
  200. dq1=dq + 7.85390E-02 ;
  201. dq1=abs dq1 ;
  202. mess 'DQ1=' DQ1 ;
  203. si ( dq1 > err2 ) ; erreur 5 ; finsi ;
  204.  
  205. FIN ;
  206.  
  207.  
  208.  
  209.  
  210.  
  211.  
  212.  
  213.  
  214.  
  215.  
  216.  
  217.  
  218.  
  219.  
  220.  

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