Characterization of Evaporation/Condensation During Pool Boiling and Flow Boiling Mostafa Mobli University of South Carolina
Total Page:16
File Type:pdf, Size:1020Kb
University of South Carolina Scholar Commons Theses and Dissertations 2018 Characterization Of Evaporation/Condensation During Pool Boiling And Flow Boiling Mostafa Mobli University of South Carolina Follow this and additional works at: https://scholarcommons.sc.edu/etd Part of the Mechanical Engineering Commons Recommended Citation Mobli, M.(2018). Characterization Of Evaporation/Condensation During Pool Boiling And Flow Boiling. (Doctoral dissertation). Retrieved from https://scholarcommons.sc.edu/etd/4898 This Open Access Dissertation is brought to you by Scholar Commons. It has been accepted for inclusion in Theses and Dissertations by an authorized administrator of Scholar Commons. For more information, please contact [email protected]. CHARACTERIZATION OF EVAPORATION/CONDENSATION DURING POOL BOILING AND FLOW BOILING by Mostafa Mobli Bachelor of Science University of Tehran, 2012 Master of Science University of South Carolina, 2014 Submitted in Partial Fulfillment of the Requirements For the Degree of Doctor of Philosophy in Mechanical Engineering College of Engineering and Computing University of South Carolina 2018 Accepted by: Chen Li, Major Professor Jamil Khan, Committee Member Tanvir Farouk, Committee Member Yi Sun, Committee Member Cheryl L. Addy, Vice Provost and Dean of the Graduate School © Copyright by Mostafa Mobli, 2018 All Rights Reserved. ii DEDICATION To my wife Sara and my parents, for their endless love, support, and encouragement. iii ACKNOWLEDGEMENTS I would like to express my deepest appreciation to my advisor Dr. Chen Li for his guidance, patience and support. Thank you for forcing me to look at research and my work in different ways. Thank you for believing in me and helping me stand up after every fall. Your support was essential to my success. I would like to thank my committee members, Dr. Jamil Khan and Dr. Tanvir Farouk, and Dr. Yi Sun for their invaluable advice and help throughout my PhD. My gratitude also extends to my friend and colleague Dr. Mehdi Famouri. Much of my work would have not been completed without your assistance. I owe my deepest gratitude to my parents, for their endless love, support, and encouragement. No words can describe my love to you both. I would like to thank my beautiful wife Sara Nakhi for her invaluable friendship and support in difficult times and my fellow graduate students for their friendship and assistance. I will always remember the great times we had. At last, I would like to thank the members of my department, Mechanical engineering. The faculty, staff, and students made my stay in Columbia a great experience. iv ABSTRACT Present dissertation has investigated pool and flow boiling and their characteristics via numerical means. A code was developed to investigate and enhance heat transfer performance during different modes of phase change phenomena. Multiphase heat transfer has proven to be one of the most effective means of heat transfer in different industries, therefore, there have been numerous experimental and numerical studies on the subject of phase change phenomena in a wide range of conditions and setups; yet there are complex bubble dynamics and heat transfer characteristics that remain unresolved. To have a more detailed look at and a better understanding of complex characteristics of phase change phenomena, our code focused on the mostly unresolved parts of this phenomena, such as spurious currents, interface diffusion, phase change modeling, micro-layer heat transfer, conjugate heat transfer effects and interfacial heat transfer coefficient. These complexities arise mostly because of small scale of the phase change phenomena and pace of phase change heat transfer, these scaling issues makes it difficult or in some cases impossible to design a robust and comprehensive experiment, which can study different aspects of phase change heat transfer. On the numerical side, lack of exact solutions and equations to phase change can cause immense problems in modeling and numerical studies. To mention a few of these numerical difficulties one can mention bubble or droplet curvature estimation which does not have an exact mathematical v solution that can translate to a viable algorithm, or interfacial temperatures which is said to be most important factor driving phase change. To address some of these difficulties in numerical simulations we have employed volume of fluid method which benefits from global mass conservation combined with level set method which shows a more promising interface curvature estimation. Combining the two methods has proven to be a challenging task and, in some cases, not so much superior; therefore, a simplified method was employed to capture the best of the two methods. Phase change source terms was simulated based on none equilibrium conditions which states that phase change happens because of deviations of interface temperature from saturation temperature, unlike equilibrium condition which maintains the interface at saturation conditions. Using none equilibrium conditions forces a smaller grid onto simulation, which was cared by introduction of smearing factor. Other numerically challenging phenomena is micro-layer heat transfer which is mostly resolved by simplifying continuity, momentum and energy equations and deriving a set ODEs that are solved outside of main simulation algorithm. This is mostly due the fact micro-layer is mostly sub grid phenomena that cannot be seen by conventional CFD codes. We have employed a method that solves the micro-layer within the main algorithm without the need of solving those simplified set of ODEs and includes the none equilibrium interface conditions in micro-layer simulations. Present dissertation contains a comprehensive literature review on numerical and experimental studies on the subject of two phase flow temperature driven phase change and heat and mass transfer. Which is followed by a detailed description of mathematical background and code algorithm. Then we have validated our code against numerous vi experimental studies available in the literature. Eventually, interfacial heat transfer coefficient during sub cooled pool and flow boiling was studied, which has never been studied numerically and there a few experimental studies on it. vii TABLE OF CONTENTS DEDICATION ............................................................................................................................... iii ACKNOWLEDGEMENTS ........................................................................................................... iv ABSTRACT .................................................................................................................................... v LIST OF TABLES .......................................................................................................................... x LIST OF FIGURES ....................................................................................................................... xii LIST OF SYMBOLS .................................................................................................................... xvi LIST OF ABBREVIATIONS .................................................................................................... xviii CHAPTER 1 - INTRODUCTION .................................................................................................. 1 LITERATURE REVIEW ................................................................................................... 4 WALL HEAT FLUX PARTITIONING ............................................................................ 9 SURFACE CHARACTERISTICS ................................................................................... 13 FLOW BOILING STUDIES ............................................................................................ 19 SIMULATION ................................................................................................................. 22 VOLUME OF FLUIDS .................................................................................................... 27 LEVEL SET METHOD ................................................................................................... 32 LATTICE BOLTZMANN METHOD ............................................................................. 35 PHASE CHANGE MODEL ............................................................................................ 40 PROPOSED MODEL IN THE PRESENT WORK ......................................................... 41 CHAPTER 2 - GOVERNING EQUATIONS AND NUMERICAL ANALYSIS ........................ 45 GOVERNING EQUATIONS .......................................................................................... 45 NUMERICAL ANALYSIS METHODS AND MODEL IMPLEMENTATION ............ 51 SOLUTION PROCEDURE AND DISCRETIZATION SCHEMES .............................. 61 A BRIEF INTRODUCTION ON DISCRETIZATION SCHEMES ................................ 64 viii CHAPTER 3 - MODEL VALIDATION AND GRID INDEPENDENCY STUDY .................... 68 SPURIOUS CURRENTS ................................................................................................. 68 BUBBLE GROWTH RATE AND PHASE CHANGE METHOD ................................. 73 GRID INDEPENDENCY ................................................................................................ 75 CHAPTER 4 - POOL BOILING ................................................................................................... 78 CASE I - POOL BOILING IN A POOL OF R113 .........................................................