Emeric Bourasseau

Mesoscopic simulations of shock-to-detonation transition in reactive liquid high explosive

An extension of the model described in a previous work (see Maillet J. B. et al., EPL, 78 (2007) 68001) based on Dissipative Particle Dynamics is presented and applied to a liquid high explosive (HE), with thermodynamic properties mimicking those of liquid nitromethane. Large scale nonequilibrium simulations of reacting liquid HE with model kinetic under sustained shock conditions allow a better understanding of the shock-to-detonation transition in homogeneous explosives. Moreover, the propagation of the reactive wave appears discontinuous since ignition points in the shocked material can be activated by the compressive waves emitted from the onset of chemical reactions.

Molecular based equation of state for shocked liquid nitromethane

An approach is proposed to obtain the equation of state of unreactive shocked liquid nitromethane. Unlike previous major works, this equation of state is not based on extended integration schemes [P.C. Lysne, D.R. Hardesty, Fundamental equation of state of liquid nitromethane to 100 kbar, J. Chem. Phys. 59 (1973) 6512]. It does not follow the way proposed by Winey et al. [J.M. Winey, G.E. Duvall, M.D. Knudson, Y.M. Gupta, Equation of state and temperature measurements for shocked nitromethane, J. Chem. Phys. 113 (2000) 7492] where the specific heat C(v), the isothermal bulk modulus B(T) and the coefficient of thermal pressure (deltaP/deltaT)(v) are modeled as functions of temperature and volume using experimental data. In this work, we compute the complete equation of state by microscopic calculations. Indeed, by means of Monte Carlo molecular simulations, we have proposed a new force field for nitromethane that lead to a good description of shock properties [N. Desbiens, E. Bourasseau, J.-B. Maillet, Potential optimization for the calculation of shocked liquid nitromethane properties, Mol. Sim. 33 (2007) 1061; A. Hervouët, N. Desbiens, E. Bourasseau, J.-B. Maillet, Microscopic approaches to liquid nitromethane detonation properties, J. Phys. Chem. B 112 (2008) 5070]. Particularly, it has been shown that shock temperatures and second shock temperatures are accurately reproduced which is significative of the quality of the potential. Here, thermodynamic derivative properties are computed: specific heats, Grüneisen parameter, sound velocity among others, along the Hugoniot curve. This work constitutes to our knowledge the first determination of the equation of state of an unreactive shocked explosive by molecular simulations.

Surface tension of spherical drops from surface of tension

The determination of surface tension of curved interfaces is a topic that raised many controversies during the last century. Explicit liquid-vapor interface modelling (ELVI) was unable up to now to reproduce interfacial behaviors in drops due to ambiguities in the mechanical definition of the surface tension. In this work, we propose a thermodynamic approach based on the location of surface of tension and its use in the Laplace equation to extract the surface tension of spherical interfaces from ELVI modelling.

Communication: Slab thickness dependence of the surface tension: Toward a criterion of liquid sheets stability

Microscopic Monte Carlo simulations of liquid sheets of copper and tin have been performed in order to study the dependence of the surface tension on the thickness of the sheet. It results that the surface tension is constant with the thickness as long as the sheet remains in one piece. When the sheet is getting thinner, holes start to appear, and the calculated surface tension rapidly decreases with thickness until the sheet becomes totally unstable and forms a cylinder. We assume here that this decrease is not due to a confinement effect as proposed by Werth et al. [Physica A 392, 2359 (2013)] on Lennard-Jones systems, but to the appearance of holes that reduces the energy cost of the surface modification. We also show in this work that a link can be established between the stability of the sheet and the local fluctuations of the surface position, which directly depends on the value of the surface tension. Finally, we complete this study by investigating systems interacting through different forms of Lennard-Jones potentials to check if similar conclusions can be drawn.

Surface tension and long range corrections of cylindrical interfaces

The calculation of the surface tension of curved interfaces has been deeply investigated from molecular simulation during this last past decade. Recently, the thermodynamic Test-Area (TA) approach has been extended to the calculation of surface tension of curved interfaces. In the case of the cylindrical vapour-liquid interfaces of water and Lennard-Jones fluids, it was shown that the surface tension was independent of the curvature of the interface. In addition, the surface tension of the cylindrical interface is higher than that of the planar interface. Molecular simulations of cylindrical interfaces have been so far performed (i) by using a shifted potential, (ii) by means of large cutoff without periodic boundary conditions, or (iii) by ignoring the long range corrections to the surface tension due to the difficulty to estimate them. Indeed, unlike the planar interfaces there are no available operational expressions to consider the tail corrections to the surface tension of cylindrical interfaces. We propose here to develop the long range corrections of the surface tension for cylindrical interfaces by using the non-exponential TA (TA2) method. We also extend the formulation of the Mecke-Winkelmann corrections initially developed for planar surfaces to cylindrical interfaces. We complete this study by the calculation of the surface tension of cylindrical surfaces of liquid tin and copper using the embedded atom model potentials.

Constant entropy sampling and release waves of shock compressions.

We present or recall several equilibrium methods that allow one to compute isentropic processes, either during the compression or the release of the material. These methods are applied to compute the isentropic release of a shocked monoatomic liquid at high pressure and temperature. Moreover, equilibrium results of isentropic release are compared to the direct nonequilibrium simulation of the same process. We show that due to the viscosity of the liquid but also to nonequilibrium effects, the release of the system is not strictly isentropic.

Parameter Optimization for Charge Equilibration Method in Molecular Simulations

In order to build a complete potential model for classical molecular dynamic simulations of dense liquids, a new optimization method is proposed and applied to the determination of transferable parameters for charge equilibration method. To make sure to obtain parameters which only describe coulomb interactions, the minimization of the error function is performed over a data set constituted of pure electrostatic results taken from ab initio calculations. This new procedure has been applied to liquid HF and to liquid nitromethane.