Datasets examples

* fichier :  testkfpt.dgibi
opti dime 2 elem QUA4;

nx = 1000;ny = 5;

p0=0 0;
p1 = 30. 0. ;
co1= p0 d dini 0.001 dfin 1. p1;

DISCR= QUAF;
*DISCR= MACRO;

COMPLET = FAUX ;
Si (NON COMPLET);
 DISCR= LINE;
p1 = 1. 0. ;
co1= p0 d 5 p1;
FINSI;



Mco1 = chan QUAF co1;

$co1 = mode Mco1  'NAVIER_STOKES' DISCR;

co1 = doma $co1  maillage ;

YP  = 0.001;


 uet = ((coor 1 co1)*0.5+1.e-5);
 teta= (coor 1 co1);
 RO =  1.  ;
 CP =  1.  ;
 mu = 1.e-6 ;
 lb = 1.e-6;
 alpha = lb / ro*cp ;
 h =  kfpt $co1 RO mu  CP lb uet yp ;
 yplus = ro * yp * uet / mu;
 tetap = ro * cp * uet *(inve h);
 y1 = extr  (evol chpo yplus co1) 'ORDO';
 t1 = extr  (evol chpo tetap co1) 'ORDO';
 ev1 = evol 'MANU' 'Y+' y1    'TETA+' t1    ;
 pr = mu/alpha;
 Si COMPLET;
 dess ev1  TITR (chai 'Prandtl ' pr)'XBOR' 1. 1.e3 'LOGX';
 Finsi ;



 mu = 1.e-7;
  h =  kfpt $co1 RO mu  CP lb uet yp ;
 yplus = ro * yp * uet / mu;
 tetap = ro * cp * uet *(inve h);
 y1 = extr  (evol chpo yplus co1) 'ORDO';
 t1 = extr  (evol chpo tetap co1) 'ORDO';
 ev01 = evol 'MANU' 'Y+' y1    'TETA+' t1    ;
 pr = mu/alpha;


 mu = 1.e-5;
  h =  kfpt $co1 RO mu  CP lb uet yp ;
 yplus = ro * yp * uet / mu;
 tetap = ro * cp * uet *(inve h);
 y1 = extr  (evol chpo yplus co1) 'ORDO';
 t1 = extr  (evol chpo tetap co1) 'ORDO';
 ev10 = evol 'MANU' 'Y+' y1    'TETA+' t1    ;
 pr = mu/alpha;



 Si COMPLET;
 dess (ev1 et ev01 et ev10) 'XBOR' 1.  1.e3 'LOGX';
 Finsi ;


FIN ;