fix addforce and variable atom

Submitted by thwagner on Sun, 04/03/2011 - 15:56

Dear Liggghts users,

I have difficulties in defining an additional force due to a couette flow. Therefor I want to add a force depending on the atom position. I define a per atom variable "b" and I want to refer to it as "v_b" in the fix addforce.

I guess I didn't got the idea how to use per atom variables within fixes like addforce.
I attached a simple input file in which I summarized, what I have already tried.

Maybe someone can show me a working example of an addforce/setforce fix with per atom force values.

With kind regards,
Thomas Wagner.

AttachmentSize
Plain text icon in.addforce.txt1.14 KB
ckloss's picture

ckloss | Mon, 04/04/2011 - 09:05

The fix addforce in the current LIGGGHTS version does not support variables. Refer to the doc in the tarball...
What you can try (at your own risk) is to use fix addforce from a recent LAMMPS version

Christoph

nt86 | Tue, 06/28/2011 - 12:57

I also bumped into this issue when tried to use fix addforce with variables. Do you have any news, when it will be implemented in LIGGGHS?

ckloss's picture

ckloss | Tue, 06/28/2011 - 13:45

at some point over the next few (1-3m months) a new major version will be released including all LAMMPS patches

Christoph