Hostname: page-component-8448b6f56d-m8qmq Total loading time: 0 Render date: 2024-04-24T11:46:29.847Z Has data issue: false hasContentIssue false
  • English
  • Français

A Fifth-Order Combined Compact Difference Scheme for Stokes Flow on Polar Geometries

Part of: Partial differential equations, boundary value problems

Published online by Cambridge University Press:  31 January 2018

Dongdong He  and
Kejia Pan
Dongdong He*
Affiliation:
School of Science and Engineering, The Chinese University of Hong Kong, Shenzhen, 518172, China
Kejia Pan*
Affiliation:
School of Mathematics and Statistics, Central South University, Changsha 410083, China
*
*Corresponding author. Email addresses:hedongdong@cuhk.edu.cn (D.D. He), pankejia@hotmail.com (K.J. Pan)
*Corresponding author. Email addresses:hedongdong@cuhk.edu.cn (D.D. He), pankejia@hotmail.com (K.J. Pan)
Get access

Abstract

Incompressible flows with zero Reynolds number can be modeled by the Stokes equations. When numerically solving the Stokes flow in stream-vorticity formulation with high-order accuracy, it will be important to solve both the stream function and velocity components with the high-order accuracy simultaneously. In this work, we will develop a fifth-order spectral/combined compact difference (CCD) method for the Stokes equation in stream-vorticity formulation on the polar geometries, including a unit disk and an annular domain. We first use the truncated Fourier series to derive a coupled system of singular ordinary differential equations for the Fourier coefficients, then use a shifted grid to handle the coordinate singularity without pole condition. More importantly, a three-point CCD scheme is developed to solve the obtained system of differential equations. Numerical results are presented to show that the proposed spectral/CCD method can obtain all physical quantities in the Stokes flow, including the stream function and vorticity function as well as all velocity components, with fifth-order accuracy, which is much more accurate and efficient than low-order methods in the literature.

Keywords

Stokes flowcombined compact difference (CCD) schemetruncated Fourier seriesshifted gridcoordinate singularity

MSC classification

