Molecular Dynamics Simulations of the Oil-Detachment from the Hydroxylated Silica Surface: Effects of Surfactants, Electrostatic Interactions, and Water Flows on the Water Molecular Channel Formation

Experimental and Molecular Dynamics Simulation to Investigate Oil Adsorption and Detachment from Sandstone/Quartz Surface by Low-Salinity Surfactant Brines

In this study, we explore the impact of monovalent (NaCl) and divalent (CaCl2) brines, coupled with sodium dodecyl sulfate (SDS) surfactant at varying low concentrations, on the detachment and displacement of oil from sandstone rock surfaces. Employing the sessile drop method and molecular dynamics simulations, we scrutinize the behavior of the brine solutions. Our findings reveal that both low salinity and low-salinity surfactant solutions induce a gradual shift in rock wettability toward a more water-wet state. This wettability transformation is not instantaneous but evolves over time, as observed through meticulous molecular motion analyses. Through contact angle measurements and molecular dynamics simulations, we delve into the molecular motion at subpore and micropore scales on sandstone/quartz surfaces. The adsorption of surface-active agents from the oil to the oil-brine interface results in a reduced interfacial tension, significantly contributing to oil displacement. Notably, low salinity concentrations ranging from 1000 to 10,000 ppm exhibit the lowest contact angles within 30 min across all solutions. However, higher concentrations deviate from this declining trend, especially with divalent ions like Ca2+, which bridge polar molecules onto the rock surface, resulting in an increased oil-wetting state. This research unveils the intricate molecular motions involved in employing low-salinity surfactant solutions for oil detachment from surfaces. Furthermore, it provides valuable insights into the underlying forces driving oil detachment and wettability alteration.

Insights into Molecular Dynamics and Oil Extraction Behavior of the Polymeric Surfactant in a Multilayered Heterogeneous Reservoir

Interlayer heterogeneity, an inevitable and complex challenge during water flooding, seriously constrains the spread of the sweep region and oil recovery enhancement in multilayered heterogeneous reservoirs. To overcome this challenge, a novel polymeric surfactant, having an excellent performance in the reduction of interfacial tension (IFT) and the increase of viscosity of displacing fluid, is applied for enlarging the sweep resonance and increasing the oil washing efficiency. Through the molecular dynamics (MD) simulation, the molecular distributing mechanisms of the polymeric surfactant at the oil-water interface are analyzed to provide the theoretical basis for explaining the microscopic mechanism of oil extraction. To directly reflect the microscopic behavior of oil extraction, multiple transparent sand-packed models are designed to investigate the flowing behavior of different fluids and the extracted mechanisms of the remaining oil in both pore and macroscales. The multilayered heterogeneous reservoirs consisting of high-, moderate-, and low-permeability layers are fabricated to represent a heterogeneous characteristic. The recognition from the visual experiment and MD simulation can study the performance control, the extracting performance of the remaining oil, and the expression of the displacing front from different perspectives. The results from MD simulation demonstrate that the polymeric surfactant can promote the disintegration of the remaining oil and enhance its mobility. The experimental results indicate that the sweep efficiency is restricted by viscous fingering and tongue advance. Through the analysis of mathematical models, the rising mobility ratio and the location of the displacing front have a strong positive relationship with viscous fingering and tongue advance, which can reasonably explain the plugging performance of the polymeric surfactant, greatly improving the sweeping effect of the whole reservoir. Moreover, the Marangoni effect generated by the IFT gradient can induce the transformation of interfacial energy to displacement kinetic energy by the emulsification of the oil-water interface so that the remaining oil in the blind-end pore can be effectively extracted. However, by comparing data from image quantification techniques and production dynamic performance, the sweep efficiency (484%) was significantly greater than that of oil recovery (300%), demonstrating that the expanded sweep effect still plays a dominant role in the extraction of remaining oil after polymeric surfactant flooding. This study provides a novel plugging and effective washing agent that is expected to be an excellent and comprehensive method for solving the problem of low oil recovery in multilayered heterogeneous reservoirs.

Molecular Simulation of Surfactant Displacement of Residual Oil in Nanopores: Formation of Water Channels and Electrostatic Interaction

