161k views
3 votes
I am preparing to conduct shock simulation for an energetic material in LAMMPS with the help of a classical forcefield. How can I apply shock absorbing boundary condition so that that I can gather information during shock properly?

1 Answer

3 votes

Final Answer:

To apply shock absorbing boundary conditions in LAMMPS, use the command `fix deform` with the "tri" style, specifying the desired strain rate components for the shock direction, such as `fix 1 all deform 1 xy erate 0 0 0.1 units box`.

Step-by-step explanation:

To simulate shock conditions in LAMMPS with shock absorbing boundary conditions, the "fix deform" command can be employed. This command deforms the simulation box in a controlled manner. For shock simulations, the "tri" style is often preferred.

The deformation gradient components control the strain rates along different axes. For instance, to simulate a shock along the x-axis, the command `fix 1 all deform 1 xy erate 0 0 0.1 units box` is used. Here, "xy" indicates deformation along the xy plane, and the erate values control the strain rates in the x, y, and z directions, respectively. The erate values can be adjusted to achieve the desired strain rate for the shock.

It's crucial to select an appropriate erate value based on the specific shock conditions and the material properties. This value determines how rapidly the simulation box deforms, influencing the shockwave's intensity. Careful consideration of the material's response to shock and the desired outcome of the simulation is essential in determining suitable erate values. Adjustments may be needed through iterative simulations to achieve accurate and meaningful results for the energetic material under study.

User Vinodh Kumar
by
9.2k points