Dear Lammps Users:
I’m facing up with a problem. For my system, it’s
Orthoganal Cubic Box, with x,y periodic and Z shrink-wrapped
I constructed a rough wall(particles) my self in xy plane(normal to z direction) both at top and bottom.
The box is fulfilled with fluid particle.
For the initial condition I have set, It will lead to expand the box (only in z direction because x,y are periodic)
Now I need to constraint the wall movement, because the fluid inside is pushing the wall hardly that the wall just blow up.
of cause, I can apply external force onto the wall to avoid blow up. But what I want is to reach a equilibrium of the whole system, in other words, I want reach a final state that average force on wall is fluctuating about 0.
I used to do it with my in-house code, which is a bit like “Fix Spring/self” command. but I need to update the initial position of the wall particles. As I go through the manual, I find that for the command “Fix Spring/self” the initial position of each particle is set as the input value and never updated. thus cannot satisfied my need.
Can anyone give possible solution for the problem
Thanks A Lot