*commentaires opti dime 2 elem qua8 ; *taille de la maille (longueur) densite 1 ; p1=0. 1.; p2= 4. 1.; p3=4. 0.; p4=20. 0.; p5=20. 3.;p7=0. 3.; p6=20. 1.; * droite d12= droite p1 p2 'DINI' 1. 'DFIN' 0.5; d23= droite p2 p3 'DINI' 0.5 'DFIN' 0.5; d34= droite p3 p4 'DINI' 0.5 'DFIN' 0.5; d45= droite p4 p5; d57= droite p5 p7; d71= droite p7 p1; d46=droite p4 p6; d65=droite p6 p4; d62=droite p6 p2; d16=droite p1 p6; *maillage mail1= surf (d12 et d23 et d34 et d45 et d57 et d71); *trace mail1; *maillage quaf mail1f=change quaf mail1; *trace mail1f; m1ns=model mail1f navier_stokes macro; *calcul rv=eqex m1ns itma 200 opti expl zone m1ns oper ns (1./150.) inco un opti centree zone m1ns oper dfdt 1. un deltat inco un; rv.pression=eqpr m1ns zone m1ns oper pression; y71= coor 2 d71; uin=(3-y71)*(y71-1); wall= d12 et d23 et d34 et d57; rv=eqex rv clim un uimp wall 0. un vimp wall 0. un uimp d71 uin un vimp d71 0. ; *condition initiale rv.inco.un = kcht m1ns vect sommet (0. 0.); exec rv; unx=exco ux rv.inco.un; uny=exco uy rv.inco.un; *vo=unx et uny; *vect1=vect 1 vo 'UX' 'UY'; trac unx mail1; opti sauv 'back.sauv'; sauv; fin;