Efficient Algorithm on a Non-staggered Mesh for Simulating Rayleigh-Benard Convection in a Box
arXiv:nlin/0302057 · doi:10.1103/PhysRevE.68.026705
Abstract
An efficient semi-implicit second-order-accurate finite-difference method is described for studying incompressible Rayleigh-Benard 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 are solved by fast direct methods. The method is numerically stable even though all field values are placed on a single non-staggered mesh commensurate with the boundaries. The efficiency and accuracy of the method are characterized for several representative convection problems.
REVTeX, 30 pages, 5 figures
References in corpus (6)
- Spiral Defect Chaos in Large Aspect Ratio Rayleigh-Benard Convection
- Mean flow and spiral defect chaos in Rayleigh-Benard convection
- Direct Transition to Spatiotemporal Chaos in Low Prandtl Number Fluids
- Power-Law Behavior of Power Spectra in Low Prandtl Number Rayleigh-Benard Convection
- Rayleigh-Benard Convection with a Radial Ramp in Plate Separation
- Finite Size Scaling of Domain Chaos