How to set rotation angle α in SRD (Stochastic rotation dynamics) of LAMMPS?

Using SRD to simulate hydrodynamic interactions requires setting the rotation angle α.

But I couldn’t find the relevant option for rotation angle in the fix srd command in LAMMPS manual. I would like to inquire about how to change the rotation angle setting. Thanks

Hello Dynamics_rookie,

as you will read in the articles cited by the fix srd documentation, the SRD implementation of LAMMPS uses an angle of 90° around one of the “six” Cartesian axes. Changing that requires modifying the source code. See also https://journals.aps.org/pre/abstract/10.1103/PhysRevE.68.036701 which does it similarly (“model B”).

Best,
November