Dissipative particle dynamics (DPD) study of hydrocarbon–water interfacial tension (IFT)

Methane-Water Interfacial Tension in Nanopores: A Dissipative Particle Dynamics Study

Imbibition of fracturing fluid in deep shale nanopores has a significant effect on shale gas production. One of the key parameters affecting imbibition is the interfacial tension of the methane-water system. However, studies on the methane-water interfacial tension in nanopores are very limited, and obtaining the accurate value of the methane-water interfacial tension at the nanoscale is difficult and time-consuming. In this work, a dissipative particle dynamics simulation model was built to study the methane-water interfacial tension in nanopores. This model provides reliable access to methane-water interfacial tension for deep shales under high-temperature, high-pressure conditions at low computation cost. It can be easily used to compute the methane-water interfacial tension in nanopores or the confined space in wide application scenarios. A sensitivity study of methane-water interfacial tension on a variety of factors was conducted. Results demonstrate that under high-pressure conditions, the increase in pressure leads to the rise of interfacial tension. When pressure increases from 20 to 120 MPa, interfacial tension rises from 0.0275 to 0.12 N/m, which contributes to the severe imbibition of fracturing fluid in deep shales. The confinement effect was observed by investigating the influence of pore size. Interfacial tension almost remains unchanged in pores smaller than 7 nm because most of the confined space is occupied by interface layer molecules in these pores. When pore size increases from 7 to 15 nm, the confinement effect is reduced. The interfacial tension experiences a growth from 0.1155 to 0.27 mN/m. Compared with pressure and pore size, the effect of temperature on interfacial tension can be neglected during deep shale gas production.

Dissipative Particle Dynamics-Based Simulation of the Effect of Asphaltene Structure on Oil-Water Interface Properties

Asphaltenes are the main substances that stabilize emulsions during the production, processing, and transport of crude oil. The purpose of this research is to investigate the process of asphaltenes forming interfacial films at the oil-water interface by means of dissipative particle dynamics (DPD) and the effect of asphaltenes of different structures on the oil-water interface during the formation of interfacial film. It is demonstrated that the thickness of the interfacial film formed at the oil-water interface gradually increases as the asphaltene concentration rises and the amount of asphaltene adsorbed at the oil-water interface gradually multiplies. Both the number and type of heteroatoms in asphaltenes affect the interfacial behavior of asphaltenes. The interface activity of asphaltenes can be enhanced by increasing the number of heteroatoms in the asphaltene, and the type of heteroatom affects as well the interfacial activity of the asphaltene as it affects the aggregation behavior of the asphaltene in the system. As the number of asphaltene aromatic rings increases, the oil-water interfacial tension (IFT) trends down gradually, while the effect of alkyl side chains on the reduction of IFT of asphaltenes is different, and asphaltenes with medium length alkyl side chains can reduce IFT more efficiently.

The effective interfacial tensions between pure liquids and rough solids: a coarse-grained simulation study

The effective solid-liquid interfacial tension (SL-IFT) between pure liquids and rough solid surfaces is studied through coarse-grained simulations. Using the dissipative particle dynamics method, we design solid-liquid interfaces, confining a pure liquid between two explicit solid surfaces with different roughness degrees. The roughness of the solid phase is characterized by Wenzel's roughness factor and the effective SL-IFT is reported as a function of it also. Two solid-liquid systems, different from each other by their solid-liquid repulsion strength, are studied to measure the effects caused by the surface roughness on the calculation of . It is found that the roughness changes the structure of the liquid, which is observed in the first layer of liquid near the solid. These changes are responsible for the effective SL-IFT increase, as surface roughness increases. Although there is a predominance of surface roughness in the calculation of it is found that the effective SL-IFT is directly proportional to the magnitude of the solid-liquid repulsion strength. The insights provided by these simulations suggest that the increase of Wenzel's roughness factor increases the number of effective solid-liquid interactions between particles, yielding significant changes in the local values of the normal and tangential components of the pressure tensor.

Comprehensive review of the interfacial behavior of water/oil/surfactant systems using dissipative particle dynamics simulation

