Identificador: TDX:332
Autores: Pérez Pellitero, Javier
Resumen:
Parallel with the increase of computer speed, in the last decade, molecular simulation techniques have emerged as important tools to predict physical properties of systems of industrial interest. These properties are essential in the chemical and petrochemical industries in order to perform process design, optimization, simulation and process control. The actual moderate cost of powerful computers converts molecular simulation into an excellent tool to provide predictions of such properties. In particular, the predictive capability of molecular simulation techniques becomes very important when dealing with extreme conditions of temperature and pressure as well as when toxic compounds are involved in the systems to be studied due to the fact that experimentation at such extreme conditions is difficult and expensive.Consequently, alternative processes must be considered in order to obtain the required properties. Chemical and petrochemical industries have made intensive use of thermophysical models including equations of state, activity coefficients models and corresponding state theories. These predictions present the advantage of providing good approximations with minimal computational needs. However, these models are often inadequate when only a limited amount of information is available to determine the necesary parameters, or when trying to reproduce complex fluid properties such as that of molecules which exhibit hydrogen bonding, polymers, etc. In addition, there is no way for dynamical properties to be estimated in a consistent manner.In this thesis, the HR and FSS techniques are combined with the main goal of extending the application of these methodologies to the calculation of the vaporliquid equilibrium and critical point of real mixtures. Before applying the methodologies to the real mixtures of industrial interest, the LennardJones fluid has been taken as a reference model and as a preliminary step. In this case, the predictions are affected only by the omnipresent statistical errors, but not by the accuracy of the model chosen to reproduce the behavior of the real molecules or the interatomic potential used to calculate the configurational energy of the system.The simulations have been performed in the grand canonical ensemble (GCMC)using the GIBBS code. Liquidvapor coexistences curves have been obtained from HR techniques for pure fluids and binary mixtures, while critical parameters were obtained from FSS in order to close the phase envelope of the phase diagrams. In order to extend the calculations to multicomponent systems modifications to the conventional HR techniques have been introduced in order to avoid the construction of histograms and the consequent need for large memory resources. In addition an alternative methodology known as the fourth order cumulant calculation, also known as the Binder parameter, has been implemented to make the location of the critical point more straightforward. In particular, we propose the use of the fourth order cumulant calculation considering two different possibilities: either the intersection of the Binder parameter for two different system sizes or the intersection of the Binder parameter with the known value for the system universality class combined with a FSS study. The development of transferable potential models able to describe the inter and intramolecular energies of the molecules involved in the simulations constitutes an important field in the improvement of Monte Carlo techniques. In the last decade, potential models, also referred to as force fields, have been developed for a wide range of compounds. One of the most common approaches for modeling hydrocarbons and other flexible molecules is the use of the unitedatoms model, where each chemical group is represented by one LennardJones center. This scheme results in a significant reduction of the computational time as compared to allatoms models since the number of pair interactions goes as the square of the number of sites. Improvements on the standard unitedatoms model, where typically a 612 LennardJones center of force is placed on top of the most significant atom, have been proposed. For instance, the AUA model consists of a displacement of the LennardJones centers of force towards the hydrogen atoms, converting the distance of displacement into a third adjustable parameter. In this thesis we have developed AUA 4 intermolecular potentials for three different families of compounds. The family of ethers is of great importance due to their applications as solvents. The other two families, thiophenes and alkanols, play an important roles in the oil and gas industry. Thiophene due to current and future environmental restrictions and alkanols due ever higher importance and presence of biofuels in this industry.
Repositori URV