Wave-Current Interaction in Water of Finite Depth Zhenhua Huang 0 6 8
Total Page:16
File Type:pdf, Size:1020Kb
Wave-current interaction in water of finite depth by S9IJV8fl Zhenhua Huang 0 6 8 B.S., Tsinghua University (1988) AO.1ONH031JO M.E., Tsinghua University(1991) _3.1J.LSNI S.L3sInHOVSSYW Submitted to the Department of Civil and Environmental Engineering in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Environmental Fluid Mechanics at the MASSACHUSETTS INSTITUTE OF TECHNOLOGY February 2004 @ Massachusetts Institute of Technology 2004. All rights reserved. A uthor .. ................... Department of Civil and Eyronmental Engineering September 12, 2003 Certifi d by .... ...... .... Certified by. ~ ~Chiang C. Mei Donald Martha Harleman Professor of Civil & Environmental Engineering Thesis Supervisor A ccepted by ............................ Heidi M.' epI Associate Professor of Civil & Environmental Engineering Chairman, Department Committee on Graduate Students MASSACHUSETS iNST E OF TECHNOLOGY FEB 19 2004 LIBRARIES 2 Wave-current interaction in water of finite depth by Zhenhua Huang Submitted to the Department of Civil and Environmental Engineering on September 12, 2003, in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Environmental Fluid Mechanics Abstract In this thesis, the nonlinear interaction of waves and current in water of finite depth is studied. Wind is not included. In the first part, a 2D theory for the wave effect on a turbulent current over rough or smooth bottom is presented. The logarithmic profile of the basic current is modified by the waves due to an effective mean shear stress on the free surface. Surface distortion of the eddy viscosity is shown to be important for the change of the mean velocity profile by waves. Both wave-following and wave-opposing current are studied. Comparisons are made with some existing laboratory experiments. In the second part, an instability theory is presented for the initiation of Langmuir cells due to waves interacting with a turbulent current maintained by tides or by an external pressure gradient. With an infinitesimal span-wise disturbance, the free surface experiences a new mean stress, which generates new vorticity to be diffused downward, and induces further growth. Various contributions to the unstable growth of Langmuir circulation are analyzed by examining the mechanical energy budget. Evidences will be shown that the surface stress contributes significantly to instability. Both wave-following current and wave-opposing current are studied. For the wave- following current, two types of Langmuir cells can grow in time; while for the wave- opposing current, only one can grow. Effects of current strength, wave conditions, and water depth on the growth of Langmuir circulation will be studied by numerical examples. Remarks on existing laboratory experiments are made Thesis Supervisor: Chiang C. Mei Title: Donald & Martha Harleman Professor of Civil & Environmental Engineering Acknowledgments In the past seven years at MIT, I owe much to many people who have contributed considerable to this thesis directly or indirectly. First, I would like to thank my advisor, Professor Chiang C. Mei. I am lucky to be able to work with and learn from him these years. His breadth of knowledge, his patience with students, and his high standards have not only guided my research to success but also helped me to grow as a scientist. I would also like to thank Caroline, Professor Mei's wife, for her kindness. I will never forget those warm parties in their home, which provide comfort for my wife and me during these years away from home. I also thank Professor 0. S. Madsen for many fruitful discussions. I would also like to thank the other members of my thesis committee Professor H. M. Nepf and Professor T. R. Akylas for taking time to review and evaluate my work. I am also indebted to the many friends in Parsons Lab who made my graduate experience educational and pleasurable. I am grateful for the financial support by US Office of Naval Research (Grant N00014-89J-3128, directed by Dr. Thomas Swean) and US National Science Founda- tion (GrantCTS-0075713, directed by Profs. John Foss, C. F. Chen and M. Plesniak). I also thank the Ippen Fellowship for supporting me to attend two conferences in the past two years. Finally, I would like to thank my wife, Lixia, and my parents. I could never have reached this point without their understanding and constant support. To them I owe more than I can ever express. 6 Contents I Effects of surface waves on a turbulent current over a smooth or rough seabed 23 1 Introduction 25 2 Formulation 31 2.1 Governing equations and boundary conditions in dimensional form . 31 2.2 Eddy viscosity . ......................... ..... 32 2.3 Basic assumptions ................... .......... 35 2.4 Normalization by core scales ..... ....... ...... ..... 37 2.5 Basic current ... ........ ......... ........ ... 40 2.6 Length scale of wave attenuation ....... ............ 41 3 Mean and oscillatory motions in the core region 43 3.1 Order estimate of wave-induced core current .. ..... ..... .. 43 3.2 Momentum balance of mean motion ..... ..... ..... .... 44 3.3 Vortex force zDQ .... ..... ..... ..... ..... ..... 46 3.4 Mean shear stress of wave-perturbed current . .... .... .... 47 4 Bottom boundary layer 51 4.1 Order of magnitude estimate ..... ......... 51 4.2 Approximate momentum balance of the mean motion 52 4.3 Oscillatory motion: b. ... .... ..... .... 54 4.4 Wave damping rate ...... ...... ...... 55 5 Free surface 57 5.1 Unimportance of surface wave boundary layer .. 57 7 5.2 Perturbed mean Reynolds stress on the still water surface . .5 59 6 Final solution for the core 63 6.1 Formulas for mean turbulent stress and current velocity ........ 63 6.2 Friction velocities ............................. 64 6.3 Numerical procedure ........................... 66 7 Comparison with experiments 69 7.1 Past experim ents .............. ............... 69 7.2 Wave-following current over a smooth bed by Kemp & Simons .. .. 71 7.3 Wave-following current over a rough bed by Kemp & Simons ..... 73 7.4 Wave-opposing current over a rough bed by Kemp & Simons ..... 76 7.5 Wave-following/opposing currents over a rough bed by Klopman . 80 7.6 Longitudinal variation of current .................... 80 8 Concluding remarks of part I 83 II Langmuir circulation in water of finite depth 85 9 Introduction 87 10 Equations and boundary conditions for total motion 95 10.1 Norm alization ...... ...... ...... ...... ....... 95 10.2 Governing equations of total motion ..... ...... ...... .. 98 10.3 Surface boundary conditions . ..... ..... ..... ..... .. 100 10.3.1 Exact surface boundary conditions on free surface .. .... 100 10.3.2 Surface boundary conditions at z = 0 . ...... ...... 101 10.3.3 Mean surface boundary conditions at z = 0 ..... ..... 106 11 Basic state 109 11.1 Basic state for Langmuir circulation .... .... .... .... ... 109 11.2 Shear rate of basic current .. ... .... ... .... ... .... 112 8 12 Linearized equations and boundary conditions for Langmuir circu- lation L17 12.1 Spanwise Perturbations ... ..... ...... ...... ..... 117 12.2 Secondary waves ...... ....... ........ ....... 119 12.3 Linearized equations governing the initial Langmuir circulation .. 122 12.3.1 Continuity equation ....... ....... ....... .. 122 12.3.2 Equation of mean longitudinal momentum . ....... .. 122 12.3.3 Equation of mean lateral momentum . ........ .... 125 12.3.4 Equation of mean vertical momentum .. ........... 127 12.3.5 Equation of mean longitudinal vorticity ......... ... 127 12.3.6 Equation of mean vertical vorticity .......... .... 129 12.4 Linearized free-surface boundary conditions for Langmuir circualtion 129 12.5 Bottom boundary conditions for perturbed current ........ 135 12.5.1 Spanwise perturbations ......... ........... 137 12.5.2 Mean motion associated Langmuir circulation inside BWBL 138 13 Initial Energy Budget of Langmuir Circulation L43 13.1 Equation of mechanical energy ..................... 143 13.2 Growth rate of the normal mode ................. ... 150 13.3 Contribution of BWBL to the growth rate ... ........... 153 14 Generation mechanism of Langmuir circulation 157 14.1 Review of inviscid problem ............. ... .. 157 14.2 The feedback mechanism of Craik and Leibovich . ... .. 159 14.2.1 Instability criterion for inviscid limit . .. .. .. 159 14.2.2 Role of viscosity ............. .. ... .. 162 14.2.3 Feedback mechanism of CL-II ....... .. .. 163 14.3 The surface stress mechanism ......... .. .. .. 166 14.4 Circualtion generated by sidewalls . ........ ... .... 170 14.5 Sum m ary ...................... .. .. .. 173 15 Numerical Scheme 175 15.1 Normal mode .................... ... .. 175 9 15.2 R escaling ........... ............. .. .. .. .. 176 15.3 Discretization Scheme ... ......... ..... .. .. .. .. 178 15.4 Coupled surface boundary condition .... ...... ..... .... 182 15.5 A test example ............... ...... ... ... ... 184 15.6 Solution procedure .. ................. ...... .. 187 15.6.1 Typical parameters ............. .. ... ... .. 187 15.6.2 Bottom roughness ........... ..... .... ..... 189 15.6.3 Maximum wave slope ......... ..... .. ....... 190 15.6.4 Solution procedure .......... ..... .. ....... 190 15.7 Numerical accuracy ....... ......... ... .... ... 191 16 Langmuir circulation in a wave-following current 195 16.1 Unstable eigen-modes .... ......... ..... .. ... ... 195 16.1.1 Eigenvalue spectrum .. ............ ......... 197 16.1.2 Energy budget .. ............... ......... 199 16.1.3 Eigenfunctions