chtMultiRegionSimpleFoam - OpenFOAM Solver

Solver: chtMultiRegionSimpleFoam   Description

chtMultiRegionSimpleFoam is a solver designed for steady-state simulations involving buoyant, turbulent fluid flow and solid heat conduction. It is dedicated to handling conjugate heat transfer (CHT) between solid and fluid regions, making it suitable for a broad range of multiphysics problems where interaction between fluid and solid phases is a key aspect.

The solver uses the PIMPLE (merged PISO-SIMPLE) algorithm for pressure-momentum coupling. This algorithm leverages the strengths of both PISO and SIMPLE methods for pressure-velocity coupling, ensuring robustness in handling transient flows with large time steps. This approach is supplemented by under-relaxation techniques to secure convergence stability. It supports Multiple Reference Frames (MRF), porosity modeling and allows easy integration of passive scalar transport equations and source terms.

Applications of the solver span across industries where thermal management and heat transfer are crucial. This includes HVAC systems, where optimizing airflow and temperature distribution can significantly impact energy efficiency and comfort. In the realm of electronics, the solver aids in the design of cooling systems, ensuring devices operate within safe temperature limits. The energy sector also benefits from this solver, especially in analyzing heat exchanger performance.

Solver: chtMultiRegionSimpleFoam   Features

  • Steady-State
  • Compressible
  • Multi-Region Mesh (Solid/Fluid)
  • Arbitrary Fluid and Solid Regions
  • Conjugate Heat Transfer (CHT)
  • Heat Source
  • Radiation
  • Laminar and Turbulent (RANS, LES, DES)
  • Equation of State Models
  • Pressure-Based Solver
  • Passive Scalar
  • Porosity Modeling
  • Buoyancy
  • Source Term (explicit/implicit)
  • SIMPLE Algorithm
  • Solution Limiters:
    • Velocity Damping
    • Pressure Limit
    • Temperature Limit

Solver: chtMultiRegionSimpleFoam   Application


  • CPU and/or GPU Cooling (e.g. Electronics cooling)
  • Mother Board Heating Component Placement
  • Heat Sink Design
  • Fan Sizing and Placement

Lightning design

  • LED Fixtures Cooling
  • LED Spacing, PCB Choice
  • Active vs. Passive Cooling



  • Heating and Cooling Systems
  • Heat Exchangers
  • Solar Load


  • Engine Cooling
  • Windshield Condensation
  • Windshield Defrost

Solver: chtMultiRegionSimpleFoam   Heat Transfer Solvers

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: chtMultiRegionSimpleFoam   Tutorial

  • Simulation of the heat transfer between a hot and cold region in a heat exchanger.
  • Electronic cooling simulation, with a Joule heat source in a CPU and a cooling domain with flowing air. The steady-state simulation calculates conjugate heat transfer (CHT) on the CPU-air interface, considering the radiation.
  • Conjugate heat transfer (CHT) simulation of a heat transfer from a heated plate on the top surface and cooled by flowing water through a channel.

Solver: chtMultiRegionSimpleFoam   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\)]


\(p\) [\(Pa\)]

Hydrostatic Perturbation Pressure

\(p - \rho gh\) [\(Pa\)]

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

Derivative Results


\(\rho\) [\(\frac{kg}{m^{3}}\)]


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

Mach Number

\(Ma\) [\(-\)]

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\)]

Wall Heat Flux

\(\phi_q\) [\(W/m^2\)]

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\)