A hybrid finite element formulation is used to model temperature evolutions in concrete structures at early ages, taking into account the thermo-chemical process of cement hydration. The most relevant fields in the application problem, namely the temperature and the hydration degree fields, are selected for direct and independent approximation. The Galerkin formulation used qualifies as hybrid because the heat flux is also approximated on the boundary of the element. The temperature and hydration degree bases are used to enforce the thermal equilibrium condition and the hydration law, respectively, and the heat flux basis is used to enforce the thermal continuity conditions. The resulting nonlinear transient hybrid formulation leads to high...