A Molecular Dynamics code has been developed which allows including two immiscible fluids and surfactant molecules. The behaviour of this system can, thus, be observed at the microscopic level. The evolution towards separation and the formation of an interface have been monitored. It was found that the surfactant chains distribute themselves across this interface and this distribution results in a decrease of the interfacial tension. In many cases, the self-assembly of micelles can also be seen. The effect of the surfactant architecture was studied by employing Dissipative Particle Dynamics, a mesoscopic tool offered in the Materials Studio software. The adsorption of ''palm tree'' and random copolymers on a specific wall (polymer particle) in the presence of a solvent (CO2) was modelled. Values of interaction parameters were screened and the ones that resulted in adsorption were selected. Depending on the surfactant architecture, the volume fraction of the surfactant in the mixture was varied. For the case of palm tree surfactants, it was found that that bilayer structures were formed onto the polymer particle surface. This could result in particle coagulation and, eventually, in the production of a poor quality product.
The phase equilibrium properties of CO2, VDF and PVDF were studied by the Gibbs Ensemble Monte-Carlo method, a simulation technique at the molecular level. The code developed is generic and can also be used for mixtures and/or polymeric chains. The phase behaviour of CO2 was modelled by using force fields from the literature, under both normal and supercritical conditions, and yielded very good agreement with the available experimental data. Intermolecular potentials for VDF have also been developed, which describe very well the corresponding phase diagrams. Using the same method, the properties of the mixture (CO2+VDF) have been predicted under selected conditions and it was found that CO2 was completely miscible with VDF for the range of pressures and temperatures studied. The latest literature force fields were tested for the chemical species under study. The COMPASS all-purpose force field produced excellent results for scCO2. For VDF, the simulation results under PCFF were in perfect agreement with the experimental data. Both force fields yielded the correct PVT behaviour for the PVDF polymer. Summation methods for the difficult-to-handle columbic forces were also tested on the grounds of accuracy and computational cost. It was found that the group-based cut off method provided accurate results and was much faster than the conventional Ewald method.
Possible end-users of the software codes developed can be found in the polymer as well as in the chemical industry in general.