It is often useful to determine temperature and heat flux in multidimensional solid domains of arbitrary shape with inaccessible boundaries. In this study, an effective algorithm for solving boundary inverse heat conduction problems (IHCPs) is implemented: transient temperatures on inaccessible boundaries are estimated from redundant simulated measurements on accessible boundaries. A nonlinear heat equation is considered, where some of the material properties are dependent on temperature. The IHCP is reformulated as an optimization problem. The resulting functional is iteratively minimized using a conjugate gradient method together with an adjoint (dual) problem approach. The associated partial differential equations are solved using the finite-element package FEniCS. Tikhonov regularization is introduced to mitigate the ill-posedness of the IHCP. The accuracy of the implemented algorithm is assessed by comparing the solutions to the IHCP with the correct temperature values, on the inaccessible boundaries. The robustness of our method is tested by adding Gaussian noise to the initial conditions and redundant boundary data in the inverse problem formulation. A mesh independence study is performed. (C) 2018 Elsevier Ltd. All rights reserved.