The main work of this paper focuses on identifying the heat flux in inverse problem of two-dimensional nonhomogeneous parabolic equation, which has wide application in the industrial field such as steel-making and continuous casting. Firstly, the existence of the weak solution of the inverse problem is discussed. With the help of forward solution and dual equation, this paper proves the Lipchitz continuity of the cost function and derives the Lipchitz constant. Furthermore, in order to accelerate the convergence rate and reduce the running time, this paper presents a sufficient descent Levenberg–Marquard algorithm with adaptive parameter (SD-LMAP) to solve this inverse problem. At last, compared with other methods, the results of simulation experiment show that this algorithm can obviously reduce the running time and iterative number.