0
Research Papers

Molecular Dynamics Study on Explosive Boiling of Thin Liquid Argon Film on Nanostructured Surface Under Different Wetting Conditions

[+] Author and Article Information
Sheikh Mohammad Shavik

Department of Mechanical Engineering,
Bangladesh University of Engineering and Technology,
Dhaka 1000, Bangladesh
e-mail: shavik@me.buet.ac.bd

Mohammad Nasim Hasan, A. K. M. Monjur Morshed

Department of Mechanical Engineering,
Bangladesh University of Engineering and Technology,
Dhaka 1000, Bangladesh

1Corresponding author.

Contributed by the Electronic and Photonic Packaging Division of ASME for publication in the JOURNAL OF ELECTRONIC PACKAGING. Manuscript received September 20, 2015; final manuscript received November 25, 2015; published online March 10, 2016. Assoc. Editor: Toru Ikeda.

J. Electron. Packag 138(1), 010904 (Mar 10, 2016) (8 pages) Paper No: EP-15-1086; doi: 10.1115/1.4032463 History: Received September 20, 2015; Revised November 25, 2015

Molecular dynamics (MDs) simulations have been performed to investigate the boiling phenomena of thin liquid film adsorbed on a nanostructured solid surface with particular emphasis on the effect of wetting condition of the solid surface. The molecular system consists of liquid and vapor argon and solid platinum wall. The nanostructures which reside on top of the solid wall have shape of rectangular block. The solid–liquid interfacial wettability, in other words whether the solid surface is hydrophilic or hydrophobic, has been altered for different cases to examine its effect on boiling phenomena. The initial configuration of the simulation domain comprises a three-phase system (solid platinum, liquid argon, and vapor argon), which was equilibrated at 90 K. After equilibrium period, the wall temperature was suddenly increased from 90 K to 250 K which is far above the critical point of argon and this initiates rapid or explosive boiling. The spatial and temporal variation of temperature and density as well as the variation of system pressure with respect to time were closely monitored for each case. The heat flux normal to the solid surface was also calculated to illustrate the effectiveness of heat transfer for different cases of wetting conditions of solid surface. The results show that the wetting condition of surface has significant effect on explosive boiling of the thin liquid film. The surface with higher wettability (hydrophilic) provides more favorable conditions for boiling than the low-wetting surface (hydrophobic), and therefore, the liquid argon responds quickly and shifts from liquid to vapor phase faster in the case of hydrophilic surface. The heat transfer rate is also much higher in the case of hydrophilic surface.

Copyright © 2016 by ASME
Your Session has timed out. Please sign back in to continue.

References

Park, H. , Kim, D. , Grigoropoulos, C. , and Tam, A. , 1996, “ Pressure Generation and Measurement in the Rapid Vaporization of Water on a Pulsed-Laser-Heated Surface,” J. Appl. Phys., 80(7), pp. 4072–4081. [CrossRef]
Xu, X. , and Song, K. , 1999, “ Interface Kinetics During Pulsed Laser Ablation,” Appl. Phys. A, 69(7), pp. S869–S873. [CrossRef]
Dou, Y. , Zhigilei, L. V. , Winograd, N. , and Garrison, B. J. , 2001, “ Explosive Boiling of Water Films Adjacent to Heated Surfaces: A Microscopic Description,” J. Phys. Chem. A, 105(12), pp. 2748–2755. [CrossRef]
Gu, X. , and Urbassek, H. M. , 2005, “ Atomic Dynamics of Explosive Boiling of Liquid-Argon Films,” Appl. Phys. B, 81(5), pp. 675–679. [CrossRef]
Nagayama, G. , Tsuruta, T. , and Cheng, P. , 2006, “ Molecular Dynamics Simulation on Bubble Formation in a Nanochannel,” Int. J. Heat Mass Transfer, 49(23–24), pp. 4437–4443. [CrossRef]
Yamamoto, T. , and Matsumoto, M. , 2012, “ Initial Stage of Nucleate Boiling: Molecular Dynamics Investigation,” J. Therm. Sci. Technol. Jpn., 7(1), pp. 334–349. [CrossRef]
Nagayama, G. , and Cheng, P. , 2004, “ Effects of Interface Wettability on Microscale Flow by Molecular Dynamics Simulation,” Int. J. Heat Mass Transfer, 47(3), pp. 501–513. [CrossRef]
Yu, J. , and Wang, H. , 2012, “ A Molecular Dynamics Investigation on Evaporation of Thin Liquid Films,” Int. J. Heat Mass Transfer, 55(4), pp. 1218–1225. [CrossRef]
Mao, Y. , and Zhang, Y. , 2014, “ Molecular Dynamics Simulation on Rapid Boiling of Water on a Hot Copper Plate,” Appl. Therm. Eng., 62(2), pp. 607–612. [CrossRef]
Maroo, S. C. , and Chung, J. N. , 2008, “ Molecular Dynamic Simulation of Platinum Heater and Associated Nano-Scale Liquid Argon Film Evaporation and Colloidal Adsorption Characteristics,” J. Colloid Interface Sci., 328(1), pp. 134–146. [CrossRef] [PubMed]
Fu, T. , Mao, Y. , Tang, Y. , Zhang, Y. , and Yuan, W. , 2015, “ Effect of Nanostructure on Rapid Boiling of Water on a Hot Copper Plate: A Molecular Dynamics Study,” Heat Mass Transfer, pp. 1–10.
Morshed, A. K. M. M. , Paul, T. C. , and Khan, J. A. , 2011, “ Effect of Nanostructures on Evaporation and Explosive Boiling of Thin Liquid Films: A Molecular Dynamics Study,” Appl. Phys. A, 105(2), pp. 445–451. [CrossRef]
Seyf, H. R. , and Zhang, Y. , 2012, “ Molecular Dynamics Simulation of Normal and Explosive Boiling on Nanostructured Surface,” ASME J. Heat Transfer, 213(135), p. 121503.
Seyf, H. R. , and Zhang, Y. , 2013, “ Effect of Nanotextured Array of Conical Features on Explosive Boiling Over a Flat Substrate: A Nonequilibrium Molecular Dynamics Study,” Int. J. Heat Mass Transfer, 66, pp. 613–624. [CrossRef]
Wang, W. , Zhang, H. , Tian, C. , and Meng, X. , 2015, “ Numerical Experiments on Evaporation and Explosive Boiling of Ultra-Thin Liquid Argon Film on Aluminum Nanostructure Substrate,” Nanoscale Res. Lett., 10(158), pp. 1–14. [PubMed]
Kimura, T. , and Maruyama, S. , 2002, “ Molecular Dynamics Simulation of Heterogeneous Nucleation of a Liquid Droplet on a Solid Surface,” Microscale Thermophys. Eng., 6(1), pp. 3–13. [CrossRef]
Maruyama, S. , and Kimura, T. , 1999, “ A Molecular Dynamics Simulation of a Bubble Nucleation on Solid Surface,” Trans. Japan Soc. Mech. Eng. B, 65(638), pp. 3461–3467. [CrossRef]
Diaz, R. , and Guo, Z. , 2015, “ A Molecular Dynamics Study of Passive Heat Flux Enhancement in Boiling Liquid Argon,” CHT-15, ICHMT International Symposium on Advances in Computational Heat Transfer.
Hens, A. , Agarwal, R. , and Biswas, G. , 2014, “ Nanoscale Study of Boiling and Evaporation in a Liquid Ar Film on a Pt Heater Using Molecular Dynamics Simulation,” Int. J. Heat Mass Transfer, 71, pp. 303–312. [CrossRef]
Gambill, W. R. , and Lienhard, J. H. , 1989, “ An Upper Bound for the Critical Boiling Heat Flux,” ASME J. Heat Transfer, 111(3), pp. 815–818. [CrossRef]
Lennard-Jones, J. E. , and Devonshire, A. F. , 1937, “ Critical Phenomena in Gases. I,” Philos. Trans. R. Soc., A, 163(912), pp. 53–70.
Kofke, D. A. , 1993, “ Direct Evaluation of Phase Coexistence by Molecular Simulation Via Integration Along the Saturation Line,” J. Chem. Phys., 98(5), pp. 4149–4162. [CrossRef]
Plimpton, S. , 1995, “ Fast Parallel Algorithms for Short-Range Molecular Dynamics,” J. Comput. Phys., 117(1), pp. 1–19. [CrossRef]
Humphrey, W. , Dalke, A. , and Schulten, K. , 1996, “ VMD: Visual Molecular Dynamics,” J. Mol. Graphics Modell., 14(1), pp. 33–38. [CrossRef]
“LAMMPS User's Manual,” http://lammps.sandia.gov

