In this study, we demonstrated that by generating an arbitrary temperature field using the FEM and incorporating this temperature field into BDSs to account for thermophoresis, we were able to obtain computational results that exhibit satisfactory agreement with experimental validation. Our computational approach is expected to serve as a flexible and effective tool for quantitatively evaluating the thermophoretic behavior of PSNPs in liquid phases. In this validation, a temperature field was formed by laser-induced single-spot heating, making the thermophoretic behavior relatively straightforward to visualize. However, given that the depletion range of the PSNPs closely matched the experimental results, we believe that the same approach could be effectively applied to more complex temperature fields, such as those involving multiple plasmonic nanoparticles, by using FEM to calculate arbitrary temperature fields and simulating the thermophoretic behavior of nanoparticles.