This paper presents a numerical method to address function estimation problems in inverse heat transfer problems using parameter estimation approach without prior information on the functional form of the variable to be estimated. Using an inverse analysis, the functional form of a time-dependent heat transfer coefficient is estimated efficiently and accurately. The functional form of the heat transfer coefficient is assumed unknown and the inverse heat transfer problem should be treated using a function estimation approach by solving sensitivity and adjoint problems during the minimization process. Based on proposing a new sensitivity matrix, however, the functional form can be estimated in an accurate and very efficient manner using a parameter estimation approach without the need for solving the sensitivity and adjoint problems and imposing extra computational cost, mathematical complexity, and implementation efforts. In the proposed sensitivity analysis scheme, all sensitivity coefficients can be computed in only one direct problem solution at each iteration. In this inverse heat transfer problem, the body shape is irregular and meshed using a body-fitted grid generation method. The direct heat conduction problem is solved using the finite-difference method. The steepest-descent method is used as a minimization algorithm to minimize the defined objective function and the termination of the minimization process is carried out based on the discrepancy principle. A test case with three different functional forms and two different measurement errors is considered to show the accuracy and efficiency of the used inverse analysis.