In this paper, we present axisymmetric numerical simulations of shock propagation in nitromethane over an aluminum particle for post-shock pressures up to 10 GPa. We use the Mie-Gruneisen equation of state to describe both the medium and the particle. The numerical method is a finite-volume based solver on a Cartesian grid, that allows for multi-material interfaces and shocks, and uses a novel constrained reinitialization scheme to precisely preserve particle mass and volume. We compute the unsteady inviscid drag coefficient as a function of time, and show that when normalized by post-shock conditions, the maximum drag coefficient decreases with increasing post-shock pressure. We also compute the mass-averaged particle pressure and show that the observed oscillations inside the particle are on the particle-acoustic time scale. Finally, we present simplified point-particle models that can be used for macroscale simulations. In the Appendix, we extend the isothermal or isentropic assumption concerning the point-force models to non-ideal equations of state, thus justifying their use for the current problem.
Sridharan, P., Jackson, T. L., Zhang, J., Balachandar, S., & Thakur, S. (2016). Shock interaction with deformable particles using a constrained interface reinitialization scheme. Journal of Applied Physics, 119(6). https://doi.org/10.1063/1.4941687
Sridharan, P., Jackson, T. L., Zhang, J., et al., "Shock interaction with deformable particles using a constrained interface reinitialization scheme," Journal of Applied Physics 119, no. 6 (2016), https://doi.org/10.1063/1.4941687
@article{osti_1421155,
author = {Sridharan, P. and Jackson, T. L. and Zhang, J. and Balachandar, S. and Thakur, S.},
title = {Shock interaction with deformable particles using a constrained interface reinitialization scheme},
annote = {In this paper, we present axisymmetric numerical simulations of shock propagation in nitromethane over an aluminum particle for post-shock pressures up to 10 GPa. We use the Mie-Gruneisen equation of state to describe both the medium and the particle. The numerical method is a finite-volume based solver on a Cartesian grid, that allows for multi-material interfaces and shocks, and uses a novel constrained reinitialization scheme to precisely preserve particle mass and volume. We compute the unsteady inviscid drag coefficient as a function of time, and show that when normalized by post-shock conditions, the maximum drag coefficient decreases with increasing post-shock pressure. We also compute the mass-averaged particle pressure and show that the observed oscillations inside the particle are on the particle-acoustic time scale. Finally, we present simplified point-particle models that can be used for macroscale simulations. In the Appendix, we extend the isothermal or isentropic assumption concerning the point-force models to non-ideal equations of state, thus justifying their use for the current problem.},
doi = {10.1063/1.4941687},
url = {https://www.osti.gov/biblio/1421155},
journal = {Journal of Applied Physics},
issn = {ISSN 0021-8979},
number = {6},
volume = {119},
place = {United States},
publisher = {American Institute of Physics},
year = {2016},
month = {02}}
Proceedings of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences, Vol. 459, Issue 2031https://doi.org/10.1098/rspa.2002.1045