SUMMARYAn implicit multigrid-driven algorithm for two-dimensional incompressible laminar viscous ows has been coupled with a solution adaptation method and a mesh movement method for boundary movement. Time-dependent calculations are performed implicitly by regarding each time step as a steady-state problem in pseudo-time. The method of artiÿcial compressibility is used to solve the ow equations. The solution mesh adaptation method performs local mesh reÿnement using an incremental Delaunay algorithm and mesh coarsening by means of edge collapse. Mesh movement is achieved by modeling the computational domain as an elastic solid and solving the equilibrium equations for the stress ÿeld. The solution adaptation method has been validated by comparison with experimental results and other computational results for low Reynolds number ow over a shedding circular cylinder. Preliminary validation of the mesh movement method has been demonstrated by a comparison with experimental results of an oscillating airfoil and with computational results for an oscillating cylinder.