-
Notifications
You must be signed in to change notification settings - Fork 4
2C. Particle type modifiers
--(***) All the commands listed here are not considered if system exportation from a previous simulation is done
Once we have our system defined with the desired shape and topography, we can change locally the constituent atoms or particles. Very similar commands to the shown ones in the previous section do the job. A region of the system sees its type changed into a final atom type:
CHANGE_AB_PLANE_TYPE (int)pos_cell_a (int)pos_cell_b (int)pos_cell_c (int)length_a (int)length_b (text)final_atom_type
CHANGE_AC_PLANE_TYPE (int)pos_cell_a (int)pos_cell_b (int)pos_cell_c (int)length_a (int)length_c (text)final_atom_type
CHANGE_BC_PLANE_TYPE (int)pos_cell_a (int)pos_cell_b (int)pos_cell_c (int)length_b (int)length_c (text)final_atom_type
CHANGE_CUBE_TYPE (double)x (double)y (double)z (double)side (text)final_atom_type
CHANGE_XY_DISLOCATION_TYPE (double)x (double)y (double)radius (text)final_atom_type
(optional line) FROM_Z_TO_Z (double)bot_z (double)top_z
(optional line) ANGLE_XZ_ANGLE_YZ (double)angle_xz (double)angle_yz
CHANGE_XZ_DISLOCATION_TYPE (double)x (double)z (double)radius (text)final_atom_type
(optional line) FROM_Y_TO_Y (double)bot_y (double)top_y
(optional line) ANGLE_XY_ANGLE_ZY (double)angle_xz (double)angle_yz
CHANGE_YZ_DISLOCATION_TYPE (double)y (double)z (double)radius (text)final_atom_type
(optional line) FROM_X_TO_X (double)bot_x (double)top_x
(optional line) ANGLE_YX_ANGLE_ZX (double)angle_yx (double)angle_zx
CHANGE_PLANE_TYPE (double)A (double)B (double)C (double)D (double)distance (text)final_atom_type
CHANGE_SPHERE_TYPE (double)x (double)y (double)z (double)radius (text)final_atom_type
CHANGE_ELLIPSOID_TYPE (double)x (double)y (double)z (double)radiusx (double)radiusy (double)radiusz (text)final_atom_type
CHANGE_GENERAL_ELLIPSOID_TYPE (double)A (double)B (double)C (double)D (double)E (double)F (double)G (double)H (double)J (double)K (text)final_atom_type
-
The parameters of this last command
CHANGE_GENERAL_ELLIPSOID_TYPEindicates the coeficients of the ellipsoid general equation Ax^2+By^2+Cz^2+Dxy+Exz+Fyz+Gx+Hy+Jz+K<1. -
The parameters of the
CHANGE_PLANE_TYPEcommand indicates the coeficients of the plane general equation Ax+By+Cz+D=0. -
the
pos_cellgoes from 0 toDIMENSION_A-1,DIMENSION_B-1 andDIMENSION_C-1 respectively.
All the previous commands can be used selectively targeting atoms or particles of certain type. They are identical to the previous ones, but they have the TO_TYPE usarname and an addition input parameter with the target type:
CHANGE_AB_PLANE_TYPE_TO_TYPE (text)atom_type (int)pos_cell_a (int)pos_cell_b (int)pos_cell_c (int)length_a (int)length_b (text)final_atom_type
CHANGE_AC_PLANE_TYPE_TO_TYPE (text)atom_type (int)pos_cell_a (int)pos_cell_b (int)pos_cell_c (int)length_b (int)length_c (text)final_atom_type
CHANGE_BC_PLANE_TYPE_TO_TYPE (text)atom_type (int)pos_cell_a (int)pos_cell_b (int)pos_cell_c (int)length_a (int)length_c (text)final_atom_type
CHANGE_CUBE_TYPE_TO_TYPE (text)atom_type (double)x (double)y (double)z (double)side (text)final_atom_type
CHANGE_XY_DISLOCATION_TYPE_TO_TYPE (text)atom_type (double)x (double)y (double)radius (text)final_atom_type
(optional line) FROM_Z_TO_Z (double)bot_z (double)top_z
(optional line) ANGLE_XZ_ANGLE_YZ (double)angle_xz (double)angle_yz
CHANGE_XZ_DISLOCATION_TYPE_TO_TYPE (text)atom_type (double)x (double)z (double)radius (text)final_atom_type
(optional line) FROM_Y_TO_Y (double)bot_y (double)top_y
(optional line) ANGLE_XY_ANGLE_ZY (double)angle_xz (double)angle_yz
CHANGE_YZ_DISLOCATION_TYPE_TO_TYPE (text)atom_type (double)y (double)z (double)radius (text)final_atom_type
(optional line) FROM_X_TO_X (double)bot_x (double)top_x
(optional line) ANGLE_YX_ANGLE_ZX (double)angle_yx (double)angle_zx
CHANGE_PLANE_TYPE_TO_TYPE (text)atom_type (double)A (double)B (double)C (double)D (double)distance (text)final_atom_type
CHANGE_SPHERE_TYPE_TO_TYPE (text)atom_type (double)x (double)y (double)z (double)radius (text)final_atom_type
CHANGE_ELLIPSOID_TYPE_TO_TYPE (text)atom_type (double)x (double)y (double)z (double)radiusx (double)radiusy (double)radiusz (text)final_atom_type
CHANGE_GENERAL_ELLIPSOID_TYPE_TO_TYPE (text)atom_type (double)A (double)B (double)C (double)D (double)E (double)F (double)G (double)H (double)J (double)K (text)final_atom_type
At this point we have our system perfectly defined. In the next step we are going to define the dissolution events that may happen in the system.
Whenever you use KIMERA, please cite "Martin, P., Gaitero, J. J., Dolado, J. S., & Manzano, H. (2020). KIMERA: A Kinetic Montecarlo Code for Mineral Dissolution. Minerals, 10(9), 825" https://www.mdpi.com/2075-163X/10/9/825/htm