(541e) Mesoscale Modeling of Thermo-Mechanical Response: Constant-Energy Dissipative Particle Dynamics Simulations

Authors: 
Brennan, J. K. - Presenter, U.S. Army Research Laboratory
Lisal, M. - Presenter, Academy of Sciences of the Czech Republic


Nanocomposites such as energetic materials (EMs) are inherently heterogeneous, containing several species types of disparate shapes and sizes. EMs are often comprised of micro-sized crystallites (or grains) held together by a small amount of polymer binder, and in some instances contain metal nanoparticles. Accordingly, mechanical stimulation (e.g., shock or shear) of these materials can incite responses over a wide range of spatial and temporal scales. For example, localized regions of elevated thermal energy, or hot spots, often occur near microscale defects in the material such as voids and grain boundaries; however the energy transfer processes are atomistically governed. While modeling and simulation of such phenomena can offer fundamental insight, it is computationally challenging. Computing power and resources continue to grow, but there still exists much behavior that cannot be studied using atomistic simulation techniques. Moreover, field-based mesoscale modeling techniques have fidelity limitations with respect to gaining a fundamental understanding.

To address these issues, we recently embarked on a study of the thermo-mechanical response of crystalline nanocomposites using a particle-based mesoscale modeling technique, namely, the constant-energy Dissipative Particle Dynamics method (DPDe) [1, 2]. A constant energy approach is required since materials subjected to mechanical stimuli will inevitably respond via energy exchange and transfer processes. In contrast to the standard Dissipative Particle Dynamics method [3, 4] that conserves only momentum, the DPDe method conserves both momentum and energy. The DPDe method assigns an internal energy to each particle, allowing particles to exchange both momentum and thermal energy. This particle internal energy is included as a separate equation of motion along with the equations of motion for the particle's position and momentum. The atomic degrees-of-freedom are explicitly resolved or coarse-grained into the particle internal energy, and included in the simulation in an averaged way.

A number of modeling challenges exist for nanocomposites, including the development of an accurate mesoscale model that can capture the following known thermo-mechanical responses: (i) phase transitions; (ii) structural rearrangements (plastic deformation, shear banding); and (iii) chemical decomposition reactions. The model parameters should be physics-based so that direct links to real material properties or higher resolution modeling are possible. The standard DPD model [5] likely cannot capture such features, fortunately though, promising alternatives exist [6-8]. Technical challenges also exist for such studies. Since while numerical integrator algorithms work quite satisfactorily for DPD [5] and DPDe [8] simulations under normal conditions, at extreme conditions such as high densities, much shorter integration timesteps are required to maintain numerical stability. In order to maintain a reasonable timestep that allows for the simulation of mesoscale events, higher quality integrators will be explored.

In this work, we will present the status of this ongoing project, including progress on nanocomposite model development, mesoscale simulation of phase transitions, and improved integrator algorithms. To date, progress has been encouraging, where we find the DPDe method to be a viable tool for simulating the thermo-mechanical response of crystalline nanocomposites to mechanical stimuli.

References

[1] Bonet Avalos, J., and A.D. Mackie, Europhysics Letters, 40, 141 (1997).

[2] Español, P., Europhysics Letters, 40, 631 (1997).

[3] Hoogerbrugge, P.J., and J.M.V.A. Koelman, Europhys. Lett., 19, 155 (1992).

[4] Koelman, J.M.V.A., and P.J. Hoogerbrugge, Europhys. Lett., 21, 363 (1993).

[5] Groot, R.D., and P.B. Warren, J. Chem. Phys., 107, 4423 (1997).

[6] Pagonabarraga, I., and D. Frenkel, Molec. Simul., 25, 167 (2000).

[7] Pagonabarraga, I., and D. Frenkel, J. Chem. Phys., 115, 5015 (2001).

[8] Groot, R.D., and S.D. Stoyanov, Phys. Rev. E, 78, 051403(1-9) (2008).

[9] Mackie, A.D., J. Bonet Avalos and V. Navas, Phys. Chem. Chem. Phys., 1, 2039 (1999).