A monolithic solver based on a diffuse-interface immersed-boundary (IB) approach for conjugate-heat-transfer (CHT) problems is presented. The IB strategy assumes that the solid which is "immersed" into the computational grid is occupied by a "virtual" fluid to facilitate construction of "unified" governing equations that are solved everywhere in the domain. A unified momentum equation is devised using the solid volume fraction that reduces to the Navier-Stokes equation outside of the solid and to the no-slip boundary condition inside of it. The "unified" energy equation is constructed in an analogous fashion reducing to a convective-diffusive equation in the fluid domain and a fully diffusive equation in the solid domain with different thermal conductivities (or diffusivities) for both domains. The resulting equations are solved in both domains simultaneously using a hybrid staggered and nonstaggered finite-volume (FV) framework for incompressible flows. The second-order accurate IB-FV solver is employed to carry out investigations for CHT problems in natural and forced convective regimes. Numerical studies for different fluid-to-solid conductivity ratios show that the monolithic IB-CHT solver is a fast, simple, and accurate framework for simulations of CHT problems for Boussinesq flows. © 2019 American Physical Society.