let me refer to an article titled "Static and dynamic behavior of water droplet on solid surfaces with pillar-type nanostructures from molecular dynamics simulation"
this is similar to my case except that I want to add heat transfer to the simulation.
I was thinking that if it was possible to add some amount to the surface atoms' velocity as oscillating velocity while it is not rescaled by fix nvt , it would be easy to implement. but the only problem is that the "velocity set v_vx v_vy v_vz sum yes" command does not take a time-dependent value.
I apologize if I could not explain my meaning clearly and strongly appreciate your help.