1
|
Salarnia M, Toghraie D, Fazilati MA, Mehmandoust B, Pirmoradian M. The effects of different nanoparticles on physical and thermal properties of water in a copper oscillating heat pipe via molecular dynamics simulation. J Taiwan Inst Chem Eng 2023. [DOI: 10.1016/j.jtice.2023.104721] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
|
2
|
Surblys D, Müller-Plathe F, Ohara T. Computing the Work of Solid-Liquid Adhesion in Systems with Damped Coulomb Interactions via Molecular Dynamics: Approaches and Insights. J Phys Chem A 2022; 126:5506-5516. [PMID: 35929812 PMCID: PMC9393893 DOI: 10.1021/acs.jpca.2c03934] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
Recently, the dry-surface method [Langmuir2016, 31, 8335−8345] has been developed
to compute the work of adhesion of solid–liquid and other interfaces
using molecular dynamics via thermodynamic integration. Unfortunately,
when long-range Coulombic interactions are present in the interface,
a special treatment is required, such as solving additional Poisson
equations, which is usually not implemented in generic molecular dynamics
software, or as fixing some groups of atoms in place, which is undesirable
most of the time. In this work, we replace the long-range Coulombic
interactions with damped Coulomb interactions, and explore several
thermal integration paths. We demonstrate that regardless of the integration
path, the same work of adhesion values are obtained as long as the
path is reversible, but the numerical efficiency differs vastly. Simple
scaling of the interactions is most efficient, requiring as little
as 8 sampling points, followed by changing the Coulomb damping parameter,
while modifying the Coulomb interaction cutoff length performs worst.
We also demonstrate that switching long-range Coulombic interactions
to damped ones results in a higher work of adhesion by about 10 mJ/m2 because of slightly different liquid molecule orientation
at the solid–liquid interface, and this value is mostly unchanged
for surfaces with substantially different Coulombic interactions at
the solid–liquid interface. Finally, even though it is possible
to split the work of adhesion into van der Waals and Coulomb components,
it is known that the specific per-component values are highly dependent
on the integration path. We obtain an extreme case, which demonstrates
that caution should be taken even when restricting to qualitative
comparison.
Collapse
Affiliation(s)
- Donatas Surblys
- Institute of Fluid Science, Tohoku University, 2-1-1 Katahira, Aoba-ku, Sendai, 980-8577, Japan
| | - Florian Müller-Plathe
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, D-64287, Germany
| | - Taku Ohara
- Institute of Fluid Science, Tohoku University, 2-1-1 Katahira, Aoba-ku, Sendai, 980-8577, Japan
| |
Collapse
|
3
|
Heterogeneous nucleation of argon vapor on the nanostructure surface with molecular dynamics simulation. J Mol Graph Model 2020; 100:107674. [PMID: 32750651 DOI: 10.1016/j.jmgm.2020.107674] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2020] [Revised: 06/23/2020] [Accepted: 06/26/2020] [Indexed: 11/21/2022]
Abstract
The study of vapor condensation on surface has important engineering significance. The condensation nucleation process of vapor on substrates with different pillar structures is studied through molecular dynamic simulation. The condensation droplets on the pillar with various heights and solid fractions exhibit Wenzel state, Cassie state and the transformation from Wenzel state to Cassie state. The results show that the condensation efficiencies are correlated with the state of droplet and it is explained from the perspective of heat transfer. For the Wenzel state, the droplet fills the gap of pillar and the form of the heat conduction change with the growth of cluster. In the initial of condensation droplets, the heat conduction is similar with various heights of pillar. As the condensation droplets grow, the efficiency of heat conduction enhances with the increasing of height of pillar. For the Cassie state, the form of heat conduction is perpetual during the condensation process with the thermal resistance of the droplet dominated due to the droplets suspended on the pillar. The efficiency of heat transfer is insensitive to the height of pillar. The form of heat conduction for the transformation state transforms from Wenzel state to Cassie state leading to the reduction of condensation rate. The droplet formed in the Wenzel state has the higher transfer efficiency than the Cassie one.
Collapse
|
4
|
Yang P, Zeng Q, Dong K, Zhu H, Yu A. A New Interaction Force Model of Gold Nanorods Derived by Molecular Dynamics Simulation. NANOMATERIALS 2020; 10:nano10071293. [PMID: 32630171 PMCID: PMC7408490 DOI: 10.3390/nano10071293] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/10/2020] [Revised: 06/26/2020] [Accepted: 06/29/2020] [Indexed: 11/18/2022]
Abstract
Interactions between nanoparticles is one of the key factors governing their assembly for ordered structures. Understanding such interactions between non-spherical nanoparticles and developing a quantitative force model are critical to achieving the ordered structures for various applications. In the present study, the non-contact interactions of two identical gold nanorods (AuNRs) with different aspect ratios have been studied by molecular dynamics simulation. A new interaction potential and force model for two nanorods approaching side-by-side has been proposed as a function of particle surface separation and their relative orientation. In addition, the interaction potentials of two nanorods approaching in other typical orientation configurations (i.e., crossed, head-to-head and head-to-side) have also been investigated.
Collapse
Affiliation(s)
- Pan Yang
- School of Engineering, Western Sydney University, Locked Bag 1797, Penrith, NSW 2751, Australia; (P.Y.); (K.D.); (H.Z)
| | - Qinghua Zeng
- School of Engineering, Western Sydney University, Locked Bag 1797, Penrith, NSW 2751, Australia; (P.Y.); (K.D.); (H.Z)
- Centre for Infrastructure Engineering, Western Sydney University, Locked Bag 1797, Penrith, NSW 2751, Australia
- Correspondence:
| | - Kejun Dong
- School of Engineering, Western Sydney University, Locked Bag 1797, Penrith, NSW 2751, Australia; (P.Y.); (K.D.); (H.Z)
- Centre for Infrastructure Engineering, Western Sydney University, Locked Bag 1797, Penrith, NSW 2751, Australia
| | - Haiping Zhu
- School of Engineering, Western Sydney University, Locked Bag 1797, Penrith, NSW 2751, Australia; (P.Y.); (K.D.); (H.Z)
| | - Aibing Yu
- Department of Chemical Engineering, Monash University, Clayton, VIC 3800, Australia;
| |
Collapse
|
5
|
Chapman J, Ramprasad R. Predicting the dynamic behavior of the mechanical properties of platinum with machine learning. J Chem Phys 2020; 152:224709. [DOI: 10.1063/5.0008955] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- James Chapman
- Department of Materials Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Rampi Ramprasad
- Department of Materials Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| |
Collapse
|
6
|
Montazeri K, Hao S, Abdolhosseini Qomi MJ, Won Y. Molecular Dynamics Investigation of Liquid and Vapor Interactions Near an Evaporating Interface: A Theoretical Genetics Perspective. ADVANCED THEORY AND SIMULATIONS 2020. [DOI: 10.1002/adts.202000017] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Kimia Montazeri
- Department of Mechanical and Aerospace EngineeringUniversity of California Irvine 5200 Engineering Hall Irvine Irvine CA 92697‐2700 USA
| | - Shuai Hao
- Department of Mechanical and Aerospace EngineeringUniversity of California Irvine 5200 Engineering Hall Irvine Irvine CA 92697‐2700 USA
| | - Mohammad Javad Abdolhosseini Qomi
- Department of Civil and Environmental EngineeringUniversity of California Irvine 5200 Engineering Hall Irvine Irvine CA 92697‐2700 USA
| | - Yoonjin Won
- Department of Mechanical and Aerospace EngineeringUniversity of California Irvine 5200 Engineering Hall Irvine Irvine CA 92697‐2700 USA
| |
Collapse
|
7
|
|
8
|
Surblys D, Matsubara H, Kikugawa G, Ohara T. Application of atomic stress to compute heat flux via molecular dynamics for systems with many-body interactions. Phys Rev E 2019; 99:051301. [PMID: 31212446 DOI: 10.1103/physreve.99.051301] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2019] [Indexed: 06/09/2023]
Abstract
Although the computation of heat flux and thermal conductivity either via Fourier's law or the Green-Kubo relation has become a common task in molecular dynamics simulation, contributions of three-body and larger many-body interactions have always proved problematic to compute. In recent years, due to the success when applying to pressure tensor computation, atomic stress approximation has been widely used to calculate heat flux, where the lammps molecular dynamics package is the most prominent propagator. We demonstrated that the atomic stress approximation, while adequate for obtaining pressure, produces erroneous results in the case of heat flux when applied to systems with many-body interactions, such as angle, torsion, or improper potentials. This also produces incorrect thermal conductivity values. To remedy this deficiency, by starting from a strict formulation of heat flux with many-body interactions, we reworked the atomic stress definition which resulted in only a simple modification. We modified the lammps package accordingly to demonstrate that the new atomic stress approximation produces excellent results close to that of a rigid formulation.
Collapse
Affiliation(s)
- Donatas Surblys
- Institute of Fluid Science, Tohoku University, 2-1-1 Katahira, Aoba-ku, Sendai 980-8577, Japan
| | - Hiroki Matsubara
- Institute of Fluid Science, Tohoku University, 2-1-1 Katahira, Aoba-ku, Sendai 980-8577, Japan
| | - Gota Kikugawa
- Institute of Fluid Science, Tohoku University, 2-1-1 Katahira, Aoba-ku, Sendai 980-8577, Japan
| | - Taku Ohara
- Institute of Fluid Science, Tohoku University, 2-1-1 Katahira, Aoba-ku, Sendai 980-8577, Japan
| |
Collapse
|
9
|
Liu L, Stetsyuk V, Kubiak KJ, Yap YF, Goharzadeh A, Chai JC. Nanoparticles for convective heat transfer enhancement: heat transfer coefficient and the effects of particle size and zeta potential. CHEM ENG COMMUN 2018. [DOI: 10.1080/00986445.2018.1525364] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2022]
Affiliation(s)
- Lande Liu
- School of Applied Sciences, University of Huddersfield, Huddersfield, UK
| | - Viacheslav Stetsyuk
- School of Computing and Engineering, University of Huddersfield, Huddersfield, UK
| | - Krzystof J. Kubiak
- School of Computing and Engineering, University of Huddersfield, Huddersfield, UK
| | - Yit Fatt Yap
- Department of Mechanical Engineering, Petroleum Institute, Khalifa University of Science and Technology, Abu Dhabi, United Arab Emirates
| | - Afshin Goharzadeh
- Department of Mechanical Engineering, Petroleum Institute, Khalifa University of Science and Technology, Abu Dhabi, United Arab Emirates
| | - John C. Chai
- School of Computing and Engineering, University of Huddersfield, Huddersfield, UK
| |
Collapse
|
10
|
Ge S, Gu Y, Chen M. A molecular dynamics simulation on the convective heat transfer in nanochannels. Mol Phys 2014. [DOI: 10.1080/00268976.2014.970593] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
11
|
Wang Y, Keblinski P, Chen Z. Viscosity calculation of a nanoparticle suspension confined in nanochannels. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2012; 86:036313. [PMID: 23031019 DOI: 10.1103/physreve.86.036313] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/29/2012] [Revised: 07/23/2012] [Indexed: 06/01/2023]
Abstract
The kinetic properties of the pressure-driven Poiseuille flow in nanochannels with and without nanoparticles were studied with a nonequilibrium molecular dynamics simulation. To allow the fluid to dissipate heat, the boundary was kept at a constant temperature. Pure fluid simulations were taken as references and also used to study the fluid-wall interfacial interaction effects. The viscosity profiles of the fluid were calculated on the basis of velocity profiles and known applied shear stress. We present the relationship between the viscosity increase and particle loading. The role of channel wall-fluid wetting properties on the flow and viscosity was also investigated.
Collapse
Affiliation(s)
- Yan Wang
- School of Energy and Environment, Southeast University, Nanjing 210096, China.
| | | | | |
Collapse
|
12
|
Spijker P, Markvoort AJ, Nedea SV, Hilbers PAJ. Computation of accommodation coefficients and the use of velocity correlation profiles in molecular dynamics simulations. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 81:011203. [PMID: 20365360 DOI: 10.1103/physreve.81.011203] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/30/2009] [Indexed: 05/29/2023]
Abstract
For understanding the behavior of a gas close to a channel wall it is important to model the gas-wall interactions as detailed as possible. When using molecular dynamics simulations these interactions can be modeled explicitly, but the computations are time consuming. Replacing the explicit wall with a wall model reduces the computational time but the same characteristics should still remain. Elaborate wall models, such as the Maxwell-Yamamoto model or the Cercignani-Lampis model need a phenomenological parameter (the accommodation coefficient) for the description of the gas-wall interaction as an input. Therefore, computing these accommodation coefficients in a reliable way is very important. In this paper, two systems (platinum walls with either argon or xenon gas confined between them) are investigated and are used for comparison of the accommodation coefficients for the wall models and the explicit molecular dynamics simulations. Velocity correlations between incoming and outgoing particles colliding with the wall have been used to compare explicit simulations and wall models even further. Furthermore, based on these velocity correlations, a method to compute the accommodation coefficients is presented, and these newly computed accommodation coefficients are used to show improved correlation behavior for the wall models.
Collapse
Affiliation(s)
- Peter Spijker
- Department of Biomedical, Eindhoven University of Technology, Eindhoven, The Netherlands.
| | | | | | | |
Collapse
|
13
|
Maroo SC, Chung J. Molecular dynamic simulation of platinum heater and associated nano-scale liquid argon film evaporation and colloidal adsorption characteristics. J Colloid Interface Sci 2008; 328:134-46. [PMID: 18829047 DOI: 10.1016/j.jcis.2008.09.018] [Citation(s) in RCA: 69] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2008] [Revised: 08/16/2008] [Accepted: 09/07/2008] [Indexed: 11/26/2022]
|
14
|
McGaughey A, Kaviany M. Phonon Transport in Molecular Dynamics Simulations: Formulation and Thermal Conductivity Prediction. ADVANCES IN HEAT TRANSFER 2006. [DOI: 10.1016/s0065-2717(06)39002-8] [Citation(s) in RCA: 152] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/03/2022]
|
15
|
Markvoort AJ, Hilbers PAJ, Nedea SV. Molecular dynamics study of the influence of wall-gas interactions on heat flow in nanochannels. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2005; 71:066702. [PMID: 16089906 DOI: 10.1103/physreve.71.066702] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/24/2004] [Revised: 03/09/2005] [Indexed: 05/03/2023]
Abstract
Especially at the nanometer scale interfaces play an important role. The effect of the wettability on the solid-liquid interface has already been studied with molecular dynamics. In this paper we study the dependence of wetting on the solid-gas interface for different density gases and investigate the influence of wetting on the heat transport properties over such an interface using molecular dynamics. Subsequently we show how the flow profile of a gas flowing along a surface also depends on this wettability. These simulations show that wettability increases the conductivity of a solid to a stationary gas and decreases the flow velocity near the interface for a gas flow. These two effects influence the cooling of a solid achieved by a cold gas flowing along its surface in opposite ways. However, we show that a higher wettability has a positive net effect on the cooling, explaining experimental results that showed an increased heat cooling effect of hydrophilic over hydrophobic microchannels.
Collapse
Affiliation(s)
- A J Markvoort
- Department of Biomedical Engineering, Technische Universiteit Eindhoven, Postbus 513, 5600 MB Eindhoven, The Netherlands.
| | | | | |
Collapse
|