Simulation of incompressible viscous flow using finite element method

¹ Department of Master of Civil Engineering, Universitas Tarumanagara, Jl. Letjen S. Parman No. 1, Jakarta, Indonesia
² Department of Civil Engineering, Pelita Harapan University, Tangerang, Indonesia
³ Department of Master of Civil Engineering, Pelita Harapan University, Tangerang, Indonesia
⁴ Department of Mechanical and Control Engineering, University of Electro-Communications, Tokyo, Japan

^* Corresponding author: watip@ft.untar.ac.id

Abstract

This study focused on simulating incompressible viscous flow using the finite element method. This study used velocity and pressure as unknowns known as primitive variable formulations. Simulation of incompressible fluid flow poses numerical challenges due to the presence of nonlinear convective terms in Navier-Stokes equations and the incompressible nature of the fluid. If the connection between velocities and pressure is not discretized correctly, the stable and convergent velocities might be gained, but the obtained pressure will be oscillatory. To avoid these difficulties, continuous quadratic and additional cubic bubble functions will be used for the velocity field and linear functions for the pressure field. This kind of discretization satisfies the Ladyzhenskaya-Babuška-Brezzi (LBB) stability condition. Two cases of different Reynolds numbers were used to test the formulation's effectiveness. In the case of Reynolds number 0.12, no vortices were formed, suggesting that the flow is primarily governed by fluid friction, and fluid inertia has minimal effect. In the case of Reynolds number 120, the vortex formation, which is known as Von Kármán vortex street, appeared. These results concluded that the formulation using the finite element method is correct.