Chưa xác định

WELL-BALANCED ALGORITHM AND HEIGHT FUNCTION METHOD FOR DYNAMIC CONTACT ANGLE IN TWO-PHASE SYSTEMS

Năm XB 2024 Tạp chí / Hội thảo Multiphase Science and Technology Volume 36 (3) DOI / Link https://doi.org/10.1615/multscientechn.2024051265 ↗

Tác giả

Tóm tắt

The well-balanced algorithm combined with dynamic contact angle was well studied in the literature but was never implemented with the pressure-implicit with splitting of operators (PISO) algorithm in a collocated grid commonly used in an incompressible, transient simulation. This article presents a well-balanced algorithm for PISO schemes coupling with the height function method for curvature estimation. The dynamic contact angle model from Kistler and Cox is also integrated to improve the modelling of the curvature at the wall boundary. In collocated finite volume schemes, the well-balanced PISO algorithm is developed by modifying both the calculation of the gradients in the momentum equation and the Rhie and Chow algorithm. This new gradient calculation method ensures that surface tension force and pressure gradient are identically discretized at the same location. The Rhie and Chow algorithm is also modified by adding the surface tension force to balance the pressure forces. The stationary droplet case in two-dimensions is presented first to validate the proposed methodology. The well-balanced algorithm coupling with the height function method shows its benefits by damping spurious currents by two to three orders of magnitude. The 3D surface-driven flow and water-spreading droplets are then simulated; the results show that the new scheme coupled with dynamic contact angle model outperforms the unbalanced scheme of the smooth void fraction method for theoretical and experimental comparisons.

Tài liệu tham khảo

[1] Abadie, T., Aubin, J., and Legendre, D., On the Combined Effects of Surface Tension Force Calculation and Interface Advection on Spurious Currents within Volume of Fluid and Level Set Frameworks, J. Comput. Phys., vol. 297, pp. 611-636, 2015. DOI: 10.1016/j.jcp.2015.04.054

[2] Afkhami, S. and Bussmann, M., Height Functions for Applying Contact Angles to 2D VOF Simulations, Int. J. Numer. Methods Fluids, vol. 57, no. 4, pp. 453-472, 2008. DOI: 10.1002/fld.1651

[3] Afkhami, S. and Bussmann, M., Height Functions for Applying Contact Angles to 3D VOF Simulations, Int. J. Numer. Methods Fluids, vol. 61, no. 8, pp. 827-847, 2009. DOI: 10.1002/fld.1974

[4] Afkhami, S., Zaleski, S., and Bussmann, M., A Mesh-Dependent Model for Applying Dynamic Contact Angles to VOF Simulations, J. Comput. Phys., vol. 228, no. 15, pp. 5370-5389, 2009. DOI: 10.1016/j. jcp.2009.04.027

[5] Aulisa, E., Manservisi, S., Scardovelli, R., and Zaleski, S., Interface Reconstruction with Least-Squares Fit and Split Advection in Three-Dimensional Cartesian Geometry, J. Comput. Phys., vol. 225, no. 2, pp. 2301-2319, 2007. DOI: 10.1016/j.jcp.2007.03.015

[6] Benoit, M., Etudes Physique, Experimentale et Numerique des Mecanismes de Base Intervenant dans les Écoulements Diphasiques en Micro-Fluidique, PhD, Université de Provence, 2003.

[7] Brackbill, J., Kothe, D., and Zemach, C., A Continuum Method for Modeling Surface Tension, J. Comput. Phys., vol. 100, no. 2, pp. 335-354, 1992. DOI: 10.1016/0021-9991(92)90240-Y

[8] Bussmann, M., Mostaghimi, J., and Chandra, S., On a Three-Dimensional Volume Tracking Model of Droplet Impact, Phys. Fluids, vol. 11, no. 6, pp. 1406-1417, 1999. DOI: 10.1063/1.870005

[9] CONVERGE CFD Software, accessed December 7, 2022, from https://convergecfd.com/, 2022.

[10] Cox, R.G., The Dynamics of the Spreading of Liquids on a Solid Surface. Part 1. Viscous Flow, J. Fluid Mech., vol. 168, no. 1, p. 169, 1986. DOI: 10.1017/S0022112086000332

[11] Cummins, S.J., Francois, M.M., and Kothe, D.B., Estimating Curvature from Volume Fractions, Comput. Struct., vol. 83, nos. 6-7, pp. 425-434, 2005. DOI: 10.1016/j.compstruc.2004.08.017

