Accueil | Cast3M

* fichier :  ccar1.dgibi
************************************************************************
************************************************************************
**
**  --- 10 OCTOBRE 1997 ---
**
**   TEST CAVITE CARRE
**
**   teste NS(SUPG) DUDW en IMPL
**   Formulation QUADR P0 et P1
 
GRAPH = 'N'                   ;
COMPLET = FAUX ;
 
SI ( COMPLET ) ;
ds1=0.02 ; ds2=0.2 ;
ds1=0.02 ; ds2=0.1 ;
SINON        ;
err1=5.e-5;
ds1=0.04 ; ds2=0.3 ;
FINSI ;
 
 
 
 DEBPROC TEST KPRESS*MOT TYPK*MOT MACRO*MOT GRAPH*MOT ;
 option dime 2 elem TYPK ;
 
 p1= 0 0 ; p12=0.5 0. ; p2= 1 0 ;
 
 ab=p1 d dini ds1 dfin ds2 p12 d dini ds2 dfin ds1 p2 ;
 mt= ab trans dini ds1 dfin ds2 (0 0.5) trans dini ds2
 dfin ds1 (0 0.5) ;
 bc=cote 2 mt ; cd=cote 3 mt ; da=cote 4 mt ;
 AA=bc moin (0.5 0.) ;
 BB=ab plus (0. 0.5) ;
 elim (AA et BB et mt) 1.e-3 ;
 
 mt = mt chan quaf ;
 $mt=mode mt 'NAVIER_STOKES' MACRO ;
 
EPSS=1.e-8;
MU  =1. ;
RO= 400. ;
DT=1. ;
 
 
   RV= eqex   'NITER' 5 'OMEGA' 0.9
  'OPTI' 'EF' 'IMPL' 'SUPG' KPRESS
  ZONE  $mt       OPER    NS 1. 'UN' (MU/RO)    INCO 'UN'
  ZONE  $mt       OPER    DUDW EPSS       INCO 'UN'
  CLIM
  UN UIMP CD   1. UN VIMP CD 0.  UN UIMP DA   0. UN VIMP DA 0.
  UN UIMP AB   0. UN VIMP AB 0.  UN UIMP BC   0. UN VIMP BC 0. ;
 
rv.inco= table inco ;
rv.inco.un =  kcht $mt  vect sommet (1.e-5 1.e-5) ;
rv.inco.pres =  kcht $mt  scal KPRESS 1.e-5  ;
 
 exec rv ;
 
 AA= chan AA QUAF ;
 BB= chan BB QUAF ;
 $AA=mode AA 'NAVIER_STOKES' MACRO ;
 $BB=mode BB 'NAVIER_STOKES' MACRO ;
 
srti=doma $AA 'MAILLAGE'        ;
srth=doma $BB 'MAILLAGE'        ;
evolV = EVOL 'CHPO' (rv.'INCO'.'UN') UX (srti  ) ;
evolH = EVOL 'CHPO' (rv.'INCO'.'UN') UY (srth  ) ;
evx=extr evolV 'ORDO' ;
evy=extr evolV 'ABSC' ;
evolV= evol 'MANU' 'Vitesse' evx 'Hauteur' evy ;
 list evx ;
rv.'EVOLV'=evolV ;
rv.'EVOLH'=evolH ;
 
si ('EGA' graph 'O' );
TAB1=TABLE;
TAB1.'TITRE'=TABLE ;
TAB1 . 1 ='MARQ      REGU '     ;
TAB1.'TITRE' . 1 = mot 'Composante_UX '     ;
DESS evolV 'TITX' 'R (m)' 'TITY' 'V (m/s)' LEGE  TAB1 ;
 
TAB1.'TITRE' . 1 = mot 'Composante_UY '     ;
DESS evolH 'TITX' 'R (m)' 'TITY' 'V (m/s)' LEGE  TAB1 ;
 c1=vect (rv.inco.'UN') 0.3 ux uy jaune ;
 trace c1 mt ;
 
finsi ;
 
FINPROC RV ;
 
 RV= test 'CENTREP0' 'QUA8' 'QUAF' GRAPH ;
 evv= (rv.'EVOLV');
 evh= (rv.'EVOLH');