The water-oil-rock system's surfactant and electrostatic interactions are essential for removing oil droplets from rock substrates. Our work illustrates the impact of surface charge on the oil contact angle in an ideal system comprising silica, water, and dodecane; smaller contact angles are observed for more polar substrates. Modifying the polarity of the model silica surface allows for the observation of the creation of heteromolecule channels and the process of stripping crude oil while accounting for the impacts of water flow and different types of surfactant molecules. In solutions containing ionic surfactants, the injection and diffusion of water molecules between the oil layer and the silica substrate are facilitated by the disturbance of the oil molecules by the surfactant molecules. By comparing different surfactants in water flow, the characterization of water molecular channels and the stripping process of crude oil can be observed. The disruption of oil molecules by the surfactant molecules has been found to enhance the injection and diffusion of water molecules between the oil layer and the silica substrate in solutions containing ionic surfactants. The size of the contact angle and the extension of the water channel are simultaneously greatly influenced by the surfactant's molecular characteristics and the substrate's polarity. These simulation results show that several factors influence the process of water molecule channel creation that water molecules diffuse, and the detachment of oil from the silica substrate is facilitated by the migration of surfactants to the bottom of the oil molecule and the electrostatic interactions between the water molecules and the silica substrate.

Molecular Alignment-Mediated Stick-Slip Poiseuille Flow of Oil in Graphene Nanochannels

The flow behavior of oil in nanochannels has attracted extensive attention for oil transport applications. In most, if not all, of the prior theoretical simulations, oil molecules were observed to flow steadily in nanochannels under pressure gradients. In this study, non-equilibrium molecular dynamics simulations are conducted to simulate the Poiseuille flow of oil with three different hydrocarbon chain lengths in graphene nanochannels. Contrary to the conventional perception of steady flows of oil in nanochannels, we find that oil molecules with the longest hydrocarbon chain (i.e., n-dodecane) exhibit notable stick-slip flow behavior. An alternation between the high average velocity of n-dodecane in the slip motion and the low average velocity in the stick motion is observed, with a drastic, abrupt velocity jolt of up to 40 times occurring at the transition in a stick-slip motion. Further statistical analyses show that the stick-slip flow behavior of n-dodecane molecules originates from the molecular alignment change of oil near the graphene wall. The molecular alignment of n-dodecane shows different statistical distributions under stick and slip motion states, leading to significant changes of friction forces and thus notable velocity fluctuations. This work provides new insights into the Poiseuille flow behavior of oil in graphene nanochannels and may offer useful guidelines for other mass transport applications.

Oil adsorbing and transporting surfaces: a simulative determination of parameters for bionic functional textiles

Certain superhydrophobic plants, such asSalvinia molesta, are able to adsorb oil films from water surfaces and thus separate the oil from the water. There are first attempts to transfer this phenomenon to technical surfaces, but the functional principle and the influence of certain parameters are not yet fully understood. The aim of this work is to understand the interaction behavior between biological surfaces and oil, and to define design parameters for transferring the biological model to a technical textile. This will reduce the development time of a biologically inspired textile. For this purpose, the biological surface is transferred into a 2D model and the horizontal oil transport is simulated in Ansys Fluent. From these simulations, the influence of contact angle, oil viscosity and fiber spacing/diameter ratio was quantified. The simulation results were verified with transport tests on spacer fabrics and 3D prints. The values obtained serve as a starting point for the development of a bio-inspired textile for the removal of oil spills on water surfaces. Such a bio-inspired textile provides the basis for a novel method of oil-water separation that does not require the use of chemicals or energy. As a result, it offers great added value compared to existing methods.

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.

A mechanism study of sodium dodecylbenzene sulfonate on oil recovery: effect of branched chain structure

The effects of the branched structures of SDBS molecules on oil recovery are investigated by molecular dynamics method. The relative density of oil molecules shows that SDBS molecule with benzene ring located near the center of alkyl chain has the best effect on oil displacement. Dynamic trajectories show that the water and SDBS molecules gradually occupy the calcite surface and replace the oil droplets. In this process, more water molecules gather nearby the polar groups of SDBS, indicating that the polar group has a significant effect on the water infiltration and the formation of water channels. Contact angle between SDBS molecule and calcite surface indicates that compared to straight chains, the branched structure tends to spread on the calcite interface. Moreover, adsorption energies of the simulation systems further prove that as the aromatic ring is closer to the middle of the alkyl chain, the oil displacement effect is better.

