Efficient algorithm on a nonstaggered mesh for simulating Rayleigh-Bénard convection in a box.

Journal Article (Journal Article)

An efficient semi-implicit second-order-accurate finite-difference method is described for studying incompressible Rayleigh-Bénard convection in a box, with sidewalls that are periodic, thermally insulated, or thermally conducting. Operator-splitting and a projection method reduce the algorithm at each time step to the solution of four Helmholtz equations and one Poisson equation, and these are solved by fast direct methods. The method is numerically stable even though all field values are placed on a single nonstaggered mesh commensurate with the boundaries. The efficiency and accuracy of the method are characterized for several representative convection problems.

Full Text

Duke Authors

Cited Authors

  • Chiam, K-H; Lai, M-C; Greenside, HS

Published Date

  • August 29, 2003

Published In

Volume / Issue

  • 68 / 2 Pt 2

Start / End Page

  • 026705 -

PubMed ID

  • 14525146

Electronic International Standard Serial Number (EISSN)

  • 1550-2376

International Standard Serial Number (ISSN)

  • 1539-3755

Digital Object Identifier (DOI)

  • 10.1103/physreve.68.026705


  • eng