Abstract
A fourth order finite difference method is presented for the 2D unsteady viscous incompressible Boussinesq equations in vorticity-stream function formulation. The method is especially suitable for moderate to large Reynolds number flows. The momentum equation is discretized by a compact fourth order scheme with the no-slip boundary condition enforced using a local vorticity boundary condition. Fourth order long-stencil discretizations are used for the temperature transport equation with one-sided extrapolation applied near the boundary. The time stepping scheme for both equations is classical fourth order Runge–Kutta. The method is highly efficient. The main computation consists of the solution of two Poisson-like equations at each Runge–Kutta time stage for which standard FFT based fast Poisson solvers are used. An example of Lorenz flow is presented, in which the full fourth order accuracy is checked. The numerical simulation of a strong shear flow induced by a temperature jump, is resolved by two perfectly matching resolutions. Additionally, we present benchmark quality simulations of a differentially-heated cavity problem. This flow was the focus of a special session at the first MIT conference on Computational Fluid and Solid Mechanics in June 2001.
This is a preview of subscription content, log in via an institution to check access.
Similar content being viewed by others
References-
Bell, J. B., and Marcus, D. L. (1992). A Second-order projection method for variable-density flows. J. Comput. Phys. 101, 334-348
Briley, W. R. (1971). A numerical study of laminar separation bubbles using the Navier–Stokes equations. J. Fluid Mech. 47, 713-736
Christon, M. (2001). Results Summary: Special session on computational predictability of natural convection flows in enclosures, http://wotan.me.unm.edu/∼christon/mit_convection/summary/-
E, Weinan, and Liu, JG. (1996). Vorticity boundary condition for finite difference schemes. J. Comput. Phys. 124, 368-382
E, Weinan, and Liu, JG. (1996). Essentially compact schemes for unsteady viscous incompressible flows. J. Comput. Phys. 126, 122-138
E, Weinan, and Shu, Chi-Wang (1994). Small-scale structures in Boussinesq convection. Phys. Fluids 6 (1), 49-58
Gill, A. E. (1982). Atmosphere-Ocean Dynamics, Academic Press
Glowinski, R., and Pironneau, O. (1979). Numerical methods for the first biharmonic equation and for the two-dimensional Stokes problem. SIAM Rev. 21, 167-212
Henshaw, W. D., Kreiss, H. O., and Reyna, L. G. M. (1994). A fourth-order accurate difference approximation for the incompressible Navier–Stokes equations. Comput. and Fluids 23, 575-593
Hou, T. Y., and Wetton, B., Stable fourth order stream-function methods for incompressible flows with boundaries, unpublished
Johnston, H., and Krasny, R. (2001). Computational predictability of natural convection flows in enclosures: A benchmark problem. In Bathe, K. J. (ed.), Computational Fluids and Solid Mechanics (Conference Proceedings), Elsevier Science
Johnston, H., and Krasny, R. (2002). Fourth order finite difference simulation of a differentially-heated cavity. To appear in Int. J. Num. Meth. Fluids
Johnston, H., and Liu#x00AD;JG. (2002). Finite difference schemes for incompressible flow based on local pressure boundary conditions. J. Comput. Phys. 180, 120-154
Orszag, S. A., and Israeli, M. (1974). Numerical simulation of viscous incompressible flows. Ann. Rev. Fluid Mech. 6, 281-318
Quartapelle, L. (1983). Numerical Solution of the Incompressible Navier–Stokes Equations, Birkhauser, Berlin
Quéré, P. L., and Behnia, M. (1998). From onset of unsteadiness to chaos in a differentially heated cavity. J. Fluid Mech. 359, 81-107
Thom, A. (1933). The flow past circular cylinders at low speeds, Proc. Roy. Soc. A 141, 651-669
Thorpe, S. A. (1968). On standing internal gravity waves of finite amplitude. J. Fluid Mech. 32, 489-528
Thorpe, S. A. (1969). Experiments on the instability of stratified shear flows: Immiscible fluids. J. Fluid Mech. 39, 25-48
Xin, S., and Le Quéré, P. L. (2001). In Bathe, K. J. (ed.), Computational Fluids and Solid Mechanics (Conference Proceedings), Elsevier Science
Wang, C., and Liu, JG. (2002). Analysis of finite difference schemes for unsteady Navier–Stokes equations in vorticity formulation. Numer. Math. 91, 543-576
Wang, C., and Liu, JG., Fourth order convergence of compact finite difference solvers for 2D incompressible flow. Accepted for publication in Comm. in Appl. Anal.-
Wang, C., Liu, JG., and Johnston, H. Analysis of a fourth order finite difference method for incompressible Boussinesq equations. Submitted to Numer. Math.--
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Liu, JG., Wang, C. & Johnston, H. A Fourth Order Scheme for Incompressible Boussinesq Equations. Journal of Scientific Computing 18, 253–285 (2003). https://doi.org/10.1023/A:1021168924020
Issue Date:
DOI: https://doi.org/10.1023/A:1021168924020