Research Papers: Multiphase Flows

Development of a New Contact Angle Control Algorithm for Level-Set Method

[+] Author and Article Information
Mengnan Li

Department of Nuclear Engineering,
North Carolina State University,
Raleigh, NC 27695-7909
e-mail: mli21@ncsu.edu

Kaiyue Zeng

Department of Nuclear Engineering,
North Carolina State University,
Raleigh, NC 27695-7909
e-mail: kzeng2@ncsu.edu

Louis Wonnell

Department of Mechanical and
Nuclear Engineering,
Kansas State University,
Manhattan, KS 66506
e-mail: lwonnell@ksu.edu

Igor A. Bolotnov

Department of Nuclear Engineering,
North Carolina State University,
Raleigh, NC 27695-7909
e-mail: igor_bolotnov@ncsu.edu

1Corresponding author.

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received April 19, 2018; final manuscript received November 2, 2018; published online December 24, 2018. Assoc. Editor: Riccardo Mereu.

J. Fluids Eng 141(6), 061301 (Dec 24, 2018) (16 pages) Paper No: FE-18-1281; doi: 10.1115/1.4041987 History: Received April 19, 2018; Revised November 02, 2018

A contact angle control algorithm is developed and implemented in the multiphase interface tracking flow solver—phasta. The subgrid force model is introduced to control the evolving contact angle. The contact angle force is applied when the current contact angle deviates from the desired value (or range of values) and decreases to zero when it reaches the desired value. The single bubble departure simulation and the capillary flat plates simulation are performed for verification purpose. The numerical results are compared with the analytical solution with good agreement. The mesh resolution sensitivity analysis and parametric study are conducted for both simulations. Coupled with the other existing capabilities in phasta like evaporation and condensation algorithm, the contact angle control algorithm will allow us to investigate the boiling phenomenon in various conditions with lower cost (by utilizing localized mesh refinement for bubble growth region) compared to uniformly refined structured meshes and in engineering geometries.

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


Basu, N. , Warrier, G. R. , and Dhir, V. K. , 2005, “ Wall Heat Flux Partitioning During Subcooled Flow Boiling—Part 1: Model Development,” ASME J. Heat Transfer, 127(2), pp. 131–140. [CrossRef]
Fritz, W. , 1935, “ Maximum Volume of Vapor Bubbles,” Phys. Z, 36(11), pp. 379–384.
Zuber, N. , 1959, “ Hydrodynamic Aspects of Boiling Heat Transfer (Thesis),” Ramo-Wooldridge Corp.,/University of California, Los Angeles, CA, Report No. AECU-4439.
Young, T. , 1832, “ An Essay on the Cohesion of Fluids,” Proc. R. Soc. Lond., 1, pp. 171–172. [CrossRef]
Wenzel, R. , 1936, “ Surface Roughness and Contact Angle,” Ind. Eng. Chem., 28, pp. 988–994. [CrossRef]
Cassie, A. , 1948, “ Contact Angles,” Discuss. Faraday Soc., 3, pp. 11–16. [CrossRef]
Lam, C. , Wu, R. , and Li, D. , 2002, “ Study of the Advancing and Receding Contact Angles: Liquid Sorption as a Cause of Contact Angle Hysteresis,” Adv. Colloid Interface Sci., 96(1–3), pp. 169–191. [CrossRef] [PubMed]
Ramanujapu, N. , and Dhir, V. , 1999, “ Dynamics of Contact Angle During Growth and Detachment of a Vapor Bubble at a Single Nucleation Site,” University of California, Los Angeles, CA, No. 20026815. https://www.osti.gov/biblio/20026815
Sobolev, V. , Churaev, N. , and Velarde, M. , 2000, “ Surface Tension and Dynamic Contact Angle of Water in Thin Quartz Capillaries,” J. Colloid Interface Sci., 222(1), pp. 51–54. [CrossRef] [PubMed]
Kandlikar, S. G. , Kuan, W. K. , and Mukherjee, A. , 2005, “ Experimental Study of Heat Transfer in an Evaporating Meniscus on a Moving Heated Surface,” ASME J. Heat Transfer, 127(3), pp. 244–252. [CrossRef]
Fukai, J. , Shiiba, Y. , and Yamamoto, T. , 1995, “ Wetting Effects on the Spreading of a Liquid Droplet Colliding With a Flat Surface: Experiment and Modeling,” Phys. Fluids, 7(2), pp. 236–247. [CrossRef]
Bussmann, M. , Mostaghimi, J. , and Chandra, S. , 1999, “ On a Three-Dimensional Volume Tracking Model of Droplet Impact,” Phys. Fluids, 11(6), pp. 1406–1417. [CrossRef]
Abarajith, H. , and Dhir, V. , 2002, “ A Numerical Study of the Effect of Contact Angle on the Dynamics of a Single Bubble During Pool Boiling,” ASME Paper No. IMECE2002-33876.
Son, G. , and Hur, N. , 2005, “ A Level Set Formulation for Incompressible Two-Phase Flows on Nonorthogonal Grids,” Numer. Heat Transfer, Part B: Fundam., 48(3), pp. 303–316. [CrossRef]
Mukherjee, A. , and Kandlikar, S. G. , 2007, “ Numerical Study of Single Bubbles With Dynamic Contact Angle During Nucleate Pool Boiling,” Int. J. Heat Mass Transfer, 50(1–2), pp. 127–138. [CrossRef]
Sato, Y. , and Ničeno, B. , 2012, “ A New Contact Line Treatment for a Conservative Level Set Method,” J. Comput. Phys., 231(10), pp. 3887–3895. [CrossRef]
Li, M. , and Bolotnov, I. A. , 2015, “ Interface Tracking Simulation of Boiling Phenomena—Single Bubble Verification,” ANS Annual Meeting 2015 Transactions, Washington, DC, Nov. 8–12, pp. 1515–1517. https://www.researchgate.net/publication/319777906_Interface_Tracking_Simulation_of_Boiling_Phenomena_-Single_Bubble_Verification
Whiting, C. H. , and Jansen, K. E. , 2001, “ A Stabilized Finite Element Method for the Incompressible Navier-Stokes Equations Using a Hierarchical Basis,” Int. J. Numer. Methods Fluids, 35(1), pp. 93–116. [CrossRef]
Fang, J. , Cambareri, J. J. , and Brown, C. S. , 2018, “ Direct Numerical Simulation of Reactor Two-Phase Flows Enabled by High-Performance Computing,” Nucl. Eng. Des., 330, pp. 409–419. [CrossRef]
Fang, J. , Rasquin, M. , and Bolotnov, I. A. , 2017, “ Interface Tracking Simulations of Bubbly Flows in PWR Relevant Geometries,” Nucl. Eng. Des., 312, pp. 205–213. [CrossRef]
Wilcox, D. C. , 1998, Turbulence Modeling for CFD, DCW Industries, La Cañada, CA.
Brackbill, J. U. , Kothe, D. B. , and Zemach, C. , 1992, “ A Continuum Method for Modeling Surface Tension,” J. Comput. Phys., 100(2), pp. 335–354. [CrossRef]
Sussman, M. , Almgren, A. S. , and Bell, J. B. , 1999, “ An Adaptive Level Set Approach for Incompressible Two-Phase Flows,” J. Comput. Phys., 148(1), pp. 81–124. [CrossRef]
Sussman, M. , and Fatemi, E. , 1999, “ An Efficient, Interface-Preserving Level Set Re-Distancing Algorithm and Its Application to Interfacial Incompressible Fluid Flow,” SIAM J. Sci. Comput., 20(4), pp. 1165–1191. [CrossRef]
Sussman, M. , Fatemi, E. , and Smereka, P. , 1998, “ An Improved Level Set Method for Incompressible Two-Phase Flows,” Comput. Fluids, 27(5), pp. 663–680. [CrossRef]
Sussman, M. , and Smereka, P. , 1997, “ Axisymmetric Free Boundary Problems,” J. Fluid Mech., 341, pp. 269–294. [CrossRef]
Sethian, A. J. , 1999, Level Set Methods and Fast Marching Methods, Cambridge University Press, Cambridge, UK.
Mishra, A. , and Bolotnov, A. I. , 2013, “ Contact Angle Control Algorithm Development for Level-Set Interface Tracking Method,” Transactions of 2013 ANS Annual Meeting, Atlanta, GA, June 16–20.
Celik, I. B. , Ghia, U. , and Roache, P. J. , 2008, “ Procedure for Estimation and Reporting of Uncertainty Due to Discretization in CFD Applications,” ASME J. Fluids Eng., 130(7), p. 078001. [CrossRef]
Laplace, P. S. , Bowditch, N. , and Bowditch, N. I. , 1829, Mécanique Céleste, Meccanica.
Young, T. , 1855, Miscellaneous Works of the Late Thomas Young, Vol. 2.


