Natural Convection CFD Analysis
Natural Convection CFD: Theoretical Foundations
Fundamentals of Natural Convection
Professor, natural convection is when flow occurs without an external power source, right? What is the driving force?
It's buoyancy. The density of the fluid changes with temperature, and in a gravitational field, this density difference creates a pressure imbalance, which generates flow. Fluid in the high-temperature region becomes less dense and rises, while fluid in the low-temperature region descends. This circulation is natural convection.
What is the indicator that represents the strength of natural convection?
The Rayleigh number $Ra$ is the fundamental indicator.
Here, $g$ is gravitational acceleration, $\beta$ is the volumetric expansion coefficient, $\Delta T$ is the temperature difference, $L$ is the characteristic length, $\nu$ is the kinematic viscosity, $\alpha$ is the thermal diffusivity. It is the product of the Grashof number $Gr$ and the Prandtl number $Pr$.
The flow regime is divided by $Ra$. For a vertical flat plate, $Ra < 10^9$ is laminar, $Ra > 10^9$ is turbulent. In an enclosed cavity, transition to unsteady flow occurs around $Ra > 10^8$.
Boussinesq Approximation
What exactly is the assumption of the Boussinesq approximation?
It's an approximation applicable when density fluctuations are sufficiently small compared to a reference density ($\beta \Delta T \ll 1$). It considers density variation only in the buoyancy term of the momentum equation, while treating density as constant in the continuity equation.
For air, a safe application range is typically $\Delta T < 30$°C. For water, since the volumetric expansion coefficient has temperature dependence, a polynomial density model should be used for wide temperature ranges.
The World of "Turbulent Natural Convection" Beyond Rayleigh Number 10⁹
The flow regime in natural convection changes dramatically with the Rayleigh number (Ra = Gr × Pr). For Ra < 10⁸, stable laminar convection cells form, and the correlation for the Nusselt number is relatively simple. However, when Ra > 10⁹, turbulence begins, and the prediction accuracy of the Nu number plummets. This scale corresponds to "building exterior walls heated by sunlight" or "oil cooling for large transformers"—conditions commonly encountered in practice. In turbulent natural convection, plume-like structures near walls randomly generate and vanish, often causing steady-state analyses to fail to converge. In such cases, switching to a strategy of taking time averages from unsteady calculations is advisable.
Computational Methods for Natural Convection CFD
Turbulence Modeling for Natural Convection
I've heard that selecting a turbulence model for natural convection is difficult.
Compared to forced convection, natural convection has lower turbulence intensity and a wider transition region. The standard k-ε model tends to produce excessive turbulent diffusion and overestimate the Nu number. The recommended order is as follows.
| Rayleigh Number Range | Recommended Model | Notes |
|---|---|---|
| $Ra < 10^9$ | Laminar (Laminar Model) | Turbulence model not required |
| $10^9 < Ra < 10^{12}$ | SST k-ω + Low-Re Wall Treatment | Wall first layer $y^+ < 1$ mandatory |
| $Ra > 10^{12}$ | SST k-ω or Realizable k-ε | Enhanced Wall Treatment |
Can't wall functions be used for natural convection?
They should generally be avoided. The velocity and temperature profiles near walls in natural convection differ from the logarithmic law of forced convection, making standard wall functions less applicable. The Low-Re approach with $y^+ \approx 1$ is recommended. Fluent's Enhanced Wall Treatment automatically switches based on $y^+$, but for natural convection, ensuring $y^+ < 1$ is the best practice.
Mesh Design
What should I be careful about when designing a mesh for natural convection?
Estimate the thickness of the thermal and velocity boundary layers and place a sufficient number of cells within each. The boundary layer thickness for natural convection can be approximated by
For example, for a vertical flat plate with $L = 0.1$m and $Ra = 10^9$, $\delta_T \approx 0.6$mm. At least 10–15 cell layers need to be placed within this thin boundary layer.
Is it okay to make the region outside the boundary layer coarse?
The core region can be relatively coarse. However, abrupt cell size changes should be avoided; a guideline is to keep the volume ratio of adjacent cells below 1.2. STAR-CCM+'s Prism Layer Mesher or Fluent's Inflation Layer can automatically refine the region near walls.
Boussinesq Approximation in Natural Convection CFD—When It Works and When It Breaks Down
The standard "Boussinesq approximation" in natural convection analysis is a technique that linearly approximates density variation only in the buoyancy term (ρ ≈ ρ₀(1-βΔT)) and treats density as constant in other terms. It stabilizes calculations and makes convergence easier, but errors become non-negligible when the temperature difference ΔT is large. As a rule of thumb, the application limit is "βΔT < 0.1 (approximately up to 10–20°C difference)." For systems with temperature differences of several hundred degrees, such as furnace combustion or solar thermal collectors, a "non-Boussinesq (variable density)" model that treats density as a full function of temperature is essential. This requires either considering low-Mach-number compressibility in a pressure-based Navier-Stokes solver or directly incorporating a density-dependent equation of state.
Natural Convection CFD in Practice
Enclosed Cavity Benchmark
Are there any benchmarks available for validating natural convection CFD?
The most famous is the differentially heated rectangular cavity problem by de Vahl Davis (1983). The left wall is hot, the right wall is cold, the top and bottom walls are adiabatic, and reference solutions for Nu number and flow field are provided for $Ra = 10^3$ to $10^6$. For CFD code validation, it is standard to confirm that the wall-averaged Nu number matches de Vahl Davis's values within 1%.
What are the specific values?
| $Ra$ | Average $Nu$ (de Vahl Davis) |
|---|---|
| $10^3$ | 1.118 |
| $10^4$ | 2.243 |
| $10^5$ | 4.519 |
| $10^6$ | 8.800 |
If CFD results deviate from these by more than 2%, there is likely a problem with the settings or mesh.
Natural Air Cooling Design for Electronic Devices
What kind of CFD design is done for fanless natural air cooling?
For natural convection CFD inside an enclosed enclosure, the procedure is: (1) Set the heat generation of each component as a volumetric source, (2) Set natural convection + radiation boundary conditions on the outer surface of the enclosure walls, (3) Solve for the air inside the enclosure as the fluid. In many cases, radiation contributes 30–50% of the total heat dissipation, so combining a Surface-to-Surface radiation model (Fluent's S2S, STAR-CCM+'s Surface-to-Surface Radiation) is essential.
Modeling the heat conduction of a circuit board is tough, right?
A PCB (printed circuit board) has a laminated structure of copper layers and glass epoxy layers, with thermal conductivity differing by more than 10 times between the in-plane and thickness directions. Tools specialized for electronic thermal design like Ansys Icepak have automatic anisotropic modeling functions for PCBs. With general-purpose CFD solvers, you need to manually set orthotropic thermal conductivity.
What about enclosures with ventilation holes?
When the enclosure has intake and exhaust openings, the flow is no longer purely natural convection inside—it transitions toward "mixed convection" driven by both density differences and pressure gradients from outside air. In such cases, the computational domain should include both the interior and a sufficiently large exterior domain. Typical practice: define the exterior domain as at least 2–5 times the enclosure size, with appropriate boundary conditions (atmospheric temperature and pressure).