Télécharger la plaquette Cast3M


fr en

Télécharger burgerpsi.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  burgerpsi.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4.  
  5. NX = 30 ; NY = 30 ;
  6. NITER = 200 ;
  7. CFL = 0.8 ;
  8. GRAPH = 'N' ;
  9. opti isov suli ;
  10.  
  11. ***********************************************************
  12. * EQUATION DE CONVECTION NON-LINEAIRE Ut + div (F(U)) = 0 *
  13. * RESOLUE SOUS FORME NON-CONSERVATIVE                     *
  14. *       dU/dt + lambda . nabla U = 0 avec lambda = dF/dU  *
  15. * AVEC L'OPTION PSI (POSITIVE STREAMWISE INVARIANT)       *
  16. ***********************************************************
  17.  
  18. ********************************
  19. * PROCEDURE POUR TRACER COUPES *
  20. ********************************
  21.  
  22. DEBPROC TRCE ;
  23. ARGU M*'MAILLAGE' ;
  24. nb = nbel M ;
  25. R = (M poin 1) droi 1 (M poin 2) ;
  26. i=2 ; itma = nb-2 ;
  27. repeter bloc1 itma ;
  28. R = R et ((M poin i) droi 1 (M poin (i+1))) ;
  29. i=i+1 ;
  30. fin bloc1 ;
  31. FINPROC R ;
  32.  
  33. ***********************************************
  34. * PROCEDURE POUR CALCULER CHAMP DE VITESSE ET *
  35. * TESTER LA CONVERGENCE                       *
  36. ***********************************************
  37.  
  38. DEBPROC CALCUL ;
  39. ARGU RVX*'TABLE' ;
  40.  
  41. RV = RVX.'EQEX' ;
  42. CN = RV.INCO.'CN' ;
  43.  
  44. DD = RV.PASDETPS.'NUPASDT' ;
  45. NN = DD/5 ;
  46. LO = (DD-(5*NN)) EGA 0 ;
  47. SI ( LO ) ;
  48.  
  49. ERR = KOPS (RV.INCO.'CN') - (RV.INCO.'CNM1') ;
  50.  
  51. ELI = MAXI ERR 'ABS' ;
  52. ELI = (LOG (ELI + 1.0E-20))/(LOG 10.) ;
  53. MESSAGE 'ITER ' RV.PASDETPS.'NUPASDT' '   ERREUR LINF ' ELI ;
  54. IT = PROG RV.PASDETPS.'NUPASDT' ;
  55. ER = PROG ELI ;
  56. RV.INCO.'IT' = (RV.INCO.'IT') ET IT ;
  57. RV.INCO.'ER' = (RV.INCO.'ER') ET ER ;
  58. FINSI ;
  59.  
  60. VY = KCHT $DOMTOT 'SCAL' 'SOMMET' 'COMP' 'UY' 1.0 ; 
  61. VX = NOMC 'UX' CN ;
  62. VV = KCHT $DOMTOT 'VECT' 'SOMMET' 'COMP' (MOTS 'UX' 'UY')
  63.         (VX ET VY) ;
  64. RV.INCO.'VITESSE' = VV ;
  65.  
  66. RV.INCO.'CNM1' = KCHT $DOMTOT 'SCAL' 'SOMMET' (RV.INCO.'CN') ;
  67. as2 ama1  = 'KOPS' 'MATRIK' ;
  68. FINPROC as2 ama1 ;
  69.  
  70. ****************************************************
  71.  
  73. opti elem tri3 ;
  74. titre 'Equation de Burger' ;
  75.  
  76. *************
  77. * DIFFUSION *
  78. *************
  79. DIF = 1.0E-10 ;
  80.  
  81. ************
  82. * MAILLAGE *
  83. ************
  84.  
  85. A1 = 0.0 0.0 ;
  86. A2 = 1.0 0.0 ;
  87. A3 = 1.0 1.0 ;
  88. A4 = 0.0 1.0 ;
  89.  
  90. FBAS = A1 'DROI' NX A2 ;
  91. FDRO = A2 'DROI' NY A3 ;
  92. FHAU = A3 'DROI' NX A4 ;
  93. FGAU = A4 'DROI' NY A1 ;
  94.  
  95. PA1 = FBAS POIN 'PROC' A1 ;
  96. PA2 = FBAS POIN 'PROC' A2 ;
  97. P12 = CHAN POI1 FBAS ;
  98. MA1 = P12 ELEM 'CONTENANT' PA1 ;
  99. MA2 = P12 ELEM 'CONTENANT' PA2 ;
  100.  
  101. DOMTOT = 'DALL' FBAS FDRO FHAU FGAU 'PLAN' ;
  102.  
  103. *********************************
  104. * CREATION MODELE NAVIER-STOKES *
  105. *********************************
  106.  
  107. XDOMTOT = CHAN DOMTOT QUAF ;
  108. XFBAS   = CHAN FBAS   QUAF ;
  109. ELIM (XDOMTOT ET XFBAS) 1.E-3 ;
  110.  
  111. $DOMTOT = MODE XDOMTOT 'NAVIER_STOKES' LINE ;
  112. $FBAS   = MODE XFBAS   'NAVIER_STOKES' LINE ;
  113.  
  114. **************************************
  115. * INITIALISATION DU CHAMP DE VITESSE *
  116. **************************************
  117.  
  118. VX = KCHT $DOMTOT 'SCAL' 'SOMMET' 'COMP' 'UX' 0.0 ;
  119. VY = KCHT $DOMTOT 'SCAL' 'SOMMET' 'COMP' 'UY' 1.0 ;
  120. CHVIT=KCHT $DOMTOT 'VECT' 'SOMMET' 'COMP' (MOTS 'UX' 'UY') (VX ET VY) ;
  121.  
  122. ************************************
  123. * PROFIL DE LA SOLUTION A L'ENTREE *
  124. ************************************
  125.  
  126. XBAS = COOR 1 FBAS ;
  127. XBAS = NOMC 'UX' XBAS ;
  128. SOLUTION = 1.5 - (2.0*XBAS) ;
  129. CHP1 = KCHT $FBAS 'SCAL' 'SOMMET' 'COMP' 'UX' SOLUTION ;
  130. CHP1 = NOMC 'CN' CHP1 ;
  131.  
  132. ***********************
  133. * TABLE DE RESOLUTION *
  134. ***********************
  135.  
  136. RV = EQEX $DOMTOT 'ITMA' NITER 'ALFA' CFL
  137.         'ZONE' $DOMTOT 
  138.         'OPER' CALCUL 
  139.         'OPTI' 'PSI'
  140.         'OPER' TSCAL DIF 'VITESSE' 0.0 'INCO' 'CN'
  141.         'OPTI' 'CENTREE'
  142.         'OPER' 'DFDT' 1.   'CN' 'DELTAT'    'INCO' 'CN'
  143.         'CLIM' 'CN' 'TIMP' FBAS CHP1 
  144.         'CLIM' 'CN' 'TIMP' FGAU 1.5
  145.         'CLIM' 'CN' 'TIMP' FDRO (-0.5)  
  146.         'CLIM' 'CN' 'TIMP' MA1 (-1.5) 
  147.         'CLIM' 'CN' 'TIMP' MA2 (0.5) ;
  148.  
  149. RV.INCO = TABLE 'INCO' ;
  150. RV.INCO.'CN' = KCHT $DOMTOT 'SCAL' 'SOMMET' 0. ;
  151. RV.INCO.'VITESSE'= CHVIT ;
  152.  
  153. RV.INCO.'CNM1' = KCHT $DOMTOT 'SCAL' 'SOMMET' 1. ;
  154. RV.INCO.'IT' = PROG 1 ;
  155. RV.INCO.'ER' = PROG 0. ;
  156.  
  157. EXEC RV ;
  158.  
  159. *************************
  160. * ANALYSE DES RESULTATS *
  161. *************************
  162.  
  163. MESSAGE 'MAX = ' (MAXI RV.INCO.'CN') 'MAXTHEORIQUE = 1.5 ' ;
  164. MESSAGE 'MIN = ' (MINI RV.INCO.'CN') 'MINTHEORIQUE = -0.5 ' ;
  165.  
  166. SI ( (MINI RV.INCO.'ER') > -4.0 ) ;
  167.             ERREUR 5 ;
  168. FINSI ;
  169. SI ( (MAXI RV.INCO.'CN') > 1.50001D0 ) ;
  170.             ERREUR 5 ;
  171. FINSI ;
  172. SI ( (MINI RV.INCO.'CN') < -0.50001D0 ) ;
  173.             ERREUR 5 ;
  174. FINSI ;
  175.  
  176. SI ( 'EGA' graph 'O') ;
  177. TRAC DOMTOT ;
  178. TRAC DOMTOT (RV.INCO.'CN') (CONT DOMTOT) 14 ;
  179. P1 = DOMTOT POIN 'PROC' (0.0 0.7) ;
  180. P2 = DOMTOT POIN 'PROC' (1.0 0.7) ;
  181. P1P2 = DOMTOT POIN 'DROI' P1 P2 0.01 ;
  182. LI12 = TRCE P1P2 ;
  183. XX = COOR 1 LI12 ;
  184. EVOL1 = EVOL 'CHPO' (RV.INCO.'CN') SCAL LI12 ;
  185. EVOL2 = EVOL 'CHPO' XX SCAL LI12 ;
  186. LIX = 'EXTR' EVOL2 ORDO ;
  187. LIU = 'EXTR' EVOL1 ORDO ;
  188. EVOL3 = EVOL 'MANU' 'X' LIX 'U(X)' LIU ;
  189. EVOL3 = EVOL3 'COUL' VERT ;
  190.  
  191. EVOL4 = EVOL 'MANU' 'ITERATIONS' (RV.INCO.'IT') 'LOG|E|inf'
  192.            (RV.INCO.'ER') ;
  193. EVOL4 = EVOL4 'COUL' VERT ;
  194. FINSI ;
  195.  
  196. ****************************************************
  197.  
  198. opti elem qua4 ;
  199.  
  200. DOMTOT = 'DALL' FBAS FDRO FHAU FGAU 'PLAN' ;
  201.  
  202. *********************************
  203. * CREATION MODELE NAVIER-STOKES *
  204. *********************************
  205.  
  206. XDOMTOT = CHAN DOMTOT QUAF ;
  207. XFBAS   = CHAN FBAS   QUAF ;
  208. ELIM (XDOMTOT ET XFBAS) 1.E-3 ;
  209.  
  210. $DOMTOT = MODE XDOMTOT 'NAVIER_STOKES' LINE ;
  211. $FBAS   = MODE XFBAS   'NAVIER_STOKES' LINE ;
  212.  
  213. ************************************
  214. * PROFIL DE LA SOLUTION A L'ENTREE *
  215. ************************************
  216.  
  217. XBAS = COOR 1 FBAS ;
  218. XBAS = NOMC 'UX' XBAS ;
  219. SOLUTION = 1.5 - (2.0*XBAS) ;
  220. CHP1 = KCHT $FBAS 'SCAL' 'SOMMET' 'COMP' 'UX' SOLUTION ;
  221. CHP1 = NOMC 'CN' CHP1 ;
  222.  
  223. ***********************
  224. * TABLE DE RESOLUTION *
  225. ***********************
  226.  
  227. CFL = 0.2 ;
  228. RV = EQEX $DOMTOT 'ITMA' NITER 'ALFA' CFL
  229.         'ZONE' $DOMTOT 
  230.         'OPER' CALCUL 
  231.         'OPTI' 'PSI'
  232.         'OPER' TSCAL DIF 'VITESSE' 0.0 'INCO' 'CN'
  233.         'OPTI' 'CENTREE'
  234.         'OPER' 'DFDT' 1.   'CN' 'DELTAT'    'INCO' 'CN'
  235.         'CLIM' 'CN' 'TIMP' FBAS CHP1 
  236.         'CLIM' 'CN' 'TIMP' FGAU 1.5
  237.         'CLIM' 'CN' 'TIMP' FDRO (-0.5)  
  238.         'CLIM' 'CN' 'TIMP' MA1 (-1.5) 
  239.         'CLIM' 'CN' 'TIMP' MA2 (0.5) ;
  240.  
  241. RV.INCO = TABLE 'INCO' ;
  242. RV.INCO.'CN' = KCHT $DOMTOT 'SCAL' 'SOMMET' 0. ;
  243. RV.INCO.'VITESSE'= CHVIT ;
  244.  
  245. RV.INCO.'CNM1' = KCHT $DOMTOT 'SCAL' 'SOMMET' 1. ;
  246. RV.INCO.'IT' = PROG 1 ;
  247. RV.INCO.'ER' = PROG 0. ;
  248.  
  249. EXEC RV ;
  250.  
  251. *************************
  252. * ANALYSE DES RESULTATS *
  253. *************************
  254.  
  255. MESSAGE 'MAX = ' (MAXI RV.INCO.'CN') 'MAXTHEORIQUE = 1.5 ' ;
  256. MESSAGE 'MIN = ' (MINI RV.INCO.'CN') 'MINTHEORIQUE = -0.5 ' ;
  257.  
  258. SI ( (MINI RV.INCO.'ER') > -4.0 ) ;
  259.             ERREUR 5 ;
  260. FINSI ;
  261. SI ( (MAXI RV.INCO.'CN') > 1.50001D0 ) ;
  262.             ERREUR 5 ;
  263. FINSI ;
  264. SI ( (MINI RV.INCO.'CN') < -0.50001D0 ) ;
  265.             ERREUR 5 ;
  266. FINSI ;
  267.  
  268. ***********************
  269. * TRACE DES RESULTATS *
  270. ***********************
  271.  
  272. SI ( 'EGA' graph 'O') ;
  273. UEXACT = KCHT $DOMTOT 'SCAL' 'SOMMET' 0. ;
  274. XX YY = 'COOR' (DOMA $DOMTOT SOMMET) ;
  275. REPETER BLOC1 (NBNO DOMTOT) ;
  276. P1 = (DOMA $DOMTOT SOMMET) POIN &BLOC1 ;
  277. X1 = 'EXTR' XX 'SCAL' P1 ;
  278. Y1 = 'EXTR' YY 'SCAL' P1 ;
  279. D1 = Y1 - (0.5*X1/0.75) ;
  280. D2 = Y1 - (2.0*(1.0-X1)) ;
  281. D3 = Y1 - 1.0 - (2.0*(X1-1.0)) ;
  282. BO1 = ( D1 > 0.) ;
  283. BO2 = ( D2 > 0.) ;
  284. BO3 = ( D3 > 0.) ;
  285. SI ( BO1 ET BO3 ) ;
  286. U1 = 1.5 ;
  287. SINON ;
  288. SI ( (NON BO1) ET (NON BO2) );
  289. U1 = (1.5-(2.0*X1))/(1.-(2.0*Y1)) ;
  290. SINON ;
  291. U1 = -0.5 ;
  292. FINSI ;
  293. FINSI ;
  294. C1 = MANU 'CHPO' P1 1 SCAL U1 ;
  295. C2 = KCHT $DOMTOT 'SCAL' 'SOMMET' C1 ;
  296. UEXACT = (UEXACT ET C2) ;
  297. FIN BLOC1 ;
  298.  
  299.  
  300. TRACE DOMTOT ;
  301. TRAC DOMTOT (RV.INCO.'CN') (CONT DOMTOT) 14 ;
  302.  
  303. P1 = DOMTOT POIN 'PROC' (0.0 0.7) ;
  304. P2 = DOMTOT POIN 'PROC' (1.0 0.7) ;
  305. P1P2 = DOMTOT POIN 'DROI' P1 P2 0.01 ;
  306. LI12 = TRCE P1P2 ;
  307.  
  308. XX = COOR 1 LI12 ;
  309.  
  310. EVOL1 = EVOL 'CHPO' (RV.INCO.'CN') SCAL LI12 ;
  311. EVOL2 = EVOL 'CHPO' XX SCAL LI12 ;
  312.  
  313. LIX = 'EXTR' EVOL2 ORDO ;
  314. LIU = 'EXTR' EVOL1 ORDO ;
  315. EVOL33 = EVOL 'MANU' 'X' LIX 'U(X)' LIU ;
  316. EVOL33 = EVOL33 'COUL' 'TURQ' ;
  317.  
  318. EVOL44 = EVOL 'CHPO' UEXACT SCAL LI12 ;
  319. LIUEX = 'EXTR' EVOL44 ORDO ;
  320. EVOL55 = EVOL 'MANU' 'X' LIX 'UEXACT(X)' LIUEX ;
  321. EVOL55 = EVOL55 'COUL' 'ROUG' ;
  322.  
  323. TAB1 = TABLE ;
  324. TAB1.1 = 'MARQ TRIA' ;
  325. TAB1.2 = 'MARQ CARR' ;
  326. TAB1.'TITRE' = TABLE ;
  327. TAB1.'TITRE' . 1 = 'MOT' 'TRIANGLE' ;
  328. TAB1.'TITRE' . 2 = 'MOT' 'QUADRANGLE' ;
  329. TAB1.'TITRE' . 3 = 'MOT' 'EXACT' ;
  330. DESS (EVOL3 ET EVOL33 ET EVOL55) 'XBOR' 0.0 1.0 'YBOR' -1. 2.
  331.                 'TITR' 'Coupe a y=0.7' LEGE TAB1 ; 
  332.  
  333. EVOL6 = EVOL 'MANU' 'ITERATIONS' (RV.INCO.'IT') 'LOG|E|inf'
  334.            (RV.INCO.'ER') ;
  335. EVOL6 = EVOL6 COUL TURQ ;
  336. DESS (EVOL4 ET EVOL6) 'XBOR' 0. 1000. 'YBOR' -10.0 0.0
  337.  LEGE TAB1 ;
  338.  
  339. FINSI ;
  340.  
  341. FIN ;
  342.  
  343.  
  344.  
  345.  
  346.  
  347.  
  348.  
  349.  
  350.  

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