Hostname: page-component-cd9895bd7-jn8rn Total loading time: 0 Render date: 2024-12-23T12:51:24.931Z Has data issue: false hasContentIssue false
  • English
  • Français

Higher-Order Compact Scheme for the Incompressible Navier-Stokes Equations in Spherical Geometry

Published online by Cambridge University Press:  20 August 2015

T. V. S. Sekhar ,
B. Hema Sundar Raju  and
Y. V. S. S. Sanyasiraju
T. V. S. Sekhar*
Affiliation:
Department of Mathematics, Pondicherry Engineering College, Puducherry-605 014, India
B. Hema Sundar Raju*
Affiliation:
Department of Mathematics, Pondicherry Engineering College, Puducherry-605 014, India
Y. V. S. S. Sanyasiraju*
Affiliation:
Department of Mathematics, Indian Institute of Technology Madras, Chennai-600 036, India
*
Corresponding author.Email:[email protected]
Email:[email protected]
Email:[email protected]
Get access

Abstract

A higher-order compact scheme on the nine point 2-D stencil is developed for the steady stream-function vorticity form of the incompressible Navier-Stokes (N-S) equations in spherical polar coordinates, which was used earlier only for the cartesian and cylindrical geometries. The steady, incompressible, viscous and axially symmetric flow past a sphere is used as a model problem. The non-linearity in the N-S equations is handled in a comprehensive manner avoiding complications in calculations. The scheme is combined with the multigrid method to enhance the convergence rate. The solutions are obtained over a non-uniform grid generated using the transformation r=e? while maintaining a uniform grid in the computational plane. The superiority of the higher order compact scheme is clearly illustrated in comparison with upwind scheme and defect correction technique at high Reynolds numbers by taking a large domain. This is a pioneering effort, because for the first time, the fourth order accurate solutions for the problem of viscous flow past a sphere are presented here. The drag coefficient and surface pressures are calculated and compared with available experimental and theoretical results. It is observed that these values simulated over coarser grids using the present scheme are more accurate when compared to other conventional schemes. It has also been observed that the flow separation initially occurred at Re = 21.

Keywords