[12] Evrard, F., Denner, F., and van Wachem, B., Estimation of Curvature from Volume Fractions Using Parabolic Reconstruction on Two-Dimensional Unstructured Meshes, J. Comput. Phys., vol. 351, pp. 271-294, 2017. DOI: 10.1016/j.jcp.2017.09.034

[13] Evrard, F., Denner, F., and van Wachem, B., Height-Function Curvature Estimation with Arbitrary Order on Non-Uniform Cartesian Grids, J. Comput. Phys.: X, vol. 7, p. 100060, 2020. DOI: 10.1016/j. jcpx.2020.100060

[14] Francois, M.M., Cummins, S.J., Dendy, E.D., Kothe, D.B., Sicilian, J.M., and Williams, M.W., A Balanced-Force Algorithm for Continuous and Sharp Interfacial Surface Tension Models within a Volume Tracking Framework, J. Comput. Phys., vol. 213, no. 1, pp. 141-173, 2006. DOI: 10.1016/j.jcp.2005.08.004

[15] Gohl, J., Mark, A., Sasic, S., and Edelvik, F., An Immersed Boundary Based Dynamic Contact Angle Framework for Handling Complex Surfaces of Mixed Wettabilities, Int. J. Multiphase Flow, vol. 109, pp. 164-177, 2018. DOI: 10.1016/j.ijmultiphaseflow.2018.08.001

[16] Han, S., Yang, R., Li, C., and Yang, L., The Wettability and Numerical Model of Different Silicon Microstructural Surfaces, Appl. Sci., vol. 9, no. 3, p. 566, 2019. DOI: 10.3390/app9030566

[17] Han, T.-Y., Zhang, J., Tan, H., and Ni, M.-J., A Consistent and Parallelized Height Function Based Scheme for Applying Contact Angle to 3D Volume-of-Fluid Simulations, J. Comput. Phys., vol. 433, p. 110190, 2021. DOI: 10.1016/j.jcp.2021.110190

[18] Hoffman, R.L., A Study of the Advancing Interface. I. Interface Shape in Liquid-Gas Systems, J. Colloid Interface Sci., vol. 50, no. 2, pp. 228-241, 1975. DOI: 10.1016/0021-9797(75)90225-8

[19] Jiang, M. and Zhou, B., Droplet Behaviors on Inclined Surfaces with Dynamic Contact Angle, Int. J. Hydrogen Energy, vol. 45, no. 54, pp. 29848-29860, 2020. DOI: 10.1016/j.ijhydene.2019.07.173

[20] Jiang, M., Zhou, B., and Wang, X., Comparisons and Validations of Contact Angle Models, Int. J. Hydrogen Energy, vol. 43, no. 12, pp. 6364-6378, 2018. DOI: 10.1016/j.ijhydene.2018.02.016

[21] Kistler, S.F., Hydrodynamics of Wetting, Wettability 6, pp. 251-310, 1993.

[22] Kothe, D.B. and Mjolsness, R.C., RIPPLE - A New Model for Incompressible Flows with Free Surfaces, AIAA J., vol. 30, no. 11, pp. 2694-2700, 1992. DOI: 10.2514/3.11286

[23] Le, T.T.H. and van Nguyen, C., Numerical Study of Partial Dam-Break Flow with Arbitrary Dam Gate Location using VOF Method, Appl. Sci., vol. 12, no. 8, p. 3884, 2022. DOI: 10.3390/app12083884

[24] Legendre, D. and Maglio, M., Comparison between Numerical Models for the Simulation of Moving Contact Lines, Comput. Fluids, vol. 113, pp. 2-13, 2015.

[25] Malgarinos, I., Nikolopoulos, N., Marengo, M., Antonini, C., and Gavaises, M., VOF Simulations of the Contact Angle Dynamics during the Drop Spreading: Standard Models and a New Wetting Force Model, Adv. Colloid Interface Sci., vol. 212, pp. 1-20, 2014. DOI: 10.1016/j.cis.2014.07.004

[26] Meier, M., Yadigaroglu, G., and Smith, B.L., A Novel Technique for Including Surface Tension in PLIC-VOF

[27] Methods, Eur. J. Mech. B/Fluids, vol. 21, no. 1, pp. 61-73, 2002. DOI: 10.1016/S0997-7546(01)01161-X

