Télécharger la plaquette Cast3M


fr en

Télécharger test_fimp_dual2DQ.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier : test_fimp_dual2DQ.dgibi
  2. *Ce test vérifie l'égalité discrète Div U - q = 0
  3. *ou q est la discrétisation via l'opérateur FIMP de
  4. *la divergence continue.
  5. * U = (u,v,w)  avec u=x**2. , v=y**2.
  6. * q=2(x+y)
  7. *Ce cas test vérifie l'option INCOD KPRES de FIMP
  8. *avec KPRES = CENTREP1,CENTRE,MSOMMET
  9. *et DISCR = QUAF,MACRO,CENTRE
  10. *
  11. * Cas 2D quadrangle
  12. *
  13. opti dime 2 elem QUA4;
  14.  
  15. opti trace 'PSC';
  16.  
  17. GRAPH = VRAI ;
  18. GRAPH = FAUX ;
  19.  
  20. L=1. ; H = 1.;  nb = 20;
  21. *L=1. ; H = 1.;  nb = 5;
  22. *L=1. ; H = 1.; nb = 2;
  23.  
  24. p0=0. 0.;
  25. p1=L  0.;
  26. p2=0. H ;
  27.  
  28. l1= p0 d nb p1;
  29. mt = l1 trans nb p2;
  30. cnt = cont mt ;
  31. *opti elem TRI3;
  32. *mt= surf cnt 'PLAN';
  33. Mmt= chan mt 'QUAF';
  34.  
  35. lmaxerr=prog;
  36. REPETER BLOCKAS 7 ;
  37.  
  38. Si (EGA &BLOCKAS 5);
  39. DISCR='LINE';KPRES='MSOMMET';
  40. Finsi;
  41.  
  42. Si (EGA &BLOCKAS 6);
  43. DISCR='MACRO';KPRES='MSOMMET';
  44. Finsi;
  45.  
  46. Si (EGA &BLOCKAS 7);
  47. DISCR='QUAF';KPRES='MSOMMET';
  48. Finsi;
  49.  
  50. Si (EGA &BLOCKAS 1);
  51. DISCR='QUAF';KPRES='CENTREP1';
  52. Finsi;
  53.  
  54. Si (EGA &BLOCKAS 2);
  55. DISCR='MACRO';KPRES='CENTREP1';
  56. Finsi;
  57.  
  58. Si (EGA &BLOCKAS 3);
  59. DISCR='MACRO';KPRES='CENTRE';
  60. Finsi;
  61.  
  62. Si (EGA &BLOCKAS 4);
  63. DISCR='LINE';KPRES='CENTRE';
  64. Finsi;
  65.  
  68. mess ' ====================================';
  69. mess ' ====================================';
  70. mess ' DISCR=' DISCR ' KPRES=' KPRES;
  71.  
  72. $mt = model Mmt 'NAVIER_STOKES' DISCR;
  73. doma $mt 'IMPR';
  74. at= doma $mt table ;
  75. mt = doma $mt maillage ;
  76. dg = doma $mt 'XXDIAGSI';
  77.  
  78.  
  79. Si(EGA KPRES 'CENTREP1');
  80. mp1=doma $mt 'CENTREP1';
  81. mp1 = mp1 coul rouge;
  82. mtp=doma $mt 'ELTP1NC';
  83. Si GRAPH;
  84. trace (at.'QUAF' et mp1);
  85. Finsi;
  86. dgp=doma $mt 'XXCTREP1';
  87. mess (somt dgp);
  88. Finsi;
  89.  
  90. Si(EGA KPRES 'CENTRE');
  91. dgp=doma $mt 'XXVOLUM';
  92. mtp=doma $mt maillage ;
  93. mp1=doma $mt 'CENTRE';
  94. mp1 = mp1 coul rouge;
  95. Si(EGA DISCR 'MACRO');
  96. mtp=at.'MELSTB ';
  97. le=lect 1 pas 4 ((nbel mtp) - 3);
  98. mtp=elem mtp le;
  99. Finsi;
  100. Finsi;
  101.  
  102. Si(EGA KPRES 'MSOMMET');
  103. dgp= doma $mt 'XXMSOMME';
  104. mess ' Différence  entre XXINTG et XXDIAGSI' ;
  105. mess ' maxi (dgp - dg)=' (maxi (dgp - dg));
  106. mtp=doma $mt 'MMAIL' ;
  107. mp1=doma $mt 'MSOMMET';
  108. Finsi;
  109.  
  110. lt = chai '  ' DISCR ' ' KPRES ;
  111.  
  112. rv= eqex $mt
  113. opti 'EF' 'CENTREE' 'INCOD' KPRES
  114. 'ZONE' $mt 'OPER' 'FIMP' 'DSRC' INCO 'PRES'
  115. opti 'EF' 'IMPL' KPRES
  116. 'ZONE' $mt 'OPER' 'KMAB' 1. 1.e-2 INCO 'UN' 'PRES'
  117. opti 'EF' 'CENTREE' 
  118. 'ZONE' $mt 'OPER' 'LAPN' 'DSRC' INCO 'TN'
  119. opti 'EF' 'CENTREE' 'INCOD' 'SOMMET'
  120. 'ZONE' $mt 'OPER' 'FIMP' 'DSRC' INCO 'TN';
  121.  ;
  122.  
  123. rv.'INCO'=table 'INCO';
  124. rv.'INCO'.'PRES'=kcht $mt scal KPRES 1.;
  125. rv.'INCO'.'UN'=kcht $mt vect sommet (1. 1.);
  126. rv.'INCO'.'TN'=kcht $mt scal sommet 1.;
  127.  
  128. cx=1. ; cy=1.;
  129. cx=1. ; cy=0.;
  130. x= coor 1 mt;
  131. y= coor 2 mt;
  132. xp= coor 1 mp1;
  133. yp= coor 2 mp1;
  134. x2= (x*x*cx) CHAN 'ATTRIBUT' 'NATURE' 'DIFFUS';
  135. y2= (y*y*cy) CHAN 'ATTRIBUT' 'NATURE' 'DIFFUS';
  136.  
  137. dvu=2*((x*cx)+(y*cy));
  138. dvup=2*((xp*cx)+(yp*cy));
  139.  
  140. rv.'INCO'.'DSRC'=kcht $mt scal sommet dvu;
  141. *rv.'INCO'.'DSRC'=kcht $mt scal sommet 1. ;
  142.  
  143. UN= (nomc x2 '1UN') et (nomc y2 '2UN');
  144.  
  145. *-------------------- CUN -----------------------------------
  146.  cs C  = kmab (rv.'2KMAB');
  147.  
  148.  
  149.  cun = kops C '*' un ;
  150.  
  151.  cund=(nomc cun 'SCAL') *(inve dgp);
  152.  cuns = kcht $mt scal KPRES (nomc cun 'SCAL') ;
  153.  cuns=elno $mt cuns KPRES ;
  154.  cund = nomc cund 'SCAL';
  155.  dcunp = somt (dgp*(cund - dvup)*(cund - dvup));
  156.  mess 'Erreur L2 cund - dvup =' dcunp;
  157.  
  158.  Si Graph;
  159.  trace cuns mt TITR (chai ' CUNS (ELNO) ' lt);
  160.  
  161.  Si((EGA KPRES 'MSOMMET') ou (EGA KPRES 'CENTREP1')
  162.  ou ((EGA KPRES 'CENTRE') et (EGA DISCR 'MACRO'))
  163.  );
  164.  trace dvup mtp TITR (chai ' DVUP' lt);
  165.  trace cun mtp TITR (chai ' CUN' lt);
  166.  trace cund mtp TITR (chai ' D-1 *  CUN' lt);
  167.  trace (dvup - cund) mtp
  168.  TITR (chai ' (dvup - (D-1 * cun)) ' lt);
  169.  Finsi ;
  170.  Finsi ;
  171. *-------------------- CUN -----------------------------------
  172.  
  173. *-------------------- MSF -----------------------------------
  174.  msf ma = fimp (rv.'1FIMP');
  175.  msfd=(nomc msf 'SCAL') *(inve dgp);
  176.  dmsf = somt (dgp*(msfd - dvup)*(msfd - dvup));
  177.  mess 'Erreur L2 msfd - dvup =' dmsf;
  178.  
  179.  msfs = elno $mt msf kpres;
  180.  dcusf = nomc (cun - msf) 'SCAL';
  181.  dcusf = somt (dgp*dcusf*dcusf);
  182.  mess 'Erreur L2 (cun - msf) =' dcusf;
  183.  
  184.  Si Graph;
  185.  trace msfs mt TITR (chai ' MSFs (elno MSF) '  lt);
  186.  
  187.  Si((EGA KPRES 'MSOMMET') ou (EGA KPRES 'CENTREP1')
  188.  ou ((EGA KPRES 'CENTRE') et (EGA DISCR 'MACRO'))
  189.  );
  190.  trace (cun - msf) mtp TITR (chai ' Différence entre CUN et MSF ' lt);
  191.  trace msf mtp TITR (chai ' MSF' lt);
  192.  Finsi ;
  193.  Finsi ;
  194. *-------------------- MSF -----------------------------------
  195.  
  196. *-------------------- MSG -----------------------------------
  197.  msg ma = fimp (rv.'4FIMP');
  198.  msg = nomc msg 'SCAL';
  199.  msgd=msg *(inve dg);
  200.  Si Graph;
  201.  trace msg mt TITR (chai ' MSG (FIMP normal) ' lt);
  202.  trace msgd mt TITR (chai ' D-1 * MSG (FIMP normal) ' lt);
  203.  trace (dvu - msgd) mt
  204.   TITR (chai ' (dvu - (D-1 * MSG)) (FIMP normal) ' lt);
  205.  Finsi;
  206.  dmsgd = somt (dg*(msgd - dvu)*(msgd - dvu));
  207.  mess 'Erreur L2 msgd - dvu =' dmsgd;
  208. *-------------------- MSG -----------------------------------
  209.  
  210.  ms ma = lapn (rv.'3LAPN');
  211.  
  214.  mess (chai lt '  Erreurs L2 : ')
  215.  ' dcunp=' dcunp ' dmsf=' dmsf ' dcusf=' dcusf ' dmsgd='  dmsgd ;
  216.  err = (MAXI (prog dcunp dcusf dmsgd));
  217.  lmaxerr = lmaxerr et (prog err);
  218.  Mess (chai lt '  MAXI des erreurs = ') err;
  221.  
  222.  
  223. FIN BLOCKAS;
  224.  
  225.  
  226.  list lmaxerr ;
  227.  er=maxi lmaxerr;
  228.  si ( er > 2.E-3); erreur 5 ; finsi ;
  229.  
  230.  FIN ;
  231.  
  232.  
  233.  
  234.  
  235.  

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