Testing various constitutive equations for debris flow modelling
|
|
- Julius Floyd
- 5 years ago
- Views:
Transcription
1 Hydrology, Water Resources and Ecology in Headwaters (Proceedings of the HeadWater'98 Conference held at Meran/Merano, Italy, April 1998). IAHS Publ. no. 48, Testing various constitutive equations for debris flow modelling THILO KOCH Swiss Federal Institute of Snow and Landscape Research, CH-8903 Birmensdorf, Switzerland Abstract On the basis of a curvilinear coordinate system and incorporation of a number of flow resistance laws, a one-dimensional model has been developed to simulate debris flows in mountain torrent channels. The model simulations are compared to debris flow events in the Kamikamihori Valley in Japan. The results suggest that the models with a Newtonian turbulent and a Voellmy fluid flow resistance, as well as a model with a new empirical flow resistance approach, are better at reproducing the observed velocities and flow-heights than the models with a Bingham, Newtonian laminar or dilatant grain-shearing flow resistance. The material constant for the dynamic part in the flow resistance law, the channel width and the total volume were found to be the most sensitive parameters for the flow simulation. INTRODUCTION Debris flows form one type of natural hazard often occurring in mountain regions. A debris flow is a rapidly flowing mixture of grains of various sizes and water initiated on a steep slope, usually travelling along torrents and coming to rest at a debris flow fan. A debris flow implies fluid-to-fluid, fluid-to-solid and solid-to-solid interactions, which are generally not well understood and not considered in most of the existing approaches to model debris flows. Most of the models developed have been restricted to the one-dimensional case, with mainly three types of constitutive equations: Bingham flow, Newtonian turbulent flow, or dilatant inertial. A more recent overview of constitutive laws to describe the rheology of debris flows is given by Jan & Shen (1993). The objective of this study is to evaluate a number of different model approaches by applying them to documented events and comparing simulated and observed flow behaviour. SIMPLE DEBRIS FLOW MODELS The debris flow mixture is described as a one-dimensional incompressible continuum. The shallow water approximation, averaging different constitutive equations and the mass and momentum balance laws, has been used. The pressuredependent deformation is assumed to comply with a Mohr-Coulomb plastic yield criterion. The model is formulated for a curvilinear coordinate system, following the torrent topography, where s is the downward coordinate. The depth-averaged
2 50 Thilo Koch equations of mass and momentum are described in a dimensionless form for a Coulomb-friction flow in Savage & Hutter (1991). A formulation in a dimensional form is presented by Hungr (1995) for several constitutive equations. I have followed Hungr (1995) to a larger extent. The numerical procedure used here is a 1-D version of a Lagrangian central finite difference scheme, completely described in Hutter & Koch (1991). The mass is directly connected with n cells (Fig. 1). The index i corresponds to the cell boundary points and j to the cell centres; k counts for the time. The calculations start with an initial geometry and velocity distribution for the mass using the geometric data (bed slope and width) of the torrent profile. The boundaries with heights Hf' x (at time k - 1) are connected with the depth-averaged velocities v*~ l/ (at time k- Vz, second-order time accuracy). As a the first step, we obtain the new coordinates s* from the previous coordinates sf~ ] and velocity v*" 1/ **=**-'+v,*- At where At is the time step. Second, the continuity equation is solved according to (D H B"+> V () where H) is an average cell height, ds* =sf +l -sf the cell length, B'. the local width of the cell boundary and F = V\ = constant, the volume of cell j. The resulting force, F.*, connected to boundary point i is composed of the tangential components of weight G*, the resisting force Tf and the pressure-distribution dependent component P t k F," = G: + p; with rpk G* =y,dsf//*,* since* (3) (4) Fig. 1 Definition of mesh cell notation for the numerical scheme used in the 1-D model.
3 Testing various constitutive equations for debrisflowmodelling 51 where y, = bulk unit weight, d sf = length, constructed from the neighbouring cells with length ds*, and a* = bed slope angle. The pressure term P k s acllpa: «sdh, -Y,M ds i-dsïh'b* sin a* (5) is different from Hungr (1995) and is in this implementation calculated by a Mohr- Coulomb plastic yield criterion (Hutter & Koch, 1991). The earth pressure coefficient with. \ act I pass IK. dv, K acl for ->0 ds-, K" ass îox -<0 ds, (6) act I pass K ss = 1 + ^/1-cos <j>/cos 5 cos -1 (V) Here, S = bed friction angle and (j> = internal angle of friction. This law has been successfully tested for dry granular flows, e.g. Koch et al. (1994). I expect it to be useful in debris flows when the granular particles stay in contact with each other. The resisting forces 7]* can be mainly written in the general, depth-averaged form ^=y,bfd,f[ Kf(vff + x, + (x c ); (8) with (t c )* =Hf cosaf tan ; (9) where = material constant, i y = cohesion and x c = coulomb yield strength. (3 and A, are exponents characterizing the material laws together with the existence of % y and In the case of a Bingham fluid the following cubic equation has to be solved in rpk. T} k Hf "m 3u,5>; It + {T?) 3 (10) where x B = Bingham yield stress and u. B = Bingham viscosity. The following equation is an explicit and, without f, an often-used formulation of v; v = CH K S B +f(x v,t c ) (11) Equation (11) is equivalent to equation (8) under the assumption that we have a
4 5 Thilo Koch stationary flow and neglecting P (usually P < < G,T ); it is also helpful in understanding how v is influenced by flow height H and slope 5. As a product of 0 (in equation (8)), we obtain the material constant C ; K and & are exponents and products of P and A. respectively, while f can contain terms of x }, and/or x c. Equation (11) is written in a non-discrete form because it is not used in the simulation routine. Table 1 lists five of the tested constitutive equations by presenting the values of P, X, K, 9 in equation (8) or (11) and marking the existence of x y and x c with a cross. The equivalent material constants are listed in 0. As previously mentioned, a Bingham laminar flow cannot be described with equation (8). Under the assumption that the Bingham yield stress is negligible (here x c -> 0), a flow resistance law for a Newtonian laminar fluid (Table 1) is obtained. The new velocity v* +l/ is obtained from the old velocity v*~" with gravity g, by solving the equation of motion: * + l/ = *-l/ + F k At r - ^ ( 1 ) The steps are repeated until the debris flow comes to rest or reaches the end of the path. Table 1 Comparison of flow resistance laws described in equations (8) or (11). Equations (8) (11) -(8), (11)- Approaches P X C K & T t c Newtonian laminar Dilatant inertial Newtonian turbulent Voellmy fluid KI f-i n C 1 K n- 1 % n l C KSB The approaches are described by using the values of the exponents p and X, K and S, respectively and, indicated with a cross, whether x y and/or x c are included. The material constants, equivalent to (equation (8)) are listed: \i= viscosity, n Manning's roughness coefficient, 4= lumped coefficient accounting for grain and concentration properties in granular flow, C = Chézy roughness coefficient, and K = coefficient from empirical investigations. SIMULATIONS AND FIELD OBSERVATIONS Comprehensive comparisons between model simulations and observations are made for debris flow events in the Kamikamihori Valley in Japan (Rickenmann & Koch, 1997). The longitudinal and bed-slope-angle profile of the Kamikamihori Valley used in the calculations is shown in Fig. (check dams were smoothed and the step-height sum reduced by 19 m). In the calculations, the model channel has a constant width of 10 m down to a fan apex at X = 1.9 km. From there, it widens to 15 m at X =.0 km; further downstream there is a widening of 10 m per 100 m longitudinal distance. This is a rough approximation of the known cross-sections.
5 Testing various constitutive equations for debris flow modelling 53 - Kamikamihori valley Slope angle [ ] 30 Middle-reach observation site CO 0- (D 0) =3 CQ_ 10 CD Zi i r Fig. Longitudinal and bed-slope profile of the Kamikamihori Valley (Japan). The arrow indicates where flow hydrograph had been measured. The model parameters have been adjusted to best reproduce the measured data on flow velocity along the flow path. For calculating the earth pressure coefficient, I set 5 = 4 and <b = 10, where 8 is found in equations with a XQ term (Table 1) to produce a stopping position of the wave front close to the measured value in the field and reasonable -values; <j> is assumed as there is only little influence of <j> on the results for long and fiat waves. The events of 3 August 1976 (Okuda et al., 1980), case A, and of 1 September 1979 (Suwa & Okuda, 1981), case B, are used here. Average debris wave-front velocities were obtained for segments along the flow path measured between signal wires tightened perpendicular to the torrent channel. Because local velocity variations are not known, I assumed an uncertainty margin of ±0% with regard to the measured values. For case A, the initial input hydrograph at X 0 = 0 km has a maximum flow depth H Q of.5 m, a length, L 0 of 70 m and a constant inflow velocity, v 0 of 6.5 m s" 1. The total volume, M 0 (sediment and water), corresponds to 6500 m 3. For each of the five approaches listed in Table 1, one simulation is shown in Fig. 3. The simulated velocities are reasonably close to the observed ones for the turbulent, the Voellmy and the KI models. For the Bingham and the dilatant models no good agreement could be obtained for the velocity development along the whole flow path. Simulations with the Newtonian laminar model show results similar to the Bingham model. While using values in a reasonable range, the first term in the brackets of equation (8) dominates the flow resistance for higher velocities, while the yield stress terms (x y,x c ) dominate the stopping process. In the case of the Voellmy and KI models, terms x y,x c produce a stopping wave in the range of the observed runout distance for the chosen parameters. For case B, the conditions of the input hydrograph are: X Q = 0.8 km, H 0 = 4.5 m, L 0 = 640 m, v 0 = 5 m s" 1, and M Q = m 3. Figure 4 shows the results of the same approaches, but with different material parameters. Again, the simulated velocities of the turbulent, the Voellmy and the KI models show an acceptable agreement with the measured ones, using almost the same flow resistance parameters as for case A. With the Bingham model, the velocity development can be reproduced somewhat better than for case A; however, the simulations started at a smaller bed
6 54 Thilo Koch turbulent: n = 0.15 Voellmy: 5 = 3.5, C? = 10 ms- - Bingham: T B= 1000 Pa, u^= 800 Pa s - dilatant: 5 = 31 rrrf- 5 s-i V observed (+/- 0%), 3 Aug 1976 y Fig. 3 Simulated and observed velocities for event A in the Kamikamihori Valley. The arrows indicate the range of observed runout distances. slope (Fig. ) than for case A. This and further experience with other events seem to confirm that the Bingham and the dilatant models cannot reproduce the flow behaviour for strongly varying slope angles. SIMULATION SENSITIVITY TO BOUNDARY AND INITIAL CONDITIONS The initial conditions and the exact cross-sections of the observed Kamikamihori events are poorly known. Therefore additional calculations have been made considering variations in the initial and boundary conditions. Figure 5 compares front velocities (case A event) for a turbulent model (reference case I: n = 0.15; turbulent model, see equation (11), with K = 0.67 (Table 1)) under the following conditions (original value in brackets): Case II: Input hydrograph: Max. depth H Q = 5.5 [.5] m, length L 0 = 135 [70] Case III: Case IV: Case V: Case VI: m. Volume: V 0 = 9735 [6500] m 3, H 0 unchanged, L 0 = 400 [70] m. Channel width: 8 [10] m, depth H 0 = 3 [.5] m, V 0 unchanged. More detailed channel slope profile, different channel widening on fan: at X = 160 m, widening to 0 m at X = 170 and at X = 1860 m additional widening of 6 m per 100 m longitudinal distance. Using a hydraulic radius calculated for a rectangular channel profile. Inflow velocity: v 0 = 5 [6.5] m s' 1. Case VII Figure 5(a) shows that using a different inflow depth but a constant volume (case II) strongly affects the front velocity during the initial phase. For X > 0.9 km, the front velocities of case II and the reference simulation are nearly identical. Using a larger volume (case III) or a smaller channel width (case IV) produces higher velocities along the whole flow path. Both cases result in higher flow depths, which yield the
7 Testing various constitutive equations for debrisflowmodelling I 8- f 6- > e 4- u_ - / i-. 11 i turbulent: n = 0.17 Vœllmy: ô=4.4, C = 10 m s- Bingham: T B= 800 Pa, 1^= 400 Pa s dilatant,=4 rrr 05 s- 1 Kl: K = m-33 S3 T = 400 Pa, 5 = 5 y \ ' "". _ V observed (+/-0'/o), 1 Sept 1979 ^M U "J 1 I I A A Fig. 4 Simulated and observed velocities for event B in the Kamikamihori Valley. The arrows indicate the range of observed runout distances. higher velocities. Figure 5(b) shows that using a hydraulic radius (case VI) or a smaller inflow velocity (case VII) obviously affects only the initial phase. Using a more detailed slope angle profile and a different channel widening on the fan (case V) produces local velocity variations which might be an improvement. A special effect is observed in the Bingham model: when there is a higher initial flow depth the wave becomes elongated, resulting in smaller flow heights downstream and thus a shorter travel distance. DISCUSSION Simulations for case A, case B and events not presented here (Saas Valley in Switzerland, lahars from Pine Creek and Muddy River at Mount St Helens, USA) lead to complementary results. Considering the front velocity along the flow path at steep variable slopes, there is better agreement between simulations and field observations for the Newtonian turbulent, the Voellmy and the KI models than for the dilatant, the Bingham and the Newtonian laminar models. Referring to Table 1, it can be seen that the better agreement is associated with flow resistance laws having smaller exponents K and &<0.5 in equation (11). Further evaluations, including models not presented here, confirm the preference of resistance laws close to turbulent clear water flow. For this reasons, models with smaller exponents seem to be more appropriate for simulations in torrents on steep and variable slopes. For less steep and less variable slopes, the differences in the results of the tested models is less significant. The Bingham model as well as the dilatant model may possibly also be suited to simulations of the flow behaviour on fans, where the bed slopes do not vary over too big a range. Simulations are generally more sensitive to input conditions, in particular to initial flow depth, for those flow resistance laws having
8 56 Thilo Koch a)? 1? fc 10.rT> n elo > c o M Cases I - IV: turbulent: n=0.15 Case I (reference simulation) Case II (new hydrograph) Case III (bigger volume) Case IV (smaller channel width) V observed (+/- 0%), 3 Aug 1976 T T (b) >>1U Cases I, V - VII: turbulent: n=0.15 Case I (reference simulation) Case V (detailed slope angle) Case VI (hydraulic radius) Case VII (smaller inflow velocity) W-S, V observed (+/- 0%), 3 Aug 1976 i 1 1 r Fig. 5 Simulated and observed front velocities for case A (debris flow events in the Kamikamihori Valley). All simulations are done with the turbulent model but for different initial and boundary conditions, (a) cases I through IV, (b) cases V through VII. higher exponents K, i.e. the Bingham, the laminar and the dilatant models. On the other hand, for the turbulent, the Voellmy and the KI models, variations in some initial conditions are damped relatively rapidly with increasing flow length. From a practical point of view, when performing simulations of debris flow events in nature, knowledge on initial and boundary conditions is only rough. Sensitivity examinations show relatively little influence of initial inflow velocity and height on the downstream flow behaviour for the preferred models. This also applies when considering side-wall effects (using a hydraulic radius) or a very detailed slope profile. Beside selecting an appropriate material constant, 0, the most accuracy is required for estimating the total volume and the channel width. These are the dominant factors controlling the dynamic behaviour of the tested models. The selection of the bed friction angle, S, and the internal angle of friction,, (equation (7)), also have only little influence on the flow behaviour. On the other hand, S in the Coulomb yield strength (equation (9)) is essential for the runout distance, which is very sensitive to values tan 5, close to tan a, in the deposition zone. 8 appears in the Voellmy, the KI and the Bingham Coulomb viscous flow. The Newtonian, dilatant and turbulent models do not simulate the stoppage. Connecting the turbulent model with a Coulomb yield strength solves this problem and produces a flow behaviour close to the Voellmy model. Acknowledgement I thank my colleague Dieter Rickenmann for the helpful discussions.
9 Testing various constitutive equations for debrisflowmodelling 57 REFERENCES Hungr, O. (1995) A model for the runout analysis of rapid flow slides, debris flows, and avalanches. Can. Geotech. J. 3, Hutter, K. & Koch, T. (1991) Motion of a granular avalanche in an exponentially curved chute: experiments and theoretical predictions. Phil. Trans. Roy. Soc. Lond. A334, Jan, C. D. & Shen, H. W. (1993) A review of debris flow analysis. In: Proc. XXVIAHR Congress, Technical Session B, vol. Ill, 5-3. Koch, T., Grève, R. & Hutter, K. (1994) Unconfined flow of granular avalanches along a partly curved surface, Part II: Experiments and numerical computations. Proc. Roy. Soc. Lond. A44S, Okuda, S., Suwa, EL, Okunishi, K., Yokoyama, K. & Nakano, M. (1980) Observations on the motion of a debris flow and its geomorphological effects. Z. Geomorph. N.F., Suppl. Bd 35, Rickenmann, D. & Koch, T. (1997) Comparison of debris flow modelling approaches. In: Debris-Flow Hazards Mitigation (Proc. Int. ASCE Conf., San Francisco), Savage, S. B & Hutter, K. (1991) The dynamics of avalanches of granular material from initiation to runout, Part I: Analysis. Acta Mechanica 86, Suwa, H. & Okuda, S. (1981) Topographical change caused by debris flow in Karnikamihori valley, Northern Japan Alps. Trans. Jap. Geomorph. Un. -,
POLITECNICO DI TORINO
POLITECNICO DI TORINO Whatever is the numerical approach to the study of rock avalanche evolution, obtained results depend on the choice of the value that is assigned to the characteristic parameters of
More informationFrequency of Rock Avalanches
Frequency of Rock Avalanches Reference Location Area (1 3 km 2 ) Geology Study Period Number of Cases Frequency 1 (cases/year /1, km 2 ) Return period (years per 1, km 2 Abele (1974) North calcareous Alps
More informationDry granular flows: gas, liquid or solid?
Dry granular flows: gas, liquid or solid? Figure 1: Forterre & Pouliquen, Annu. Rev. Fluid Mechanics, 2008 1 Characterizing size and size distribution Grains are not uniform (size, shape, ) Statistical
More informationLecture 3: Fundamentals of Fluid Flow: fluid properties and types; Boundary layer structure; unidirectional flows
GEOL 440 Sedimentology and stratigraphy: processes, environments and deposits Lecture 3: Fundamentals of Fluid Flow: fluid properties and types; Boundary layer structure; unidirectional flows Why study
More informationA discussion on the velocity of debris flow
Erosion, Debris Flows and Environment in Mountain Regions (Proceedings of the Chengdu Symposium, July 1992). IAHS Publ. no. 209, 1992. 369 A discussion on the velocity of debris flow ZHU PENGCHENG Institute
More informationPhysical modelling of sediment transport in mountain torrents upstream of open check dams
Physical modelling of sediment transport in mountain torrents upstream of open check dams Authors: Sebastian SCHWINDT Dr. Mário J. FRANCA Check dam in the region of Trent (Italy) Paper Code: EGU2015-6166
More informationLatest Trends in Engineering Mechanics, Structures, Engineering Geology
Numerical modeling of a debris flow event occurred in Campania region, southern Italy: consideration on the rheological model parameters on the run-out ANNA MARIA PELLEGRINO Department of Engineering University
More informationSIMULATION OF A 2D GRANULAR COLUMN COLLAPSE ON A RIGID BED
1 SIMULATION OF A 2D GRANULAR COLUMN COLLAPSE ON A RIGID BED WITH LATERAL FRICTIONAL EFFECTS High slope results and comparison with experimental data Nathan Martin1, Ioan Ionescu2, Anne Mangeney1,3 François
More informationSPH Molecules - a model of granular materials
SPH Molecules - a model of granular materials Tatiana Capone DITS, Univeristy of Roma (la Sapienza) Roma, Italy Jules Kajtar School of Mathematical Sciences Monash University Vic. 3800, Australia Joe Monaghan
More informationDebris flows + sediment transport in torrents. Debris flows and sediment. transport in steep catchments
International Workshop "Erosion, Transport Debris flows and sediment transport in steep catchments Dieter Rickenmann WSL - Swiss Federal Research Institute, Mountain Hydrology and Torrents, Birmensdorf,
More informationSwiss Avalanche-Dynamics Procedures for Dense Flow Avalanches
Swiss Avalanche-Dynamics Procedures for Dense Flow Avalanches H. Gubler 1 Swiss Avalanche-Dynamics Procedures for Dense Flow Avalanches H.Gubler AlpuG Richtstattweg2 CH-7270 Davos Platz, Switzerland Introduction
More informationQUASI-THREE DIMENSIONAL TWO-PHASE DEBRIS FLOW MODEL ACOUNTING FOR BOULDER TRANSPORT
DOI: 10.4408/IJEGE.2011-03.B-051 QUASI-THREE DIMENSIONAL TWO-PHASE DEBRIS FLOW MODEL ACOUNTING FOR BOULDER TRANSPORT C.E. MARTINEZ (*), F. MIRALLES-WILHELM (**) & R. GARCIA-MARTINEZ (***) (*) Department
More informationDevelopment of Kanako, a wide use 1-D and 2-D debris flow simulator equipped with GUI
Monitoring, Simulation, Prevention and Remediation of Dense Debris Flows II 49 Development of Kanako, a wide use 1-D and 2-D debris flow simulator equipped with GUI K. Nakatani, T. Wada, Y. Satofuka &
More informationLandslides & Debris Flows
T.#Perron# #12.001# #Landslides#&#Debris#Flows# 1# Landslides & Debris Flows Many geologic processes, including those shaping the land surface, are slowacting, involving feedbacks that operate over many
More informationDebris flow modelling accounting for large boulder transport
Monitoring, Simulation, Prevention and Remediation of Dense and Debris lows III 9 Debris flow modelling accounting for large boulder transport C. Martinez 1,. Miralles-Wilhelm 1 & R. Garcia-Martinez 1
More informationFluid Mechanics Prof. S.K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur
Fluid Mechanics Prof. S.K. Som Department of Mechanical Engineering Indian Institute of Technology, Kharagpur Lecture - 42 Flows with a Free Surface Part II Good morning. I welcome you to this session
More informationNumerical analysis of snow avalanche mechanics and of its interaction with structures
Numerical analysis of snow avalanche mechanics and of its interaction with structures Eloïse Bovet 1, Bernardino Chiaia 1, Luigi Preziosi 2 1 Department of Structural and Geotechnical Engineering, Politecnico
More informationExperimentally determined distribution of granular-flow characteristics in collisional bed load transport
Experimentally determined distribution of granular-flow characteristics in collisional bed load transport Václav Matoušek 1,*, Štěpán Zrostlík 1, Luigi Fraccarollo 2, Anna Prati 2, and Michele Larcher
More informationSome considerations on the application of the FLO-2D model for debris flow hazard assessment
Monitoring, Simulation, Prevention and Remediation of Dense and Debris Flows 159 Some considerations on the application of the FLO-2D model for debris flow hazard assessment V. D Agostino 1 & P. R. Tecca
More informationCHAPTER 2- BACKGROUND. INVESTIGATIONS OF COMPOSITE ROUGHNESS COEFFICIENT IN A RIVER WITH LOW FLOW
2. Background 2.1 Introduction The estimation of resistant coefficient and hence discharge capacity in a channel or river is one of the fundamental problems facing river engineers. When applying Manning
More informationIMPLEMENTATION OF NON-NEWTONIAN RHEOLOGY FOR GRANULAR FLOW SIMULATION
9. National Conference on Computational Mechanics MekIT 17 B. Skallerud and H I Andersson (Eds) IMPLEMENTATION OF NON-NEWTONIAN RHEOLOGY FOR GRANULAR FLOW SIMULATION Petter Fornes 1,, Hans Bihs 1 and Steinar
More informationMass Wasting. Requirements for Mass Wasting. Slope Stability. Geol 104: mass wasting
Mass Wasting Movement of earth materials downslope, driven by Gravitational Forces. Landslides - general term for rock or soil movement. In U.S., on average, mass wasting causes 1 to 2 billion dollars
More informationch-01.qxd 8/4/04 2:33 PM Page 1 Part 1 Basic Principles of Open Channel Flows
ch-01.qxd 8/4/04 2:33 PM Page 1 Part 1 Basic Principles of Open Channel Flows ch-01.qxd 8/4/04 2:33 PM Page 3 Introduction 1 Summary The introduction chapter reviews briefly the basic fluid properties
More informationMulti-step hazard assessment of debris flows in an Alpine region
Multi-step hazard assessment of debris flows in an Alpine region R. Genevois Department of Geosciences, University of Padova, Padova, Italy P.R. Tecca CNR-IRPI, Padova, Italy M. Floris, C. Squarzoni &
More informationHydraulics for Urban Storm Drainage
Urban Hydraulics Hydraulics for Urban Storm Drainage Learning objectives: understanding of basic concepts of fluid flow and how to analyze conduit flows, free surface flows. to analyze, hydrostatic pressure
More informationFE Fluids Review March 23, 2012 Steve Burian (Civil & Environmental Engineering)
Topic: Fluid Properties 1. If 6 m 3 of oil weighs 47 kn, calculate its specific weight, density, and specific gravity. 2. 10.0 L of an incompressible liquid exert a force of 20 N at the earth s surface.
More informationOpen Channel Flow Part 2. Ch 10 Young, notes, handouts
Open Channel Flow Part 2 Ch 10 Young, notes, handouts Uniform Channel Flow Many situations have a good approximation d(v,y,q)/dx=0 Uniform flow Look at extended Bernoulli equation Friction slope exactly
More informationMeasurements of dense snow avalanche basal shear to normal stress ratios (S/N)
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 34, L07501, doi:10.1029/26gl028670, 27 Measurements of dense snow avalanche basal shear to normal stress ratios (S/N) K. Platzer, 1 P. Bartelt,
More informationI NTERREG I VC alcotra Marseille, January 27 th École Polytechnique Fédérale de Lausanne Laboratory for Rock Mechanics (LMR)
POLITECNICO DI TORINO I NTERREG I VC alcotra 27 23 Numerical modelling w ith a continuum mechanics approach: updatings (2) Claudio Scavia, () Vincent Labiouse, (2) Marina Pirulli, () Claire Sauthier, (2,)
More informationDebris flow: categories, characteristics, hazard assessment, mitigation measures. Hariklia D. SKILODIMOU, George D. BATHRELLOS
Debris flow: categories, characteristics, hazard assessment, mitigation measures Hariklia D. SKILODIMOU, George D. BATHRELLOS Natural hazards: physical phenomena, active in geological time capable of producing
More informationPrediction of landslide-induced debris flow hydrograph: the Atsumari debris flow disaster in Japan
Monitoring, Simulation, Prevention and Remediation of Dense and Debris Flows 27 Prediction of landslide-induced debris flow hydrograph: the Atsumari debris flow disaster in Japan H. Takaoka 1, H. Hashimoto
More informationA STUDY ON DEBRIS FLOW DEPOSITION BY THE ARRANGEMENT OF SABO DAM
Annual Journal of Hydraulic Engineering, JSCE, Vol.57, 2013, February A STUDY ON DEBRIS FLOW DEPOSITION BY THE ARRANGEMENT OF SABO DAM Namgyun Kim 1, Hajime NAKAGAWA 2, Kenji KAWAIKE 3, and Hao ZHANG 4
More information5. MODELING OF NON-STRATIFIED MIXTURE FLOWS (Pseudo-homogeneous flows)
5. MODELING OF NON-STRATIFIED MIXTURE FLOWS (Pseudo-homogeneous flows) Uniform (or almost uniform) distribution of transported solids across a pipeline cross section is characteristic of pseudo-homogeneous
More informationFactors affecting confluence scour
& Wang (eds) River Sedimentation 1999., Balkema, Rotterdam. ISBN 9 9 3. 17 19 Factors affecting confluence scour R. B. Rezaur & A. W. Jayawardena. Department of Civil Engineering, The University of Hong
More informationSegregation Effects in Granular Collapses
Segregation Effects in Granular Collapses Joshua S. Caplan, Saran Tunyasuvunakool, Stuart B. Dalziel & Nathalie M. Vriend Department of Applied Mathematics and Theoretical Physics University of Cambridge,
More information1. The Properties of Fluids
1. The Properties of Fluids [This material relates predominantly to modules ELP034, ELP035] 1.1 Fluids 1.1 Fluids 1.2 Newton s Law of Viscosity 1.3 Fluids Vs Solids 1.4 Liquids Vs Gases 1.5 Causes of viscosity
More informationAdvanced Hydraulics Prof. Dr. Suresh A. Kartha Department of Civil Engineering Indian Institute of Technology, Guwahati
Advanced Hydraulics Prof. Dr. Suresh A. Kartha Department of Civil Engineering Indian Institute of Technology, Guwahati Module - 2 Uniform Flow Lecture - 1 Introduction to Uniform Flow Good morning everyone,
More informationNumerical Study of Relationship Between Landslide Geometry and Run-out Distance of Landslide Mass
Numerical Study of Relationship Between Landslide Geometry and Run-out Distance of Landslide Mass Muneyoshi Numada Research Associate, Institute of Industrial Science, The University of Tokyo, Japan Kazuo
More informationComparison of three landslide runout models on the Turnoff Creek rock avalanche, British Columbia
Landslide Processes: From geomorphologic mapping to dnamic modelling. A tribute to Prof. Dr. Theo van Asch CERG, UU & US - Strasbourg, 6-7 Februar 2009 Comparison of three landslide runout models on the
More informationMODELING OF RUNOUT LENGTH OF HIGH-SPEED GRANULAR MASSES
DOI: 10.4408/IJEGE.2013-06.B-12 MODELING OF RUNOUT LENGTH OF HIGH-SPEED GRANULAR MASSES Francesco FEDERICO & Chiara CESALI University of Rome Tor Vergata - Rome, Italy ABSTRACT The power balance of a high
More informationVARIATION OF MANNING S ROUGHNESS COEFFICIENT WITH SEEPAGE IN SAND-BED CHANNEL *Satish Patel 1 and Bimlesh Kumar 2
International Journal of Science, Environment and Technology, Vol. 5, No 6, 2016, 3678 3685 ISSN 2278-3687 (O) 2277-663X (P) VARIATION OF MANNING S ROUGHNESS COEFFICIENT WITH SEEPAGE IN SAND-BED CHANNEL
More informationthree-dimensional topographic data
Annals ofglaciology 18 1993 International Glaciological Society Analysis of powder-snow avalanches using three-dimensional topographic data YUSUKE FUKUSHIMA AND NORIO HAYAKAWA Department of Civil Engineering,
More informationENGINEERING GEOLOGY. Landslides by U.S. Region. WV Has 1st or 2nd Highest Landslide Damage Per Capita $ /person/yr UTAH
Essay Questions for 2007 http://www.geo.wvu.edu/~kite/2007test1ques321.html New Topic: SLOPES Colluvial - Related to Hill-Slope or Mountain-Slope Processes, Especially Those Driven Directly by Gravity
More informationMODELING FLUID FLOW IN OPEN CHANNEL WITH HORSESHOE CROSS SECTION
July. 2. Vol. 7. No. 2 MODELING FLUID FLOW IN OPEN CHANNEL WITH HORSESHOE CROSS SECTION 1 J. JOMBA, 2 D.M.THEURI, 2 E. MWENDA, 2 C. CHOMBA ABSTRACT Flow in a closed conduit is regarded as open channel
More informationCOMPARISON OF TRANSPORT AND FRICTION OF MONO- SIZED AND TWO-SPECIES SEDIMENT IN UPPER PLANE BED REGIME
ISBN 978-83-927084-8-3 ISSN 0867-7964 COMPARISON OF TRANSPORT AND FRICTION OF MONO- SIZED AND TWO-SPECIES SEDIMENT IN UPPER PLANE BED REGIME Štěpán Zrostlík, Vojtěch Bareš, Jan Krupička, Tomáš Picek, Václav
More informationEFFECT OF GRAIN DENSITY ON PLANE BED FRICTION. Václav Matoušek, Vojtěch Bareš, Jan Krupička, Tomáš Picek, Štěpán Zrostlík
ISBN 978-83-927084-8-3 ISSN 0867-7964 EFFECT OF GRAIN DENSITY ON PLANE BED FRICTION Václav Matoušek, Vojtěch Bareš, Jan Krupička, Tomáš Picek, Štěpán Zrostlík Czech Technical University in Prague, Faculty
More informationFlood Routing by the Non-Linear Muskingum Model: Conservation of Mass and Momentum
Archives of Hydro-Engineering and Environmental Mechanics Vol. 56 (29), No. 3 4, pp. 121 137 IBW PAN, ISSN 1231 3726 Flood Routing by the Non-Linear Muskingum Model: Conservation of Mass and Momentum Dariusz
More informationDebris flow hazard and mitigation works in Fiames slope (Dolomites, Italy)
Monitoring, Simulation, Prevention and Remediation of Dense and Debris Flows 15 Debris flow hazard and mitigation works in Fiames slope (Dolomites, Italy) 1 2 2 P. R. Tecca, C. Armento & R. Genevois 1
More informationTechnical Memorandum. To: From: Copies: Date: 10/19/2017. Subject: Project No.: Greg Laird, Courtney Moore. Kevin Pilgrim and Travis Stroth
Technical Memorandum To: From: Greg Laird, Courtney Moore Kevin Pilgrim and Travis Stroth 5777 Central Avenue Suite 228 Boulder, CO 80301 www.otak.com Copies: [Electronic submittal] Date: 10/19/2017 Subject:
More informationGEOMORPHIC CHANGES OF A LANDSLIDE DAM BY OVERTOPPING EROSION
DOI: 10.4408/IJEGE.2011-03.B-087 GEOMORPHIC CHANGES OF A LANDSLIDE DAM BY OVERTOPPING EROSION K. Yoshino (*), T. Uchida (*), T. shimizu (*) & K. Tamura (*) (*) Incorporated Administrative Agency, Public
More informationFlow and Bed Topography in a 180 Curved Channel
Flow and Bed Topography in a 180 Curved Channel Jae Wook Jung 1, Sei Eui Yoon 2 Abstract The characteristics of flow and bed topography has been analyzed by changing the bed materials in a 180-degree,
More informationDAM-BREAK FLOWS OF DRY GRANULAR MATERIAL ON GENTLE SLOPES
DOI: 10.4408/IJEGE.2011-03.B-056 DAM-BREAK FLOWS OF DRY GRANULAR MATERIAL ON GENTLE SLOPES Luca SARNO (*), Maria Nicolina PAPA (**) & Riccardo MARTINO (*) (*) University of Naples Federico II - Dept. of
More informationLecture 6: Flow regimes fluid-like
Granular Flows 1 Lecture 6: Flow regimes fluid-like Quasi-static granular flows have plasticity laws, gaseous granular flows have kinetic theory -- how to model fluid-like flows? Intermediate, dense regime:
More informationDry and wet granular flows. Diego Berzi
Dry and wet granular flows Diego Berzi Outline 2 What? Why? How? When? Who? Where? Then? What? Granular flows many solid moving particles 3 particle mass is large (at least 10 20 molecular masses). Hence,
More informationTHE HYDRAULIC PERFORMANCE OF ORIENTED SPUR DIKE IMPLEMENTATION IN OPEN CHANNEL
Tenth International Water Technology Conference, IWTC10 2006, Alexandria, Egypt 281 THE HYDRAULIC PERFORMANCE OF ORIENTED SPUR DIKE IMPLEMENTATION IN OPEN CHANNEL Karima Attia 1 and Gamal El Saied 2 1
More information6. Basic basic equations I ( )
6. Basic basic equations I (4.2-4.4) Steady and uniform flows, streamline, streamtube One-, two-, and three-dimensional flow Laminar and turbulent flow Reynolds number System and control volume Continuity
More informations and FE X. A. Flow measurement B. properties C. statics D. impulse, and momentum equations E. Pipe and other internal flow 7% of FE Morning Session I
Fundamentals of Engineering (FE) Exam General Section Steven Burian Civil & Environmental Engineering October 26, 2010 s and FE X. A. Flow measurement B. properties C. statics D. impulse, and momentum
More informationESTIMATION OF MORPHOLOGICAL IMPACT OF GROYNE LENGTHENING I. RÁTKY, ÉVA RÁTKY
ESTIMATION OF MORPHOLOGICAL IMPACT OF GROYNE LENGTHENING I. RÁTKY, ÉVA RÁTKY Abstract. Hydraulic-morphological calculations in open channel flows still cause problems for modellers, partially because of
More informationNPTEL Quiz Hydraulics
Introduction NPTEL Quiz Hydraulics 1. An ideal fluid is a. One which obeys Newton s law of viscosity b. Frictionless and incompressible c. Very viscous d. Frictionless and compressible 2. The unit of kinematic
More informationAn-Najah National University Civil Engineering Department. Fluid Mechanics. Chapter 1. General Introduction
1 An-Najah National University Civil Engineering Department Fluid Mechanics Chapter 1 General Introduction 2 What is Fluid Mechanics? Mechanics deals with the behavior of both stationary and moving bodies
More informationAbstract. 1 Introduction
One-dimensional unsteady flow computation in channels with floodplains D. Bousmar, R. Scherer & Y. Zech Civil Engineering Dept., Universite Catholique de Louvain, Place du Levant, 1, B-1348 Louvain-la-Neuve,
More informationOPEN CHANNEL FLOW. One-dimensional - neglect vertical and lateral variations in velocity. In other words, Q v = (1) A. Figure 1. One-dimensional Flow
OPEN CHANNEL FLOW Page 1 OPEN CHANNEL FLOW Open Channel Flow (OCF) is flow with one boundary exposed to atmospheric pressure. The flow is not pressurized and occurs because of gravity. Flow Classification
More informationHelical Coil Flow: a Case Study
Excerpt from the Proceedings of the COMSOL Conference 2009 Milan Helical Coil Flow: a Case Study Marco Cozzini Renewable Energies and Environmental Technologies (REET) Research Unit, Fondazione Bruno Kessler
More informationAvalanches. Avalanche s
Avalanches Avalanche s Avalanches were first imagined as giant snowballs which increased in size from accretion of underlying snow What are avalanches? They are flows which move under the influence of
More information9 Stresses. 9.1 Stress states in silos
9 Stresses The knowledge of the stresses prevailing in bulk solids, especially when being stored in bins and silos, is extremely important when considering the following topics: Silo design for flow (e.g.,
More informationThe SedAlp Project: WP6: INTERACTION WITH STRUCTURES
The SedAlp Project: WP6: INTERACTION WITH STRUCTURES 6 th International conference Water in the Alps Herrenchiemsee, 12 th of October, 2016 Jošt Sodnik Contents Structure of SedAlp project (www.sedalp.eu)
More informationIntroduction to Marine Hydrodynamics
1896 1920 1987 2006 Introduction to Marine Hydrodynamics (NA235) Department of Naval Architecture and Ocean Engineering School of Naval Architecture, Ocean & Civil Engineering First Assignment The first
More informationHow do landscape materials get from mountain tops to valley floors?
How do landscape materials get from mountain tops to valley floors? The processes that move materials into stream, creeks, and rivers are collectively called mass movements or mass wasting. This includes
More informationRapid Mass Movements System RAMMS
Rapid Mass Movements System RAMMS Yves Bühler, Marc Christen, Perry Bartelt, Christoph Graf, Werner Gerber, Brian McArdell Swiss Federal Institute for Forest, Snow and Landscape Research WSL WSL Institute
More informationModelling of dispersed, multicomponent, multiphase flows in resource industries. Section 3: Examples of analyses conducted for Newtonian fluids
Modelling of dispersed, multicomponent, multiphase flows in resource industries Section 3: Examples of analyses conducted for Newtonian fluids Globex Julmester 017 Lecture # 04 July 017 Agenda Lecture
More informationIMPLICIT NUMERICAL SCHEME FOR REGULATING UNSTEADY FLOW IN OPEN CHANNEL Mohamed. T. Shamaa 1, and Hmida M. Karkuri 2
IMPLICIT NUMERICAL SCHEME FOR REGULATING UNSTEADY FLOW IN OPEN CHANNEL Mohamed. T. Shamaa 1, and Hmida M. Karkuri 2 1 Associated Professor, Irrigation and Hydraulic Department, College of Technical Engineering,
More informationMEASUREMENT OF 3D FLOW FIELD IN A 90 BEND WITH ULTRASONIC DOPPLER VELOCITY PROFILER
MEASUREMENT OF 3D FLOW FIELD IN A 90 BEND WITH ULTRASONIC DOPPLER VELOCITY PROFILER Daniel S. Hersberger 1 1 Research associate, Laboratory of Hydraulic Constructions (LCH), Swiss Federal Institute of
More informationSpreading of kaolin and sand mixtures on a horizontal plane: physical experiments and SPH numerical modelling
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 00 (2016) 000 000 www.elsevier.com/locate/procedia 1st International Conference on the Material Point Method, MPM 2017 Spreading
More information1. Introduction, fluid properties (1.1, 2.8, 4.1, and handouts)
1. Introduction, fluid properties (1.1, 2.8, 4.1, and handouts) Introduction, general information Course overview Fluids as a continuum Density Compressibility Viscosity Exercises: A1 Fluid mechanics Fluid
More informationGRANULAR FLOWS AND NUMERICAL MODELLING OF LANDSLIDES
DAMOCLES DEBRISFALL ASSESSMENT IN MOUNTAIN CATCHMENTS FOR LOCAL END-USERS Contract No EVG1 - CT-1999-00007 GRANULAR FLOWS AND NUMERICAL MODELLING OF LANDSLIDES Giovanni B. Crosta 1, Francesco Calvetti
More informationFluvial Dynamics. M. I. Bursik ublearns.buffalo.edu October 26, Home Page. Title Page. Contents. Page 1 of 18. Go Back. Full Screen. Close.
Page 1 of 18 Fluvial Dynamics M. I. Bursik ublearns.buffalo.edu October 26, 2008 1. Fluvial Dynamics We want to understand a little of the basic physics of water flow and particle transport, as so much
More informationLecture Note for Open Channel Hydraulics
Chapter -one Introduction to Open Channel Hydraulics 1.1 Definitions Simply stated, Open channel flow is a flow of liquid in a conduit with free space. Open channel flow is particularly applied to understand
More informationStudy of Pile Interval of Landslide Restraint Piles by Centrifuge Test and FEM Analysis
Disaster Mitigation of Debris Flows, Slope Failures and Landslides 113 Study of Pile Interval of Landslide Restraint Piles by Centrifuge Test and FEM Analysis Yasuo Ishii, 1) Hisashi Tanaka, 1) Kazunori
More informationPrinciples of Convection
Principles of Convection Point Conduction & convection are similar both require the presence of a material medium. But convection requires the presence of fluid motion. Heat transfer through the: Solid
More informationScaled laboratory experiments on the evolution of fluidised avalanches.
Scaled laboratory experiments on the evolution of fluidised avalanches. J. Koegl 1,, A. Graf 1, L. Rammer, K. Kleemayr 1, M.A. Kern 1, P. Gauer 1, G. Kapeller and M. Aufleger 1 Federal Research and Training
More information16 Rainfall on a Slope
Rainfall on a Slope 16-1 16 Rainfall on a Slope 16.1 Problem Statement In this example, the stability of a generic slope is analyzed for two successive rainfall events of increasing intensity and decreasing
More informationIntroduction and Fundamental Concepts (Lectures 1-7)
Introduction and Fundamental Concepts (Lectures -7) Q. Choose the crect answer (i) A fluid is a substance that (a) has the same shear stress at a point regardless of its motion (b) is practicall incompressible
More informationSoft Bodies. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies
Soft-Body Physics Soft Bodies Realistic objects are not purely rigid. Good approximation for hard ones. approximation breaks when objects break, or deform. Generalization: soft (deformable) bodies Deformed
More informationQuantifying the hazard of catastrophic rock avalanches
IRASMOS Symposium, Davos, Switzerland, 15-16 May 2008 Quantifying the hazard of catastrophic rock avalanches Acknowledging John J. Clague Ken Hewitt Reginald L. Hermanns Isaac J. Larsen Alexander L. Strom
More informationSand Ripple Dynamics on the Inner Shelf
Sand Ripple Dynamics on the Inner Shelf Donald N. Slinn Department of Civil and Coastal Engineering, University of Florida Gainesville, FL 32611-6590, Phone: (352) 392-9537 x 1431 Fax: (352) 392-3466 E-mail:
More informationOpen Channel Flow K Subramanya Solution Manual
We have made it easy for you to find a PDF Ebooks without any digging. And by having access to our ebooks online or by storing it on your computer, you have convenient answers with open channel flow k
More informationCFD Simulation in Helical Coiled Tubing
Journal of Applied Science and Engineering, Vol. 19, No. 3, pp. 267 272 (2016) DOI: 10.6180/jase.2016.19.3.04 CFD Simulation in Helical Coiled Tubing Z. Y. Zhu Department of Petroleum Engineering, China
More informationNUMERICAL SIMULATION OF EROSION PROCESSES ON CROSSBAR BLOCK RAMPS
E-proceedings of the 36 th IAHR World Congress NUMERICAL SIMULATION OF EROSION PROCESSES ON CROSSBAR BLOCK RAMPS MARIO OERTEL (1), JAN P. BALMES (2), DANIEL B. BUNG (3) (1) Hydraulic Engineering Section,
More informationFLUID MECHANICS PROF. DR. METİN GÜNER COMPILER
FLUID MECHANICS PROF. DR. METİN GÜNER COMPILER ANKARA UNIVERSITY FACULTY OF AGRICULTURE DEPARTMENT OF AGRICULTURAL MACHINERY AND TECHNOLOGIES ENGINEERING 1 5. FLOW IN PIPES 5.1.3. Pressure and Shear Stress
More informationLECTURE 2: FLUID MECHANICS
LECTURE 2: FLUID MECHANICS Introduction Conservation of mass and momentum General types of flow Laminar vs. turbulent flow Shear Stress Reach-average shear stress Bed roughness and reach average flow velocity
More informationRivers T. Perron
1 Rivers T. Perron 12.001 After our discussions of large-scale topography, how we represent topography in maps, and how topography interacts with geologic structures, you should be frothing at the mouth
More information150A Review Session 2/13/2014 Fluid Statics. Pressure acts in all directions, normal to the surrounding surfaces
Fluid Statics Pressure acts in all directions, normal to the surrounding surfaces or Whenever a pressure difference is the driving force, use gauge pressure o Bernoulli equation o Momentum balance with
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationDynamic Avalanche Modeling in Natural Terrain
Dynamic Avalanche Modeling in Natural Terrain J.T. Fischer,a,b, J. Kowalski a, S.P. Pudasaini b, S. A. Miller b a WSL Institute for Snow and Avalanche Research SLF, Davos Dorf, Switzerland b Department
More informationErosion Rate is a Function of Erodibility and Excess Shear Stress = k ( o - c ) From Relation between Shear Stress and Erosion We Calculate c and
Equilibrium, Shear Stress, Stream Power and Trends of Vertical Adjustment Andrew Simon USDA-ARS, Oxford, MS asimon@msa-oxford.ars.usda.gov Non-Cohesive versus Cohesive Materials Non-cohesive: sands and
More informationFLOW IN CONDUITS. Shear stress distribution across a pipe section. Chapter 10
Chapter 10 Shear stress distribution across a pipe section FLOW IN CONDUITS For steady, uniform flow, the momentum balance in s for the fluid cylinder yields Fluid Mechanics, Spring Term 2010 Velocity
More informationSedimentation Scour Model Gengsheng Wei, James Brethour, Markus Grünzner and Jeff Burnham August 2014; Revised October 2014
Flow Science Report 03-14 Sedimentation Scour Model Gengsheng Wei, James Brethour, Markus Grünzner and Jeff Burnham August 2014; Revised October 2014 1. Introduction The three-dimensional sediment scour
More informationWATER INJECTION DREDGING by L.C. van Rijn
WATER INJECTION DREDGING by L.C. van Rijn (info@leovanrijn-sediment.com) Description of method Almost all harbour basins suffer from the problem of siltation of sediments. Usually, the deposited materials
More informationNumerical Hydraulics
ETH Zurich, Fall 2017 Numerical Hydraulics Assignment 2 Numerical solution of shallow water wave propagation (www.surfertoday.com) 1 Introduction 1.1 Equations Understanding the propagation of shallow
More informationProf. B V S Viswanadham, Department of Civil Engineering, IIT Bombay
19 Module 5: Lecture -1 on Stability of Slopes Contents Stability analysis of a slope and finding critical slip surface; Sudden Draw down condition, effective stress and total stress analysis; Seismic
More information