SI ( NON COMPLET ) ;
lrr=prog 1.00000E-10 -2.22749E-02 -3.80702E-02 -5.75663E-02 -7.35660E-02
 -9.15793E-02 -9.26008E-02 -7.81172E-02 -5.32730E-02 -1.67005E-02
 -6.87561E-03 -1.41800E-03  5.91220E-02   .15384       .39620
 .65237       .99999 ;
 
 lr=extr (rv.'EVOLV') 'ABSC' ;
 ER=SOMM( abs (lr - lrr) ) ;
 mess ' Ecart sur CENTREP0 QUA8 QUADR ' er ;
*si ( er > err1 ) ; erreur 5 ; finsi ;
FINSI ;
 
 
 RV= test 'CENTREP1' 'QUA8' 'QUAF' GRAPH ;
 evv= evv et (rv.'EVOLV');
 evh= evh et (rv.'EVOLH');
SI ( NON COMPLET ) ;
lrr=prog
  1.00000E-10 -1.87211E-02 -3.55770E-02 -6.17477E-02 -8.77244E-02
  -.13311      -.17850      -.17244     -8.83844E-02 -2.12847E-02
  5.57624E-02   .10414       .17627       .26503       .47630  
   .70881       .99999 ;
 
 lr=extr (rv.'EVOLV') 'ABSC' ;
 ER=SOMM( abs (lr - lrr) ) ;
 mess ' Ecart sur CENTREP1 QUA8 QUADR ' er ;
*si ( er > err1 ) ; erreur 5 ; finsi ;
FINSI ;
 
 RV= test 'CENTREP0' 'TRI6' 'QUAF' GRAPH ;
 evv=evv et (rv.'EVOLV');
 evh= evh et (rv.'EVOLH');
SI ( NON COMPLET ) ;
lrr=prog
1.00000E-10 -1.56447E-02 -2.79912E-02 -4.64291E-02 -6.41763E-02
 -8.99376E-02  -.11453      -.10648     -6.93449E-02 -2.83652E-02
 4.67842E-03  2.71908E-02  8.75278E-02   .18690       .41224
 .67173       .99999 ;
 
 lr=extr (rv.'EVOLV') 'ABSC' ;
 ER=SOMM( abs (lr - lrr) ) ;
 mess ' Ecart sur CENTREP0 TRI6 QUADR ' er ;
*si ( er > err1 ) ; erreur 5 ; finsi ;
FINSI ;
 
 RV= test 'CENTREP1' 'TRI6' 'QUAF' GRAPH ;
 evv= evv et (rv.'EVOLV');
 evh= evh et (rv.'EVOLH');
SI ( NON COMPLET ) ;
lrr=prog
 1.00000E-10 -2.34228E-02 -4.66537E-02 -7.48354E-02  -.10855
 -.14910      -.19298      -.17228     -9.81535E-02 -1.20163E-02
7.45734E-02   .11687       .19620       .27172       .48463
.70677       .99999 ;
 
 lr=extr (rv.'EVOLV') 'ABSC' ;
 ER=SOMM( abs (lr - lrr) ) ;
 mess ' Ecart sur CENTREP1 TRI6 QUADR ' er ;
*si ( er > err1 ) ; erreur 5 ; finsi ;
FINSI ;
 
 
si ('EGA' graph 'O' );
titre 'Coupe oy à x=1/2 ' ;
TAB1=TABLE;
TAB1.'TITRE'= table ;
TAB1 . 1 ='TIRR MARQ CROI REGU ';
TAB1.'TITRE' . 1 ='P0-Q2';
TAB1 . 2 ='TIRM MARQ PLUS REGU ';
TAB1.'TITRE' . 2 ='P1-Q2 ';
TAB1 . 3 ='TIRR MARQ LOSA REGU ';
TAB1.'TITRE' . 3 ='P0-P2+bulle';
TAB1 . 4 ='TIRM MARQ CARR REGU ';
TAB1.'TITRE' . 4 ='P1-P2+bulle';
DESS EVV 'TITX' 'Oy'    LEGE  TAB1 ;
titre 'Coupe ox à y=1/2 ' ;
DESS EVH 'TITX' 'Ox'    LEGE  TAB1 ;
FINSI ;
 
 
FIN ;