76D0535Q3565N06Fourth order compact schemeNavier-stokes equationsspherical polar coordinatesdrag coefficient
Type
Research Article
Information
Communications in Computational Physics , Volume 11 , Issue 1 , January 2012 , pp. 99 - 113
Copyright
Copyright © Global Science Press Limited 2012

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]Ghia, U., Ghia, K. N. and Shin, C. T., High-Re solutions for incompressible flow using Navier-Stokes equations and a multigrid method, J. Comput. Phys., 48 (1982), 387–411.Google Scholar
[2] G. H Juncu and Mihail, R., Numerical solution of the steady incompressible Navier-Stokes equations for the flow past a sphere by a multigrid defectcorrection technique, Int. J. Numer. Meth. Fluids, 11 (1990), 379–395.Google Scholar
[3]Sekhar, T. V. S., Sivakumar, R., Ravi, T. V. R. Kumar and Subbarayudu, K., High Reynolds number incompressible MHD flow under low Rm approximation, Int. J. Nonlinear Mech., 43 (2008), 231–240.CrossRefGoogle Scholar
[4]Baranyi, L., Computation of unsteady momentum and heat transfer from a fixed circular cylinder in laminar flow, J. Comput. Appl. Mech., 4 (2003), 13–25.Google Scholar
[5]Gupta, M. M., High accuracy solutions of incompressible Navier-Stokes equations, J. Comput. Phys., 93 (1991), 343–359.CrossRefGoogle Scholar
[6]Li, M., Tang, T. and Fornberg, B., A compact fourth-order finite difference scheme for the steady incompressible Navier-Stokes equations, Int. J. Numer. Meth. Fluids, 20 (1995), 1137–1151.Google Scholar
[7]Kalita, J. C., Dalal, D. C. and Dass, A. K., Fully compact higher order computation of steady-state natural convection in a square cavity, Phys. Rev. E, 64 (2001), 066703.CrossRefGoogle Scholar
[8]Spotz, W. F. and Carey, G. F., Higher order compact scheme for the steady stream-function vorticity equations, Int. J. Numer. Meth. Eng., 38 (1995), 3497–3512.Google Scholar
[9]Erturk, E. and Gokcol, C., Fourth-order compact formulation of Navier-Stokes equations and driven cavity flow athigh Reynolds numbers, Int. J. Numer. Meth. Fluids, 50(2006), 421–436.Google Scholar
[10]Iyengar, S. R. K. and Manohar, R. P., High order difference methods for heat equation in polar cylindrical coordinates, J. Comput. Phys., 72 (1988), 425–438.Google Scholar
[11]Jain, M. K., Jain, R. K. and Krishna, M., A fourth-order difference scheme for quasilinear Poisson equation in polar coordinates, Commun. Numer. Methods Eng., 10 (1994), 791–797.Google Scholar
[12]Lai, M. C., A simple compact fourth-order Poisson solver on polar geometry, J. Comput. Phys., 182 (2002), 337–345.Google Scholar
[13]Sanyasiraju, Y. V. S. S. and Manjula, V., Flow past an impulsively started circular cylinder using a higher-order semicompact scheme, Phys. Rev. E, 72 (2005), 016709.Google Scholar
[14]Sanyasiraju, Y. V. S. S. and Manjula, V., Fourth-order semi compact scheme for flow past a rotating and translating cylinder, J. Sci. Comput., 30 (2007), 389–407.CrossRefGoogle Scholar
[15]Kalita, J. C. and Ray, R. K., A transformation-free HOC scheme for incompressible viscous flows past an impulsively started circular cylinder, J. Comput. Phys., 228 (2009), 5207–5236.Google Scholar
[16]Ray, R. K. and Kalita, J. C., A transformation-free HOC scheme for incompressible viscous flows on nonuniform polar grids, Int. J. Numer. Fluids, 62 (2010), 683–708.Google Scholar
[17]Briley, W. R., A numerical study of laminar seperation bubbles using Navier-stokes equations, J. Fluid Mech., 47 (1971), 713–736.Google Scholar
[18]Sekhar, T. V. S., Kumar, T. V. R. R. and Kumar, H., MHD flow past a sphere at low and moderate Reynolds numbers, Comput. Mech., 31 (2003), 437–444.Google Scholar
[19]Wesseling, P., Report NA-37, Delft University of Technology, The Netherlands, 1980.Google Scholar
[20]Goldstein, S., Concerning some solutions of the boundary layer equations in hydrodynamics, Proc. Cambridge Philos. Soc., 26 (1930), 1–30.Google Scholar
[21]Proudman, I. and Pearson, J. R. A., Expansions at small Reynolds numbers for the flow past a sphere and a circular cylinder, J. Fluid Mech., 2 (1956), 237–262.Google Scholar
[22]Chester, W., Breach, D. R. and Proudman, I., On the flow past a sphere at low Reynolds number, J. Fluid Mech., 37 (1969), 751–760.CrossRefGoogle Scholar
[23]Dennis, S. C. R. and Walker, J. D. A., Calculations of the steady flow past a sphere at low and moderate Reynolds numbers, J. Fluid Mech., 48 (1971), 771–789.Google Scholar
[24]Chang, E. J. and Maxey, M. R., Unsteady flow about a sphere at low to moderate Reynolds numbers, J. Fluid Mech., 227 (1994), 347–379.Google Scholar
[25]Le Clair, B. P., Hamielec, A. E. and Pruppracher, H. R., A numerical study of the drag on a sphere at intermediate Reynolds and Peclet numbers, J. Atmosph. Sci., 27 (1970), 308–315.Google Scholar
[26]Fornberg, B., Steady viscous flow past a sphere at high Reynolds numbers, J. Fluid Mech., 190 (1988), 471–489.Google Scholar
[27]Feng, Z. G. and Michaelides, E. E., A numerical study on the transient heat transfer from a sphere at high Reynolds and Peclet numbers, Int. J. Heat Mass Trans., 43 (2000), 219–229.Google Scholar
[28]Atefi, G. H., Niazmand, H. and Meigourpoory, M. R., Numerical analysis of 3-D flow past a stationary sphere with slip condition at low and moderate Reynolds numbers, J. Disper. Sci. Technol., 28 (2007), 591–602.Google Scholar
[29]Roos, F. W. and Willmarth, W. W., Some experimental results on sphere and disk drag, AIAA J., 9 (1971), 285–290.Google Scholar
[30]Clift, R. and Grace Bubbles, J. R., Drops and Particles, New York, Academic Press, 1978.Google Scholar
[31]Liao, S. J., An analytic approximation of the drag coefficient for the viscous flow past a sphere, Int. J. Nonlinear Mech., 37 (2002), 1–18.Google Scholar
[32]Kawaguti, M., Rep. Inst. Sci. Tokyo, 4 (1950), 154.Google Scholar
[33]Lister, M., Ph.D Thesis, London 1953.Google Scholar
[34]Dennis, S. C. R. and Walker, M. S., Forced convection from heated spheres, Aero. Res. Counc., 26 (1964), 105.Google Scholar
[35]Hamielec, A. E., Hoffman, T. W. and Ross, L. L., Numerical solutions of the Navier-Stokes equation for flow past spheres part I: viscous flow around spheres with finite radial mass efflux, AIChE J., 13 (1950), 212–219.Google Scholar
[36]Taneda, S., Experimental investigation of the wake behind a sphere at low Reynolds numbers, J. Phys. Soc. Japan, 11 (1956), 1104–1108.Google Scholar
[37]Zou, J.-F., Ren, A.-L. and Deng, J., Numerical investigations of wake and force for flow past a sphere, Acta Aerodynamica Sinica, 03 (2004).Google Scholar
[38]Pruppacher, H. R., Le Clair, B. P. and Hamielec, A. E., Some relations between drag and flow pattern of viscous flow past a sphere and a cylinder at low and intermediate Reynolds numbers, J. Fluid Mech., 44 (1970), 781–799.Google Scholar
[39]Lee, S., A numerical study of the unsteady wake behind a sphere in a uniform flow at moderate Reynolds numbers, Comput. Fluids, 29 (2000), 639–667.Google Scholar