THREE DIMENSIONAL MOBILE BED DYNAMICS FOR SEDIMENT TRANSPORT MODELING
DISSERTATION
Presented in Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy in the
Graduate School of The Ohio State University
By
Sean O’Neil, B.S., M.S.
*****
The Ohio State University
2002
Dissertation Committee: Approved by
Professor Keith W. Bedford, Adviser Professor Carolyn J. Merry Adviser Professor Diane L. Foster Civil Engineering Graduate Program c Copyright by
Sean O’Neil
2002 ABSTRACT
The transport and fate of suspended sediments continues to be critical to the understand- ing of environmental water quality issues within surface waters. Many contaminants of environmental concern within marine and freshwater systems are hydrophobic, thus read- ily adsorbed to bed material or suspended particles. Additionally, management strategies for evaluating and remediating the effects of dredging operations or marine construction, as well as legacy pollution from military and industrial processes requires knowledge of sediment-water interactions. The dynamic properties within the bed, the bed-water column inter-exchange and the transport properties of the flowing water is a multi-scale nonlin- ear problem for which the mobile bed dynamics with consolidation (MBDC) model was formulated.
A new continuum-based consolidation model for a saturated sediment bed has been developed and verified on a stand-alone basis. The model solves the one-dimensional, vertical, nonlinear Gibson equation describing finite-strain, primary consolidation for satu- rated fine sediments. The consolidation problem is a moving boundary value problem, and has been coupled with a mobile bed model that solves for bed level variations and grain size fraction(s) in time within a thin layer at the bed surface. The MBDC model represents the first attempt to unify bed exchange and accounting mechanisms with vertically varying bed properties under a single mechanistic framework.
ii A suspended sediment transport solver, with parameterizations for noncohesive grain size settling velocity, erosion and depositions source sink terms has been extended to in- clude parameterizations for cohesive grain sizes. Further, the consolidation model has been integrated into the mobile bed modeling framework. The new fine-grained sediment trans- port model, MBDC, was configured to simulate the flow, transport and bottom evolution within an expansion channel serving as an idealized conical estuary. MBDC model results, are compared with model results from literature, demonstrating qualitative agreement and model efficacy. The MBDC model approach, though requiring more site specific data for auxiliary parameterizations, yields a more complete physical and dynamic description of bed sediment transport processes.
iii To my dearest friend and wife Chen Hui.
You have inspired me to be more than I thought I could ever be.
iv ACKNOWLEDGMENTS
I would like to acknowledge the enthusiastic help from my friends and colleagues at the
Great Lakes Forecasting System Laboratory, Dr. Philip Chu, Dr. David Welsh, Takis and
Vasso Velisariou, Guo Yong, and especially Dr. Jennifer Shore and Heather Smith. Past members of the GLFS Lab and the “Dirt Group”, who also helped me were Drs. David
Podber, John Kelley, James Yen, W. K. Yeo, and my brothers Dr. Jongkook Lee, Rob Van
Evra and Dr. Onyx Wai.
I would also like to mention some of the faculty and staff members who made a differ- ence during my stay at OSU including Dr. Robert Sykes, Dr. Bill Wolfe, Dr. Vince Ricca,
Dr. Ellen MacDonald and especially Ray Hunter. I would also like to thank my Disserta- tion Reading Committee members, Dr. Carolyn Merry and Dr. Diane Foster both of whom offered much support freely and enthusiastically.
I acknowledge the help of my HydroQual colleagues and friends including Jim Hallden,
Luca Liberti, Nicholas Kim, Dr. Pravi Shrestha, Dr. Alan Blumberg and many others. I would also like to thank my very good friend David Driscoll. I would like to thank my sisters Mary and Colleen and my father Pat who never doubted me.
Most importantly, I recognize my advisor Professor Keith W. Bedford; I will never forget him for his guidance and support.
v VITA
1987 ...... B.S. Physics, University of Minnesota, Minneapolis, MN 1993 ...... M.S. Civil Engineering, The Ohio State University, Columbus, OH 1999-2001 ...... Engineer, HydroQual, Inc., Mahwah, NJ
1991-1999,2001-present ...... Graduate Research and Teaching Asso- ciate, Civil and Environmental Engineer- ing and Geodetic Science, The Ohio State University, Columbus, OH
PUBLICATIONS
Wai, O. W.-H., Y. S. Xiong, S. O’Neil and K. W. Bedford (2001). “Parameter Estima- tion for Suspended Sediment Transport Processes”, The Science of the Total Environment, 226(1-3), 49-59.
O’Neil, S. and D. P. Podber (1997). “Sediment Transport Dynamics in a Dredged Tribu- tary,” Int. Conf. Estuarine and Coastal Modeling, eds. A. Blumberg and M. Spaulding, 5, 781-791.
O’Neil, S., K. W. Bedford and D. P. Podber (1996). “Storm-Derived Bar/Sill Dynamics in a Dredged Channel,” Proc. Int. Conf. Coastal Eng., ASCE, 25, 4289-4299.
Lee, J., S. O’Neil, K. W. Bedford and R. E. Van Evra (1994). “A Bottom Boundary Layer Sediment Response to Wave Groups,” Proc. Int. Conf. Coastal Eng., American Society of Civil Engineers, 24, 1827-1837.
Wai, O. W.-H., K. W. Bedford and S. O’Neil (1994). “Principal Components Time Spectra of Suspended Sediment in Random Waves,” Coastal Dynamics ’94, ASCE, Barcelona, 296-305.
vi Bedford, K. W., O. W.-H. Wai, S. O’Neil and M. Abdelrhman (1991). “Operational Pro- cedures for Estimating Bottom Exchange Rates,” in Hydraulic Engineering, Ed. R. Shane, American Society of Civil Engineers, 465-470.
Zhang, S., D. J. S. Welsh, K. W. Bedford, P. Sadayappan and S. O’Neil (1998). “Cou- pling of Circulation, Wave and Sediment Models,” Technical Report CEWES MSRC/PET TR/98-15. The Ohio State University, 32 pp.
Bedford, K. W., S. O’Neil, R. E. Van Evra and J. Lee (1994). “Ohio State University Measurements at SUPERTANK,” in SUPERTANK Laboratory Data Collection Project, Volume 1. Eds. Nicholas C. Kraus and Jane McKee Smith. U.S. Army Corps of Engineers, Waterways Experiment Station, Technical Report CERC-94-3, pp. 152-184.
O’Neil, S. (1993). Comparison of Sediment Transport Due to Monochromatic and Spec- trally Equivalent Random Waves. MS thesis, The Ohio State University, Columbus, Ohio.
Bedford, K. W., S. O’Neil, R. E. Van Evra and J. Lee (1993). “The Ohio State University Offshore ARMS Data - Boundary Layer, Entrainment and Resuspension: Overview plus Appendix,” Project Report, U.S. Army Corps of Engineers, Vicksburg, MS, 141 pp.
FIELDS OF STUDY
Major Field: Civil Engineering
Studies in: Models in Water Resources Engineering Sediment Transport Phenomena Prof. Keith W. Bedford Coastal Engineering
Applied Mathematics/Computational Science Profs. G. Baker, E. Overman Aerospace Engineering Prof. R. Bodonyi
vii TABLE OF CONTENTS
Page
Abstract ...... ii
Dedication ...... iv
Acknowledgments ...... v
Vita...... vi
List of Tables ...... xi
List of Figures ...... xii
Chapters:
1. Introduction ...... 1
2. Sediment Consolidation - Theoretical and Numerical Models ...... 23
2.1 One-Dimensional, Large Strain, Self-Weight, Primary Consolidation . . 32 2.1.1 The Governing Equation ...... 32 2.1.2 Force Balance ...... 35 2.1.3 Material Equilibrium ...... 35 2.1.4 Governing Equation ...... 37 2.1.5 Boundary Conditions ...... 38 2.1.6 Initial Conditions ...... 41 2.2 Numerical Solution ...... 42 2.2.1 Finite Difference Method ...... 42 2.2.2 Newton’s Method Solution ...... 44 2.2.3 Boundary Condition Implementation ...... 46 2.2.4 Stresses, Pressures and Settlement ...... 47
viii 2.2.5 Stability, Consistency and Convergence ...... 48 2.3 Summary ...... 49
3. Sediment Consolidation - Model Applications ...... 51
3.1 Uniform, Single Layer Consolidation ...... 53 3.1.1 Townsend - Scenario A ...... 53 3.1.2 Cargill - Craney Island Dredged Fill Material ...... 57 3.1.3 Been and Sills - Combwich on Somerset Clay ...... 60 3.2 Time-varying Sediment Loading ...... 63 3.2.1 Multi-Deposited Sediment Loading ...... 64 3.2.2 Multi-Eroded Sediment (Un-)Loading ...... 66 3.2.3 Sequential Loading and Unloading of Sediment Layers ..... 68 3.3 Summary ...... 73
4. Three-Dimensional Sediment Transport ...... 75
4.1 Coupling Water Column Transport and Mobile Bed ...... 78 4.1.1 Hydrodynamics ...... 79 4.1.2 Bed Mass Conservation ...... 80 4.1.3 Suspended Sediment Transport ...... 84 4.2 Boundary Conditions ...... 85 4.3 Numerical Solution ...... 86 4.3.1 Additional Considerations ...... 88 4.4 Outline of Modifications to the Model ...... 88 4.4.1 Coupling Water Column, Bed Load and Bed Evolution ...... 91 4.5 Summary ...... 96
5. Source/Sink Terms, Auxiliary Relations and Other Parameterizations ...... 97
5.1 Sediment Grain Sizes ...... 99 5.2 Settling Velocity ...... 100 5.2.1 Non-cohesive Settling ...... 101 5.2.2 Cohesive Settling ...... 104 5.3 Deposition ...... 108 5.3.1 Non-cohesive Deposition ...... 109 5.3.2 Cohesive Deposition ...... 111 5.3.3 Deposition Computation ...... 112 5.4 Erosion ...... 112 5.4.1 Non-cohesive Sediment Erosion ...... 112 5.4.2 Cohesive Sediment Erosion ...... 122 5.4.3 Bed Shear Strength Modeling ...... 126
ix 5.4.4 Unconsolidation, Fluidization and Liquefaction ...... 134 5.4.5 Erosion Computation ...... 138 5.5 Non-homogeneous Mixtures ...... 139 5.6 The Active Layer ...... 141 5.7 Other Important Parameterizations ...... 143 5.7.1 Near-bed Concentration ...... 144 5.7.2 Diffusion Coefficients ...... 145 5.7.3 Bed Roughness ...... 146 5.8 Summary ...... 147
6. Modeling Applications ...... 149
6.1 Transport in an Expansion Region ...... 149 6.1.1 Hydrodynamics of the Expansion Region ...... 152 6.1.2 Single Grain Size Sediment Transport ...... 160 6.1.3 Two Grain Size Sediment Transport ...... 171 6.1.4 Three Grain Size Sediment Transport ...... 187 6.2 Summary ...... 216
7. Conclusions ...... 219
7.1 Summary ...... 219 7.2 Conclusions ...... 220 7.3 Future Work ...... 223
Appendices:
A. Third Order Polynomial Interpolation ...... 225
B. Sediment-Water Relations ...... 229
B.1 Sediment-Water Mixture Density ...... 231
Bibliography ...... 235
x LIST OF TABLES
Table Page
1.1 Sediment transport models and their characteristics: model acronym, hy- drodynamic forcing (external, NF-no feedback, F-feedback, FE-finite el- ement, FD-finite difference), and type of sediment that can be modeled (NC-non-cohesive, C-cohesive)...... 17
1.2 Sediment transport models and their characteristics: water column trans- port of suspended sediment (AD-advection-diffusion, NA-not appropriate) and bedload flux formulation...... 17
1.3 Sediment transport models and their characteristics: erosion formulation (AP-an Ariathurai-Partheniades excess stress formula), deposition and bed model formulations ...... 19
2.1 Consistent sets of input units for consolidation model...... 49
5.1 Fine-grained sediment sizes and classes (after Vanoni, 1975)...... 99
xi LIST OF FIGURES
Figure Page 2.1 Lagrangian coordinates ( ) and convective coordinates ( )...... 33
2.2 Soil volume elements at various times. The left-most volume represents the Lagrangian coordinate, the middle volume the convective coordinate and the right-most volume is the material coordinate...... 34
3.1 Constitutive relations for Townsend’s Scenario A material (Townsend and McVay, 1990)...... 54
3.2 Density profiles in time - Townsend’s Scenario A material...... 55
3.3 Contour of height with density and time - Townsend’s Scenario A material. 56
3.4 Constitutive relations for the Craney Island dredged fill material (Cargill, 1982)...... 57
3.5 Density profiles in time - Craney Island dredged fill material...... 58
3.6 Contour of height with density and time - Craney Island dredged fill material. 59
3.7 Constitutive relations for the Combwich on Somerset clay material (Been and Sills, 1981)...... 61
3.8 Density profiles in time - Combwich on Somerset clay material...... 63
3.9 Contour of height with density and time - Combwich on Somerset clay material...... 64
3.10 Density profiles in time. Multiply deposited layers for the Been and Sills (1981) case. Base of 10 cm with the addition of 10 cm, followed by an additional deposition of 12 cm...... 65
xii 3.11 Contours height with density and time. Multi-deposited layers for the Been and Sills (1981) case. Base of 10 cm with the addition of 10 cm, followed by an additional deposition of 12 cm...... 67
3.12 Density profiles in time. Multiply eroded layers for the Been and Sills (1981) case. Base of 20 cm with the erosion of 5 cm, followed by an additional erosion of 8 cm...... 68
3.13 Variation in soil layer height with time. Multiply eroded layers for the Been and Sills (1981) case. Base of 20 cm with the erosion of 5 cm, followed by an additional erosion of 8 cm...... 69
3.14 Variation in soil layer height with time. Sequential loading for the Been and Sills (1981) case. Base of 10 cm with the addition of 10 cm, followed by erosion of 8 cm...... 70
3.15 Density profiles with time. Sequential loading for the Been and Sills (1981) case. Base of 10 cm with the addition of 10 cm, followed by erosion of 8 cm. 71
3.16 Variation in soil layer height with time. Sequential loading for the Been and Sills (1981) case. Base of 20 cm with the erosion of 5 cm, followed by deposition of 8 cm...... 72
3.17 Density profiles with time. Sequential loading for the Been and Sills (1981) case. Base of 20 cm with the erosion of 5 cm, followed by deposition of 8 cm...... 73
4.1 Block diagram of MBDC operations. The smain module controls the suspended, bedload, mobile bed and consolidation...... 90
4.2 Block diagram of MBDC operations. The initialization module is called from the smain routine. Blue blocks represent newly added modules. Red blocks have been altered from the original mobile bed routines. . . . . 91
4.3 Block diagram of MBDC operations. The sedcom module controls the interaction of suspended and bed sediment operations with source terms. Blue blocks represent newly added modules. Red blocks have been altered from the original mobile bed routines...... 93
xiii 4.4 Block diagram of MBDC operations. The bedsed module controls the
bed computations and the iterative loop for solving the bed level and coupled system as well as the call to the consolidation solver. Blue blocks represent newly added modules. Red blocks have been altered from the original mobile bed routines...... 94
5.1 Comparison of continuous and piecewise-continuous formulas for deter- mining settling velocity...... 103
5.2 Comparison of single continuous and piecewise-continuous formulas for determining the critical Shields parameter...... 118
6.1 Horizontal segment of expansion channel grid. Upstream is left, down- stream is right...... 150
6.2 Elevation at run time of 5 minutes in the expansion channel...... 154
6.3 Elevation at run time of 15 minutes in the expansion channel...... 155
6.4 Temperature at run time of 5 minutes in the expansion channel. The veloc- ity scale for each panel is local and is 1 m/s...... 157
6.5 Temperature at run time of 25 minutes in the expansion channel. The ve- locity scale for each panel is local and is 1 m/s...... 158
6.6 Concentration of 726 m particles at a run time of 15 minutes in the ex- pansion channel. The velocity scale for each panel is local and is 1 m/s. . . 161
6.7 Suspended sediment concentration (mg/L) for the single grain case (726