NOTICE: Due to a lapse in annual appropriations, most of this website is not being updated. Learn more.
Form submissions will still be accepted but will not receive responses at this time. Sections of this site for programs using non-appropriated funds (such as NVLAP) or those that are excepted from the shutdown (such as CHIPS and NVD) will continue to be updated.
An official website of the United States government
Here’s how you know
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
Secure .gov websites use HTTPS
A lock (
) or https:// means you’ve safely connected to the .gov website. Share sensitive information only on official, secure websites.
Boussinesq Algorithm for Enclosed Buoyant Convection in Two Dimensions.
Published
Author(s)
Ronald G. Rehm, Hai C. Tang, Howard R. Baum, James S. Sims, Daniel M. Corley
Abstract
Approximate equations for a Boussinesq model with viscous dissipation and thermal conduction describing buoyant convection driven by a heat source in rectangular enclosure are derived. The finite-difference algorithm for computing transient solutions in two dimensions to these equations is presented. The algorithm allows the enclosure fluid to be stratified in a direction parallel to the enclosure walls initially, or for gravity to have an arbitrary direction relative to the enclosure (but with no initial stratification). Computational results of transient, two-dimensional buoyant convection for very high resolution are presented. The hydrodynamics is directly based on the time-dependent Navier-Stokes equations; the model is valid in the Boussinesq approximation. No turbulence model or other empirical parameters are introduced. There is no inflow or outflow at boundaries; this assumption, although rather restrictive, allows the mathematical problem to be properly formulated so that no other empiricism is introduced by specification of the algorithmic boundary conditions. A finite-difference scheme second-order in space and first-order in time is used to integrate the evolution equations, and an elliptic solver is used to solve the pressure equation. The algorithms have been verified by comparisons with exact solutions to the equations in simple, special cases, and predictions of the overall model when the viscosity and thermal conductivity are zero have been compared with experimental results. The use of Lagrangian particle tracking allows one to visualize the flow patterns.
Rehm, R.
, Tang, H.
, Baum, H.
, Sims, J.
and Corley, D.
(1991),
Boussinesq Algorithm for Enclosed Buoyant Convection in Two Dimensions., NIST Interagency/Internal Report (NISTIR), National Institute of Standards and Technology, Gaithersburg, MD, [online], https://doi.org/10.6028/NIST.IR.4540
(Accessed October 14, 2025)