Effects of Surfactant Charge and Molecular Structure on Wettability Alteration of Calcite: Insights from Molecular Dynamics Simulations

Wettability alteration of oil-wet calcite by surfactants was studied by means of molecular dynamics (MD) simulations. The simulations use the recently developed model for positively charged calcite surface, whose oil-wet state originates from binding of negatively charged carboxylate molecules contained in the oil, consistently with the bulk of the available experimental data. The ability to alter the surface wettability, which can be directly quantified by the release of the surface-bound carboxylates, is tested for nine different surfactants of all charge types-cationic, anionic, nonionic, and zwitterionic-and compared to that of brine. It was found that only the cationic surfactants are able to detach the organic carboxylates more efficiently than brine, while the neutral and anionic surfactants do not seem to have any measurable effect on the wettability. The outperformance of the cationic surfactants is generally consistent with the majority of previously published experimental observations. The data also point toward a consistently better performance of single-tailed cationic surfactants over the two-tailed structure. Molecular mechanism of the wettability alteration by different types of surfactants is discussed, along with the implications of the results for the design of new surfactant formulations for the enhanced oil recovery.

Are all-atom any better than united-atom force fields for the description of liquid properties of alkanes?

Alkanes are a fundamental part in empirical force fields (FF) not only due to their technological relevance, but also due to the prevalence of alkane moieties in organic molecules, e.g., compounds containing a saturated carbon chain. Therefore, a good description of alkane interactions is crucial for determining the quality of a FF. In this study, the performance of 12 empirical force fields (FF) was evaluated in the context of reproducing liquid properties of alkanes. More specifically, n-octane was chosen as a reference compound since it is a liquid in a broad temperature range and it has numerous experimental data for thermodynamic, transport, and structural properties, as well as for their temperature dependencies. A normalized root-mean-square deviation (NRMSD) analysis was used to rank the force fields in their ability to reproduce the experimental data. Five out of the six best force fields considered were united-atom models. The GROMOS force field showed the smallest deviation in terms of NRMSD, followed by TRAPPE-EH, NERD, CHARMM-UA, TRAPPE-UA, and OPLS-UA. This overall better performance of the united-atom force fields indicates that complexity does not always bring quality.

Interfacial and molecular interactions between fractions of heavy oil and surfactants in porous media: Comprehensive review

The oil production by the natural energy in oil reservoirs is decreasing gradually. Only 25-30% of the world's reservoirs can be produced naturally, and different methods are employed to recover the remaining oil. The use of surfactants is one of the promising methods for unlocking the residual oil after natural depletion. In such a method, one of the main challenges is to study how surfactant, oil, and water interact and how porous media affect these interactions. Molecular dynamics (MD) simulation provides an opportunity to gain insights into this challenge. MD simulation can be used to study interactions between surfactant, oil, and water statically and dynamically in porous media. This paper presents a comprehensive review of interactions between surfactants and fractions of oil/heavy oil, including asphaltene, resin, aromatics, and saturates. Also, it explains the probable mechanisms of oil detachment from reservoir rock in the presence of surfactants. A thorough grasp of molecular interactions between surface-active agents and different fractions of oil helps us to develop successful surfactant-based oil recovery methods.

Molecular Dynamics Simulations and Density Functional Theory on Unraveling Photoresponsive Behavior of Wormlike Micelles Constructed by 12-2-12·2Br- and trans-ortho-Methoxy Cinnamate

