One important issue in understanding and modeling of turbulent heat transfer is the behavior of fluctuating temperature close to the wall. Common engineering computational approach assumes constant heat flux boundary condition on heated walls. In the present paper constant heat flux boundary condition was assumed and effects of temperature fluctuations are investigated using large eddy simulations (LES) approach. A series of large eddy simulations for two geometries is performed: First, forced convection in channels and second, forced convection over a backward facing step. LES simulation data is statistically analyzed and compared with results of direct numerical simulations (DNS) from the literature which apply three cases of heat flux boundary conditions: 1. ideal heat flux boundary condition, 2. non-ideal heat flux boundary condition, 3. conjugate heat transfer boundary condition. For low Prandtl number flows LES results show that, despite very good agreement for velocities and mean temperature, predictions of temperature fluctuations may have strong deficiencies if simplified boundary conditions are applied.