Télécharger la plaquette Cast3M


fr en

Télécharger consmasse.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  consmasse.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. COMPLET = FAUX ;
  5. SI ( COMPLET ) ;
  6.       NX = 40 ; NY = 20 ;
  7.       NITER = 800 ;
  8.       CFL = 0.9 ;
  9.       TOLCONV = -5.0 ;
  10.       TOLDEBI = 8.0E-4 ;
  11.       TOLSOLU = 0.8 ;
  12. SINON ;
  13.       NX = 8 ; NY = 4 ;
  14.       NITER = 200 ;
  15.       CFL = 0.9 ;
  16.       TOLCONV = -4.0 ;
  17.       TOLDEBI = 3.0E-3 ;
  18.       TOLSOLU = 1.7 ;
  19. FINSI ;
  20. GRAPH = 'N' ;
  21.  
  22. ******************************************************************
  23. * TESTE LA CONSERVATION DE LA MASSE POUR L'EQUATION SOUS FORME   *
  24. * CONSERVATIVE : dC/dt + div ( U C ) = 0 AVEC U CHAMP DE VITESSE *
  25. * A DIVERGENCE NON NULLE                                         *
  26. * COMPARAISON AVEC SOLUTION ANALYTIQUE EXACTE                    *         
  27. * ON TESTE SUCCESSIVEMENT MAILLAGE QUA4 ET MAILLAGE TRI3         *
  28. ******************************************************************
  29.  
  30. DEBPROC CALCUL ;
  31. ARGU RVX*'TABLE' ;
  32.  
  33. RV = RVX.'EQEX' ;
  34. CN = RV.INCO.'CN' ;
  35.  
  36. DD = RV.PASDETPS.'NUPASDT' ;
  37. NN = DD/5 ;
  38. LO = (DD-(5*NN)) EGA 0 ;
  39. SI ( LO ) ;
  40.  
  41. ERR = KOPS (RV.INCO.'CN') - (RV.INCO.'CNM1') ;
  42. ELI = MAXI ERR 'ABS' ;
  43. ELI = (LOG (ELI + 1.0E-20))/(LOG 10.) ;
  44. MESSAGE 'ITER ' RV.PASDETPS.'NUPASDT' '   ERREUR LINF ' ELI ;
  45. IT = PROG RV.PASDETPS.'NUPASDT' ;
  46. ER = PROG ELI ;
  47. RV.INCO.'IT' = (RV.INCO.'IT') ET IT ;
  48. RV.INCO.'ER' = (RV.INCO.'ER') ET ER ;
  49. FINSI ;
  50.  
  51. RV.INCO.'CNM1' = KCHT $DOMTOT 'SCAL' 'SOMMET' (RV.INCO.'CN') ;
  52.  
  53. as2 ama1  = 'KOPS' 'MATRIK' ;
  54. FINPROC as2 ama1 ;
  55.  
  57. opti elem qua4 ;
  58. opti isov suli ;
  59. titre 'Conservation de la masse' ;
  60.  
  61. *************
  62. * DIFFUSION *
  63. *************
  64. DIF = 1.0E-20 ;
  65.  
  66. ************
  67. * MAILLAGE *
  68. ************
  69.  
  70. A1 = 0.0 0.0 ;
  71. A2 = 2.0 0.0 ;
  72. A3 = 2.0 1.0 ;
  73. A4 = 0.0 1.0 ;
  74.  
  75. FBAS = A1 'DROI' NX A2 ;
  76. FDRO = A2 'DROI' NY A3 ;
  77. FHAU = A3 'DROI' NX A4 ;
  78. FGAU = A4 'DROI' NY A1 ;
  79.  
  80. PA1 = FBAS POIN 'PROC' A1 ;
  81. PA2 = FBAS POIN 'PROC' A2 ;
  82. P12 = CHAN POI1 FBAS ;
  83. MA1 = P12 ELEM 'CONTENANT' PA1 ;
  84. MA2 = P12 ELEM 'CONTENANT' PA2 ;
  85.  
  86. DOMTOT = 'DALL' FBAS FDRO FHAU FGAU 'PLAN' ;
  87.  
  88. ************************************
  89. * CREATION DU MODELE NAVIER_STOKES *
  90. ************************************
  91.  
  92. MDOMTOT = CHAN DOMTOT QUAF ;
  93. $DOMTOT = MODE MDOMTOT 'NAVIER_STOKES' LINE ; DOMA $DOMTOT 'IMPR' ;
  94. DOMTOT = DOMA $DOMTOT MAILLAGE ;
  95. CNT = 'CONT' DOMTOT ;
  96.  
  97. FDRO = CHAN FDRO QUAF ;
  98. FGAU = CHAN FGAU QUAF ;
  99. ELIM (MDOMTOT ET FDRO ET FGAU ) 1.E-3 ;
  100.  
  101. $FDRO = MODE FDRO 'NAVIER_STOKES' LINE ;
  102. $FGAU = MODE FGAU 'NAVIER_STOKES' LINE ;
  103. FDRO = DOMA $FDRO 'MAILLAGE' ;
  104. FGAU = DOMA $FGAU 'MAILLAGE' ;
  105.  
  106. **************************************
  107. * INITIALISATION DU CHAMP DE VITESSE *
  108. **************************************
  109.  
  110. alpha = 1. ;
  111. beta =  2. ;
  112. gamma = 1. ;
  113.  
  114. XX = COOR 1 (DOMA $DOMTOT MAILLAGE) ;
  115. YY = COOR 2 (DOMA $DOMTOT MAILLAGE) ;
  116. XX = KCHT $DOMTOT 'SCAL' 'SOMMET' XX ;
  117. YY = KCHT $DOMTOT 'SCAL' 'SOMMET' YY ;
  118. U1 = KOPS XX '*' alpha ;
  119. U2 = KOPS YY '*' beta ;
  120. U3 = gamma ;
  121. U = KOPS ( KOPS U1 '+' U2) '+' U3 ;
  122. CEXACT = KOPS ( KOPS U2 '+' U3 ) '/' U ;
  123. U = NOMC 'UX' U 'NATU' 'DISCRET' ;
  124. VX = KCHT $DOMTOT 'SCAL' 'SOMMET' 'COMP' 'UX' U ;
  125. VY = KCHT $DOMTOT 'SCAL' 'SOMMET' 'COMP' 'UY' 0.0 ;
  126. CHVIT=KCHT $DOMTOT 'VECT' 'SOMMET' 'COMP' (MOTS 'UX' 'UY') (VX ET VY);
  127.  
  128. **********************************************
  129. * TABLE EQEX CONTENANT L'EQUATION A RESOUDRE *
  130. **********************************************
  131.  
  132. RV = EQEX $DOMTOT 'ITMA' NITER 'ALFA' CFL 
  133.         'ZONE' $DOMTOT 
  134.         'OPER' CALCUL 
  135.         'OPTI' 'CONS' 'SUPG' 
  136.         'OPER' TSCAL DIF 'VITESSE' 0.0 'CN' 'INCO' 'CN'
  137.         'OPTI' 'CENTREE'
  138.         'OPER' DFDT  1. 'CN' 'DELTAT'       'INCO' 'CN'
  139.         'CLIM' 'CN' 'TIMP' FGAU 1.0 ;
  140.  
  141. RV.INCO = TABLE 'INCO' ;
  142. RV.INCO.'CN' = KCHT $DOMTOT 'SCAL' 'SOMMET' 0. ;
  143. RV.INCO.'VITESSE'= CHVIT ;
  144. RV.INCO.'CNM1' = KCHT $DOMTOT 'SCAL' 'SOMMET' 1. ;
  145. RV.INCO.'IT' = PROG 1 ;
  146. RV.INCO.'ER' = PROG 0. ;
  147. RV.INCO.'NBIT' = 0 ;
  148.  
  149. EXEC RV ;
  150.  
  151. *************************************************************
  152. * ANALYSE DES RESULTATS :                                   *
  153. * ON TESTE SUR LA CONVERGENCE VERS LA SOLUTION STATIONNAIRE *
  154. *          SUR LA PRECISION ENTRE DEBITS ENTRANT/SORTANT    *
  155. *          SUR LA PRECISION SOLUTIONS NUMERIQUE/EXACTE      *
  156. *************************************************************
  157.  
  158. GX = KOPS VX '*' RV.INCO.'CN' ;
  159. GY = KOPS VY '*' RV.INCO.'CN' ;
  160. GX = NOMC 'UX' GX 'NATU' 'DISCRET' ;
  161. GY = NOMC 'UY' GY 'NATU' 'DISCRET' ;
  162. G = KCHT $DOMTOT VECT SOMMET (GX ET GY) ;
  163. QENTREE = DBIT G $FGAU ;
  164. QSORTIE = DBIT G $FDRO ;
  165. ERQ = (QENTREE+QSORTIE)/QENTREE ;
  166. ERQ = ABS(ERQ)*100. ;
  167. MESSAGE 'ERREUR RELATIVE ENTRE DEBIT ENTRANT ET SORTANT ' 
  168.          ERQ ;
  169.  
  170. ERR = KOPS CEXACT '-' RV.INCO.'CN' ;
  171. ERR = KOPS ERR '/' CEXACT ;
  172. ERR = KOPS ERR '*' 100. ;
  173. ERR = ABS ERR ;
  174. MESSAGE 'ERREUR RELATIVE ENTRE SOL. NUMERIQUE ET EXACTE' 
  175.         (MAXI ERR) ;
  176.  
  177. SI ( (MINI RV.INCO.'ER') > TOLCONV ) ;
  178.             ERREUR 5 ;
  179. FINSI ;
  180. SI ( ERQ > TOLDEBI ) ;
  181.             ERREUR 5 ;
  182. FINSI ;
  183. SI ( (MAXI ERR) > TOLSOLU ) ;
  184.             ERREUR 5 ;
  185. FINSI ;
  186.  
  187. ***********************
  188. * TRACE DES RESULTATS *
  189. ***********************
  190.  
  191. SI ( 'EGA' graph 'O') ;
  192.  
  193. VV = VECT RV.INCO.'VITESSE' 0.1 UX UY JAUNE ;
  194. TRACE DOMTOT 'TITR' 'MAILLAGE' ;
  195. TRACE VV DOMTOT CNT 'TITR' 'CHAMP DE VITESSE' ;
  196. TRAC DOMTOT CEXACT CNT 'TITR' 'SOLUTION EXACTE' ;
  197. TRAC DOMTOT (RV.INCO.'CN') CNT 'TITR' 'SOLUTION NUMERIQUE' ;
  198.  
  199. LISO = PROG 0. 2. 4. 6. 8. 10. 12. 14. 16. 20. ;
  200. TRACE ERR DOMTOT (CONT DOMTOT) LISO 'TITR' 'ERREUR RELATIVE' ;
  201.  
  202. EVOL4 = EVOL 'MANU' 'ITERATIONS' (RV.INCO.'IT') 'LOG|E|inf'
  203.            (RV.INCO.'ER') ;
  204. DESS EVOL4 'XBOR' 0. 2000. 'YBOR' -12.0 0.0 
  205.         'TITR' 'HISTOIRE DE CONVERGENCE' ;
  206. FINSI ;
  207.  
  208. ************************************
  209. * ON TESTE MAINTENANT SUR DES TRI3 *
  210. ************************************
  211.  
  212. opti elem tri3 ;
  213. DOMTOT = 'DALL' FBAS FDRO FHAU FGAU 'PLAN' ;
  214.  
  215. ************************************
  216. * CREATION DU MODELE NAVIER_STOKES *
  217. ************************************
  218.  
  219. MDOMTOT = CHAN DOMTOT QUAF ;
  220. $DOMTOT = MODE MDOMTOT 'NAVIER_STOKES' LINE ; DOMA $DOMTOT 'IMPR' ;
  221. DOMTOT = DOMA $DOMTOT MAILLAGE ;
  222. CNT = 'CONT' DOMTOT ;
  223.  
  224. FDRO = CHAN FDRO QUAF ;
  225. FGAU = CHAN FGAU QUAF ;
  226. ELIM (MDOMTOT ET FDRO ET FGAU ) 1.E-3 ;
  227.  
  228. $FDRO = MODE FDRO 'NAVIER_STOKES' LINE ;
  229. $FGAU = MODE FGAU 'NAVIER_STOKES' LINE ;
  230.  
  231. FDRO = DOMA $FDRO 'MAILLAGE' ;
  232. FGAU = DOMA $FGAU 'MAILLAGE' ;
  233.  
  234. **************************************
  235. * INITIALISATION DU CHAMP DE VITESSE *
  236. **************************************
  237.  
  238. alpha = 1. ;
  239. beta =  2. ;
  240. gamma = 1. ;
  241.  
  242. XX = COOR 1 (DOMA $DOMTOT MAILLAGE) ;
  243. YY = COOR 2 (DOMA $DOMTOT MAILLAGE) ;
  244. XX = KCHT $DOMTOT 'SCAL' 'SOMMET' XX ;
  245. YY = KCHT $DOMTOT 'SCAL' 'SOMMET' YY ;
  246. U1 = KOPS XX '*' alpha ;
  247. U2 = KOPS YY '*' beta ;
  248. U3 = gamma ;
  249. U = KOPS ( KOPS U1 '+' U2) '+' U3 ;
  250. CEXACT = KOPS ( KOPS U2 '+' U3 ) '/' U ;
  251. U = NOMC 'UX' U 'NATU' 'DISCRET' ;
  252. VX = KCHT $DOMTOT 'SCAL' 'SOMMET' 'COMP' 'UX' U ;
  253. VY = KCHT $DOMTOT 'SCAL' 'SOMMET' 'COMP' 'UY' 0.0 ;
  254. CHVIT=KCHT $DOMTOT 'VECT' 'SOMMET' 'COMP' (MOTS 'UX' 'UY') (VX ET VY);
  255.  
  256. **********************************************
  257. * TABLE EQEX CONTENANT L'EQUATION A RESOUDRE *
  258. **********************************************
  259.  
  260. RV = EQEX $DOMTOT 'ITMA' NITER 'ALFA' CFL 
  261.         'ZONE' $DOMTOT 
  262.         'OPER' CALCUL 
  263.         'OPTI' 'CONS' 'SUPG' 
  264.         'OPER' TSCAL DIF 'VITESSE' 0.0 'CN' 'INCO' 'CN'
  265.         'OPTI' 'CENTREE'
  266.         'OPER' DFDT  1. 'CN' 'DELTAT'       'INCO' 'CN'
  267.         'CLIM' 'CN' 'TIMP' FGAU 1.0 ;
  268.  
  269. RV.INCO = TABLE 'INCO' ;
  270. RV.INCO.'CN' = KCHT $DOMTOT 'SCAL' 'SOMMET' 0. ;
  271. RV.INCO.'VITESSE'= CHVIT ;
  272. RV.INCO.'CNM1' = KCHT $DOMTOT 'SCAL' 'SOMMET' 1. ;
  273. RV.INCO.'IT' = PROG 1 ;
  274. RV.INCO.'ER' = PROG 0. ;
  275. RV.INCO.'NBIT' = 0 ;
  276.  
  277. EXEC RV ;
  278.  
  279. *************************************************************
  280. * ANALYSE DES RESULTATS :                                   *
  281. * ON TESTE SUR LA CONVERGENCE VERS LA SOLUTION STATIONNAIRE *
  282. *          SUR LA PRECISION ENTRE DEBITS ENTRANT/SORTANT    *
  283. *          SUR LA PRECISION SOLUTIONS NUMERIQUE/EXACTE      *
  284. *************************************************************
  285.  
  286. GX = KOPS VX '*' RV.INCO.'CN' ;
  287. GY = KOPS VY '*' RV.INCO.'CN' ;
  288. GX = NOMC 'UX' GX 'NATU' 'DISCRET' ;
  289. GY = NOMC 'UY' GY 'NATU' 'DISCRET' ;
  290. G = KCHT $DOMTOT VECT SOMMET (GX ET GY) ;
  291. QENTREE = DBIT G $FGAU ;
  292. QSORTIE = DBIT G $FDRO ;
  293. ERQ = (QENTREE+QSORTIE)/QENTREE ;
  294. ERQ = ABS(ERQ)*100. ;
  295. MESSAGE 'ERREUR RELATIVE ENTRE DEBIT ENTRANT ET SORTANT ' 
  296.          ERQ ;
  297.  
  298. ERR = KOPS CEXACT '-' RV.INCO.'CN' ;
  299. ERR = KOPS ERR '/' CEXACT ;
  300. ERR = KOPS ERR '*' 100. ;
  301. ERR = ABS ERR ;
  302. MESSAGE 'ERREUR RELATIVE ENTRE SOL. NUMERIQUE ET EXACTE' 
  303.         (MAXI ERR) ;
  304.  
  305. SI ( (MINI RV.INCO.'ER') > TOLCONV ) ;
  306.             ERREUR 5 ;
  307. FINSI ;
  308. SI ( ERQ > TOLDEBI ) ;
  309.             ERREUR 5 ;
  310. FINSI ;
  311. SI ( (MAXI ERR) > TOLSOLU ) ;
  312.             ERREUR 5 ;
  313. FINSI ;
  314.  
  315. ***********************
  316. * TRACE DES RESULTATS *
  317. ***********************
  318.  
  319. SI ( 'EGA' graph 'O') ;
  320.  
  321. VV = VECT RV.INCO.'VITESSE' 0.1 UX UY JAUNE ;
  322. TRACE DOMTOT 'TITR' 'MAILLAGE' ;
  323. TRACE VV DOMTOT CNT 'TITR' 'CHAMP DE VITESSE' ;
  324. TRAC DOMTOT CEXACT CNT 'TITR' 'SOLUTION EXACTE' ;
  325. TRAC DOMTOT (RV.INCO.'CN') CNT 'TITR' 'SOLUTION NUMERIQUE' ;
  326.  
  327. LISO = PROG 0. 2. 4. 6. 8. 10. 12. 14. 16. 20. ;
  328. TRACE ERR DOMTOT (CONT DOMTOT) LISO 'TITR' 'ERREUR RELATIVE' ;
  329.  
  330. EVOL4 = EVOL 'MANU' 'ITERATIONS' (RV.INCO.'IT') 'LOG|E|inf'
  331.            (RV.INCO.'ER') ;
  332. DESS EVOL4 'XBOR' 0. 2000. 'YBOR' -12.0 0.0 
  333.         'TITR' 'HISTOIRE DE CONVERGENCE' ;
  334. FINSI ;
  335.  
  336. FIN ;
  337.  
  338.  
  339.  
  340.  
  341.  
  342.  
  343.  

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