Secondary: 65N06: Finite difference methods
Type
Research Article
Information
East Asian Journal on Applied Mathematics , Volume 7 , Issue 4 , November 2017 , pp. 714 - 727
Copyright
Copyright © Global-Science Press 2017 

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1] Chu, P. and Fan, C., A three-point combined compact difference scheme, J. Comput. Phys. 140, 370–399 (1998).Google Scholar
[2] Chen, B.Y., He, D.D. and Pan, K.J., A linearized high-order combined compact difference scheme for multi-dimensional coupled Burgers’ equations, Numer. Math.-Theory. Me. In press.Google Scholar
[3] Denni, S.C.R. and Ng, M., Nguyen, P., Numerical solution for the steady motion of a viscous fluid inside a circular boundary using integral conditions, J. Comput. Phys. 108, 142–152 (1993).Google Scholar
[4] Ehrlich, L.W. and Gupta, M.M., Some difference schemes for the biharmonic equation, SIAM J. Numer. Anal. 12, 773–789 (1975).Google Scholar
[5] Greengard, L. and Kropinski, M.C., An integral equation approach to the incompressible Navier-Stokes equations in two dimensions, SIAM J. Sci. Comput. 20, 318–336 (1998).Google Scholar
[6] Gao, G.H. and Sun, H.W., Three-point combined compact alternating direction implicit difference schemes for two-dimensional time-fractional advection-diffusion equations, Commun. Comput. Phys. 17, 487–509 (2015).Google Scholar
[7] Gao, G.H. and Sun, H.W., Three-point combined compact difference schemes for time-fractional advection-diffusion equations with smooth solutions, J. Comput. Phys. 298, 520–538 (2015).Google Scholar
[8] He, D.D., An unconditionally stable spatial sixth-order CCD-ADI method for the two-dimensional linear telegraph equation, Numer. Algorithms 72, 1103–1117 (2016).Google Scholar
[9] He, D.D. and Pan, K.J., An unconditionally stable linearized CCD-ADI method for generalized nonlinear Schrödinger equations with variable coefficients in two and three dimensions, Comput. Math. Appl. 73, 2360–2374 (2017).Google Scholar
[10] Hanse, M.O.L. and Shen, W.Z., Vorticity-velocity formulation of the 3D Navier-Stokes equations in cylindrical co-ordinates, Int. J. Numer. Meth. Fl. 41, 29–45 (2003).Google Scholar
[11] Huang, W. and Tang, T., Pseudospectral solutions for steady motion of a viscous fluid inside a circular boundary, Appl. Numer. Math. 33, 167–173 (2000).Google Scholar
[12] Ito, K. and Qiao, Z.H., A high order compact MAC finite difference scheme for the Stokes equations: Augmented variable approach, J. Comput. Phys. 227, 8177–8190 (2008).Google Scholar
[13] Kohr, M. and Pop, I., Viscous incompressible flow for low Reynolds numbers, WIT, Southampton, Boston, 2004.Google Scholar
[14] Karageorghis, A. and Tang, T., A spectral domain decomposition approach for steady Navier-Stokes problems in circular geometries, Comput. Fluids 25, 541–549 (1996).Google Scholar
[15] Lai, M.C., A simple compact fourth-order poisson solver on polar geometry, J. Comput. Phys. 182, 337–345 (2002).Google Scholar
[16] Lai, M.C., Fourth-order finite difference scheme for the incompressible Navier-Stokes equations in a disk, Int. J. Numer. Meth. Fl. 42, 909–922 (2003).Google Scholar
[17] Lai, M.C., Lin, W.W. and Wang, W., A fast spectral/difference method without pole conditions for Poisson-type equations in cylindrical and spherical geometries, IMA J. Numer. Anal. 22, 537–548 (2002).Google Scholar
[18] Lai, M.C. and Liu, H.C., Fast direct solver for the biharmonic equation on a disk and its application to incompressible flows, Appl. Math. Comput. 164, 679–695 (2005).Google Scholar
[19] Lopez, J.M., Marques, F. and Shen, J., An efficient spectral-projection method for the Navier-Stokes equations in cylindrical geometries: II. three-dimensional cases, J. Comput. Phys. 176, 384–401 (2002).Google Scholar
[20] Lopez, J.M. and Shen, J., An efficient spectral-projection method for the Navier-Stokes equations in cylindrical geometries: I. axisymmetric cases, J. Comput. Phys. 139, 308–326 (1998).Google Scholar
[21] Lai, M.C. and Wang, W.C., Fast direct solvers for Poisson equation on 2D polar and spherical geometries, Numer. Meth. Part. D. E. 18, 56–68 (2002).Google Scholar
[22] Mahesh, K., A family of high order finite difference schemes with good spectral resolution, J. Comput. Phys. 145, 332–358 (1998).Google Scholar
[23] Mohseni, K. and Colonius, T., Numerical treatment of polar coordinate singularities, J. Comput. Phys. 157, 787–795 (2000).Google Scholar
[24] Nygård, F. and Andersson, H.I., On pragmatic parallelization of a serial Navier-Stokes solver in cylindrical coordinates, Int. J. Numer. Method. H. 22, 503–511 (2012).Google Scholar
[25] Nihei, T. and Ishii, K., A fast solver of the shallow water equations on a sphere using a combined compact difference scheme, J. Comput. Phys. 187, 639–659 (2003).Google Scholar
[26] Purcell, E.M., Life at low Reynolds number, Am. J. Phys. 45, 3–11 (1977).Google Scholar
[27] Pandit, S.K. and Karmakar, H., An efficient implicit compact streamfunction velocity formulation of two dimensional flows, J. Sci. Comput. 68, 653–688 (2016).Google Scholar
[28] Pulicani, J.P. and Ouazzani, J., A Fourier-Chebyshev pseudospectral method for solving steady 3-D Navier-Stokes and heat equations in cylindrical cavities, Comput. Fluids 20, 93–109 (1991).Google Scholar
[29] Sha, W., Nakabayashi, K. and Ueda, H., An accurate second-order approximation factorization method for time-dependent incompressible Navier-Stokes equations in spherical polar coordinates, J. Comput. Phys. 142, 47–66 (1998).Google Scholar
[30] Sun, H.W. and Li, L.Z., A CCD-ADI method for unsteady convection-diffusion equations, Comput. Phys. Commun. 185, 790–797 (2014).Google Scholar
[31] Lee, S.T., Liu, J. and Sun, H.W., Combined compact difference scheme for linear second-order partial differential equations with mixed derivative, J. Comput. Appl. Math. 264, 23–37 (2014).Google Scholar
[32] Sengupta, T.K., Lakshmanan, V. and Vijay, V., A new combined stable and dispersion relation preserving compact scheme for non-periodic problems, J. Comput. Phys. 228, 3048–3071 (2009).Google Scholar
[33] Sengupta, T.K., Vijay, V. and Bhaumik, S., Further improvement and analysis of CCD scheme: dissipation discretization and de-aliasing properties, J. Comput. Phys. 228, 6150–6168 (2009).Google Scholar
[34] Torres, D.J. and Coutsias, E.A., Pseudospectral solution of the two-dimensional Navier-Stokes equations in a disk, SIAM J. Sci. Comput. 21, 378–403 (1999).Google Scholar
[35] Tian, Z.F. and Yu, P.X., An efficient compact difference scheme for solving the streamfunction formulation of the incompressible Navier-Stokes equations, J. Comput. Phys. 230, 6404–6419 (2011).Google Scholar
[36] Verzicco, R. and Orlandi, P., A finite-difference scheme for three-dimensional incompressible flows in cylindrical coordinates, J. Comput. Phys. 123, 402–414 (1996).Google Scholar
[37] Yu, P.X. and Tian, Z.F., A compact scheme for the streamfunction-velocity formulation of the 2D steady incompressible Navier-Stokes equations in polar coordinaes, J. Sci. Comput. 56, 165–189 (2013).Google Scholar
[38] Zielinski, A.P., On trial functions applied in the generalized Trefftz method, Adv. Eng. Software 24, 147–155 (1995).Google Scholar