Development of a phase change model for volume-of-fluid method in OpenFOAM

Document Type : Full Lenght Research Article

Authors

Faculty of Mechanical Engineering, Babol Noshirvani University of Technology, Babol, Iran

Abstract

In this present study, volume of fluid method in OpenFOAM open source CFD package will be extended to consider phase change phenomena with modified model due to condensation and boiling processes. This model is suitable for the case in which both unsaturated phase and saturated phase are present and for beginning boiling and condensation process needn't initial interface. Both phases (liquid-vapor) are incompressible and immiscible. Interface between two phases is tracked with color function volume of fluid (CF-VOF) method. Surface Tension is taken into consideration by Continuous Surface Force (CSF) model. Pressure-Velocity coupling will be solved with PISO algorithm in the collocated grid. The accuracy of this phase-change model is verified by two evaporation problems (a one-dimensional Stefan problem and a two-dimensional film boiling problem) and two condensation problem (a one-dimensional Stefan problem and Filmwise condensation). The simulation results of this model show good agreement with the classical analytical or numerical results, proving its accuracy and feasibility.

Keywords

Main Subjects


[1] D.L. Youngs, Time-dependent multi-material flow with large fluid distortion, Numerical methods for fluid dynamics, 24 (1982) 273-285.
[2] C.W. Hirt, B.D. Nichols, Volume of fluid (VOF) method for the dynamics of free boundaries, Journal of Computational Physics, 39 (1981) 201-225.
[3] M. Huang, L. Wu, B. Chen, A piecewise linear interface-capturing volume-of-fluid method based on unstructured grids, Numerical Heat Transfer, Part B: Fundamentals, 61 (2012) 412-437.
[4] S. Osher, R.P. Fedkiw, Level set methods: an overview and some recent results, Journal of Computational Physics, 169 (2001) 463-502.
[5] P. Wang, H. Sun, P.Y. Wong, H. Fukuda, T. Ando, Modeling of droplet-based processing for the production of high-performance particulate materials using the level set method, Numerical Heat Transfer, Part A: Applications, 61 (2012) 401-416.
[6] M. Sussman, P. Smereka, S. Osher, A level set approach for computing solutions to incompressible two-phase flow, Journal of Computational physics, 114 (1994) 146-159.
[7] S.W. Welch, J. Wilson, A volume of fluid based method for fluid flows with phase change, Journal of computational physics, 160 (2000) 662-682.
[8] M.W. Akhtar, S.J. Kleis, Boiling flow simulations on adaptive octree grids, International Journal of Multiphase Flow, 53 (2013) 88-99.
[9] S. Hardt, F. Wondra, Evaporation model for interfacial flows based on a continuum-field representation of the source terms, Journal of Computational Physics, 227 (2008) 5871-5895.
[10] R.W. Schrage, A theoretical study of interphase mass transfer, Columbia University Press, 1953.
[11] D. Juric, G. Tryggvason, Computations of boiling flows, International Journal of Multiphase Flow, 24 (1998) 387-410.
[12] B.A. Nichita, J.R. Thome, A level set method and a heat transfer model implemented into FLUENT for modeling of microscale two phase flows, in:  AVT-178 Specialists" Meeting on System Level Thermal Management for Enhanced Platform Efficiency, 2010.
[13] W. Mao, Numerical Simulation of Vapor–liquid Phase Change Heat Transfer and Micromixing in Microfluidic Systems, Master's thesis, GuangZhou Institute of Energy Conversion Chinese Academy of Sciences, China, (2009).
[14] D. Guo, D. Sun, Z. Li, W. Tao, Phase change heat transfer simulation for boiling bubbles arising from a vapor film by the VOSET method, Numerical Heat Transfer, Part A: Applications, 59 (2011) 857-881.
[15] H.G. Weller, G. Tabor, H. Jasak, C. Fureby, A tensorial approach to computational continuum mechanics using object-oriented techniques, Computers in physics, 12 (1998) 620-631.
[16] J.U. Brackbill, D.B. Kothe, C. Zemach, A Continuum Method for Modeling Surface Tension, J. Comput. Phys., 100 (1992) 335-354.
[17] J. Klostermann, K. Schaake, R. Schwarze, Numerical simulation of a single rising bubble by VOF with surface compression, International Journal for Numerical Methods in Fluids, 71 (2013) 960-982.
[18] A. Albadawi, D. Donoghue, A. Robinson, D. Murray, Y. Delauré, Influence of surface tension implementation in volume of fluid and coupled volume of fluid with level set methods for bubble growth and detachment, International Journal of Multiphase Flow, 53 (2013) 11-28.
[19] N. Samkhaniani, A. Gharehbaghi, Z. Ahmadi, Numerical simulation of reaction injection molding with polyurethane foam, Journal of Cellular Plastics, 49 (2013) 405-421.
[20] W.H. Lee, Pressure iteration scheme for two-phase flow modeling, IN" MULTIPHASE TRANSPORT: FUNDAMENTALS, REACTOR SAFETY, APPLICATIONS". (1980) 407-432.
[21] M. Bahreini, A. Ramiar, A.A. Ranjbar, Numerical simulation of bubble behavior in subcooled flow boiling under velocity and temperature gradient, Nuclear Engineering and Design, 293 (2015) 238-248.
[22] M. Bahreini, A. Ramiar, A. Ranjbar, Numerical simulation of subcooled flow boiling under conjugate heat transfer and microgravity condition in a vertical mini channel, Applied Thermal Engineering, 113 (2017) 170-185.
[23] H.G. Weller, G. TaboraI, H. Jasak, C. Fureby, A tensorial approach to computational continuum mechanics using object-oriented techniques, Comput. Phys., 12 (1998) 620-632.
[24] C. Rhie, W. Chow, Numerical study of the turbulent flow past an airfoil with trailing edge separation, AIAA journal, 21 (1983) 1525-1532.
[25] H. Jasak, Error analysis and estimation for the finite volume method with applications to fluid flows, in:  Imperial College, Vol. PhD Thesis, University of London, 1996.
[26] H. Weller, A new approach to VOF-based interface capturing methods for incompressible and compressible flow, OpenCFD Ltd., Report TR/HGW/04, (2008).
[27] B. Van Leer, Towards the ultimate conservative difference scheme. II. Monotonicity and conservation combined in a second-order scheme, Journal of Computational Physics, 14 (1974) 361-370.
[28] M. Darwish, F. Moukalled, TVD schemes for unstructured grids, International Journal of Heat and Mass Transfer, 46 (2003) 599-611.
[29] B. Shu, F. Dammel, P. Stephan, Phase change model for two-phase fluid flow based on the volume of fluid method, ICHMT DIGITAL LIBRARY ONLINE, 13 (2008).
[30] G. Son, V. Dhir, Numerical simulation of film boiling near critical pressures with a level set method, Journal of Heat Transfer, 120 (1998) 183-192.
[31] P. Berenson, Film-boiling heat transfer from a horizontal surface, Journal of Heat Transfer, 83 (1961) 351-356.
[32] A. Gavrish, V.G. Rifert, Investigation of transition mode phenomenon from dropwise to filmwise steam condensation, Heat Transfer Research, 30 (1999).
[33] V. Semenov, N. Nikitin, Condensation heat transfer on noncircular pipes in stationary vapor, Heat Transfer Research, 39 (2008).
[34] Z. Liu, B. Sunden, J. Yuan, VOF modeling and analysis of filmwise condensation between vertical parallel plates, Heat Transfer Research, 43 (2012).
[35] S.M. Ghiaasiaan, Convective heat and mass transfer, Cambridge University Press, 2011.
[36] Z. Liu, B. Sunden, J. Yuan, Numerical Modelling of Condensation of Multiple Bubbles in Subcooled Flow Boiling With VOF Method, in:  ASME 2013 Heat Transfer Summer Conference collocated with the ASME 2013 7th International Conference on Energy Sustainability and the ASME 2013 11th International Conference on Fuel Cell Science, Engineering and Technology, American Society of Mechanical Engineers, 2013, pp. V001T003A014-V001T003A014.
[37] C. Aghanajafi, K. Hesampour, Heat transfer analysis of a condensate flow by VOF method, Journal of Fusion Energy, 25 (2006) 219-223.
[38] R. Marek, J. Straub, Analysis of the evaporation coefficient and the condensation coefficient of water, International Journal of Heat and Mass Transfer, 44 (2001) 39-53.