Hi All,
In a thread from May 2018 a similar question was discussed under the thread “Helium flow” but I did not get a clear enough answer from reading that thread.
I am new to LAMMPS but have done extensive MD modeling and have even written my own QMD code in the past.
I want to simulate Lennard-Jones or similar gaseous flow through various shaped channels by introducing new atoms, with a specified average net velocity, in the upstream side of a simulation cell. On the downstream side, I want the old atoms (that have passed through already) to be either eliminated to be replaced by “new” upstream atoms. I worry that atoms that have passed through the channel will have a hysteresis (turbulence?) that will affect their subsequent passages through the channel if I were to use periodic boundary conditions (PBC). Equivalently, if PBC are used, the velocities would have to be reset to avoid hysteresis.
Is there a mechanism in the current LAMMPS to allow this kind of control?
If not, could I create an external module “add-on” to do this?
If not either, is there a NEMD code that could do this?
Any advice would be very helpful.
Derrick E. Boucher, Ph.D.
Associate Professor of Physics and Physics Coordinator
Florida Gulf Coast University
10501 FGCU Blvd. South
Fort Myers, FL 33965
U.S.A.
Phone: (239) 590-7170 Office: Whitaker Hall 255
----------------- Required legal notice for Florida State employees----------------------------
Florida has a very broad public records law. As a result, any written communication created or received by Florida Gulf Coast University employees is subject to disclosure to the public and the media, upon request, unless otherwise exempt. Under Florida law, e-mail addresses are public records. If you do not want your email address released in response to a public records request, do not send electronic mail to this entity. Instead, contact this office by phone or in writing.