TITLE 'test' VARIABLES vx vy p Xm = MOVE(x) Ym = MOVE(y) SELECT ngrid=25 DEFINITIONS visc=1000 delp=10000 dens=1000 g=9.81 v=vector( vx, vy) vm=magnitude( v) R0 = 1 natp=0 cont=dx(vx)+dy(vy) EQUATIONS vx: dx( p)- visc*div( grad( vx))=0 vy: dy( p)- visc*div( grad( vy))+dens*g=0 p: div( grad( p))=1000*(dx(vx)+dy(vy)) Xm: div(grad(Xm)) = 0 Ym: div(grad(Ym)) = 0 BOUNDARIES REGION 1 mesh_spacing=1 START(-2,-7) value( vx)=0 natural( vy)=0 natural( p)=natp VELOCITY(Xm)=0 VELOCITY(Ym)=0 LINE TO (6,-7) natural( vx)=0 value( vy)=0 value( p)=0 LINE TO (6,3) value( vx)=0 natural( vy)=0 natural( p)=natp LINE TO (-2,3) value( vx)=300 value( vy)=0 value( p)=delp LINE TO CLOSE REGION 2 mesh_spacing=1 visc=3000 dens=100 g=9.81*(100-1000)/1000 START 'ring' (R0,0) VELOCITY(Xm) = vx VELOCITY(Ym) = vy ARC(CENTER=0,0) ANGLE=360 TO CLOSE TIME 0 TO 0.01 PLOTS FOR T = 0.001 BY 0.001 TO 0.01 contour(vx) vector(v) norm contour(p) contour(vy) contour(div( v)) painted END