TITLE 'test' VARIABLES vx(10) Xm = MOVE(x) DEFINITIONS visc=1000 R0 = 1 EQUATIONS vx: visc*div( grad( vx))=0 Xm: dt(Xm) = vx BOUNDARIES REGION 1 mesh_spacing=1 START(-5,-7) value( vx)=300*URAMP(t,t-0.01) LINE TO (6,-7) natural( vx)=0 LINE TO (6,3) value( vx)=-50*URAMP(t,t-0.01) LINE TO (-5,3) natural( vx)=0 LINE TO CLOSE REGION 2 mesh_spacing=1 visc=3000 START 'ring' (R0,-2) VELOCITY(Xm) = vx ARC(CENTER=0,-2) ANGLE=360 TO CLOSE TIME 0 TO 0.05 PLOTS FOR T = 0.001 BY 0.001 TO 0.05 contour(vx) END