Meso- and molecular-scale modeling to provide new insights into interfacial and structural properties of hydrocarbon/water/surfactant systems

Review of the Interfacial Structure and Properties of Surfactants in Petroleum Production and Geological Storage Systems from a Molecular Scale Perspective

Surfactants play a crucial role in tertiary oil recovery by reducing the interfacial tension between immiscible phases, altering surface wettability, and improving foam film stability. Oil reservoirs have high temperatures and high pressures, making it difficult and hazardous to conduct lab experiments. In this context, molecular dynamics (MD) simulation is a valuable tool for complementing experiments. It can effectively study the microscopic behaviors (such as diffusion, adsorption, and aggregation) of the surfactant molecules in the pore fluids and predict the thermodynamics and kinetics of these systems with a high degree of accuracy. MD simulation also overcomes the limitations of traditional experiments, which often lack the necessary temporal-spatial resolution. Comparing simulated results with experimental data can provide a comprehensive explanation from a microscopic standpoint. This article reviews the state-of-the-art MD simulations of surfactant adsorption and resulting interfacial properties at gas/oil-water interfaces. Initially, the article discusses interfacial properties and methods for evaluating surfactant-formed monolayers, considering variations in interfacial concentration, molecular structure of the surfactants, and synergistic effect of surfactant mixtures. Then, it covers methods for characterizing microstructure at various interfaces and the evolution process of the monolayers' packing state as a function of interfacial concentration and the surfactants' molecular structure. Next, it examines the interactions between surfactants and the aqueous phase, focusing on headgroup solvation and counterion condensation. Finally, it analyzes the influence of hydrophobic phase molecular composition on interactions between surfactants and the hydrophobic phase. This review deepened our understanding of the micro-level mechanisms of oil displacement by surfactants and is beneficial for screening and designing surfactants for oil field applications.

Effects of concentration and chain length of the sequence copolymer on interfacial properties of homopolymers/sequence copolymers ternary blends: A DPD simulation study

The effect of the concentration and chain length of the copolymer AB with sequence length τ = 8 on the interfacial properties of the ternary mixtures A₁₀/AB/B₁₀ are investigated by the dissipative particle dynamics (DPD) simulations. It is found that: i) As the copolymer concentration varies from 0.05 to 0.15, increasing the copolymer enrichment at the center of the interface enlarges the interface width ω and reduces the interfacial tension. However, as the concentration of the sequence copolymers further increases to 0.2, because the interface has formed micelles and the micellization could lower the efficiency of copolymers as a compatibilizer, the interfacial tension exhibits a slightly increase; ii) elevating the copolymer chain length, the copolymer volumes vary from a cylinder shape to a pancake shape. The blends of the copolymer with chain length N_cp = 24 exhibit a wider interfacial width w and a lower interfacial tension γ, which indicates that the sequenced copolymer N_cp = 24 exhibits a better performance as the compatibilizers. This study illustrates the correlations between the reduction in interfacial tension produced by the sequence copolymers and their molecular parameters, which guide a rational design of an efficient compatibilizer.

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.

DPD Study on the Interfacial Properties of PEO/PEO-PPO-PEO/PPO Ternary Blends: Effects of Pluronic Structure and Concentration

Using the method of dissipative particle dynamics (DPD) simulations, we investigated the interfacial properties of PEO/PEO-PPO-PEO/PPO ternary blends composed of the Pluronics L64(EO₁₃PO₃₀EO₁₃), F68(EO₇₆PO₂₉EO₇₆), F88(EO₁₀₄PO₃₉EO₁₀₄), or F127(EO₁₀₆PO₇₀EO₁₀₆) triblock copolymers. Our simulations show that: (i) The interfacial tensions (γ) of the ternary blends obey the relationship γF68 < γL64 < γF88 < γF127, which indicates that triblock copolymer F68 is most effective in reducing the interfacial tension, compared to L64, F88, and F127; (ii) For the blends of PEO/L64/PPO and the F64 copolymer concentration ranging from c_cp = 0.2 to 0.4, the interface exhibits a saturation state, which results in the aggregation and micelle formation of F64 copolymers added to the blends, and a lowered efficiency of the L64 copolymers as a compatibilizer, thus, the interfacial tension decreases slightly; (iii) For the blends of PEO/F68/PPO, elevating the Pluronic copolymer concentration can promote Pluronic copolymer enrichment at the interfaces without forming the micelles, which reduces the interfacial tension significantly. The interfacial properties of the blends contained the PEO-PPO-PEO triblock copolymer compatibilizers are, thus, controlled by the triblock copolymer structure and the concentration. This work provides important insights into the use of the PEO-PPO-PEO triblock copolymer as compatibilizers in the PEO and PPO homopolymer blend systems.

