Datasets examples

* fichier :  prodt.dgibi
* Teste la procedure PRODT production d'energie turbulente
*

opti dime 3 elem cub8 ;

DISCR = 'MACRO' ;

a = 1.13 ;
b = 2.16 ;
n1= 2;
n2= 3;
n3= 5;
P0 = 0 0 0; 
p1 = a 0 0;
p2 = a b 0;
p3 = 0 b 0;
cmt1= p0 d n1 p1 d n2 p2 d n1 p3 d n2 p0;
mt1 = surf 'PLAN' cmt1 ;
mt= mt1 volu n3 'TRAN' (0 0 3.51);
Mmt= chan mt QUAF ;
$mt= mode Mmt 'NAVIER_STOKES' DISCR;
mt=doma $mt maillage ;
x=coor 1 mt ;
y=coor 2 mt ;
z=coor 3 mt ;

u1= x + (2.*y) + (3.*z);
u2= (4.*x) + (5.*y) + (6.*z);
u3= (7.*x) + (8.*y) + (9.*z);
u=(nomc 'UX' u1) et (nomc 'UY' u2) et (nomc 'UZ' u3);
t=(8*x) + (16*y) + (32*z) ;


** Cas ou dt/dz >0  est de meme signe que gb
** Stratification instable => G=0.
gb=(1 1 1);
P=prodt u $mt gb t ;
mi= mini P ;
ma= maxi P ;
d=abs ( mi + ma - 1092.) ;
mess mi ma 'Solution : 546  d=' d ;
Si(d > 1.e-10);
erreur 5 ;
Finsi ;


** Cas ou dt/dz >0  est de signe contraire a gb (<0)
** Stratification stable => (P + G)(1 + C3*Rf)
** Rf=-G /(P+G)
gb=(1 1 1)*(-1.);
P=prodt u $mt gb t ;
mi= mini P ;
ma= maxi P ;
d=abs ( mi + ma - 1060.) ;
mess mi ma 'Solution : 530  d=' d ;
Si(d > 1.e-10);
erreur 5 ;
Finsi ;
FIN ;