Features of glottal flow are visualized by numerical simulation, using a Navier-Stokes solver for two-dimensional, incompressible flow operating with the MAC method. Results include pressure distributions leading to the forces acting on the vocal folds. Using a vocal fold model with translational and rotational mechanical resonators the boundaries corresponding to the vocal folds are allowed to move over the computation mesh to simulate dynamic behaviour.