Experimental and Numerical Study of Taylor-Couette Flow" (2015)
Total Page:16
File Type:pdf, Size:1020Kb
Iowa State University Capstones, Theses and Graduate Theses and Dissertations Dissertations 2015 Experimental and numerical study of Taylor- Couette flow Haoyu Wang Iowa State University Follow this and additional works at: https://lib.dr.iastate.edu/etd Part of the Mechanical Engineering Commons Recommended Citation Wang, Haoyu, "Experimental and numerical study of Taylor-Couette flow" (2015). Graduate Theses and Dissertations. 14462. https://lib.dr.iastate.edu/etd/14462 This Dissertation is brought to you for free and open access by the Iowa State University Capstones, Theses and Dissertations at Iowa State University Digital Repository. It has been accepted for inclusion in Graduate Theses and Dissertations by an authorized administrator of Iowa State University Digital Repository. For more information, please contact [email protected]. Experimental and numerical study of Taylor-Couette flow by Haoyu Wang A dissertation submitted to the graduate faculty in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY Major: Mechanical Engineering Program of Study Committee: Michael G. Olsen, Major Professor Terry Meyer Shankar Subramaniam James Christian Hill Richard Dennis Vigil Iowa State University Ames, Iowa 2015 Copyright c Haoyu Wang, 2015. All rights reserved. ii TABLE OF CONTENTS LIST OF TABLES . iv LIST OF FIGURES . v ACKNOWLEDGEMENTS . xv ABSTRACT . xvi CHAPTER 1. INTRODUCTION . 1 General Introduction . .1 Dissertation Organization . .8 CHAPTER 2. POWER SPECTRUM ANALYSIS IN TAYLOR-COUETTE FLOW ......................................... 9 Abstract . .9 Introduction . 10 Experiment Apparatus and Procedure . 13 Data analysis . 17 Results and Discussion . 20 Summary and Conclusion . 28 CHAPTER 3. NUMERICAL INVESTIGATION OF TAYLOR-COUETTE FLOW WITH k − " MODEL . 43 Abstract . 43 Introduction . 44 Experimental Section . 46 CFD Methodology . 47 iii Results and Discussion . 52 Summary and Conclusion . 60 CHAPTER 4. NUMERICAL INVESTIGATION OF TAYLOR-COUETTE FLOW WITH k − ! MODEL . 93 Abstract . 93 Introduction . 94 Experimental Section . 96 CFD Methodology . 97 Results and Discussion . 102 Summary and Conclusion . 109 CHAPTER 5. CONCLUSION AND FUTURE DIRECTIONS . 136 Summary and Discussion . 136 Future Directions . 141 BIBLIOGRAPHY . 143 iv LIST OF TABLES Table 1.1 Numerical simulation works of Taylor-Couette flow . .7 v LIST OF FIGURES Figure 2.1 Example flow regimes of Taylor-Couette flow: (a) TVF (b) WVF (c) MWVF (d) TTVF (Fenstermacher et al., 1979) . 30 Figure 2.2 Apparatus of experiment: Cylinder Length: 43.18cm; Inside Cylin- der Diameter: 6.985cm; Outside Cylinder Diameter: 9.525cm .. 30 Figure 2.3 Schematic of experiment setup . 31 Figure 2.4 Example of PIV velocity fields . 31 Figure 2.5 PIV velocity fields for R=6..................... 32 Figure 2.6 Power Spectral Density for R=6 . 32 Figure 2.7 PIV velocity fields for R=11 . 33 Figure 2.8 Power Spectral Density for R=11 . 33 Figure 2.9 PIV velocity fields for R=15 . 34 Figure 2.10 Power Spectral Density for R=15 . 34 Figure 2.11 PIV velocity fields for R=17 . 35 Figure 2.12 Power Spectral Density for R=17 . 35 Figure 2.13 Contour plot of axial velocity for R=6 . 36 Figure 2.14 Contour plot of axial velocity for R=11 . 36 Figure 2.15 Contour plot of axial velocity for R=15 . 37 Figure 2.16 Contour plot of axial velocity for R=17 . 37 Figure 2.17 First Peak in PSD of R from 6 through 17 . 38 Figure 2.18 Second Peak in PSD of R from 6 through 17 . 38 Figure 2.19 PIV velocity fields for R=19 . 39 vi Figure 2.20 Power Spectral Density for R=19 . 39 Figure 2.21 PIV velocity fields for R=21 . 40 Figure 2.22 Power Spectral Density for R=21 . 40 Figure 2.23 PIV velocity fields for R=24 . 41 Figure 2.24 Power Spectral Density for R=24 . 41 Figure 2.25 PIV velocity fields for R=30 . 42 Figure 2.26 Power Spectral Density for R=30 . 42 Figure 3.1 3D Mesh generated for RANS simulation . 62 Figure 3.2 Grid resolution comparison of 2D: (a) Axial Velocity (mm=s), (b) Radial Velocity (mm=s), (c) Azimuthal Velocity (mm=s).... 63 Figure 3.3 Grid resolution comparison of 3D: (a) Axial Velocity (mm=s), (b) Radial Velocity (mm=s), (c) Azimuthal Velocity (mm=s).... 64 Figure 3.4 Example of simulation mean velocity field. (a) Full mean velocity profile of 2D simulation, (b) Mean velocity profile measurement domain. The axial coordinate is defined such that 0 location corresponds to the outflow boundary. 65 Figure 3.5 2D Velocity Profile and TKE Comparison at R = 60: (a) Ax- ial Velocity (mm=s), (b) Radial Velocity (mm=s), (c) Azimuthal Velocity (mm=s), (d) Normalized TKE . 65 Figure 3.6 2D Velocity Profile and TKE Comparison at R = 80: (a) Ax- ial Velocity (mm=s), (b) Radial Velocity (mm=s), (c) Azimuthal Velocity (mm=s), (d) Normalized TKE . 66 Figure 3.7 2D Velocity Profile and TKE Comparison at R = 130: (a) Ax- ial Velocity (mm=s), (b) Radial Velocity (mm=s), (c) Azimuthal Velocity (mm=s), (d) Normalized TKE . 66 vii Figure 3.8 2D Velocity Profile and TKE Comparison at R = 160: (a) Ax- ial Velocity (mm=s), (b) Radial Velocity (mm=s), (c) Azimuthal Velocity (mm=s), (d) Normalized TKE . 67 Figure 3.9 2D Velocity Profile and TKE Comparison at R = 190: (a) Ax- ial Velocity (mm=s), (b) Radial Velocity (mm=s), (c) Azimuthal Velocity (mm=s), (d) Normalized TKE . 67 Figure 3.10 2D Velocity Profile and TKE Comparison near wall (1/8 gap dis- tance from cylinder wall) at R = 60: (a) Axial Velocity (mm=s), (b) Radial Velocity (mm=s), (c) Azimuthal Velocity (mm=s), (d) Normalized TKE . 68 Figure 3.11 2D Velocity Profile and TKE Comparison near wall (1/8 gap dis- tance from cylinder wall) at R = 80: (a) Axial Velocity (mm=s), (b) Radial Velocity (mm=s), (c) Azimuthal Velocity (mm=s), (d) Normalized TKE . 69 Figure 3.12 2D Velocity Profile and TKE Comparison near wall (1/8 gap dis- tance from cylinder wall) at R = 130: (a) Axial Velocity (mm=s), (b) Radial Velocity (mm=s), (c) Azimuthal Velocity (mm=s), (d) Normalized TKE . 69 Figure 3.13 2D Velocity Profile and TKE Comparison near wall (1/8 gap dis- tance from cylinder wall) at R = 160: (a) Axial Velocity (mm=s), (b) Radial Velocity (mm=s), (c) Azimuthal Velocity (mm=s), (d) Normalized TKE . 70 Figure 3.14 2D Velocity Profile and TKE Comparison near wall (1/8 gap dis- tance from cylinder wall) at R = 190: (a) Axial Velocity (mm=s), (b) Radial Velocity (mm=s), (c) Azimuthal Velocity (mm=s), (d) Normalized TKE . 70 viii Figure 3.15 3D Velocity Profile and TKE Comparison at R = 60: (a) Ax- ial Velocity (mm=s), (b) Radial Velocity (mm=s), (c) Azimuthal Velocity (mm=s), (d) Normalized TKE . 71 Figure 3.16 Normalized Turbulence Dissipation Rate Comparison between PIV and Simulation at R=60 . 71 Figure 3.17 Normalized Turbulence Dissipation Rate Comparison between PIV and Simulation at R=80 . 72 Figure 3.18 Normalized Turbulence Dissipation Rate Comparison between PIV and Simulation at R=130 . 72 Figure 3.19 Normalized Turbulence Dissipation Rate Comparison between PIV and Simulation at R=160 . 73 Figure 3.20 Normalized Turbulence Dissipation Rate Comparison between PIV and Simulation at R=190 . 73 Figure 3.21 Normalized Turbulence Dissipation Rate Comparison near wall (1/8 gap distance from cylinder wall) between PIV and Simulation at R=60 ............................... 74 Figure 3.22 Normalized Turbulence Dissipation Rate Comparison near wall (1/8 gap distance from cylinder wall) between PIV and Simulation at R=80 ............................... 74 Figure 3.23 Normalized Turbulence Dissipation Rate Comparison near wall (1/8 gap distance from cylinder wall) between PIV and Simulation at R=130 . 75 Figure 3.24 Normalized Turbulence Dissipation Rate Comparison near wall (1/8 gap distance from cylinder wall) between PIV and Simulation at R=160 . 75 ix Figure 3.25 Normalized Turbulence Dissipation Rate Comparison near wall (1/8 gap distance from cylinder wall) between PIV and Simulation at R=190 . 76 Figure 3.26 Normalized Normal stress (u'u') of PIV data at different Reynolds numbers . 77 Figure 3.27 Normalized Normal stress (v'v') of PIV data at different Reynolds numbers . 77 Figure 3.28 Normalized Normal stress (w'w') of PIV data at different Reynolds numbers . 78 Figure 3.29 Normalized Shear stress (u'v') of PIV data at different Reynolds numbers . 78 Figure 3.30 Normalized Shear stress (u'w') of PIV data at different Reynolds numbers . 79 Figure 3.31 Normalized Shear stress (v'w') of PIV data at different Reynolds numbers . 79 Figure 3.32 Normalized Normal stress (u'u') near wall (1/8 gap distance from cylinder wall) of PIV data at different Reynolds numbers . 80 Figure 3.33 Normalized Normal stress (v'v') near wall (1/8 gap distance from cylinder wall) of PIV data at different Reynolds numbers . 80 Figure 3.34 Normalized Normal stress (w'w') near wall (1/8 gap distance from cylinder wall) of PIV data at different Reynolds numbers . 81 Figure 3.35 Normalized Shear stress (u'v') near wall (1/8 gap distance from cylinder wall) of PIV data at different Reynolds numbers . 81 Figure 3.36 Normalized Shear stress (u'w') near wall (1/8 gap distance from cylinder wall) of PIV data at different Reynolds numbers . 82 Figure 3.37 Normalized Shear stress (v'w') near wall (1/8 gap distance from cylinder wall) of PIV data at different Reynolds numbers . 82 x Figure 3.38 2D RANS Power Spectral Density for R=17 .