// Resolution de l'equation // de Stokes dans un carre // avec penalisation //Definition du maillage mesh Th=square(30,30,[-1+2.*x,-1+2.*y]); //Definition de l'espace des elements finis P1 associe //au maillage Th fespace Vh(Th,[P1,P1]); fespace Wh(Th,P1); //u et v sont des elements de Vh Vh [ux,uy],[ux1,uy1],[vx,vy]; Wh p,q,p1,deltap; real epsilon=0.00001; //Definition du probleme variationnel problem stokes([ux,uy,p],[vx,vy,q],solver=LU)= int2d(Th)( dx(ux)*dx(vx)+dy(ux)*dy(vx) +dx(uy)*dx(vy)+dy(uy)*dy(vy) +dx(p)*vx+dy(p)*vy +(dx(ux)+dy(uy))*q +epsilon*(dx(p)*dx(q)+dy(p)*dy(q)) ) +on(2,3,4,ux=0.,uy=0.)+on(1,ux=1.,uy=0.) +on(2,p=0); ; //Resolution du probleme stokes; //On affiche le resultat plot([ux,uy],wait=1); plot(p,wait=1);