An efficient new method has been developed to incorporate the effects of heat transfer in a liquid pool into models of heat conduction with solidification. The procedure has been added into the commercial package Abaqus  as a user-defined subroutine (UMATHT). Computational results of fluid flow and heat transfer in a liquid domain can be characterized by the heat flux crossing the boundary representing the solidification front, or liquidus temperature. This superheat flux can be incorporated into an uncoupled transient simulation of heat transfer phenomena in the mushy and solid regions by enhancing latent heat. The new method has been validated and compared to semianalytical solutions and two other numerical methods on simple test problems: two-dimensional, steady-state ledge formation in cryolite in aluminum extraction cells, and shell thinning in continuous casting of steel. Its real power, however, is for multiphysics simulations involving complex phenomena, such as solidification stress analysis with nonlinear constitutive equations. Including the superheat flux from a thermal-fluid flow simulation of the liquid pool into the latent heat provides a very efficient and robust method for incorporating the effects of fluid flow in the liquid pool into thermal-stress problems, especially for transient problems.