Effects of Repulsion Parameter and Chain Length of Homopolymers on Interfacial Properties of An/Ax/2BxAx/2/Bm Blends: A DPD Simulation Study

We explored the effects of the repulsion parameter (aAB) and chain length (N_HA or N_HB) of homopolymers on the interfacial properties of A_n/A_x/2B_xA_x/2/B_m ternary polymeric blends using dissipative particle dynamics (DPD) simulations. Our simulations show that: (i) The ternary blends exhibit the significant segregation at the repulsion parameter (aAB = 40). (ii) Both the interfacial tension and the density of triblock copolymer at the center of the interface increase to a plateau with increasing the homopolymer chain length, which indicates that the triblock copolymers with shorter chain length exhibit better performance as the compatibilizers for stabilizing the blends. (iii) For the case of N_HA = 4 (chain length of homopolymers A_n) and N_HB (chain length of homopolymers B_m) ranging from 16 to 64, the blends exhibit larger interfacial widths with a weakened correlation between bead A_n and B_m of homopolymers, which indicates that the triblock copolymer compatibilizers (A_x/2B_xA_x/2) show better performance in reducing the interfacial tension. The effectiveness of triblock copolymer compatibilizers is, thus, controlled by the regulation of repulsion parameters and the homopolymer chain length. This work raises important considerations concerning the use of the triblock copolymer as compatibilizers in the immiscible homopolymer blend systems.

Dissipative Particle Dynamics Study on Interfacial Properties of Symmetric Ternary Polymeric Blends

We investigated the interfacial properties of symmetric ternary A_n/A_mB_m/B_n and A_n/A_m/2B_mA_m/2/B_n polymeric blends by means of dissipative particle dynamics (DPD) simulations. We systematically analyzed the effects of composition, chain length, and concentration of the copolymers on the interfacial tensions, interfacial widths, and the structures of each polymer component in the blends. Our simulations show that: (i) the efficiency of the copolymers in reducing the interfacial tension is highly dependent on their compositions. The triblock copolymers are more effective in reducing the interfacial tension compared to that of the diblock copolymers at the same chain length and concentration; (ii) the interfacial tension of the blends increases with increases in the triblock copolymer chain length, which indicates that the triblock copolymers with a shorter chain length exhibit a better performance as the compatibilizers compared to that of their counterparts with longer chain lengths; and (iii) elevating the triblock copolymer concentration can promote copolymer enrichment at the center of the interface, which enlarges the width of the phase interfaces and reduces the interfacial tension. These findings illustrate the correlations between the efficiency of copolymer compatibilizers and their detailed molecular parameters.

Computational investigation of a switchable emulsion stabilized by the mixture of a surfactant and tertiary amine

Molecular dynamics simulations were performed to investigate the CO2-responsiveness of an oil-in-water (O/W) emulsion stabilized by sodium oleate (NaOA) with a tertiary amine additive, named pentamethyl diethylenetriamine (PMA). The simulated results were in accordance with the experimental observations. That is, the surfactant NaOA itself can stabilize dodecane/water emulsions in aqueous solution, while the CO2-reponsiveness was strongly related to the added PMA. The electroneutral PMA molecules preferred to be located in the core region of the droplets. Thus, under the same conditions, the size of the droplet containing PMA is predictably larger than that without PMA. The increased extent of the charged surfactant headgroups distribution can increase the electrostatic repulsion between the droplets in the emulsion solution, which is the important reason why a much more stable emulsion is obtained by adding PMA. When PMA molecules were protonated to PMA2+ by bubbling CO2, they migrated from the interior to the surface of the droplets under electrostatic attraction, forming ion pairs with OA-. The binding between PMA2+ and OA- made the distribution of the surfactants very concentrated on the droplet surface, leading to large hydrophobic areas exposed to water. Besides, the hydration interactions of OA- headgroups decreased because they were covered by PMA2+. The calculated potential of mean force (PMF) confirmed that the electrostatic repulsion between droplets was crucial for the emulsion stabilization.