Charlles R.A. Abreu

Influence of particle shape on the packing and on the segregation of spherocylinders via Monte Carlo simulations

Knowledge of the properties of granular materials is important for efficient and safe design of industrial equipment. In this work, the Monte Carlo method is used for simulating granular systems of spherocylindrical particles. After presenting an overview of such method and of overlap detection in systems of hard spherocylinders, the application of the method for granular systems is discussed. Then, porosities, calculated for simulated monodispersed beds, are presented as functions of the particle elongation. Next, results for vibration-induced segregation of binary mixtures of spherocylinders with identical volume and density, but with different elongations, are evaluated, showing the influence of the particle shape on this phenomenon. Finally, effects of size and shape on such segregation are contrasted using results with simultaneous variation of particle volume and elongation.

Transport properties of carbon dioxide and methane from molecular dynamics simulations

Transport properties of carbon dioxide and methane are predicted for temperatures between (273.15 and 573.15) K and pressures up to 800 MPa by molecular dynamics simulations. Viscosities and thermal conductivities were obtained through the Green-Kubo formalism, whereas the Einstein relation was used to provide self-diffusion coefficient estimates. The differences in property predictions due to the force field nature and parametrization were investigated by the comparison of seven different CO2 models (two single-site models, three rigid three-site models, and two fully flexible three-site models) and three different CH4 models (two single-site models and one fully flexible five-site model). The simulation results show good agreement with experimental data, except for thermal conductivities at low densities. The molecular structure and force field parameters play an important role in the accuracy of the simulations, which is within the experimental deviations reported for viscosities and self-diffusion coefficients considering the most accurate CO2 and CH4 models studied. On the other hand, the molecular flexibility does not seem to improve accuracy, since the explicit account of vibrational and bending degrees of freedom in the CO2 flexible models leads to slightly less accurate results. Nonetheless, the use of a correctional term to account for vibrational modes in rigid models generally improves estimations of thermal conductivity values. At extreme densities, the caging effect observed with single-site representations of the molecules restrains mobility and leads to an unphysical overestimation of viscosities and, conversely, to the underestimation of self-diffusion coefficients. This result may help to better understand the limits of applicability of such force fields concerning structural and transport properties of dense systems.

A Monte Carlo simulation of the packing and segregation of spheres in cylinders

In this work, the Monte Carlo method (MC) was extended to simulate the packing and segregation of particles subjected to a gravitational field and confined inside rigid walls. The method was used in systems containing spheres inside cylinders. The calculation of void fraction profiles in both the axial and radial directions was formulated, and some results are presented. In agreement with experimental data, the simulations show that the packed beds present structural ordering near the cylindrical walls up to a distance of about 4 particle diameters. The simulations also indicate that the presence of the cylindrical wall does not seem to have a strong effect on the gravitational segregation phenomenon.

A phase stability analysis of the combinatorial term of the UNIQUAC model

The conditions for the onset of liquid-phase instability in the combinatorial contribution of UNIQUAC are discussed. A peculiar and non-physical aspect of this predicted instability is that it is independent of temperature. It is recognized that by renormalizing the r and q parameters, this type of phase transition can be prevented.

Molecular dynamics with rigid bodies: Alternative formulation and assessment of its limitations when employed to simulate liquid water

Sets of atoms collectively behaving as rigid bodies are often used in molecular dynamics to model entire molecules or parts thereof. This is a coarse-graining strategy that eliminates degrees of freedom and supposedly admits larger time steps without abandoning the atomistic character of a model. In this paper, we rely on a particular factorization of the rotation matrix to simplify the mechanical formulation of systems containing rigid bodies. We then propose a new derivation for the exact solution of torque-free rotations, which are employed as part of a symplectic numerical integration scheme for rigid-body dynamics. We also review methods for calculating pressure in systems of rigid bodies with pairwise-additive potentials and periodic boundary conditions. Finally, simulations of liquid phases, with special focus on water, are employed to analyze the numerical aspects of the proposed methodology. Our results show that energy drift is avoided for time step sizes up to 5 fs, but only if a proper smoothing is applied to the interatomic potentials. Despite this, the effects of discretization errors are relevant, even for smaller time steps. These errors induce, for instance, a systematic failure of the expected equipartition of kinetic energy between translational and rotational degrees of freedom.

MONTE CARLO SIMULATIONS OF THE ADSORPTION OF DIMERS ON STRUCTURED HETEROGENEOUS SURFACES

The effect of surface topography upon the adsorption of dimer molecules is analyzed by means of grand canonical ensemble Monte Carlo simulations. Heterogeneous surfaces were assumed to consist of a square lattice containing active sites with two different energies. These were distributed in three different configurations: a random distribution of isolated sites; a random distribution of grains with four high-energy sites; and a random distribution of grains with nine high-energy sites. For the random distribution of isolated sites, the results are in good agreement with the molecular simulations performed by Nitta et al. (1997). In general, the comparison with theoretical models shows that the Nitta et al. (1984) isotherm presents good predictions of dimer adsorption both on homogeneous and heterogeneous surfaces with sites having small differences in characteristic energies. The molecular simulation results also show that the energy topology of the solid surfaces plays an important role in the adsorption of dimers on solids with large differences in site energies. For these cases, the Nitta et al. model does not describe well the data on dimer adsorption on random heterogeneous surfaces (grains with one acid site), but does describe reasonably well the adsorption of dimers on more patchwise heterogeneous surfaces (grains with nine acid sites).