A comprehensive understanding of interfacial behavior in water/oil/surfactant systems is critical to evaluating the performance of emulsions in various industries, specifically in the oil and gas industry. To gain fundamental knowledge regarding this interfacial behavior, atomistic methods, e.g., molecular dynamics (MD) simulation, can be employed; however, MD simulation cannot handle phenomena that require more than a million atoms. The coarse-grained mesoscale methods were introduced to resolve this issue. One of the most effective mesoscale coarse-grained approaches for simulating colloidal systems is dissipative particle dynamics (DPD), which bridges the gap between macroscopic time and length scales and molecular-scale simulation. This work reviews the fundamentals of DPD simulation and its progress on colloids and interface systems, especially surfactant/water/oil mixtures. The effects of temperature, salt content, a water/oil ratio, a shear rate, and a type of surfactant on the interfacial behavior in water/oil/surfactant systems using DPD simulation are evaluated. In addition, the obtained results are also investigated through the lens of the chemistry of surfactants and emulsions. The outcome of this comprehensive review demonstrates the importance of DPD simulation in various processes with a focus on the colloidal and interfacial behavior of surfactants at water-oil interfaces.

Coarse Grained Modeling of Multiphase Flows with Surfactants

Coarse-grained modeling methods allow simulations at larger scales than molecular dynamics, making it feasible to simulate multifluid systems. It is, however, critical to use model parameters that represent the fluid properties with fidelity under both equilibrium and dynamic conditions. In this work, dissipative particle dynamics (DPD) methods were used to simulate the flow of oil and water in a narrow slit under Poiseuille and Couette flow conditions. Large surfactant molecules were also included in the computations. A systematic methodology is presented to determine the DPD parameters necessary for ensuring that the boundary conditions were obeyed, that the oil and water viscosities were represented correctly, and that the velocity profile for the multifluid system agreed with the theoretical expectations. Surfactant molecules were introduced at the oil–water interface (sodium dodecylsulfate and octaethylene glycol monododecyl ether) to determine the effects of surface-active molecules on the two-phase flow. A critical shear rate was found for Poiseuille flow, beyond which the surfactants desorbed to form the interface forming micelles and destabilize the interface, and the surfactant-covered interface remained stable under Couette flow even at high shear rates.

Effect of sequence distribution of block copolymers on the interfacial properties of ternary mixtures: a dissipative particle dynamics simulation

In this paper, the dissipative particle dynamics (DPD) simulations method is used to study the effect of sequence distribution of block copolymers on the interfacial properties between immiscible homopolymers. Five block copolymers with the same composition but different sequence lengths are utilized for simulation. The sequence distribution is varied from the alternating copolymer to the symmetric diblock copolymer. Our simulations show that the efficiency of the block copolymer in reducing the interfacial tension is highly dependent on both the degree of penetration of the copolymer chain into the homopolymer phase and the number of copolymers at the interface per area. The linear block copolymers AB with the sequence length of τ = 8 could both sufficiently extend into the homopolymer phases and exhibit a larger number of copolymers at the interface per area. Thereby the copolymer with the sequence length τ = 8 is more effective in reducing the interfacial tension compared to that of diblock copolymers and the alternating copolymers at the same concentration. This work offers useful tips for copolymer compatibilizer selection at the immiscible homopolymer mixture interfaces.

Experimental and Theoretical Investigation of Glycol-Based Hydrogels through Waterflooding Processes in Oil Reservoirs Using Molecular Dynamics and Dissipative Particle Dynamics Simulation

Enhanced oil processing aims to retrieve petroleum fluids from depleted reservoirs after traditional processing. Hydrogels and polymeric macromolecules are considered effective displacing agents in oil reservoirs. In the current work, the authors used hydrophilic hydrogels based on poly(ethylene glycol)/poly(propylene glycol) (PEG/PPG) surfmers for oil displacement processes. Statistical modeling of the rheological properties at 80 °C for the two hydrogels indicates that the viscosity-shearing profile obeys the power-law model. Also, shear stress scanning follows the Herschel-Bulkley and the Bingham plastic models. The two hydrogels exhibit an initial yield stress owing to the formation of a three-dimensional (3D) structure at zero shearings. Furthermore, PEG and PPG hydrogels can retain the viscosity after a shear rate of 64.68 S-1. On the scale of surface activity, the two hydrogels exhibit higher surface areas (A m) of 0.1088 and 0.1058 nm2 and lower surface excess concentrations (Γm) of 1.529 and 1.567 × 1010 mol/cm2, respectively. A molecular dynamics (MD) simulation was conducted to explore the Flory-Huggins chi parameter, the solubility parameter, and the cohesive energy density. The results indicate a negative magnitude of chi parameter (χ ij ) for water and salt, which indicates that the two hydrogels have a good tendency toward saline formation water in the underground petroleum reservoir. Furthermore, the dissipative particle dynamics (DPD) was performed on a mesoscale to investigate the interfacial tension, the radius of gyration, the concentration profile, and the radial distribution function. The increased radius of gyration (R g) confirms that the two hydrogels are more overextended and can align perpendicularly toward the water/oil boundary. Experimental displacement was operated on a linear sandpack model using different slug concentrations. The oil recovery factor, the water-cut, and the differential pressure data during the flooding process were estimated as a function of the injected pore volume. The obtained results show that the oil recovery factor reaches 72 and 88% in the cases of PEG and PPG hydrogels at 80 °C with concentrations of 1.0 and 1.5 g/L, which reveals that both hydrogels are effective enhanced oil recovery (EOR) agents for the depleted reservoirs. This study establishes a new route that employs MD and DPD simulation in the field of enhanced oil recovery and the petroleum industry.