Photoresponsive systems with controllable self-assembly morphologies and adjustable rheological properties have attracted widespread interest by researchers in the past few years. Among them, the photoresponsive systems consisting of ortho-methoxycinnamic (OMCA) and Gemini surfactants are endowed with rich self-assemblies with different states and in different scales including spherical micelles, wormlike micelles, vesicles, aqueous two-phase system (ATPS), etc. All these self-assemblies display excellent photoresponsive behavior. However, the mechanism of these photoresponsive behaviors has not been unraveled systematically so far. In this study, molecular dynamics (MD) simulations, density functional theory (DFT) calculations, transmission electron microscopy, and rheology are employed to investigate the photoresponsive behaviors of wormlike micelles caused by photoisomerization of trans-OMCA in 12-2-12·2Br^-/trans-OMCA solutions and to unravel the underlying mechanisms of these photoresponsive behaviors. The experimental results show that 12-2-12·2Br^-/trans-OMCA micelles display photoresponsiveness after UV-light irradiation, with the transformation of micellar morphologies from wormlike micelle to spherical micelles. In MD simulations, certain micelle morphologies in experiments and the specific packing between 12-2-12·2Br^-/OMCA were successfully captured. The larger three-dimensional structure and steric hindrance of cis-OMCA disturb the interior structure of micelles. The stronger hydrophilicity of cis-OMCA induces the escape of cis-OMCA from the interval of micelles to the solution. The energy results prove that trans-OMCA associates more strongly with 12-2-12·2Br^- than cis-OMCA. These causes lead to the fission and repacking of wormlike micelles.

Evidence for water ridges at oil-water interfaces: implications for ion transport

Understanding ion transport across interfaces is of fundamental importance in many processes such as liquid-liquid extraction, phase transfer catalysis, enhanced oil recovery and emulsion stabilisation. However, the factors that control ion transport across interfaces are poorly known due to a lack of knowledge of structural changes at interfaces. We studied here the effects of ionic concentration and external force on the transport of ions across the decane-water interface using classical molecular dynamics simulations. The results show that the evolution of interfacial structures during ion transfer across the interface is controlled by hydrogen bonding and ionic interactions at the interface. We also identified a new mode of ion transfer across the interface at low ionic concentrations, involving a 'water ridge', rather that the classical 'water finger'. In the water ridge mode, hydrogen bonds are not broken due to low ion levels, and the water ridge induces gradual interface deformation. Whereas, at high ionic concentrations, hydrogen bonds are broken by the strong ion electrostatic repulsion, thus inducing the formation of a water finger. We also found that the variation of the Gibbs free energy during ion transfer is directly relevant to the ionic concentration. The water ridge at low ionic concentrations, which displaces more water molecules towards the decane phase, induces less free energy variation than the water finger at high ionic concentrations.

Design, synthesis, biological evaluation and molecular dynamics studies of 4-thiazolinone derivatives as protein tyrosine phosphatase 1B (PTP1B) inhibitors

Protein tyrosine phosphatase 1B (PTP1B) is a key negative regulator of insulin signaling pathway, and more and more studies have shown that it is a potential target for the treatment of type 2 diabetes mellitus (T2DM). In this study, 17 new 4-thiazolinone derivatives were designed and synthesized as novel PTP1B inhibitors, and ADMET prediction confirmed that these compounds were to be drug-like. In vitro enzyme activity experiments were performed on these compounds, and it was found that a plurality of compounds had good inhibitory activity and high selectivity against PTP1B protein. Among them, compound 7p exhibited the best inhibitory activity with an IC₅₀ of 0.92 μM. The binding mode of compound 7p and PTP1B protein was explored, revealing the reason for its high efficiency. In addition, molecular dynamics simulations for the PTP1B^WT and PTP1B^comp#7p systems revealed the effects of compound 7p on PTP1B protein at the molecular level. In summary, the study reported for the first time that 4-thiazolinone derivatives as a novel PTP1B inhibitor had good inhibitory activity and selectivity for the treatment of T2DM, providing more options for the development of PTP1B inhibitors. AbbreviationsBBBblood-brain barrierCDC25Bcell division cycle 25 homolog BCYP2D6Cytochrome P450 2D6 bindingDCCMdynamic cross-correlation mapDSDiscovery StudioH bondhydrogen bondHIAhuman intestinal absorptionLARleukocyte antigen-related phosphataseMDmolecular dynamicsMEG-2maternal-effect germ-cell defective 2MM-PBSAmolecular mechanics Poisson Boltzmann surface area)PCAprincipal component analysisPDBProtein Data BankpNPPp-nitrophenyl phosphatePPBplasma protein bindingPTP1Bprotein tyrosine phosphotase 1BRMSDroot mean square deviationRMSFroot mean square fluctuationSHP-1src homologous phosphatase-1SHP-2src homologous phosphatase-2SPCsingle-point chargeTCPTPT cell protein tyrosine phosphataseT2DMType 2 diabetes mellitusVDWvan der WaalsCommunicated by Ramaswamy H. Sarma.