Abstract
Various characteristics of a thin liquid film in its vapor-phase are investigated using the molecular dynamics technique. Local distributions of the temperature, density, normal and tangential pressure components, and stress are calculated for various film thicknesses and temperature levels. Distributions of local stresses change considerably with respect to film thicknesses, and interfacial regions on both sides of the film start to overlap with each other as the film becomes thinner. Integration of the local stresses, i.e., the surface tension, however, does not vary much regardless of the interfacial overlap. The minimum thickness of a liquid film before rupturing is estimated with respect to the calculation domain sizes and is compared with a simple theoretical relation.
Similar content being viewed by others
Abbreviations
- di :
-
Interface thickness
- F:
-
Force
- kb :
-
Boltzmann constant
- L:
-
Simulation domain size
- m:
-
Molecular mass
- n:
-
Number density
- N:
-
Total number of molecules
- ?:
-
Pressure
- r:
-
Molecular position
- rc :
-
Cutoff radius
- rij :
-
Inter-distance between moleculesi andJ
- t:
-
Time
- T:
-
Temperature
- υ:
-
Velocity
- z:
-
Direction normal to the film
- Zo :
-
Parameter for a fitting function
- γ:
-
Surface tension
- ε:
-
Energy parameter
- σ:
-
Length parameter
- Φ:
-
Potential function
- *:
-
Dimensionless
- f:
-
Film
- I:
-
Intermolecular
- K:
-
Kinetic
- n:
-
Normal
- sl:
-
Slab
- t:
-
Tangential
- x,y,z:
-
Directions in rectangular coordinate
References
Abramson, A. R. and Tien, C. L., 1999, “Recent Developments in Microscale Thermophysical Engineering,”Microscale Thermophysical Engineering, Vol. 3, pp. 229–244.
Allen, M. P. and Tildesley, D. L., 1987,Computer Simulation of Liquids, Oxford University Press, New York.
Chapela, G. A., Saville, G., Thompson, S. M. and Rowlinson, J. S., 1977, “Computer Simulation of a Gas-Liquid Surface,” J. Chem. Soc. Faraday Trans. II, Vol. 73, p. 1133.
Haile, J. M., 1992,Molecular Dynamics Simulation, John Wiley & Sons, pp. 260-267.
Hwang, C. C., Hsieh, J. Y., Chang, K. H. and Liao, J. J., 1988, “A Study of Rupture Process of Thin Liquid Films by a Molecular Dynamics Simulation,” Physica A 256, pp. 333–341.
Israelachvili, J., 1992,Intermoleeular and Surface Forces, Academic, London, pp. 176–212.
Majumdar, A. and Mezic, I., 1998, “Stability Regimes of Thin Liquid Films,”Microscale Thermophysical Engineering, Vol. 2, pp. 203–213.
Nijmeijer, M. J. P.. Bakker, A. F., Bruin, C., and Sikkenk, J. H., 1988, “A Molecular Dynamics Simulation of the Lennard-Jones Liquid-Vapor Interface,”Journal of Chemical Physics. Vol. 89, pp. 3789–3792.
Rowlinson, J. S. and Widom, B., 1982,Molecular Theory of Capillary, Clarendon. Oxford, pp. 69–128.
Swope, W. C., Anderson, H. C., Berens, P. H. and Wilson, K. R., 1982, “A Computer Simulation Method for the Calculation of Equilibrium Constants for the Formation of Physical Clusters of Molecules: Application to Small Water Clusters,”Journal of Chemical Physics, Vol. 76, pp. 637–649.
Vargaftik, N. B., 1975,Table on the Thermophysical Properties of Liquids and Gases, Hemisphere, Washington, D. C. P. 543.
Weng, J. G., Park, S. H. and Tien, C. L., 2000a, “Interfacial Ambiguities in Microdroplets and Microbubbles,”Microscale Thermophysical Engineering, Vol. 4, pp. 83–87.
Weng, J. G., Park, S. H., Lukes, J. R. and Tien, C. L., 2000b, “Molecular Dynamics Investigation of Thickness Effect on Liquid Films,”Journal of Chemical Physics, Vol. 113, pp. 5917–5923.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Lee, J., Park, S., Kwon, O. et al. Characterization of thin liquid films using molecular dynamics simulation. KSME International Journal 16, 1477–1484 (2002). https://doi.org/10.1007/BF02985141
Received:
Revised:
Issue Date:
DOI: https://doi.org/10.1007/BF02985141