Evaluations of Molecular Dynamics Methods for Thermodiffusion in Binary Mixtures

Author(s):  
Seyedeh H. Mozaffari ◽  
Seshasai Srinivasan ◽  
M. Ziad Saghir

The objective of this paper is to investigate the behavior of two well-known boundary-driven molecular dynamics (MD) approaches, namely, reverse nonequilibrium molecular dynamics (RNEMD) and heat exchange algorithm (HEX), as well as introducing a modified HEX model (MHEX) that is more accurate and computationally efficient to simulate the mass and heat transfer mechanism. For this investigation, the following binary mixtures were considered: one equimolar mixture of argon (Ar) and krypton (Kr), one nonequimolar liquid mixture of hexane (nC6) and decane (nC10), and three nonequimolar mixtures of pentane (nC5) and decane. In estimating the Thermodiffusion factor in these mixtures using the three methods, it was found that consistent with the findings in the literature, RNEMD predictions have the largest error with respect to the experimental data. Whereas, the MHEX method proposed in this work is the most accurate, marginally outperforming the HEX method. Most importantly, the computational efficiency of MHEX method is the highest, about 7% faster than the HEX method. This makes it more suitable for integration with multiscale computational models to simulate thermodiffusion in a large system such as an oil reservoir.

2021 ◽  
Author(s):  
Seyedeh Hoda Mozaffari

Thermodiffusion phenomenon in fluid mixtures has been investigated by several scientists in theoretical as well as experimental fields for decades. Nevertheless, due to shortcomings of both methods, interest in searching for alternative approaches to shed some light on molecular scale of the phenomenon has spurred. The objective of this thesis is to develop an accurate molecular dynamics (MD) algorithm that can predict thermodiffusive separation in binary and ternary fluid mixtures. More importantly, the proposed algorithm should be computationally efficient in order to be suitable for integration into multi-scale computational models to simulate thermodiffusion in a large system such as an oil reservoir. In developing such an effective and efficient computational tool, this thesis introduces a modified heat exchange algorithms, wherein, a new mechanism is introduced to rescale velocities which curbs the energy loss in the system and at the same time minimizes the computational time. The performance of the new algorithm in studying Soret effect for binary and ternary mixtures has been compared with other non-equilibrium molecular dynamics (NEMD) models including regular heat exchange algorithm (HEX) and reverse non-equilibrium molecular dynamics (RNEMD). Different types of binary mixtures were studied including one equimolar mixture of argon (Ar)-krypton (Kr) above its triple point, non-equimolar normal alkane mixtures of hexane (nC6)-decane (nC10) as well as hexane (nC6)-dodecane (nC12) for six compositions, three non-equimolar mixtures of pentane (nC5) decane (nC10) at atmospheric temperature and pressure. Additionally, the new algorithm was validated for different ternary mixtures including ternary normal alkanes methane (nC1)-butane (nC4)- dodecane (nC12) for three compositions, and one composition of different types of alkane mixture of 1,2,3,4-tetrahydronaphthalene (THN)-dodecane (nC12)-sobutylbenzene (IBB). The new algorithm demonstrates a significant improvement in reducing the energy loss by nearly 32%. Additionally, the new algorithm is about 7-9% more computationally efficient than the regular HEX for medium and large systems. In terms of direction of thermodiffusive segregations in binary mixtures, in agreement with the experimental data, the new algorithm shows that the heavier component moves towards the cold region whereas the lighter component accumulates near the hot zone. Additionally, the strength of segregation process diminishes as the concentration of heavy component in the mixture increases. The new algorithm improved the prediction of thermodiffusion factor in binary mixtures by 24% in binary mixtures. With respect to the ternary mixtures, similarly to binary mixtures the heaviest and lightest component in the mixture move towards, cold and hot zones, respectively. While the intermediate component shows the least tendency to segregate. In terms of the strength of Soret effect, the new algorithm is about 17% more accurate than the regular HEX algorithm with respect to experimental data.


2021 ◽  
Author(s):  
Seyedeh Hoda Mozaffari

Thermodiffusion phenomenon in fluid mixtures has been investigated by several scientists in theoretical as well as experimental fields for decades. Nevertheless, due to shortcomings of both methods, interest in searching for alternative approaches to shed some light on molecular scale of the phenomenon has spurred. The objective of this thesis is to develop an accurate molecular dynamics (MD) algorithm that can predict thermodiffusive separation in binary and ternary fluid mixtures. More importantly, the proposed algorithm should be computationally efficient in order to be suitable for integration into multi-scale computational models to simulate thermodiffusion in a large system such as an oil reservoir. In developing such an effective and efficient computational tool, this thesis introduces a modified heat exchange algorithms, wherein, a new mechanism is introduced to rescale velocities which curbs the energy loss in the system and at the same time minimizes the computational time. The performance of the new algorithm in studying Soret effect for binary and ternary mixtures has been compared with other non-equilibrium molecular dynamics (NEMD) models including regular heat exchange algorithm (HEX) and reverse non-equilibrium molecular dynamics (RNEMD). Different types of binary mixtures were studied including one equimolar mixture of argon (Ar)-krypton (Kr) above its triple point, non-equimolar normal alkane mixtures of hexane (nC6)-decane (nC10) as well as hexane (nC6)-dodecane (nC12) for six compositions, three non-equimolar mixtures of pentane (nC5) decane (nC10) at atmospheric temperature and pressure. Additionally, the new algorithm was validated for different ternary mixtures including ternary normal alkanes methane (nC1)-butane (nC4)- dodecane (nC12) for three compositions, and one composition of different types of alkane mixture of 1,2,3,4-tetrahydronaphthalene (THN)-dodecane (nC12)-sobutylbenzene (IBB). The new algorithm demonstrates a significant improvement in reducing the energy loss by nearly 32%. Additionally, the new algorithm is about 7-9% more computationally efficient than the regular HEX for medium and large systems. In terms of direction of thermodiffusive segregations in binary mixtures, in agreement with the experimental data, the new algorithm shows that the heavier component moves towards the cold region whereas the lighter component accumulates near the hot zone. Additionally, the strength of segregation process diminishes as the concentration of heavy component in the mixture increases. The new algorithm improved the prediction of thermodiffusion factor in binary mixtures by 24% in binary mixtures. With respect to the ternary mixtures, similarly to binary mixtures the heaviest and lightest component in the mixture move towards, cold and hot zones, respectively. While the intermediate component shows the least tendency to segregate. In terms of the strength of Soret effect, the new algorithm is about 17% more accurate than the regular HEX algorithm with respect to experimental data.


Sign in / Sign up

Export Citation Format

Share Document