My problem is related with rotate option of region command. I would appreciate if you state where I am doing something wrong.
I am trying to create 6 cylindrical regions which are connected to each other forming a polygonal ring. Orientations of the cylindrical regions are different so their longitudinal axes are not on the same plane.
For the generation of each rotated region, first I assume a region along the x axis starting at the same point and having the same length. Let’s say this vector along the x axis is A, and the vector that I want to generate is B. Then, I calculate the rotation axis as the cross product of A and B, i.e. cross(A,B), and normalize it to get a unit vector. Angle between two vectors is also calculated conventionally. So, I want to rotate the vector A to B. But, it results in cylindrical regions that are disconnected. I am using the windows version (27 Nov 2013-ICMS) of Lammps.
Attached you can see my input script for Lammps.
den.in (2.12 KB)