The p-version of finite element methods for steady viscoelastic flows between eccentric ro- tating cylinders

摘要

The steady-state, two-dimensional creeping flow of an upper-convected Maxwell (UCM) fluid between two eccentric cylinders, with the inner one rotating, is Computed by p-version of EVSS/SUPG finite element method. The solutions converge as The polynomial order of the approximation increases. An upper Deborah number (De) Limit attains 30 (p ≤ 5). With De increasing, The boundary layers form and develop in The stress, which match closely with those predicted by asymptotic analysis. The results Show that numerical oscillations is caused by the boundary layers of stress and can be re- Duced by increasing the polynomial order of the approximation.