Compression and Ordering of Microgels in Monolayers Formed at Liquid-Liquid Interfaces: Computer Simulation Studies

Monolayers of polymer microgels adsorbed at the liquid interfaces were studied by dissipative particle dynamics simulations. The results demonstrated that the compressibility of the monolayers can be widely tuned by varying the cross-linking density of the microgels and their (in)compatibility with the immiscible liquids. In particular, the compression of the monolayers (increase of 2D concentration of the microgels) leads to the decrease of their lateral size. Herewith, the shape of the individual soft particles gradually changes from oblate (diluted 2D system) to nearly spherical (compressed monolayer). The polymer concentration profiles plotted along the normal to the interface reveal a nonmonotonous shape with a sharp maximum at the interface. This is a consequence of the shielding effect: saturation of the interface by monomer units of the subchains is driven by minimization of unfavorable contacts between the immiscible liquids and is opposed by elasticity of the network. The decrease of the interfacial tension upon concentration (compression) of the monolayer is quantified. It has been demonstrated that the interfacial tension significantly differs if the solubility of the polymer chains of the microgel network in the liquids changes. These results correlate well with experimental data. The examination of the microgels' crystalline ordering in monolayers demonstrated a nonmonotonous dependency on the compression degree (microgel concentration). Finally, the worsening of the solvent quality leads to the collapse of the microgels in monolayer and nonmonotonous behavior of the interfacial tension.

Effect of Aggregation and Adsorption Behavior on the Flow Resistance of Surfactant Fluid on Smooth and Rough Surfaces: A Many-Body Dissipative Particle Dynamics Study

To study the effect of surfactant on the resistance of wall-bound flow, the adsorption and aggregation behaviors of surfactant fluid on both smooth and groove-patterned rough surface are investigated through many-body dissipative particle dynamics (MDPD) simulation. The MDPD models of surfactants were carefully parametrized and have been validated to be able to simulate the aggregation and adsorption behavior of surfactants. The simulation results show that the surfactant in laminar flow can only increase the flow resistance on the smooth surface. On the rough surface, surfactant with strong adsorption performance on the channel wall shows a drag reduction effect at moderate concentration. The surfactant with weak adsorption properties can enhance the flow resistance, which is even more significant than that of those surfactants with no adsorption capacity. Although heating (high temperature) can generally reduce the viscosity and flow resistance of surfactant fluid, it would cause a poor drag reduction efficiency. It may arise from the destruction of the adsorption layer and the interruption of the fluid/boundary interface. Surfactant adsorption can tune the roughness of the fluid boundary on either the smooth or rough surface in a different manner, which turns out to be highly correlated to the change in flow resistance. Compared with the adsorption layer, surfactant in the bulk fluid makes a greater contribution to enhancing the flow resistance as the concentration rises. This study is expected to be helpful in guiding the application of surfactants on the micro- and nanoscale such as lab-on-a-chip nanodevices and EOR in a low-permeability porous medium.

Amphiphilic microgels adsorbed at oil-water interfaces as mixers of two immiscible liquids

