Eddy Interaction Model for Turbulent Suspension in Reynolds-Averaged

Total Page:16

File Type:pdf, Size:1020Kb

Eddy Interaction Model for Turbulent Suspension in Reynolds-Averaged Advances in Water Resources 111 (2018) 435–451 Contents lists available at ScienceDirect Advances in Water Resources journal homepage: www.elsevier.com/locate/advwatres Eddy interaction model for turbulent suspension in Reynolds-averaged T Euler–Lagrange simulations of steady sheet flow ⁎ Zhen Cheng ,1,a, Julien Chauchatb, Tian-Jian Hsua, Joseph Calantonic a Civil and Environmental Engineering, University of Delaware, Newark, DE 19716, USA b Laboratory of Geophysical and Industrial Flows (LEGI), BP 53, Grenoble Cedex 9 38041, France c Marine Geosciences Division, U.S. Naval Research Laboratory, Stennis Space Center, MS 39529, USA ARTICLE INFO ABSTRACT Keywords: A Reynolds-averaged Euler–Lagrange sediment transport model (CFDEM-EIM) was developed for steady sheet Euler–Lagrange model flow, where the inter-granular interactions were resolved and the flow turbulence was modeled with a low Eddy interaction model Reynolds number corrected kω− turbulence closure modified for two-phase flows. To model the effect of Turbulent suspension turbulence on the sediment suspension, the interaction between the turbulent eddies and particles was simulated Steady sheet flow with an eddy interaction model (EIM). The EIM was first calibrated with measurements from dilute suspension Rouse profile experiments. We demonstrated that the eddy-interaction model was able to reproduce the well-known Rouse Sediment transport rate profile for suspended sediment concentration. The model results were found to be sensitive to the choice of the coefficient, C0, associated with the turbulence-sediment interaction time. A value C0 = 3 was suggested to match the measured concentration in the dilute suspension. The calibrated CFDEM-EIM was used to model a steady sheet flow experiment of lightweight coarse particles and yielded reasonable agreements with measured velo- city, concentration and turbulence kinetic energy profiles. Further numerical experiments for sheet flow sug- gested that when C0 was decreased to C0 < 3, the simulation under-predicted the amount of suspended sediment in the dilute region and the Schmidt number is over-predicted (Sc > 1.0). Additional simulations for a range of Shields parameters between 0.3 and 1.2 confirmed that CFDEM-EIM was capable of predicting sediment transport rates similar to empirical formulations. Based on the analysis of sediment transport rate and transport layer thickness, the EIM and the resulting suspended load were shown to be important when the fall parameter is less than 1.25. 1. Introduction The conventional single-phase-based sediment transport models assume the dynamics of transport can be subjectively separated into bedload Studying sediment transport in rivers and coastal regions is critical and suspended load (e.g., van Rijn, 1984a; 1984b). While the suspended to understand the fluvial geomorphology, loss of wetland, and beach load are directly resolved, the bedload are parameterized by empirical erosion. For example, significant engineering efforts were devoted to formulations. Several laboratory measurements of sheet flow with the control the river discharge and sediment budget to reduce the loss of full profile of sediment transport flux and net transport rate indicated Louisiana wetland (Allison et al., 2012; Mossa, 1996). In the Indian that the split of bedload and suspended load may be too simple because River inlet, significant erosion of the north beach is mitigated through sediment entrainment/deposition is a continuous and highly dynamic proper beach nourishment that interacts with littoral drift process near the mobile bed (e.g., O’Donoghue and Wright, 2004; Revil- (Keshtpoor et al., 2013). The characteristics of sediment transport vary Baudard et al., 2015). In sheet flow, the two prevailing mechanisms significantly with sediment properties and flow conditions, and it is driving the sediment transport are inter-granular interactions and tur- widely believed that sheet flow plays a dominant role in nearshore bulent suspension (Jenkins and Hanes, 1998; Revil-Baudard et al., beach erosion and riverine sediment delivery, especially during storm 2015). In order to model the full profile of sediment transport, both and flood conditions, respectively. mechanisms must be taken into account. In the past decade, many Sheet flow is an intense sediment transport mode, in which a thick Eulerian two-phase flow models have been developed for sheet flow layer of concentrated sediment is mobilized above the quasi-static bed. transport in steady (Jenkins and Hanes, 1998; Longo, 2005; Revil- ⁎ Corresponding author. E-mail address: [email protected] (Z. Cheng). 1 Now at Applied Ocean Physics & Engineering, Woods Hole Oceanographic Institution, Woods Hole, MA 02543, USA. https://doi.org/10.1016/j.advwatres.2017.11.019 Received 15 April 2017; Received in revised form 31 October 2017; Accepted 19 November 2017 Available online 21 November 2017 0309-1708/ © 2017 Elsevier Ltd. All rights reserved. Z. Cheng et al. Advances in Water Resources 111 (2018) 435–451 Baudard and Chauchat, 2013) and oscillatory flows (Dong and Zhang, a Lagrangian approach. 2002; Hsu et al., 2004), Amoudry et al.,(Chen et al., 2011; Cheng et al., In the stochastic Lagrange model for particle dispersion, the tur- 2017a; Liu and Sato, 2006). By solving the mass and momentum bulent agitation to the sediment particles are considered either through equations of fluid phase and sediment phase with appropriate closures a random-walk model (RWM) or an eddy interaction model (EIM). In for interphase momentum transfer, turbulence, and intergranular the RWM framework, the strength of particle velocity fluctuations are stresses, these models are able to resolve the entire profiles of sediment typically assumed to be similar to the fluid turbulence, and a series of transport without the assumption of bedload and suspended load. random velocity fluctuations are directly added to the particle velo- In the continuum description of the sediment phase, the assumption cities. While the Lagrange model with RWM is successful in studying of uniform particle properties and spherical particle shapes are usually the particle dispersion in mixing layer (Coimbra et al., 1998) and dilute adopted. To better capture the polydisperse nature of sediment trans- suspension (Shi and Yu, 2015), the assumption of estimating the par- port and irregular particle shapes, the Lagrangian approach for the ticle velocity fluctuations based on the fluid turbulence is crucial, and particle phase, namely the Discrete Element Method (DEM, Cundall and many researchers found that the correlation between the particle and Strack, 1979; Maurin et al., 2015; Sun and Xiao, 2016a) is superior to fluid fluctuations are highly dependent on the particle Stokes number, the Eulerian approach because individual particle properties may be St= tpl/ t (Balachandar and Eaton, 2010), where tp is the particle re- uniquely specified (Calantoni et al., 2004; Fukuoka et al., 2014; Harada sponse time, and tl is the characteristic time scale of energetic eddies. and Gotoh, 2008). One of the main challenges in modeling sheet flow For the particles with very small inertia (St ≪ 1), they can closely follow arise from the various length scales involved in inter-granular interac- the eddy motion. However, if St ≫ 1, the particle trajectory is hardly tions and sediment-turbulence interactions. To resolve the flow turbu- affected by the fluid eddy motion. Due to the particle inertia effect, it lence and turbulence-sediment interactions in sheet flow, the compu- was found that the fluid turbulent intensity needs to be enhanced for tational domain needs to be sufficiently large to resolve the largest medium to coarse particles (Shi and Yu, 2015). This problem can be eddies, while the grid resolution should be small enough to resolve the largely remedied by the EIM (Matida et al., 2004), where the fluid energy containing turbulent eddies. This constrain becomes even more velocity fluctuations associated with the fluid turbulence are added challenging in the Euler–Lagrange modeling framework. Large domains through the particle-sediment interaction force, i.e., the drag force. This require both a large number of grid points to resolve a sufficient approach incorporates the particle inertia effect naturally and it is ap- amount of turbulence energy cascade (i.e., large-eddy simulation) and a plicable for a wide range of sediment properties. Graham (1996) found large number of particles in a given simulation (e.g., Finn et al., 2016). that the dispersion of inertial particles may be correctly represented It is well-established that in sheet flow, the transport layer thickness with a suitable choice of maximum interaction time and length scales scales with the grain size and the Shields parameter (Wilson, 1987), with the eddies. This model was later improved by using a randomly suggesting that a common sheet flow layer thickness must be about sampled eddy interaction time, in which more realistic turbulent scales several tens of grain diameters. To simulate the largest eddies and their become possible, and the enhanced dispersion of high-inertia particles subsequent cascade, the domain lengths in the two horizontal directions are captured. In the previous studies of particle dispersion (e.g., Shi and must be proportional to the boundary layer thickness, which is usually Yu, 2015), the turbulent intensity is either prescribed from the em- about several tens of centimeters. For a bed layer thickness of 50 grains pirical formula, or modeled using clear fluid turbulence closure without with
Recommended publications
  • Particle Size Distribution Controls the Threshold Between Net Sediment
    PUBLICATIONS Geophysical Research Letters RESEARCH LETTER Particle Size Distribution Controls the Threshold Between 10.1002/2017GL076489 Net Sediment Erosion and Deposition in Suspended Key Points: Load Dominated Flows • New flow power sediment entrainment model scales with flow R. M. Dorrell1 , L. A. Amy2, J. Peakall3, and W. D. McCaffrey3 depth instead of particle size and incorporates capacity and 1Energy and Environment Institute, University of Hull, Hull, UK, 2School of Earth Sciences, University College Dublin, Dublin, competence limits 3 • Suspended load, at the threshold Ireland, School of Earth and Environment, University of Leeds, Leeds, UK between net erosion and deposition, is shown to be controlled by the particle size distribution Abstract The central problem of describing most environmental and industrial flows is predicting when • Flow power model offers an order of material is entrained into, or deposited from, suspension. The threshold between erosional and magnitude, or better, improvement in fl predicting the erosional-depositional depositional ow has previously been modeled in terms of the volumetric amount of material transported in threshold suspension. Here a new model of the threshold is proposed, which incorporates (i) volumetric and particle size limits on a flow’s ability to transport material in suspension, (ii) particle size distribution effects, and (iii) a fi Supporting Information: new particle entrainment function, where erosion is de ned in terms of the power used to lift mass from the • Table S1 bed. While current suspended load transport models commonly use a single characteristic particle size, • Supporting Information SI the model developed herein demonstrates that particle size distribution is a critical control on the threshold between erosional and depositional flow.
    [Show full text]
  • Delft Hydraulics
    Prepared for: DG Rijkswaterstaat, Rijksinstituut voor Kust en Zee | RIKZ Description of TRANSPOR2004 and Implementation in Delft3D-ONLINE FINAL REPORT Report November 2004 Z3748.00 WL | delft hydraulics Prepared for: DG Rijkswaterstaat, Rijksinstituut voor Kust en Zee | RIKZ Description of TRANSPOR2004 and Implementation in Delft3D-ONLINE FINAL REPORT L.C. van Rijn, D.J.R. Walstra and M. van Ormondt Report Z3748.10 WL | delft hydraulics CLIENT: DG Rijkswaterstaat Rijks-Instituut voor Kust en Zee | RIKZ TITLE: Description of TRANSPOR2004 and Implementation in Delft3D-ONLINE, FINAL REPORT ABSTRACT: In 2003 much effort has been spent in the improvement of the DELFT3D-ONLINE model based on the engineering sand transport formulations of the TRANSPOR2000 model (TR2000). This work has been described in Delft Hydraulics Report Z3624 by Van Rijn and Walstra (2003). However, the engineering sand transport model TR2000 has recently been updated into the TR2004 model within the EU-SANDPIT project. The most important improvements involve the refinement of the predictors for the bed roughness and the suspended sediment size. Up to now these parameters had to be specified by the user of the models. As a consequence of the use of predictors for bed roughness and suspended sediment size, it was necessary to recalibrate the reference concentration of the suspended sediment concentration profile. The formulations (including the newly derived formulations of TR2004) implemented in this 3D-model are described in detail. The implementation of TR2004 in Delft3D-ONLINE is part of an update of Delft3D which involves among others: the extension of the model to be run in profile mode, the synchronisation of the roughness formulations and inclusion of two breaker delay concepts.
    [Show full text]
  • Formulae of Sediment Transport in Unsteady Flows (Part 2) Shu-Qing Yang
    Chapter Formulae of Sediment Transport in Unsteady Flows (Part 2) Shu-Qing Yang Abstract Sediment transport (ST) in unsteady flows is a complex phenomenon that the existing formulae are often invalid to predict. Almost all existing ST formulae assume that sediment transport can be fully determined by parameters in streamwise direc- tion without parameters in vertical direction. Different from this assumption, this paper highlights the importance of vertical motion and the vertical velocity is suggested to represent the vertical motion. A connection between unsteadiness and vertical velocity is established. New formulae in unsteady flows have been developed from inception of sediment motion, sediment discharge to suspension’s Rouse num- ber. It is found that upward vertical velocity plays an important role for sediment transport, its temporal and spatial alternations are responsible for the phase lag phenomenon and bedform formation. Reasonable agreement between the measured and the proposed conceptual model was achieved. Keywords: unsteady flow, shields number, rouse number, sediment transport, vertical velocity 1. Introdution Sediment transport is the movement of solid particles driven by fluid like water or wind in rivers, lakes, reservoirs, coastal waters. Generally, in the real world the flow is unsteady like flood waves, tidal waves and wind waves, because steady and uniform flows are very rare in reality. Even so, it is understandable that sediment transport is first observed under well controlled conditions in laboratory, and then the data are collected to calibrate the models. These formulae are further examined using field data by assuming the laboratory flow conditions (generally steady and uniform flows) can be extended to rivers and coastal waters (generally unsteady and non-uniform).
    [Show full text]
  • Two-Phase Modeling of Granular Sediment for Sheet Flows and Scour
    TWO-PHASE MODELING OF GRANULAR SEDIMENT FOR SHEET FLOWS AND SCOUR. A Dissertation Presented to the Faculty of the Graduate School of Cornell University in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy by Laurent Olivier Amoudry January 2008 c 2008 Laurent Olivier Amoudry ALL RIGHTS RESERVED TWO-PHASE MODELING OF GRANULAR SEDIMENT FOR SHEET FLOWS AND SCOUR. Laurent Olivier Amoudry, Ph.D. Cornell University 2008 Even though sediment transport has been studied extensively in the past decades, not all physical processes involved are yet well understood and repre- sented. This results in a modeling deficiency in that few models include complete and detailed descriptions of the necessary physical processes and in that mod- els that do usually focus on the specific case of sheet flows. We seek to address this modeling issue by developing a model that would describe appropriately the physics and would not only focus on sheet flows. To that end, we employ a two-phase approach, for which concentration- weighted averaged equations of motion are solved for a sediment and a fluid phase. The two phases are assumed to only interact through drag forces. The correlations between fluctuating quantities are modeled using the turbulent viscosity and the gradient diffusion hypotheses. The fluid turbulent stresses are calculated using a modified k−ε model that accounts for the two-way particle-turbulence interaction, and the sediment stresses are calculated using a collisional granular flow theory. This approach is used to study three different problems: dilute flow modeling, sheet flows and scouring. In dilute models sediment stresses are neglected.
    [Show full text]
  • A Two-Phase Flow Model of Sediment Transport: Transition from Bedload To
    J. Fluid Mech. (2014), vol. 755, pp. 561–581. c Cambridge University Press 2014 561 doi:10.1017/jfm.2014.422 A two-phase flow model of sediment transport: transition from bedload to suspended load Filippo Chiodi1, Philippe Claudin1 and Bruno Andreotti1, † 1Laboratoire de Physique et Mécanique des Milieux Hétérogènes, (PMMH UMR 7636 ESPCI – CNRS – University Paris Diderot – University P. M. Curie) 10 rue Vauquelin, 75005 Paris, France (Received 7 December 2013; revised 2 May 2014; accepted 21 July 2014) The transport of dense particles by a turbulent flow depends on two dimensionless numbers. Depending on the ratio of the shear velocity of the flow to the settling velocity of the particles (or the Rouse number), sediment transport takes place in a thin layer localized at the surface of the sediment bed (bedload) or over the whole water depth (suspended load). Moreover, depending on the sedimentation Reynolds number, the bedload layer is embedded in the viscous sublayer or is larger. We propose here a two-phase flow model able to describe both viscous and turbulent shear flows. Particle migration is described as resulting from normal stresses, but is limited by turbulent mixing and shear-induced diffusion of particles. Using this framework, we theoretically investigate the transition between bedload and suspended load. Key words: alluvial dynamics, multiphase and particle-laden flows, particle/fluid flow 1. Introduction When a sedimentary bed is sheared by a water flow of sufficient strength, the particles are entrained into motion. In a homogeneous and steady situation, the fluid flow can be characterized by a unique quantity: the shear velocity u .
    [Show full text]
  • Suspended Sediment Concentration Profiles in the Bottom Boundary Layer
    W&M ScholarWorks Dissertations, Theses, and Masters Projects Theses, Dissertations, & Master Projects 1993 Suspended sediment concentration profiles in the bottom boundary layer Jingping Xu College of William and Mary - Virginia Institute of Marine Science Follow this and additional works at: https://scholarworks.wm.edu/etd Part of the Geology Commons, and the Oceanography Commons Recommended Citation Xu, Jingping, "Suspended sediment concentration profiles in the bottom boundary layer" (1993). Dissertations, Theses, and Masters Projects. Paper 1539616914. https://dx.doi.org/doi:10.25773/v5-55tq-zj21 This Dissertation is brought to you for free and open access by the Theses, Dissertations, & Master Projects at W&M ScholarWorks. It has been accepted for inclusion in Dissertations, Theses, and Masters Projects by an authorized administrator of W&M ScholarWorks. For more information, please contact [email protected]. INFORMATION TO USERS This manuscript has been reproduced from the microfilm master. UMI films the text directly from the original or copy submitted. Thus, some thesis and dissertation copies are in typewriter face, while others may be from any type of computer printer. The quality of this reproduction is dependent upon the quality of the copy submitted. Broken or indistinct print, colored or poor quality illustrations and photographs, print bleedthrough, substandard margins, and improper alignment can adversely afreet reproduction. In the unlikely event that the author did not send UMI a complete manuscript and there are missing pages, these will be noted. Also, if unauthorized copyright material had to be removed, a note will indicate the deletion. Oversize materials (e.g., maps, drawings, charts) are reproduced by sectioning the original, beginning at the upper left-hand corner and continuing from left to right in equal sections with small overlaps.
    [Show full text]
  • Design of a Sediment Mitigation System for Conowingo Dam Presenters: Sheri Gravette Kevin Cazenas Said Masoud Rayhan Ain Sediment Plume
    Design of a Sediment Mitigation System for Conowingo Dam Presenters: Sheri Gravette Kevin Cazenas Said Masoud Rayhan Ain Sediment Plume Sponsor: Faculty Advisor: Lower George Donohue Susquehanna Riverkeeper Agenda • Context • Stakeholders • Problem/Need Statement • Mission Requirements • Design Alternatives • Technical Approach • Preliminary Results • Project Management 2 Chesapeake Bay and The Susquehanna River • Chesapeake Bay is the largest estuary in the United States • 3 largest tributaries of the Bay are the Susquehanna, Potomac and James rivers – Provide more than 80% of the Bay’s freshwater • Susquehanna River is the Bay’s largest tributary – Provides nearly 50% of freshwater to the Bay – Flows from NY to PA to MD Map of the Chesapeake Bay Watershed Source: The PA Dept. of Environmental Protection 3 Lower Susquehanna River and Water Quality • Flows through Pennsylvania and Maryland • Quality of water is vital to the bay’s health – Improvement in water quality thus far can be attributed to US Army Corp. of Engineers • Provides power for turbines in hydroelectric plants and clean water to people • Contains 4 Dams: York Haven, Safe Harbor, Holtwood, Conowingo – York Haven, Safe Harbor and Holtwood have reached steady state - dam has completely silted up and is no longer able to retain sediment; dams are at maximum capacity Map of Conowingo Reservoir Source: US Army Corps of Engineers, (2013) 4 Conowingo Dam • Constructed in 1928 • Southernmost Dam of the Lower Susquehanna sediment • Location of Conowingo Hydroelectric Station – Mainly provides power to Philadelphia, PA – A black start power source – Provides 1.6 billion kWh annually • Traps sediment and nutrients from reaching the Chesapeake Bay – Water quality is closely related to sediment Conowingo Dam deposition Source: J.
    [Show full text]
  • Chapter 10 Movement of Sediment by Water Flows
    CHAPTER 10 MOVEMENT OF SEDIMENT BY WATER FLOWS INTRODUCTION 1 A simple flume experiment on sediment movement by a unidirectional current of water in a flume serves to introduce the material in this chapter. Place a layer of sediment in the flume, level it to have a planar surface, and establish a uniform flow at a certain depth and velocity. Gradually, in steps, increase the strength of the flow beyond the condition for incipient movement. The magnitude of the flow strength relative to what is required for incipient movement of the bed sediment is conventionally called the flow intensity, and is usually taken to be the ratio τo/τoc (or, what is the same, u*/u*c), where the subscript c denotes the threshold (“critical”) condition. 2 At first the particles move as bed load, by hopping, rolling, and/or sliding. Particle movement is neither continuous nor uniform over the bed: brief gusts or pulses of movement affect groups of particles locally, and seemingly randomly, on the bed. Particles move a short distance, stop, and then move again. Even when they are moving, they are generally not moving as fast as the fluid near the bed surface. 3 As the flow becomes stronger, some of the particles moving near the bed are lifted upward by upward-moving turbulent eddies and travel for more or less long distances downstream as suspended load. The stronger the flow and/or the finer the sediment, the greater is the concentration of suspended sediment, the higher it can travel in the flow, and the longer it moves downstream before returning to the bed.
    [Show full text]
  • Chapter 29 Sediment Transport
    Chapter 29 Sediment transport 29.1 Introduction The following seven test cases have been developed to illustrate the various options implemented in the sediment transport module bedload Compares the different formulations for bed load transport by im- posed currents (bottom stress) values (see Section 7.4). totload Compares the different formulations for total load transport by im- posed current (bottom stress) values (see Section 7.5). wavload Compares formulations for bed and total load transport with wave effects using imposed values for currents and wave parameters (wave height, period and direction). sedvprof Simulates the diffusion and settling of sediments in a water column by a bottom stress imposed via an external pressure gradient (surface slope). sedhprof This test case, taken from Hjelmfelt & Lenau (1970), simulates the transition between a non-eroding and an eroding bed. This is modeled by applying a boundary condition with a concentration of zero at the upstream boundary. Note that the velocity field is supplied as an initial condition, rather than calculated by COHERENS. Both 2-D and 3-D simulations are performed and different model settings are applied to simulate the convergence to an equilibrium state. seddens Simulates the evolution of an initial sediment distribution in a closed channel, consisting of vertically uniform concentrations in a sec- tion at the left of the channel and zero concentrations to the right. 793 794 CHAPTER 29. SEDIMENT TRANSPORT The sediment is allowed to settle due to is own weight. A current is generated by the horizontal concentration gradient. thacker This test case, after Thacker (1981), simulates the oscillations in a rotating parabolic basin with moving boundaries (shores lines).
    [Show full text]
  • Numerical Simulation of the Inception of Channel Meandering
    Earth Surface Processes and Landforms EarthSimulation Surf. Process. of the inceptionLandforms of30 channel, 1093–1110 meandering (2005) 1093 Published online in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/esp.1264 Numerical simulation of the inception of channel meandering Jennifer G. Duan1* and Pierre Y. Julien2 1 Division of Hydrologic Sciences, Desert Research Institute, 755 E. Flamingo Road, Las Vegas, Nevada 89119, USA 2 Engineering Research Centre, Colorado State University, Fort Collins, CO 80523, USA *Correspondence to: J. G. Duan, Abstract Division of Hydrologic Sciences, Desert Research Institute, The inception of channel meandering is the result of the complex interaction between flow, University and College System of bed sediment, and bank material. A depth-averaged two-dimensional hydrodynamic model Nevada, 755 E. Flamingo Road, is developed to simulate the inception and development of channel meandering processes. Las Vegas, NV 89119, USA. The sediment transport model calculates both bedload and suspended load assuming equi- E-mail: [email protected] librium sediment transport. Bank erosion consists of two interactive processes: basal erosion and bank failure. Basal erosion is calculated from a newly derived equation for the entrain- ment of sediment particles by hydrodynamic forces. The mass conservation equation, where basal erosion and bank failure are considered source terms, was solved to obtain the rate of bank erosion. The parallel bank failure model was tested with the laboratory experiments of Friedkin on the initiation and evolution processes of non-cohesive meandering channels. The model replicates the downstream translation and lateral extension of meandering loops reasonably well. Plots of meandering planforms illustrate the evolution of sand bars and redistribution of flow momentum in meandering channels.
    [Show full text]
  • Universit`A Degli Studi Di Trieste Large Eddy Simulation of Suspended
    Universita` degli Studi di Trieste XXVII Ciclo del Dottorato di Ricerca in Environmental and Industrial Fluid Mechanics Large Eddy Simulation of Suspended Sediment in Turbulent Open-Channel Flow Settore scientifico-disciplinare: ICAR/01 Ph.D. Student: Maryam Dallali Supervisor and Director: Prof. Vincenzo Armenio Co-Supervisor: Prof. Ana Maria Ferreira da Silva Co-Supervisor: Prof. Andrew Pollard Academic Year 2014-2015 ii Abstract This study aims to shed light on the dark sides of sediment transport and meandering channel flow patterns. For this purpose, the Large Eddy Simula- tion methodology is used to investigate sediment transport in open-channel flow and also to predict the meandering channel flow pattern in two sine- generated channels. In the turbulent flows the transported suspended sediments is usually assumed in a passive form. In other word, the feedback effect of sediment suspension on near-wall momentum transport is neglectable. The first phase of this research is devoted to quantify the just mentioned effect in the open channel flow. A single phase approach is used to gain better description of sediment-turbulence interaction. The effect of suspended particles on flow properties is presented through a buoyancy term in momentum equations. The initial analysis of concentration profiles shows that the combina- tion of large near-wall concentration and small deposition term, enhances the suspension concentration in case of small particles. Considering the buoyancy effect (two-way coupling model), the suspension concentration is reduced compared to the passive case. In addition, the buoyancy effects on the velocity field are also discussed. The comparison of streamwise mean velocity profiles reveals that the interaction of suspended particles with the underlying flow suppresses the vertical turbulent fluctuations and conse- quently reduces the von K`arm`anconstant (κ).
    [Show full text]
  • Numerical Study of Particle Bed Scour by Vortices Daniel S
    University of Vermont ScholarWorks @ UVM Graduate College Dissertations and Theses Dissertations and Theses 2014 Numerical study of particle bed scour by vortices Daniel S. Hagan University of Vermont Follow this and additional works at: https://scholarworks.uvm.edu/graddis Part of the Mechanical Engineering Commons Recommended Citation Hagan, Daniel S., "Numerical study of particle bed scour by vortices" (2014). Graduate College Dissertations and Theses. 267. https://scholarworks.uvm.edu/graddis/267 This Thesis is brought to you for free and open access by the Dissertations and Theses at ScholarWorks @ UVM. It has been accepted for inclusion in Graduate College Dissertations and Theses by an authorized administrator of ScholarWorks @ UVM. For more information, please contact [email protected]. NUMERICAL STUDY OF PARTICLE BED SCOUR BY VORTICES A Thesis Presented by Dan Hagan to The Faculty of the Graduate College of The University of Vermont In Partial Fullfillment of the Requirements for the Degree of Master of Science Specializing in Mechanical Engineering October, 2014 Accepted by the Faculty of the Graduate College, The University of Vermont, in partial fulfillment of the requirements for the degree of Master of Science, specializing in Mechanical Engineering. Thesis Examination Committee: Advisor Yves Dubief, Ph.D. Advisor Mandar Dewoolkar, Ph.D. Darren Hitt, Ph.D. Chairperson George Pinder, Ph.D. Dean, Graduate College Cynthia J. Forehand, Ph.D. Date: July 2, 2014 Abstract Scouring is the process of soil or sediment erosion due to flowing water, which can lead to bed degradation and compromised transportation infrastructure. In the decade before 2000, over half of the 500 bridge failures in the United States were caused by flooding or scouring.
    [Show full text]