Flume Studies of Gravel Bed Surface Response to Flowing
Total Page:16
File Type:pdf, Size:1020Kb
FLUME STUDIES OF GRAVEL BED SURFACE RESPONSE TO FLOWING WATER By JOHN FREDRIC WOLCOTT B.Sc, B.F.A., University of Washington, 1982 M.Sc, University of British Columbia, 1984 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY in THE FACULTY OF GRADUATE STUDIES Department of Geography We accept this thesis as conforming to the required standard THE UNIVERSITY OF BRITISH COLUMBIA January, 1990 ® John Fredric Wolcott, 1990 08 In presenting this thesis in partial fulfilment of the requirements for an advanced degree at the University of British Columbia, I agree that the Library shall make it freely available for reference and study. I further agree that permission for extensive copying of this thesis for scholarly purposes may be granted by the head of my department or by his or her representatives. It is understood that copying or publication of this thesis for financial gain shall not be allowed without my written permission. Department of Geography The University of British Columbia Vancouver, Canada Date Ffthrnary 7. DE-6 (2/88) ii ABSTRACT Almost all sediment transport equations incorporate the Shields parameter, which is a ratio of the total boundary shear stress as a driving force and the particle weight as a resisting force. Shields (1936) equated particle resistance to entrainment with particle weight, which is proportional to particle diameter, or bed texture. The present work analyses the particle resistance term in the Shields parameter. As the bed material adjusts to a given flow condition, bed stability increases. The arrangement of particles into more stable configurations is here termed geometric structure, and includes the formation of pebble clusters, and imbrication. After an initial surface coarsening, here termed textural structure, particle resistance to movement is a function primarily of geometric structure. The Shields number for entrainment is thus a measure of particle resistance due to both types of bed structure rather than the conventional notion of particle resistance due to particle weight. The response of a mobile bed surface composed of < 8 mm diameter gravels to flowing water was explored in a 6 meter by 0.5 meter flume using four different slopes and various water depths. Corrected bed shear stresses varied between 0.05 and 2.79 Pa. Step increases in discharge with a constant slope caused the bed surface to develop a structure which was more stable at the end of a run than at the beginning. Under these conditions, the Shields number for incipient motion was found to vary between 0.001 and 0.066. This variability can be explained by the degree of geometric structure present. Previous studies, including Shields' work (1936), have implicitly included the effects of geometric structure on incipient motion. iii Surface coarsening develops with very low flows, but subsequent coarsening in higher flows is minor, with less than 5% increase in median diameter following a 50% increase in bed shear stress. Calculations of Manning's n based on depth, slope, and velocity measurements show an increase in flow resistance as structure develops. The development of a coarse surface layer appears to be limited by flow characteristics near the bed which are in turn modified by the development of structure. Measurements of the area occupied by the largest stones show that they do not cover more than 14% of the surface during maximum coarsening. Froude scaling of the flume data indicates that the time necessary for development of maximum strength is on the order of a month for natural rivers under steady flow conditions. This suggests that gravel river beds are rarely in equilibrium with natural flow conditions. iv TABLE OF CONTENTS Abstract ii Table of Contents iv List of Tables vi List of Figures vii Acknowledgements ix CHAPTER ONE: INTRODUCTION 1 Prologue 1 Field Studies , 2 Flume Experiments 4 Shields' Work 6 Other Flume Studies 9 Objectives of This Study 14 CHAPTER TWO: METHODOLOGY 17 Justification of Approach 17 Models.... 18 A Note on Units 21 Experimental Apparatus 22 Experimental Procedure 28 Summary of Runs 34 CHAPTER THREE: OBSERVATIONS 41 Series One Runs 41 A Qualitative Model 54 CHAPTER FOUR: BED STABILITY AND INCIPIENT MOTION 57 An Overview 57 The Effect of Bed Structure Development 60 Bed Structure and Incipient Motion 63 V The Shields Parameter Re-Examined 71 Bed Structure and Manning's n 77 Summary 78 CHAPTER FIVE: LAG TIME 80 CHAPTER SIX: LIMITS TO COARSENING 88 CHAPTER SEVEN: CONCLUSIONS 96 REFERENCES 103 APPENDIX 108 VI LIST OF TABLES 2.1 Velocity Profiles Across the Test Section 25 2.2 Hydraulic Parameters of the Experiments 35 4.1 Measured and Calculated Sediment Transport Values for Later Runs 64 4.2 Shields Parameter Based on Various D50 68 4.3 The D50 and the Amounts of the +4 mm and +5.66 mm Material Captured in Series 3 Runs 69 5.1 Variations in Transport Rates through Time for Series 3 Runs 82 5.2 Variations in Transport Rates through Time for Series 4 Runs 83 5.3 r^ for Linearized Exponential and Power functions of Sediment Transport through Time 86 6.1 Weights and Percentage of the Total Area of the + 5.6 mm Size Class for Several Surface Samples 91 A.l Weights and Areas of Various Size Classes 110 A.2 Ratio of A and B Axes and Areas Occupied by weight Ill vii LIST OF FIGURES 1.1 Shields Curve, after Shields (1936) 8 1.2 Illustration of a particle cluster and imbrication 12 2.1 Sketch of flume 23 2.2 Cumulative grain size distribution curves for sediment used in the experimental runs 29 3.1 Loose, well mixed material (top) and armored, and settled material after Run 1-3 (bottom) 42 3.2 Surface web of larger stones as seen on bed (top), and with smaller stones removed for clarity (bottom) 44 3.3 Photograph of incipient web development after Run 4-3 45 3.4 Comparison of bulk sample with surface samples after 0.02 m, 0.03 m, 0.05 m deep flows (Runs 1-1 to 1-4), and thalweg from degrading run (Run 1-5) 49 3.5 Textural differences between channel surface (above) and bar top surface (below). Flow from left to right 52 3.6 Comparison of bulk sample, sediment trapped during degrading run, and surface sample from thalweg and top of bar 53 4.1 Proposed relation between hydraulic variables and sediment supply (Parker, 1987, pers. com.) 58 4.2 Cumulative grain size distribution of trapped material during Series 4 runs 62 4.3 Decline in transport rate through time 66 4.4 Proposed model of bed structure development 73 viii 4.5 Shields parameter based on material moved compared to Shields parameter based on subsurface D50 74 4.6 Difference in Shields parameters of Figure 4.5 (which is the excess bed strength due to geometric structure) compared to Shields parameter based on subsurface D50 75 5.1 Coarsening of the surface median sizes following disturbance by three major storms, with discharges of 0.068, 0.039, and 0.056 m3/s 84 6.1 Comparison of typical surface (top, after 0.02 m flow depth, Run 1-2) and very well sorted surface (bottom, after Run 1-5) 95 ix ACKNOWLEDGEMENTS Since it is probable that this is the most read page of any thesis, I feel compelled to mention two things, before getting on with the traditional thank yous. First, in my two and half decades of doing research in two countries, which included one province and three states, the most enjoyable experience was completing the M.Sc. The most miserable was dealing with the examining committee during and following the PhD defense, and the administration. The bitterness has not yet dissipated. In fact, it is only through incredible encouragement of a small circle of very good friends that this thesis has continued to the final stage. May it rest in peace on some dark and dusty library shelf. Second, Dr. Humphrey of Cal. Tech. has pointed out the following irony. Since an examiner has asked that the centimeter/gram/second system of units be changed to SI units, this thesis has grown in length and therefore stature simply by the addition of zeros. Think about this when you try and figure out what 0.000158 m^/s means (it is 9.5 liters per minute or 2.5 gallons per minute), or recognize that if a 2.3 mm particle rolls 0.005 m, it has gone 2 grain diameters. To the examiner who edited the original acknowledgements page I apologize. I guess the thesis was dropped on the way to the library and the wrong page inserted. The success of this thesis is the result of generous help from many individuals, a few of whom are mentioned here. First and foremost, I would like to thank Dr. M. Church, one of the best supervisors I have ever met. He supported my initial explorations into this subject in spite of no clear objective on my part, other than curiosity. Later, through suggested readings and provocative questions, he guided me to the point where everything fell into place, and even I could see the way home. Dr. O. Slaymaker of UBC and Dr. T. Hicken of SFU provided moral support and helpful suggestions during Dr. Church's sabbatical leave, and I am very grateful for their encouragement during a discouraging period. Dr. P. Ward of Ward and Associates, Vancouver, was extremely helpful with pump and motor calculations and the design of fluid flow systems. In addition to my committee, I would like to thank Dr.