Grahic Jump Location
Fig. 1

The mesh design of bubble departure case

Grahic Jump Location
Fig. 2

The schematic plot of contact angle calculation

Grahic Jump Location
Fig. 3

The schematic plot of target contact angle function

Grahic Jump Location
Fig. 4

The decision-making procedure for contact angle control modeling

Grahic Jump Location
Fig. 5

The schematic of contact angle control force application

Grahic Jump Location
Fig. 6

The schematic plot of F1 function

Grahic Jump Location
Fig. 7

Force balance analysis for the single bubble on the wall: (a) bubble division for force analysis and (b) force applied on the bubble

Grahic Jump Location
Fig. 8

The initial condition of single bubble departure case

Grahic Jump Location
Fig. 9

The bubble departure process for different target contact angles. The surface tension is equal to 0.0729 N/m in the simulations. CA indicates the prescribed target contact angle. It is noted that only the near-departure screenshots are shown here, and the first column is not the initial contact angle of the bubble. The gray level indicates the velocity magnitude field.

Grahic Jump Location
Fig. 10

The contact angle evolution history for different target contact angles

Grahic Jump Location
Fig. 11

The contact angle evolution for critical contact angle simulation

Grahic Jump Location
Fig. 12

The contact angle evolution against time for different meshes

Grahic Jump Location
Fig. 13

The schematic of capillary effect verification case

Grahic Jump Location
Fig. 14

The domain design of the capillary effect verification case. The gray color region on the left represented the computational domain. The zoom-in picture on the right shows the boundary layer mesh design.

Grahic Jump Location
Fig. 15

The mesh design for mesh resolution sensitivity study. The zoom-in figures show only part of the capillary domain. The elements parallel to the wall surface are boundary layer elements.

Grahic Jump Location
Fig. 16

The simulation results of capillary effect case (the coarsest mesh is on the left side while the finest one is on the right side). The target contact angle in this case is 45deg.

Grahic Jump Location
Fig. 17

The evolution of liquid rise and observed contact angle for different mesh sizes

Grahic Jump Location
Fig. 18

Screenshots of liquid rise between the flat plates for different contact angles

Grahic Jump Location
Fig. 19

The evolution of liquid rise and observed contact angle simulation for different contact angles



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