A boundary element method for steady viscous fluid flow using penalty function formulation

✍ Scribed by M. M. Grigoriev; A. V. Fafurin

• Publisher: John Wiley and Sons
• Year: 1997
• Tongue: English
• Weight: 392 KB
• Volume: 25
• Category: Article
• ISSN: 0271-2091
• DOI: 10.1002/(sici)1097-0363(19971030)25:8<907::aid-fld592>3.0.co;2-t

✦ Synopsis

A new boundary element method is presented for steady incompressible ¯ow at moderate and high Reynolds numbers. The whole domain is discretized into a number of eight-noded cells, for each of which the governing boundary integral equation is formulated exclusively in terms of velocities and tractions. The kernels used in this paper are the fundamental solutions of the linearized Navier±Stokes equations with arti®cial compressibility. Signi®cant attention is given to the numerical evaluation of the integrals over quadratic boundary elements as well as over quadratic quadrilateral volume cells in order to ensure a high accuracy level at high Reynolds numbers. As an illustration, square driven cavity ¯ows are considered for Reynolds numbers up to 1000. Numerical results demonstrate both the high convergence rate, even when using simple (direct) iterations, and the appropriate level of accuracy of the proposed method. Although the method yields a high level of accuracy in the primary vortex region, the secondary vortices are not properly resolved.