Dear Lammps users,
I am new to lammps. I am trying to model the classic example of poiseuille flow in 3D. So I tired to modify the example problem given in the Lammps package mainly for a real wall model ( velocities are not set to zero for wall particles and a berendsen thermostat is used on wall and fluid particles separately).
In current simulation Argon gas is confined in two parallel Platinum plates. Now the problem I am facing while running the code is the very high temperature in the wall atoms with a time step less than 0.0001 fs or atoms are lost with a time step greater than 0.0001 fs.
A part of the code related to the potentials are given below:
pair_style lj/cut 10.215
pair_coeff 1 1 3.45 0.238 10.215
pair_coeff 2 2 2.47 7.495 6
pair_coeff 1 2 2.94 0.157 9
[ 1 for argon particle, 2 for platinum particle]
How can we reduce the high velocity/energy on wall atoms due to the potential?
Thanking you for your time.
Regards
Abdul Gafoor