Télécharger hy4.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : hy4.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. *$$$ HY4
  5.  
  6.  
  7. ** Exemple HY4
  8. ** ___________
  9. **
  10. ** --- 2 FEVRIER 2002 ---
  11. **
  12. ** CANAL LONGUEUR 3 X 10. LARGEUR 1.
  13. ** test cas isotherme NS et FROT (Faisceau de tube)
  14. ** Sur le tronçon du milieu on impose une forte perte de charge
  15. ** qui a pour effet de redresser le champ de vitesse
  16. ** Les conditions limites (adhérence en paroi) conviennent pour
  17. ** la zone Poiseuille mais physiquement ne sont pas adaptée pour
  18. ** le tronçon faisceau de tube. Il faudrait imposer un glissement.
  19. ** Mais ça marche !
  20.  
  21. GRAPH=VRAI ;
  22. GRAPH=FAUX ;
  23. DISCR=QUAF ;
  24. DISCR=LINE ;
  25. KPRES=MSOMMET ;
  26.  
  27. type=mot tri6 ;
  28. type=mot QUA8 ;
  29. *mess 'Type d élément TRI6 ou QUA8 ? ' ;
  30. *obtenir TYPE*mot ;
  31.  
  32. nbe=7 ; nbv=10 ;
  33.  
  34. p1=0 0.;
  35. p2=1. 0.;
  36. entree= p1 d nbe p2 ;
  37. c1= 5. 0. ;
  38.  
  39. ikas=0 ;
  40. *mess 'ikas= 0 --> DROIT (option par defaut) ikas=1 --> COURBE ? ';
  41. *obtenir ikas*entier ;
  42.  
  43. si (EGA ikas 1) ;
  44.  
  45. q1=p1 tour c1 -90. ;
  46. q2=p2 tour c1 -90. ;
  47. pp1=p1 c c1 q1 nbv;
  48. pp2=p2 c c1 q2 nbv;
  49. sortie=entree tour c1 -90 ;
  50.  
  51. sinon ;
  52.  
  53. q1=p1 plus (0 10) ;
  54. q2=p2 plus (0 10) ;
  55. pp1=p1 d q1 nbv;
  56. pp2=p2 d q2 nbv;
  57. sortie=entree plus (0 10) ;
  58. finsi ;
  59.  
  60. pp1=inve pp1 ;
  61. sortie=inve sortie;
  62. elim 0.0001 (sortie et pp1 et pp2 et entree );
  63. cnt=entree et pp2 et sortie et pp1 ;
  64. *bell=surf cnt ;
  65. bell=daller entree pp2 sortie pp1 ;
  66. angle=0. ;
  67. *obtenir angle*flottant ;
  68.  
  69. bell=bell tour c1 angle ;
  70. entree=entree tour c1 angle ;
  71. sortie=sortie tour c1 angle ;
  72. pp1=pp1 tour c1 angle ;
  73. pp2=pp2 tour c1 angle ;
  74.  
  75. elim (pp1 et pp2 et entree et sortie et bell) 0.0001 ;
  76.  
  77.  
  78. * On met bout a bout 3 longueur de troncon
  79.  
  80. bell_2 sortie_2 pp1_2 pp2_2 = bell sortie pp1 pp2 'PLUS' (0. 10.) ;
  81. bell_3 sortie_3 pp1_3 pp2_3 = bell_2 sortie_2 pp1_2 pp2_2
  82. 'PLUS' (0. 10.) ;
  83.  
  84. bell = bell et bell_2 et bell_3;
  85. pp1 = pp1 et pp1_2 et pp1_3 ;
  86. pp2 = pp2 et pp2_2 et pp2_3 ;
  87. sortie = sortie_3 ;
  88.  
  89. elim (pp1 et pp2 et entree et sortie et bell) 0.0001 ;
  90.  
  91. Mbell = chan bell QUAF ;
  92. MbellM= chan bell_2 QUAF ;
  93. Mentree = chan entree QUAF ;
  94. Msortie = chan sortie QUAF ;
  95. Mpp1 = chan pp1 QUAF ;
  96. Mpp2 = chan pp2 QUAF ;
  97.  
  98. Elim (Mbell et Mentree et Msortie et Mpp1 et Mpp2 et MbellM) 1.e-5 ;
  99.  
  100. $bell=MODE Mbell 'NAVIER_STOKES' DISCR ;
  101. $bellM=MODE MbellM 'NAVIER_STOKES' DISCR ;
  102. $entree = mode Mentree 'NAVIER_STOKES' DISCR ;
  103. $sortie = mode Msortie 'NAVIER_STOKES' DISCR ;
  104. $pp1 = mode Mpp1 'NAVIER_STOKES' DISCR ;
  105. $pp2 = mode Mpp2 'NAVIER_STOKES' DISCR ;
  106. PP1 = DOMA $PP1 MAILLAGE ;
  107. PP2 = DOMA $PP2 MAILLAGE ;
  108. pp1p = DOMA $PP1 MMAIL ;
  109. pp2p = DOMA $PP2 MMAIL ;
  110. ENTREE = DOMA $ENTREE MAILLAGE ;
  111. SORTIE = DOMA $SORTIE MAILLAGE ;
  112. SORTIP = DOMA $SORTIE MSOMMET ;
  113.  
  114.  
  115. mu=5.E-2 ;
  116. ro=1 ;
  117. nu=mu/ro ;
  118.  
  119. ck= 1.E+5 1.E+2;
  120. cb= 2. 2. ;
  121.  
  122. DT = 1. ;
  123.  
  124. rv=eqex $bell 'DUMP' ALFA 1. ITMA 10
  125. OPTI 'EF' 'IMPL'
  126. ZONE $BELLM OPER FROT CK CB UN INCO 'UN'
  127. ZONE $BELL OPER NS 1. 'UN' NU INCO 'UN'
  128. OPTI 'EF' 'IMPL' 'CENTREE'
  129. ZONE $BELL OPER DFDT 1. 'UNM' DT INCO 'UN'
  130. CLIM
  131. UN UIMP (PP1 ET PP2 et entree) 0.
  132. UN VIMP entree 1. UN VIMP (PP1 ET PP2) 0.
  133. ;
  134.  
  135. RVP= EQEX
  136. 'OPTI' 'EF' KPRES
  137. ZONE $bell OPER KBBT (-1.) INCO 'UN' 'PRES'
  138. CLIM PRES TIMP sortip 0.
  139. ;
  140.  
  141. RV.'PROJ'= RVP ;
  142.  
  143.  
  144. rv.'INCO'=table 'INCO' ;
  145. rv.'INCO'.'UN' = kcht $bell VECT SOMMET (0 1) ;
  146. rv.'INCO'.'UNM' = kcht $bell VECT SOMMET (0 1) ;
  147. rv.inco.'PRES'= kcht $bell scal kpres 0.;
  148.  
  149. lh= (noeu 10) et (noeu 20) et (noeu 30) et (noeu 40) et
  150. (noeu 50) et (noeu 60) ;
  151. lj= (manu poi1 ((doma $bell maillage) poin proc( 0.5 0.5) ) ) ;
  152. his = khis 'UN' 1 lh 'UN' 2 (lh et lj) ;
  153. rv.'HIST'=his ;
  154. exec rv ;
  155.  
  156. *dessin (his.'TABD') (his.'1UN') ;
  157. *dessin (his.'TABD') (his.'2UN') ;
  158.  
  159. Si (EGA KPRES 'MSOMMET') ;
  160. mtp= doma $bell'MMAIL' ;
  161. p= (rv.inco.pression);
  162.  
  163. evolp1 = EVOL 'CHPO' p 'PRES' pp1p ;
  164. evolp2 = EVOL 'CHPO' p 'PRES' pp2p ;
  165.  
  166. lpr1=(extr evolp1 'ORDO') ;
  167. lpr2=(extr evolp2 'ORDO') ;
  168. list lpr1 ; list lpr2 ;
  169.  
  170. Lpr= prog
  171. 110.68 110.06 109.23 108.66 108.08
  172. 107.54 106.83 106.48 105.36 105.84
  173. 103.10 94.036 84.523 75.032 65.533
  174. 56.035 46.538 37.039 27.546 18.038
  175. 8.8027 4.5193 4.5952 4.0910 3.4839
  176. 2.8817 2.2759 1.7034 1.0727 .55958
  177. 0.00000E+00;
  178.  
  179. ERPR1=SOMM( abs (lPR - lpr1 ) ) /100. ;
  180. ERPR2=SOMM( abs (lPR - lpr2 ) ) /100. ;
  181.  
  182. Si GRAPH ;
  183. dess (evolp1 et evolp2 ) ;
  184. trace p mtp (cont mtp) ;
  185. FINSI ;
  186.  
  187. Sinon;
  188. p=elno $mt (kcht $mt scal kpres
  189. (exco (rv.inco.pression) 'PRES')) kpres;
  190. trace p mt TITRE ctit;
  191. finsi ;
  192.  
  193.  
  194. ung1 = vect 0.5 (rv.'INCO'.'UN') ux uy jaune ;
  195.  
  196. qe=dbit (rv.'INCO'.'UN') $entree ;
  197. qs=dbit (rv.'INCO'.'UN') $sortie ;
  198. dq=(abs qs )-(abs qe) ;
  199. mess (' Bilan : dq=') dq ;
  200.  
  201. Si GRAPH ;
  202. trace ung1 bell ;
  203. FINSI ;
  204.  
  205. Mess ' ERPR1=' ERPR1 ' ERPR2=' ERPR2 ;
  206. Si (ERPR1 '>' 2.E-2) ; erreur 5 ; Finsi ;
  207. Si (ERPR2 '>' 2.e-2) ; erreur 5 ; Finsi ;
  208.  
  209. FIN ;
  210.  
  211.  
  212.  
  213.  
  214.  
  215.  
  216.  
  217.  
  218.  
  219.  
  220.  
  221.  
  222.  

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