[28] Patel, H.V., Kuipers, J., and Peters, E., Computing Interface Curvature from Volume Fractions: A Hybrid Approach, Comput. Fluids, vol. 161, pp. 74-88, 2018. DOI: 10.1016/j.compfluid.2017.11.011

[29] Popinet, S., An Accurate Adaptive Solver for Surface-Tension-Driven Interfacial Flows, J. Comput. Phys., vol. 228, no. 16, pp. 5838-5866, 2009. DOI: 10.1016/j.jcp.2009.04.042

[30] Popinet, S., Numerical Models of Surface Tension, Ann. Rev. Fluid Mech., vol. 50, no. 1, pp. 49-75, 2018. DOI: 10.1146/annurev-fluid-122316-045034

[31] Raessi, M., Mostaghimi, J., and Bussmann, M., Advecting Normal Vectors: A New Method for Calculating Interface Normals and Curvatures when Modeling Two-Phase Flows, J. Comput. Phys., vol. 226, no. 1, pp. 774-797, 2007. DOI: 10.1016/j.jcp.2007.04.023

[32] Rhie, C.M. and Chow, W.L., Numerical Study of the Turbulent Flow past an Airfoil with Trailing Edge Separation, AIAA J., vol. 21, no. 11, pp. 1525-1532, 1983. DOI: 10.2514/3.8284

[33] Rider, W.J. and Kothe, D.B., Reconstructing Volume Tracking, J. Comput. Phys., vol. 141, no. 2, pp. 112-152, 1998. DOI: 10.1006/jcph.1998.5906

[34] Rioboo, R., Marengo, M., and Tropea, C., Time Evolution of Liquid Drop Impact onto Solid, Dry Surfaces, Exp. Fluids, vol. 33, no. 1, pp. 112-124, 2002. DOI: 10.1007/s00348-002-0431-x

[35] Roisman, I.V., Opfer, L., Tropea, C., Raessi, M., Mostaghimi, J., and Chandra, S., Drop Impact onto a Dry Surface: Role of the Dynamic Contact Angle, Colloids Surf. A: Physicochem. Eng. Aspects, vol. 322, nos. 1-3, pp. 183-191, 2008. DOI: 10.1016/j.colsurfa.2008.03.005

[36] Scardovelli, R. and Zaleski, S., Direct Numerical Simulation of Free-Surface and Interfacial Flow, Ann. Rev. Fluid Mech., vol. 31, no. 1, pp. 567-603, 1999. DOI: 10.1146/annurev.fluid.31.1.567

[37] Shikhmurzaev, Y.D., The Moving Contact Line on a Smooth Solid Surface, Int. J. Multiphase Flow, vol. 19, no. 4, pp. 589-610, 1993. DOI: 10.1016/0301-9322(93)90090-H

[38] Sikalo, S., Wilhelm, H.-D., Roisman, I.V., Jakirlic, S., and Tropea, C., Dynamic Contact Angle of Spreading Droplets: Experiments and Simulations, Phys. Fluids, vol. 17, no. 6, p. 62103, 2005. DOI: 10.1063/1.1928828

[39] Sui, Y. and Spelt, P.D., An Efficient Computational Model for Macroscale Simulations of Moving Contact Lines, J. Comput. Phys., vol. 242, pp. 37-52, 2013. DOI: 10.1016/j.jcp.2013.02.005

[40] Tryggvason, G., Scardovelli, R., and Zaleski, S., Direct Numerical Simulations of Gas-Liquid Multiphase Flows, Cambridge, UK: Cambridge University Press, 2011.

[41] Ukiwe, C. and Kwok, D.Y., On the Maximum Spreading Diameter of Impacting Droplets on Well-Prepared Solid Surfaces, Langmuir, vol. 21, no. 2, pp. 666-673, 2005. DOI: 10.1021/la0481288

[42] Youngs, Time-Dependent Multi-Material Flow with Large Fluid Distortion, Numerical Methods in Fluid Dynamics, Cambridge, MA: Academic Press, p. 273, 1982.

[43] Youngs, D.L., An Interface Tracking Method for a 3D Eulerian Hydrodynamics Code, Tech. Rep., AWRE, 1984.

[44] Zhang, K., Li, Y., Chen, Q., and Lin, P., Numerical Study on the Rising Motion of Bubbles near the Wall, Appl. Sci., vol. 11, no. 22, p. 10918, 2021. DOI: 10.3390/app112210918