Télécharger la plaquette Cast3M


fr en

Télécharger test_fimp_dual3DT.dgibi

Retour à la liste

Numérotation des lignes : 

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

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