𝔖 Bobbio Scriptorium
✦   LIBER   ✦

Fourth-order finite difference scheme for the incompressible Navier–Stokes equations in a disk

✍ Scribed by Ming-Chih Lai


Publisher
John Wiley and Sons
Year
2003
Tongue
English
Weight
212 KB
Volume
42
Category
Article
ISSN
0271-2091

No coin nor oath required. For personal study only.

✦ Synopsis


Abstract

We develop an efficient fourth‐order finite difference method for solving the incompressible Navier–Stokes equations in the vorticity‐stream function formulation on a disk. We use the fourth‐order Runge–Kutta method for the time integration and treat both the convection and diffusion terms explicitly. Using a uniform grid with shifting a half mesh away from the origin, we avoid placing the grid point directly at the origin; thus, no pole approximation is needed. Besides, on such grid, a fourth‐order fast direct method is used to solve the Poisson equation of the stream function. By Fourier filtering the vorticity in the azimuthal direction at each time stage, we are able to increase the time step to a reasonable size. The numerical results of the accuracy test and the simulation of a vortex dipole colliding with circular wall are presented. Copyright © 2003 John Wiley & Sons, Ltd.


📜 SIMILAR VOLUMES


A fourth-order compact finite difference
✍ Zhenfu Tian; Yongbin Ge 📂 Article 📅 2003 🏛 John Wiley and Sons 🌐 English ⚖ 256 KB 👁 1 views

## Abstract A fourth‐order compact finite difference scheme on the nine‐point 2D stencil is formulated for solving the steady‐state Navier–Stokes/Boussinesq equations for two‐dimensional, incompressible fluid flow and heat transfer using the stream function–vorticity formulation. The main feature o

On a compact mixed-order finite element
✍ Morten M. T. Wang; Tony W. H. Sheu 📂 Article 📅 1997 🏛 John Wiley and Sons 🌐 English ⚖ 134 KB 👁 2 views

Our work is an extension of the previously proposed multivariant element. We assign this re®ned element as a compact mixed-order element in the sense that use of this element offers a much smaller bandwidth. The analysis is implemented on quadratic hexahedral elements with a view to analysing a thre