Viscous incompressible flow simulation using penalty finite element

摘要

A new set of boundary conditions for velocity-divergence formulation of incompressible Navier-Stokes equations are derived. The boundary is characterized by traction due to friction and surface tension. The pressure being unknown is decoupled from the computation of velocity and is determined by post processing of the velocity field. Numerical results are presented for - classical lid-driven cavity flow-widely used by numerous authors due to its simple geometry and complicated flow behavior and squeezed flow between two parallel plates amenable to analytical solution and compared with benchmark and analytical results.