Hi Lammps users and developers,
I want to conduct a "Poiseuille Flow Between Infinite Parallel Plates"
simulation and I notice that there is an example for this. However, the
example uses fix enforce2D and is only 2 dimensional. I want to perform
this in three dimensions and any suggestions into how to go about this is
appreciated.
How do you impose a Pressure gradient? I see in the example that it uses
an addforce to the group of atoms that flow.However the top boundary has a
fix addforce also orthogonal to the flow direction. Why is this?
I plan to you use fix nvt hence I need only to thermostat the 2 axes
orthogonal to the flow so I will use compute temp/partial for this.Is this
the only thing I need to consider?
Jan-Michael Carrillo