Amphiphilic microgels adsorbed at an oil-water interface were studied by means of dissipative particle dynamics (DPD) simulations. The hydrophobic (A) and hydrophilic (B) monomer units in the polymer network are considered to be randomly distributed. Effects of the crosslinking density, interfacial tension between the liquids, their selectivity as solvents towards species A and B, and the degree of incompatibility between the A and B units on the internal microgel structure and distribution of the liquids are considered. The most important predictions are that (i) two immiscible liquids can homogeneously be mixed within the microgels and (ii) the adsorbed microgels contain a high fraction of the liquids (they are swollen at the interface). Simultaneous fulfillment of these two conditions can have a high impact on the design of new and efficient catalytic systems. In particular, such microgels can mix immiscible reactants dissolved in water and oil and trigger chemical reactions in the presence of a catalyst embedded into the microgel.

Microgels Adsorbed at Liquid-Liquid Interfaces: A Joint Numerical and Experimental Study

Soft particles display highly versatile properties with respect to hard colloids and even more so at fluid-fluid interfaces. In particular, microgels, consisting of a cross-linked polymer network, are able to deform and flatten upon adsorption at the interface due to the balance between surface tension and internal elasticity. Despite the existence of experimental results, a detailed theoretical understanding of this phenomenon is still lacking due to the absence of appropriate microscopic models. In this work, we propose an advanced modeling of microgels at a flat water/oil interface. The model builds on a realistic description of the internal polymeric architecture and single-particle properties of the microgel and is able to reproduce its experimentally observed shape at the interface. Complementing molecular dynamics simulations with in situ cryo-electron microscopy experiments and atomic force microscopy imaging after Langmuir-Blodgett deposition, we compare the morphology of the microgels for different values of the cross-linking ratios. Our model allows for a systematic microscopic investigation of soft particles at fluid interfaces, which is essential to develop predictive power for the use of microgels in a broad range of applications, including the stabilization of smart emulsions and the versatile patterning of surfaces.

Application of a Novel Coarse-Grained Soil Organic Matter Model in the Environment

Soil organic matter (SOM) is ubiquitous in the environment. Intensive efforts have been made to find effective ways to assess the interaction of SOM with contaminants since such interactions are one of the important criteria used to evaluate the migration, persistency and bioavailability of chemicals in the environment. This study aims to extend the application of coarse-grained (CG) dissipative particle dynamics (DPD) to the water/SOM system and predict contaminant mobility in the system. The CG model was based on the Vienna Soil-Organic-Matter Modeler, which can generate flexible condensed-phase models of SOM. A series of DPD simulations was performed to investigate the mobility of perfluorinated sulfonic acids (PFSAs) and hexachlorobutadiene (HCBD). The results indicated that the mobility of PFSAs decreased with increasing length in the carbon chain. In addition, HCBD and hexachlorobenzene (HCB) have similar diffusion coefficients, indicating analogous behavior in SOM. Moreover, water-containing SOM layers may reflect a more realistic situation. This work, coupling the CG method with DPD simulation, provides a new high-efficiency tool to assess the behavior of contaminants in the environment.

Simulations of Interfacial Tension of Liquid-Liquid Ternary Mixtures Using Optimized Parametrization for Coarse-Grained Models

In this work, liquid-liquid systems are studied by means of coarse-grained Monte Carlo simulations (CG-MC) and Dissipative Particle Dynamics (DPD). A methodology is proposed to reproduce liquid-liquid equilibrium (LLE) and to provide variation of interfacial tension (IFT), as a function of the solute concentration. A key step is the parametrization method based on the use of the Flory-Huggins parameter between DPD beads to calculate solute/solvent interactions. Parameters are determined using a set of experimental compositional data of LLE, following four different approaches. These approaches are evaluated, and the results obtained are compared to analyze advantages/disadvantages of each one. These methodologies have been compared through their application on six systems: water/benzene/1,4-dioxane,water/chloroform/acetone, water/benzene/acetic acid, water/benzene/2-propanol, water/hexane/acetone, and water/hexane/2-propanol. CG-MC simulations in the Gibbs (NVT) ensemble have been used to check the validity of parametrization approaches for LLE reproduction. Then, CG-MC simulations in the osmotic (μ_soluteN_solventP _zzT) ensemble were carried out considering the two liquid phases with an explicit interface. This step allows one to work at the same bulk concentrations as the experimental data by imposing the precise bulk phase compositions and predicting the interface composition. Finally, DPD simulations were used to predict IFT values for different solute concentrations. Our results on variation of IFT with solute concentration in bulk phases are in good agreement with experimental data, but some deviations can be observed for systems containing hexane molecules.

Amphiphilic nanosheet self-assembly at the water/oil interface: computer simulations

The self-assembly of amphiphilic Janus triangular-plates at the water/oil interface is simulated for the first time.