| Autor: |
Lísal M; Department of Molecular and Mesoscopic Modelling, The Czech Academy of Sciences, Institute of Chemical Process Fundamentals, Prague 165 01, Czech Republic.; Department of Physics, Faculty of Science, Jan Evangelista Purkyně University in Ústí nad Labem, Ústí n. Lab. 400 96, Czech Republic., Larentzos JP; U.S. Army Combat Capabilities Development Command (DEVCOM) Army Research Laboratory, Aberdeen Proving Ground, Maryland 21005, United States., Avalos JB; Departament d'Enginyeria Química, ETSEQ, Universitat Rovira i Virgili, Tarragona 43007, Spain., Mackie AD; Departament d'Enginyeria Química, ETSEQ, Universitat Rovira i Virgili, Tarragona 43007, Spain., Brennan JK; U.S. Army Combat Capabilities Development Command (DEVCOM) Army Research Laboratory, Aberdeen Proving Ground, Maryland 21005, United States. |
| Abstrakt: |
We present an extension of the generalized energy-conserving dissipative particle dynamics method (J. Bonet Avalos, et al., Phys Chem Chem Phys, 2019, 21 , 24891-24911) to include chemical reactivity, denoted GenDPDE-RX. GenDPDE-RX provides a means of simulating chemical reactivity at the micro- and mesoscales, while exploiting the attributes of density- and temperature-dependent many-body force fields, which include improved transferability and scalability compared to two-body pairwise models. The GenDPDE-RX formulation considers intra-particle reactivity via a coarse-grain reactor construct. Extent-of-reaction variables assigned to each coarse-grain particle monitor the temporal evolution of the prescribed reaction mechanisms and kinetics assumed to occur within the particle. Descriptions of the algorithm, equations of motion, and numerical discretization are presented, followed by verification of the GenDPDE-RX method through comparison with reaction kinetics theoretical model predictions. Demonstrations of the GenDPDE-RX method are performed using constant-volume adiabatic heating simulations of three different reaction models, including both reversible and irreversible reactions, as well as multistep reaction mechanisms. The selection of the demonstrations is intended to illustrate the flexibility and generality of the method but is inspired by real material systems that span from fluids to solids. Many-body force fields using analytical forms of the ideal gas, Lennard-Jones, and exponential-6 equations of state are used for demonstration, although application to other forms of equation of states is possible. Finally, the flexibility of the GenDPDE-RX framework is addressed with a brief discussion of other possible adaptations and extensions of the method. |