Figures

Grahic Jump Location
Fig. 1

Initial configuration of the simulation

Grahic Jump Location
Fig. 2

Solid wall configurations: (a) flat surface, (b) nanostructured surface 1, (c) nanostructured surface 2, and (d) nanostructured surface 3

Grahic Jump Location
Fig. 3

Number density profile of argon on a flat hydrophilic (εAr–PtAr–Ar = 2) surface during equilibrium period

Grahic Jump Location
Fig. 4

Temporal variation of system energy

Grahic Jump Location
Fig. 5

Temperature history of argon (Ar) and solid wall for hydrophilic surface case

Grahic Jump Location
Fig. 6

Temperature history of argon (Ar) and solid wall for hydrophobic surface case

Grahic Jump Location
Fig. 7

Pressure history of the simulation domain for hydrophilic surface case

Grahic Jump Location
Fig. 8

Pressure history of the simulation domain for hydrophobic surface case

Grahic Jump Location
Fig. 9

Snapshots from the simulation domain for hydrophilic surface case: (a) flat surface, (b) surface 1, (c) surface 2, and (d) surface 3

Grahic Jump Location
Fig. 10

Snapshots from the simulation domain for hydrophobic surface case: (a) flat surface, (b) surface 1, (c) surface 2, and (d) surface 3

Grahic Jump Location
Fig. 11

Spatial temperature distribution of argon in the case of flat hydrophilic surface at different times

Grahic Jump Location
Fig. 12

Number density profile of argon in the case of flat hydrophilic surface at different times

Grahic Jump Location
Fig. 13

Spatial temperature distribution of argon at 3 ns for hydrophilic surface case

Grahic Jump Location
Fig. 14

Number density profile of argon at 3 ns for hydrophilic surface case

Grahic Jump Location
Fig. 15

Spatial temperature distribution of argon at 3 ns for hydrophobic surface case

Grahic Jump Location
Fig. 16

Number density profile of argon at 3 ns for hydrophobic surface case

Grahic Jump Location
Fig. 17

Heat flux normal to solid wall for hydrophilic surface case

Grahic Jump Location
Fig. 18

Heat flux normal to solid wall for hydrophobic surface case

Grahic Jump Location
Fig. 19

Net evaporation number for hydrophilic surface case

Grahic Jump Location
Fig. 20

Net evaporation number for hydrophobic surface case

Grahic Jump Location
Fig. 21

Nonevaporating layer in the case of hydrophilic surface: (a) surface 2 and (b) surface 3

Grahic Jump Location
Fig. 22

Nonevaporating layer in the case of hydrophobic surface: (a) surface 2 and (b) surface 3

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In