buoyantBoussinesqSimpleFoam - OpenFOAM Solver

Solver: buoyantBoussinesqSimpleFoam   Description

buoyantBoussinesqSimpleFoam is a pressure-based solver designed for steady-state simulations of incompressible flows. It handles laminar and turbulent, single-phase flows. The unique feature of the solver is its use of the Boussinesq approximation, which simplifies the computations related to buoyancy by linearly relating the density changes in the fluid to temperature changes, based on a reference temperature. The solver is particularly accurate and efficient when the changes in density compared to a reference density are small.

The solver uses the SIMPLE (Semi-Implicit Method for Pressure-Linked Equations) algorithm for pressure-momentum coupling, augmented by under-relaxation techniques to enhance convergence. It supports Multiple Reference Frame (MRF) and porosity modeling and allows easy integration of passive scalar transport equations and source terms.

A common application of the solver is limited to incompressible flows with relatively low-temperature variation. For example, for the air, the variation should be less than 15 Celsius degrees and less than 2 Celsius degrees for water.

For HVAC (Heating, Ventilation, and Air Conditioning) applications, the solver can be used in room ventilation simulations, and chimney effect in buildings. In electronics applications, the solver can be applied to study passive cooling systems.

Solver: buoyantBoussinesqSimpleFoam   Features

  • Steady-State
  • Incompressible
  • Single-Phase
  • Buoyancy using Boussinesq Approximation
  • Heat Transfer
  • Heat Source (temperature source therm)
  • Radiation
  • Laminar and Turbulent (RANS, LES, DES)
  • Newtonian and Non-Newtonian Fluid
  • Pressure-Based Solver
  • Rotating Objects:
    • Multiple Reference Frames (MRF)
  • Passive Scalar
  • Porosity Modeling
  • Buoyancy
  • Source Term (explicit/implicit)
  • SIMPLE Algorithm
  • Solution Limiters:
    • Velocity Damping

Solver: buoyantBoussinesqSimpleFoam   Application


  • Ventilation
  • Water Heating

Solver: buoyantBoussinesqSimpleFoam   Heat Transfer Solvers Comparison

Heat Transfer Solvers In this group, we have included solvers that are designed to model: Heat Transfer, Radiation, Natural and Forced Convection, Conjugate Heat Transfer (CHT).

Heat Transfer, Single Fluid

Heat Transfer, Single Fluid - Boussinesq

Heat Transfer, Single Solid

CHT, Multiple Fluids / Solids

  • CHT - Conjugate Heat Transfer
  • MRF - Multiple Reference Frame
  • Overset - also known as Chimera Grid (Method)

Solver: buoyantBoussinesqSimpleFoam   Alternative Solvers

In this section, we propose alternative solvers from different categories, distinct from the current solver. While they may fulfill similar purposes, they diverge significantly in approach and certain features.

  • simpleFoam base version of buoyantBoussinesqSimpleFoam without buoyancy and heat transfer

Solver: buoyantBoussinesqSimpleFoam   Results Fields

This solver provides the following results fields:

  • Primary Results Fields - quantities produced by the solver as default outputs
  • Derivative Results - quantities that can be computed based on primary results and supplementary models. They are not initially produced by the solver as default outputs.

Primary Results Fields


\(U\) [\(\frac{m}{s}\)]


\(T\) [\(K\)]

Kinematic Pressure \(p/\rho\)

\(p\) [\(\frac{m^{2}}{s^{2}}\)]

Kinematic Hydrostatic Perturbation Pressure

\(p - \rho gh\) [\(\frac{m^{2}}{s^{2}}\)]

Hydrostatic Perturbation Pressure This value represents the pressure without the hydrostatic component (minus gravitational potential). Read More: Hydrostatic Pressure Effects

Kinematic Pressure It is a pressure normalized by density. To obtain pressure in Pascals [Pa], multiply kinematic pressure by the fluid’s reference density. Read More: Kinematic Fluid Properties

Derivative Results


\(P\) [\(Pa\)]


\(\omega\) [\(\frac{1}{s}\)]

Peclet Number

\(Pe\) [\(-\)]

Stream Function

\(\psi\) [\(\frac{m^2}{s}\)]

Q Criterion

\(Q\) [\(-\)]

Wall Functions (for RANS/LES turbulence)

\(y^+\) [\(-\)]

Wall Shear Stress

\(WSS\) [\(Pa\)]

Turbulent Fields (for RANS/LES turbulence)

\(k\) \(\epsilon\) \(\omega\) \(R\) \(L\) \(I\) \(\nu_t\) \(\alpha_t\)

Volumetric Stream

\(\phi\) [\(\frac{m^{3}}{s}\)]

Passive Scalar

\(scalar_i\) [\(-\)]

Forces and Torque acting on the Boundary

\(F\) [\(N\)] \(M\) [\(-\)]

Force Coefficients

\(C_l\) [\(-\)] \(C_d\) [\(-\)] \(C_m\) [\(-\)]

Average, Minimum or Maximum in Volume from any Result Field

\(Avg\) \(Min\) \(Max\)