Many metal manufacturing processes involve with near the melting processes such as solidification. Computational materials models such as phase-field models can be used to investigate these processes provided near the melting thermomechanical properties. While conducting experiments at those high temperatures is very challenging, molecular dynamic (MD) simulation offers a powerful tool to that end. Here, we present a systematic method to determine thermal linear expansion coefficient, heat capacity, and elastic constants of Fe,Cu,Ni, and Al during solidification using modified embedded-atom method (MEAM) MD simulations. MEAM potentials are developed by considering a wide range of low temperature properties and the melting point of the metal. Due to the huge fluctuations of atoms at high temperatures, the calculation of elastic constants is carried out by a universal, fast deformation–fluctuation hybrid approach which is significantly faster than direct method and results in much smaller standard deviations for the calculated data.